Prediction of pH Effect on the Octanol–Water Partition Coefficient of Ionizable Pharmaceuticals

The effects of commercial 2,4-D herbicide on game fish species: Natural lake water vs. laboratory system water

Aquatic herbicides with active ingredient 2,4-dichlorophenoxyacteic acid (2,4-D) are commonly used to control and combat aquatic non-native species that cause detrimental impacts including habitat destruction, strained resources among biota, and biodiversity loss. While many (eco)toxicology studies are performed in the laboratory under highly controlled circumstances, these studies may disregard the nuances and disorder that come with the complexity of natural aquatic ecosystems. Therefore, we conducted a series of laboratory experiments using laboratory system water, different lake waters, and different water parameters to determine the effects of ecologically relevant concentrations of 2,4-D (0.00-4.00 ppm a.e.) on the development and survival of two freshwater game species (Sander vitreus and Esox lucius). For 2,4-D exposures using different water sources, there were significant main effects of 2,4-D concentration and water source on walleye embryo and larval survival, however, there was no significant interaction between 2,4-D exposure and water source. For 2,4-D exposures and pH (5-9 pH), there were significant main effects of 2,4-D concentration and pH on walleye and northern pike embryo survival and a significant interaction between 2,4-D exposure and pH. Our results indicate that 2,4-D exposures in controlled laboratory system water can predict similar outcomes as 2,4-D exposures in natural lake water. Moreover, individual water parameters, such as pH, play a significant role in the toxicity of 2,4-D. Taken together, these results suggest that highly controlled laboratory studies are a useful tool for predicting impacts on survival of non-target fish in natural waters, but it is crucial for management agencies to consider individual water sources and specific lake water parameters in herbicide risk assessments to minimize the impacts to non-target organism.

A sustainable and highly efficient fossil-free carbon from olive stones for emerging contaminants removal from different water matrices

The olive stone is a large waste product of the olive oil extraction industry. The present study investigates developing activated carbon from olive stone waste (OSAC) to remove pharmaceuticals from water. Different temperatures and olive stone: KOH ratios were studied. The OSAC produced at 750 °C and 1:3 ratio was found to have the highest porosity and surface area and was tested in the adsorption process. Diclofenac and ciprofloxacin were selected as model contaminants. The adsorption process was optimized with regards to OSAC dosage, pH, temperature, and initial concentration of adsorbate. The OSAC was found to be effective for a wide pH range (2-11) with an optimum dosage of 1 g/L at 25 °C. The pharmaceuticals were almost completely removed in 75 min. The adsorption was endothermic and followed first-order kinetics with physical mechanisms such as electrostatic possibly being the main driver. The optimum conditions were applied to test the removal of diclofenac and ciprofloxacin in synthetic water, lake water (Lake Balaton, Hungary) and secondary wastewater for seven cycles. There was little difference between the removal of the tested water matrices highlighting the potency of OSAC as an adsorbent for pharmaceutical removal in industrial applications. The removal dropped from >99% for the first cycle to 20-30% for the seventh cycle.

COSMO-RS blind prediction of distribution coefficients and aqueous pKa values from the SAMPL8 challenge

The SAMPL8 blind prediction challenge, which addresses the acid/base dissociation constants (pKa) and the distribution coefficients (logD), was addressed by the Conductor like Screening Model for Realistic Solvation (COSMO-RS). Using the COSMOtherm implementation of COSMO-RS together with a rigorous conformational sampling, yielded logD predictions with a root mean square deviation (RMSD) of 1.36 log units over all 11 compounds and seven bi-phasic systems of the data set, which was the most accurate of all contest submissions (logD).For the SAMPL8 pKa competition, participants were asked to report the standard state free energies of all microstates, which were then used to calculate the macroscopic pKa. We have used COSMO-RS based linear free energy fit models to calculate the requested energies. The assignment of the calculated and experimental pKa values was made on the basis of the popular transitions, i.e. the transition hat was predicted by the majority of the submissions. With this assignment and a model that covers both, pKa and base pKa, we achieved an RMSD of 3.44 log units (18 pKa values of 14 molecules), which is the second place of the six ranked submissions. By changing to an assignment that is based on the experimental transition curves, the RMSD reduces to 1.65. In addition to the ranked contribution, we submitted two more data sets, one for the standard pKa model and one or the standard base pKa model of COSMOtherm. Using the experiment based assignment with the predictions of the two sets we received a RMSD of 1.42 log units (25 pKa values of 20 molecules). The deviation mainly arises from a single outlier compound, the omission of which leads to an RMSD of 0.89 log units.

Partitioning of Antioxidants in Edible Oil–Water Binary Systems and in Oil-in-Water Emulsions

In recent years, partitioning of antioxidants in oil–water two-phase systems has received great interest because of their potential in the downstream processing of biomolecules, their benefits in health, and because partition constant values between water and model organic solvents are closely related to important biological and pharmaceutical properties such as bioavailability, passive transport, membrane permeability, and metabolism. Partitioning is also of general interest in the oil industry. Edible oils such as olive oil contain a variety of bioactive components that, depending on their partition constants, end up in an aqueous phase when extracted from olive fruits. Frequently, waste waters are subsequently discarded, but their recovery would allow for obtaining extracts with antioxidant and/or biological activities, adding commercial value to the wastes and, at the same time, would allow for minimizing environmental risks. Thus, given the importance of partitioning antioxidants, in this manuscript, we review the background theory necessary to derive the relevant equations necessary to describe, quantitatively, the partitioning of antioxidants (and, in general, other drugs) and the common methods for determining their partition constants in both binary (PWOIL) and multiphasic systems composed with edible oils. We also include some discussion on the usefulness (or not) of extrapolating the widely employed octanol–water partition constant (PWOCT) values to predict PWOIL values as well as on the effects of acidity and temperature on their distributions. Finally, there is a brief section discussing the importance of partitioning in lipidic oil-in-water emulsions, where two partition constants, that between the oil-interfacial, POI, and that between aqueous-interfacial, PwI, regions, which are needed to describe the partitioning of antioxidants, and whose values cannot be predicted from the PWOIL or the PWOCT ones.

Impacts of wastewater effluents and seasonal trends on levels of antipsychotic pharmaceuticals in water and sediments from two cold-region rivers

Most pharmaceuticals are found at trace concentrations in aquatic systems, but their continuous release and potential accumulation can lead to adverse health effects in exposed organisms. Concentrations can vary temporally, driven by variations in discharges of receiving waters, sorption to sediments, and other biotic and abiotic exchange processes. The principal aim of this research was to better understand the occurrence, trends, and dynamics of pharmaceuticals in a cold-climate, riverine environment. To this end, a suite of seven representative antipsychotic pharmaceuticals was measured upstream and downstream of two wastewater treatment plants (WWTPs) in Saskatchewan, Canada, located in the South Saskatchewan River and Wascana Creek, respectively, across three seasons. Concentrations of analytes were in the ng/L range and generally greater downstream of both WWTPs compared to upstream. Some compounds, including the tricyclic antidepressant amitriptyline, which was the most abundant analyte in water and sediment from both sites and across seasons, reached low μg/L concentrations. Data collected from this research effort indicate contamination with antipsychotic pharmaceuticals, with the potential to adversely impact exposed organisms.

Synthesis, Characterization and Anticancer Efficacy Evaluation of Benzoxanthone Compounds toward Gastric Cancer SGC-7901

Three benzoxanthone derivatives were synthesized through a new photochemical strategy. The in vitro cytotoxic activity of these compounds was evaluated by 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay, and their partition coefficients (logP) were measured by shake flask method. The pK_a values of the compounds were detected by potentionmetric titration. The interaction of the compounds with calf thymus DNA (CT-DNA) was investigated by electronic absorption, luminescence spectra and viscosity. A molecular docking analysis was performed. The antitumor efficacy of the compounds was evaluated by cell apoptosis, cell cycle arrest, intracellular Ca²⁺ concentrations and reactive oxygen species (ROS) levels. The mitochondrial membrane potential was assayed using JC-1 (5,5′,6,6′-tetrachloro-1,1,3′,3′-tetraethyl-imidacarbocyanine iodide) as the fluorescence probe. The expression of Bcl-2 family protein, caspase 3 and poly ADP-ribose polymerase (PARP) was explored by western blot. The results showed that the compounds induced apoptosis through a ROS-mediated mitochondrial dysfunction pathway. This work provides an efficient approach to synthesize benzoxanthone derivatives, and is helpful for understanding the apoptotic mechanism.

Novel (Phenothiazinyl)Vinyl-Pyridinium Dyes and Their Potential Applications as Cellular Staining Agents

We report here the synthesis and structural characterization of novel cationic (phenothiazinyl)vinyl-pyridinium (PVP) dyes, together with optical (absorption/emission) properties and their potential applicability as fluorescent labels. Convective heating, ultrasound irradiation and mechanochemical synthesis were considered as alternative synthetic methodologies proficient for overcoming drawbacks such as long reaction time, nonsatisfactory yields or solvent requirements in the synthesis of novel dye (E)-1-(3-chloropropyl)-4-(2-(10-methyl-10H-phenothiazin-3-yl)vinyl)pyridin-1-ium bromide 3d and its N-alkyl-2-methylpyridinium precursor 1c. The trans geometry of the newly synthesized (E)-4-(2-(7-bromo-10-ethyl-10H-phenothiazin-3-yl)vinyl)-1-methylpyridin-1-ium iodide 3b and (E)-1-methyl-4-(2-(10-methyl-10H-phenothiazin-3-yl)vinyl)pyridin-1-ium tetrafluoroborate 3a′ was confirmed by single crystal X-ray diffraction. A negative solvatochromism of the dyes in polar solvents was highlighted by UV-Vis spectroscopy and explanatory insights were supported by molecular modeling which suggested a better stabilization of the lowest unoccupied molecular orbitals (LUMO). The photostability of the dye 3b was investigated by irradiation at 365 nm in different solvents, while the steady-state and time-resolved fluorescence properties of dye 3b and 3a′ in solid state were evaluated under one-photon excitation at 485 nm. The in vitro cytotoxicity of the new PVP dyes on B16-F10 melanoma cells was evaluated by WST-1 assay, while their intracellular localization was assessed by epi-fluorescence conventional microscopy imaging as well as one- and two-photon excited confocal fluorescence lifetime imaging microscopy (FLIM). PVP dyes displayed low cytotoxicity, good internalization inside melanoma cells and intense fluorescence emission inside the B16-F10 murine melanoma cells, making them suitable staining agents for imaging applications.

Estimation of the octanol-water distribution coefficient of acidic compounds by microemulsion electrokinetic chromatography

The feasibility of extending the determination of the lipophilicity of partially ionized acids (log D_o/w) by microemulsion electrokinetic chromatography (MEEKC) is tested. Theoretical considerations predict that a linear log D_o/w vs. log k correlation can be obtained only when the neutral and ionic forms of an acid follow the same correlation equation and the slope of the correlation is unity. In practice, since the lipophilicity of the neutral acid is much higher than that of the ionic form and the correlation slope is not very different from 1, the general linear correlation for neutral compounds can be applied across most of the ionization range of the acid. The linear correlation between log P_o/w and log k of 20 neutral solutes (calibration curve) has been established and extended to 6 acids used as models, tested across their full ionization range. log D_o/w-pH, and log k-pH profiles have been obtained for these 6 acids, and plotted log D_o/w against log k for any acid at any degree of ionization. Furthermore, the log D_o/w of the acids has been estimated from the calibration curve and log k-pH profile, and compared to values in the literature determined using reference methods such as the shake-flask one. Accurate values have been obtained using the MEEKC method when the acids are in their neutral form or partially ionized (ionization degree, α < 0.995). However, this parameter is overestimated when the acids are highly or fully ionized (α ≈ 1). Finally, in order to test the applicability of this method, we have applied the same procedure to estimate log D_o/w at pH = 7.4 (blood physiological pH) of a set of 30 additional compounds (including partially and fully ionized acids). The results at this pH follow the same trend observed in the 6 model acids, and validate the application of the method for D_o/w determination, except when α is very close to 1.

Recommendations for Improving Methods and Models for Aquatic Hazard Assessment of Ionizable Organic Chemicals

Ionizable organic chemicals (IOCs) such as organic acids and bases are an important substance class requiring aquatic hazard evaluation. Although the aquatic toxicity of IOCs is highly dependent on the water pH, many toxicity studies in the literature cannot be interpreted because pH was not reported or not kept constant during the experiment, calling for an adaptation and improvement of testing guidelines. The modulating influence of pH on toxicity is mainly caused by pH-dependent uptake and bioaccumulation of IOCs, which can be described by ion-trapping and toxicokinetic models. The internal effect concentrations of IOCs were found to be independent of the external pH because of organisms' and cells' ability to maintain a stable internal pH milieu. If the external pH is close to the internal pH, existing quantitative structure-activity relationships (QSARs) for neutral organics can be adapted by substituting the octanol-water partition coefficient by the ionization-corrected liposome-water distribution ratio as the hydrophobicity descriptor, demonstrated by modification of the target lipid model. Charged, zwitterionic and neutral species of an IOC can all contribute to observed toxicity, either through concentration-additive mixture effects or by interaction of different species, as is the case for uncoupling of mitochondrial respiration. For specifically acting IOCs, we recommend a 2-step screening procedure with ion-trapping/QSAR models used to predict the baseline toxicity, followed by adjustment using the toxic ratio derived from in vitro systems. Receptor- or plasma-binding models also show promise for elucidating IOC toxicity. The present review is intended to help demystify the ecotoxicity of IOCs and provide recommendations for their hazard and risk assessment. Environ Toxicol Chem 2020;39:269-286. © 2019 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals, Inc. on behalf of SETAC.

Concepts of Artificial Intelligence for Computer-Assisted Drug Discovery

Artificial intelligence (AI), and, in particular, deep learning as a subcategory of AI, provides opportunities for the discovery and development of innovative drugs. Various machine learning approaches have recently (re)emerged, some of which may be considered instances of domain-specific AI which have been successfully employed for drug discovery and design. This review provides a comprehensive portrayal of these machine learning techniques and of their applications in medicinal chemistry. After introducing the basic principles, alongside some application notes, of the various machine learning algorithms, the current state-of-the art of AI-assisted pharmaceutical discovery is discussed, including applications in structure- and ligand-based virtual screening, de novo drug design, physicochemical and pharmacokinetic property prediction, drug repurposing, and related aspects. Finally, several challenges and limitations of the current methods are summarized, with a view to potential future directions for AI-assisted drug discovery and design.

Measuring log Kow coefficients of neutral species of perfluoroalkyl carboxylic acids using reversed-phase high-performance liquid chromatography

Accurate measuring n-octanol/water partition coefficients (log K_ow) of perfluoroalkyl carboxylic acids (PFCAs) using experimental approach has been proven to be very difficult due to their special properties. The ionizable carboxyl groups in PFCAs make their log K_ow dependent on pH. In this study, the log K_ow values of neutral species of PFCAs (C_4≤n≤14) were measured based on reversed-phase high-performance liquid chromatography (RP-HPLC) with the mobile phase pH varying in the range of 1.09-5.00. The relationship between log K_ow and retention times was established using some reference compounds (including agrochemicals, polycyclic aromatic hydrocarbons) with known log K_ow values, and then validated with alkyl fatty acids, which have similar chemical structures as PFCAs. The apparent log K_ow (i.e., log D_ow) of the C_4-14 PFCAs were calculated based on their retention times using the established model, and they displayed a negative linear relationship with the mobile phase pH in the range of 1.09-4.00. Consequently, the log D_ow values were converted to the corresponding log K_ow values (1.05-7.19) based on the relationship of log D_ow = log K_ow + pK_a - pH. The log K_ow increased with perfluorinated carbon chain length with a greater rate for C₄ to C₅ PFCAs than for C_5-14 PFCAs.

Benchmark of different charges for prediction of the partitioning coefficient through the hydrophilic/lipophilic index

A few different theoretical methods for assigning the partial atomic charges were benchmarked for calculation of the hydrophilic/lipophilic index (HLI). The coefficients were selected to produce the best correlation of the HLI values with the experimental octanol-water partition. Different parameters were checked in calculations of partial charges to get the best performance of the HLI values obtained. Thus, four partitioning schemes (Coulson, Mulliken, Merz-Kollman, Ford-Wang) were benchmarked for calculations of atomic charges with six semiempirical methods (AM1, PM3, RM1, PM6, PM6-D3H4, PM7). Moreover, five distinct types of partial atomic charges (Mulliken, Hirshfeld, Löwdin, CHELPG, NPA), obtained at the Hartree-Fock and DFT levels of theory with three basis sets, were tested for their ability to produce the HLI values with the best correlation to experimental logP coefficients of 50 mono-charged organic anions. In the case of the semiempirical methods, the best correlation between the HLI and logP values (the correlation coefficient r = 0.9216) was obtained with the AM1 Ford-Wang parametric electrostatic potential charges. The Mulliken and Coulson charges calculated with the PM7 method can be used as an alternative to AM1, with the r values of 0.9107 and 0.8984, respectively. In the case of the DFT, the PBE/def2-TZVP natural population analysis charges produce the best correlation (r = 0.9220). Nevertheless, in spite of a marginally lower performance (r = 0.9159), the NPA charges computed at the PBE/def2-SVP level are more robust and can be regarded as the optimum choice for calculating the HLI values. Graphical abstract The hydrophilic/lipophilic index (HLI).

Prediction of partition and distribution coefficients in various solvent pairs with COSMO-RS

Performance of COSMO-RS method as a tool for partition and distribution modeling in 20 solvent pairs-composed of neutral or acidic aqueous solution and organic solvents of different polarity, ranging from alcohols to toluene and hexane-was evaluated. Experimental partition/distribution data of lignin-related and drug-like compounds (neutral, acidic, moderately basic) were used as reference. Several aspects of partition modeling were addressed: accounting for mutual saturation of aqueous and organic phases, variability of systematic prediction errors across solvent pairs, taking solute ionization into account. COSMO-RS was found to predict extraction outcome for both ligneous and drug-like compounds in various solvent pairs fairly well without any additional empirical input. The solvent-specific systematic errors were found to be moderate, despite being statistically significant, and related to the solvent hydrophobicity. Accounting for mutual solubilities of the two liquids was proven crucial in cases where water was considerably soluble in the organic solvent. The root mean square error of a priori logP prediction varied, depending mainly on the solvent pair, from 0.2 to 0.7, overall value being 0.6 log units. The accuracy was higher in case of hydrophilic than hydrophobic solvents. The logD predictions were less accurate, due to pK_a prediction being an additional source of error, and also because of the complexity of modeling the behaviour of ionic species in the two-phase system. A simple correction for partitioning of free ions was found to notably improve logD prediction accuracy in case of the most hydrophilic organic phase (butanol/water).

Lipophilicity Studies on Thiosemicarbazide Derivatives

The lipophilicity of two series of thiosemicarbazide derivatives was assessed by the RP-HPLC method with the RP-18 chromatographic column and the methanol–water mixture as the mobile phase. Distribution coefficients logP_HPLC were compared to calculated values generated by commonly used AClogP software and quantum chemical calculations. The reliability of the predictions was evaluated using the correlation matrix and PCA. For 4-benzoylthiosemicarbazides, a high correlation between theoretical and experimental logP parameters was obtained using the XlogP3 algorithm, while for 4-aryl/(cyclohexyl)thiosemicarbazides, the XlogP2 parameter was strongly correlated with the experimentally obtained logP.

Explicit consideration of spatial hydrogen bonding direction for activity coefficient prediction based on implicit solvation calculations

Directional hydrogen bonding is introduced to implicit solvation calculations for improved prediction of solvation properties and phase equilibria of associating fluids.