Development of quantitative structure-property relationship models for early ADME evaluation in drug discovery. 2. Blood-brain barrier penetration

Astragaloside IV as a Memory-Enhancing Agent: In Silico Studies with In Vivo Analysis and Post Mortem ADME-Tox Profiling in Mice

Many people around the world suffer from neurodegenerative diseases associated with cognitive impairment. As life expectancy increases, this number is steadily rising. Therefore, it is extremely important to search for new treatment strategies and to discover new substances with potential neuroprotective and/or cognition-enhancing effects. This study focuses on investigating the potential of astragaloside IV (AIV), a triterpenoid saponin with proven acetylcholinesterase (AChE)-inhibiting activity naturally occurring in the root of Astragalus mongholicus, to attenuate memory impairment. Scopolamine (SCOP), an antagonist of muscarinic cholinergic receptors, and lipopolysaccharide (LPS), a trigger of neuroinflammation, were used to impair memory processes in the passive avoidance (PA) test in mice. This memory impairment in SCOP-treated mice was attenuated by prior intraperitoneal (ip) administration of AIV at a dose of 25 mg/kg. The attenuation of memory impairment by LPS was not observed. It can therefore be assumed that AIV does not reverse memory impairment by anti-inflammatory mechanisms, although this needs to be further verified. All doses of AIV tested did not affect baseline locomotor activity in mice. In the post mortem analysis by mass spectrometry of the body tissue of the mice, the highest content of AIV was found in the kidneys, then in the spleen and liver, and the lowest in the brain.

Least Absolute Shrinkage and Selection Operator-Based Prediction of the Binding Constant of p-Sulfonatocalix[6]/[8]arenes with Alkaloids

p-Sulfonatocalix[n]arenes (SCnA) have demonstrated great potential for drug encapsulation through host-guest complexation to improve solubility, stability, and bioavailability. In this study, the solubilization effect of SCnA (n = 4, 6, 8) on 95 active compounds derived from traditional Chinese medicine (TCM) was investigated. Based on the significant solubilization effect on alkaloids, SC6A/SC8A and 76 alkaloids were selected as the host and guest, respectively, to determine the binding constant by competitive fluorescence titration. LASSO regression was adopted to investigate the mechanism of the complex of SCnA with alkaloids. The binding constant of alkaloids-SC6A and alkaloids-SC8A was related to the alkaloid alkalinity. Also, the electronegativity, polarization, first ionization potential, hydrogen bond potential, the molecular size, and shape of alkaloids are critical properties to determine alkaloids-SC6A binding constant as well as electronegativity, polarization, hydrophobicity, and the molecular size and shape of alkaloids play an important role for the alkaloids-SC8A binding constant.

Machine Learning-Assisted Computational Screening of Adhesive Molecules Derived from Dihydroxyphenyl Alanine

Marine mussels adhere to virtually any surface via 3,4-dihydroxyphenyl-L-alanines (L-DOPA), an amino acid largely contained in their foot proteins. The biofriendly, water-repellent, and strong adhesion of L-DOPA are unparalleled by any synthetic adhesive. Inspired by this, we computationally designed diverse derivatives of DOPA and studied their potential as adhesives or coating materials. We used first-principles calculations to investigate the adsorption of the DOPA derivatives on graphite. The presence of an electron-withdrawing group, such as nitrogen dioxide, strengthens the adsorption by increasing the π-π interaction between DOPA and graphite. To quantify the distribution of electron charge and to gain insights into the charge distribution at interfaces, we performed Bader charge analysis and examined charge density difference plots. We developed a quantitative structure-property relationship (QSPR) model using an artificial neural network (ANN) to predict the adsorption energy. Using the three-dimensional and quantum mechanical electrostatic potential of a molecule as a descriptor, the present quantum NN model shows promising performance as a predictive QSPR model.

Modeling the Blood-Brain Barrier Permeability of Potential Heterocyclic Drugs via Biomimetic IAM Chromatography Technique Combined with QSAR Methodology

Penetration through the blood-brain barrier (BBB) is desirable in the case of potential pharmaceuticals acting on the central nervous system (CNS), but is undesirable in the case of drug candidates acting on the peripheral nervous system because it may cause CNS side effects. Therefore, modeling of the permeability across the blood-brain barrier (i.e., the logarithm of the brain to blood concentration ratio, log BB) of potential pharmaceuticals should be performed as early as possible in the preclinical phase of drug development. Biomimetic chromatography with immobilized artificial membrane (IAM) and the quantitative structure-activity relationship (QSAR) methodology were successful in modeling the blood-brain barrier permeability of 126 drug candidates, whose experimentally-derived lipophilicity indices and computationally-derived molecular descriptors (such as molecular weight (MW), number of rotatable bonds (NRB), number of hydrogen bond donors (HBD), number of hydrogen bond acceptors (HBA), topological polar surface area (TPSA), and polarizability (α)) varied by class. The QSARs model established by multiple linear regression showed a positive effect of the lipophilicity (log kw, IAM) and molecular weight of the compound, and a negative effect of the number of hydrogen bond donors and acceptors, on the log BB values. The model has been cross-validated, and all statistics indicate that it is very good and has high predictive ability. The simplicity of the developed model, and its usefulness in screening studies of novel drug candidates that are able to cross the BBB by passive diffusion, are emphasized.

A machine learning-based quantitative model (LogBB_Pred) to predict the blood-brain barrier permeability (logBB value) of drug compounds

MOTIVATION

Efficient assessment of the blood-brain barrier (BBB) penetration ability of a drug compound is one of the major hurdles in central nervous system drug discovery since experimental methods are costly and time-consuming. To advance and elevate the success rate of neurotherapeutic drug discovery, it is essential to develop an accurate computational quantitative model to determine the absolute logBB value (a logarithmic ratio of the concentration of a drug in the brain to its concentration in the blood) of a drug candidate.

RESULTS

Here, we developed a quantitative model (LogBB_Pred) capable of predicting a logBB value of a query compound. The model achieved an R2 of 0.61 on an independent test dataset and outperformed other publicly available quantitative models. When compared with the available qualitative (classification) models that only classified whether a compound is BBB-permeable or not, our model achieved the same accuracy (0.85) with the best qualitative model and far-outperformed other qualitative models (accuracies between 0.64 and 0.70). For further evaluation, our model, quantitative models, and the qualitative models were evaluated on a real-world central nervous system drug screening library. Our model showed an accuracy of 0.97 while the other models showed an accuracy in the range of 0.29-0.83. Consequently, our model can accurately classify BBB-permeable compounds as well as predict the absolute logBB values of drug candidates.

AVAILABILITY AND IMPLEMENTATION

Web server is freely available on the web at http://ssbio.cau.ac.kr/software/logbb_pred/. The data used in this study are available to download at http://ssbio.cau.ac.kr/software/logbb_pred/dataset.zip.

Neuroprotective Properties of Oleanolic Acid-Computational-Driven Molecular Research Combined with In Vitro and In Vivo Experiments

Oleanolic acid (OA), as a ubiquitous compound in the plant kingdom, is studied for both its neuroprotective and neurotoxic properties. The mechanism of acetylcholinesterase (AChE) inhibitory potential of OA is investigated using molecular dynamic simulations (MD) and docking as well as biomimetic tests. Moreover, the in vitro SH-SY5Y human neuroblastoma cells and the in vivo zebrafish model were used. The inhibitory potential towards the AChE enzyme is examined using the TLC-bioautography assay (the IC50 value is 9.22 μM). The CH-π interactions between the central fragment of the ligand molecule and the aromatic cluster created by the His440, Phe288, Phe290, Phe330, Phe331, Tyr121, Tyr334, Trp84, and Trp279 side chains are observed. The results of the in vitro tests using the SH-SY5Y cells indicate that the viability rate is reduced to 71.5%, 61%, and 43% at the concentrations of 100 µg/mL, 300 µg/mL, and 1000 µg/mL, respectively, after 48 h of incubation, whereas cytotoxicity against the tested cell line with the IC50 value is 714.32 ± 32.40 µg/mL. The in vivo tests on the zebrafish prove that there is no difference between the control and experimental groups regarding the mortality rate and morphology (p > 0.05).

Predictions from First-Principles of Membrane Permeability to Small Molecules: How Useful Are They in Practice?

Predicting from first-principles the rate of passive permeation of small molecules across the biological membrane represents a promising strategy for screening lead compounds upstream in the drug-discovery and development pipeline. One popular avenue for the estimation of permeation rates rests on computer simulations in conjunction with the inhomogeneous solubility-diffusion model, which requires the determination of the free-energy change and position-dependent diffusivity of the substrate along the translocation pathway through the lipid bilayer. In this Perspective, we will clarify the physical meaning of the membrane permeability inferred from such computer simulations, and how theoretical predictions actually relate to what is commonly measured experimentally. We will also examine why these calculations remain both technically challenging and overly computationally expensive, which has hitherto precluded their routine use in nonacademic settings. We finally synopsize possible research directions to meet these challenges, increase the predictive power of physics-based rates of passive permeation, and, by ricochet, improve their practical usefulness.

Parsing structural fragments of thiazolidin-4-one based α-amylase inhibitors: A combined approach employing in vitro colorimetric screening and GA-MLR based QSAR modelling supported by molecular docking, molecular dynamics simulation and ADMET studies

α-Amylase (EC.3.2.1.1) is a ubiquitous digestive endoamylase. The abrupt rise in blood glucose levels due to the hydrolysis of carbohydrates by α-amylase at a faster rate is one of the main reasons for type 2 diabetes. The inhibitors prevent the action of digestive enzymes, slowing the digestion of carbs and eventually assisting in the management of postprandial hyperglycemia. In the course of developing α-amylase inhibitors, we have screened 2-aryliminothiazolidin-4-one based analogs for their in vitro α-amylase inhibitory potential and employed various in silico approaches for the detailed exploration of the bioactivity. The DNSA bioassay revealed that compounds 5c, 5e, 5h, 5j, 5m, 5o and 5t were more potent than the reference drug (IC₆₀ value = 22.94 ± 0.24 μg mL^-1). The derivative 5o with -NO₂ group at both the rings was the most potent analog with an IC₆₀ value of 19.67 ± 0.20 μg mL^-1 whereas derivative 5a with unsubstituted aromatic rings showed poor inhibitory potential with an IC₆₀ value of 33.40 ± 0.15 μg mL^-1. The reliable QSAR models were developed using the QSARINS software. The high value of R²_ext = 0.9632 for model IM-9 showed that the built model can be applied to predict the α-amylase inhibitory activity of the untested molecules. A consensus modelling approach was also employed to test the reliability and robustness of the developed QSAR models. Molecular docking and molecular dynamics were employed to validate the bioassay results by studying the conformational changes and interaction mechanisms. A step further, these compounds also exhibited good ADMET characteristics and bioavailability when tested for in silico pharmacokinetics prediction parameters.

New Ni(II) and Pd(II) complexes bearing derived sulfa drug ligands: synthesis, characterization, DFT calculations, and in silico and in vitro biological activity studies

In the present study, the synthesis of six new Ni(II) and Pd(II) complexes with three derived sulfamethoxazole drug ligands is reported. The coordination mode, geometry, and chemical formula of all the synthesized compounds have been determined by elemental analysis, mass spectrometry, emission atomic spectroscopy, conductivity measurements, magnetic susceptibility, FTIR, TGA, ¹H-NMR, electronic absorption spectroscopy, SEM-EDX along with DFT calculations. The Schiff Base ligands were found to be bidentate and coordinated to the metal ions through sulfonamidic nitrogen and oxazolic nitrogen atoms leading to a square planar geometry for palladium (II) while a distorted octahedral geometry around Nickel (II) ion was suggested. Biological applications of the new complexes including in vitro antimicrobial, antioxidant and anticancer properties were investigated. The results showed that the new metal (II) compounds exhibit remarkable antibacterial inhibition activity against both Gram-positive and Gram-negative bacteria, in addition to noticeable DPPH free radical scavenging activity. The in vitro cytotoxicity assay of the complexes against cell lines of chronic myelogenous leukaemia (K562) showed promising potential for the application of the coordination compounds in antitumor therapy. Subsequently, to evaluate the pharmaceutical potential of the metal-containing compounds, pharmacokinetics and toxicity were studied by ADMET simulations while interactions between the complexes and bacterial proteins were evaluated by molecular docking.

Combined Micellar Liquid Chromatography Technique and QSARs Modeling in Predicting the Blood-Brain Barrier Permeation of Heterocyclic Drug-like Compounds

The quantitative structure-activity relationship (QSAR) methodology was used to predict the blood-brain permeability (log BB) for 65 synthetic heterocyclic compounds tested as promising drug candidates. The compounds were characterized by different descriptors: lipophilicity, parachor, polarizability, molecular weight, number of hydrogen bond acceptors, number of rotatable bonds, and polar surface area. Lipophilic properties of the compounds were evaluated experimentally by micellar liquid chromatography (MLC). In the experiments, sodium dodecyl sulfate (SDS) as the effluent component and the ODS-2 column were used. Using multiple linear regression and leave-one-out cross-validation, we derived the statistically significant and highly predictive quantitative structure-activity relationship models. Thus, this study provides valuable information on the expected properties of the substances that can be used as a support tool in the design of new therapeutic agents.

Quantum Artificial Neural Network Approach to Derive a Highly Predictive 3D-QSAR Model for Blood-Brain Barrier Passage

A successful passage of the blood–brain barrier (BBB) is an essential prerequisite for the drug molecules designed to act on the central nervous system. The logarithm of blood–brain partitioning (LogBB) has served as an effective index of molecular BBB permeability. Using the three-dimensional (3D) distribution of the molecular electrostatic potential (ESP) as the numerical descriptor, a quantitative structure-activity relationship (QSAR) model termed AlphaQ was derived to predict the molecular LogBB values. To obtain the optimal atomic coordinates of the molecules under investigation, the pairwise 3D structural alignments were conducted in such a way to maximize the quantum mechanical cross correlation between the template and a target molecule. This alignment method has the advantage over the conventional atom-by-atom matching protocol in that the structurally diverse molecules can be analyzed as rigorously as the chemical derivatives with the same scaffold. The inaccuracy problem in the 3D structural alignment was alleviated in a large part by categorizing the molecules into the eight subsets according to the molecular weight. By applying the artificial neural network algorithm to associate the fully quantum mechanical ESP descriptors with the extensive experimental LogBB data, a highly predictive 3D-QSAR model was derived for each molecular subset with a squared correlation coefficient larger than 0.8. Due to the simplicity in model building and the high predictability, AlphaQ is anticipated to serve as an effective computational screening tool for molecular BBB permeability.

Pharmacophore model, docking, QSAR, and molecular dynamics simulation studies of substituted cyclic imides and herbal medicines as COX-2 inhibitors

Cyclooxygenase-2 (COX-2) enzyme inhibitors have not eliminated the necessity for developed drugs not only in the nonsteroidal anti-inflammatory drug (NSAIDs) area, but also in other therapeutic applications including prevention of cancer and Alzheimer's disease. A series of novel substituted cyclic imides have been reported as selective COX-2 inhibitors. To understand the structural features responsible for their activity, a 3D validated pharmacophore and quantitative structure−activity relationship (QSAR) model have been developed. The values of enrichment factor (EF), goodness of hit score (GH), area under the ROC curve (AUC), sensitivity, and specificity refer to the good ability of the pharmacophore model to identify active compounds. Multiple linear regression (MLR) produced statistically significant QSAR model with (R²_training = 0.763, R²_test = 0.96) and predictability (Q²_training = 0.66, Q²_test = 0.84). Then, using the pharmacophore and QSAR models, eight authenticated botanicals in two herbal medicines and the ZINC compounds database, were virtually screened for ligands to COX-2. The retrieved hits which also obey lipinski's rule of five (RO5) were docked in the COX-2 3D structure to investigate their binding mode and affinity. Finally, based on the docking results, nine molecules were prioritized as promising hits that could be used as leads to discover novel COX-2 inhibitors. COX-2 inhibition of most of these hits has not been reported previously. Ten-nanosecond molecular dynamics simulation (10-ns MD) was performed on the initial structure COX-2 complex with ZINC000113253375 and ZINC000043170560 resulted from the docking. Our utilization of the 3D pharmacophore model, QSAR, molecular docking, and molecular dynamics simulation trials can be a potent strategy to successfully predict activity, efficiently design drugs, and screen large numbers of new compounds as active drug candidates.

Drug properties and host factors contribute to biochemical presentation of drug-induced liver injury: a prediction model from a machine learning approach

Drug-induced liver injury (DILI) presentation varies biochemically and histologically. Certain drugs present quite consistent injury patterns, i.e., DILI signatures. In contrast, others are manifested as broader types of liver injury. The variety of DILI presentations by a single drug suggests that both drugs and host factors may contribute to the phenotype. However, factors determining the DILI types have not been yet elucidated. Identifying such factors may help to accurately predict the injury types based on drugs and host information and assist the clinical diagnosis of DILI. Using prospective DILI registry datasets, we sought to explore and validate the associations of biochemical injury types at the time of DILI recognition with comprehensive information on drug properties and host factors. Random forest models identified a set of drug properties and host factors that differentiate hepatocellular from cholestatic damage with reasonable accuracy (69-84%). A simplified logistic regression model developed for practical use, consisting of patient's age, drug's lipoaffinity, and hybridization ratio, achieved a fair prediction (68-74%), but suggested potential clinical usability, computing the likelihood of liver injury type based on two properties of drugs taken by a patient and patient's age. In summary, considering both drug and host factors in evaluating DILI risk and phenotypes open an avenue for future DILI research and aid in the refinement of causality assessment.

Prediction of the Partition Coefficient between Adipose Tissue and Blood for Environmental Chemicals: From Single QSAR Models to an Integrated Approach

The adipose tissue:blood partition coefficient is a key-endpoint to predict the pharmacokinetics of chemicals in humans and animals, since other organ:blood affinities can be estimated as a function of this parameter. We performed a search in the literature to select all the available rat in vivo data. This approach resulted into two improvements to existing models: a homogeneous definition of the endpoint and an expanded data collection. The resulting dataset was used to develop QSAR models as a function of linear and non-linear algorithms. Several applicability domain definitions were assessed and the definition corresponding to a good balance between performance and coverage was retained. We assessed the pertinence of combining single models into integrated approaches to increase the accuracy in predictions. The best integrated model outperformed the single models and it was characterized by an external mean absolute error (MAE) equal to 0.26, while preserving an adequate coverage (84 %). This performance is comparable to experimental variability and it highlights the pertinence of the integrated model.

Nonlinear Method to Predict the Distribution of Structurally Diverse Compounds between Blood and Tissue

In this work, methods for predicting the distribution of compounds between blood and tissue were investigated using nonlinear regression analysis. For the tissue/blood partition coefficient, 282 compounds and 810 activity data for seven tissues were selected. Twenty-four parameters were studied for each state of the compound. A study set with the most compounds and activity data in similar studies was established, and a model with good prediction ability was obtained. A total of 773 data points were randomly divided into two data sets: the training set contained 623 data points (n = 623, r = 0.822, s = 0.438, F = 142.2, and Q = 0.814) and the test set contained 150 data points (n = 150, r = 0.814, and s = 0.334). Furthermore, individual tissue/blood distribution coefficients were also studied in depth to obtain separate models for predicting the distribution coefficient of a drug between blood and a tissue. By applying these models, not only can the tissue/blood distribution coefficient or single tissue/blood distribution coefficient of the seven tissues and organs of the human body be predicted, but the distribution of the drug molecules in different states (neutral, cation, and anion) in the three tissue components can also be predicted.

Antiproliferative activity of cardenolides on cell line A549: structure-activity relationship analysis

Since the beginning, natural products have represented an important source of bioactive molecules for cancer treatment. Among them, cardenolides attract the attention of different research groups due to their cardiotonic and antitumor activity. The observed biological activity is closely related to their Na⁺/K⁺-ATPase inhibition potency. Currently, the discovery of new compounds against cancer is an urgent need in modern pharmaceutical research. Thus, the aim of this work is to determine the physicochemical properties and substituent effects that module the antiproliferative activity of cardenolides on the human lung cancer cell line A549. We build and curate a library with results obtained from literature; molecular descriptors were calculated in PaDEL software, and SAR/QSAR analysis was performed. The SAR results showed that cardenolides were sensitive to modifications in C and D steroidal ring and required substituent groups with the function of hydrogen bond acceptor at the C3 position. QSAR models to doubly linked-type cardenolides indicated that properties as lipoaffinity and atoms with the capacity to be hydrogen bond acceptors are involved in the increment of antiproliferative activity on A549 cell line. In contrast, the presence and position of very electro-negative atoms on the molecule decreased the antiproliferative effect on A549 cells. These results suggest that the antiproliferative capacity of cardenolides on the cell line A549 is strongly related to substituent groups on the C3 position, which must not be carbohydrate. Additionally, the steroidal rings C and D must remain without modifications.

Novel and Predictive QSAR Model for Steroidal and Nonsteroidal 5α- Reductase Type II Inhibitors

AIMS AND OBJECTIVE

In this study, a novel quantitative structure activity relationship (QSAR) model has been developed for inhibitors of human 5-alpha reductase type II, which are used to treat benign prostate hypertrophy (BPH).

METHODS

The dataset consisted of 113 compounds-mainly nonsteroidal-with known inhibitory concentration. Then 3D structures of compounds were optimized and molecular structure descriptors were calculated. The stepwise multiple linear regression was used to select descriptors encoding the inhibitory activity of the compounds. Multiple linear regression (MLR) was used to build up the linear QSAR model.

RESULTS

The results obtained revealed that the descriptors which best describe the activity were atom type electropological state, carbon type, radial distribution function (RDF), barysz matrix and molecular linear free energy relation. The suggested model could achieve satisfied square correlation coefficient of R² = 0.72, higher than of many previous studies, indicating its superiority. Rigid validation criteria were met using external data with Q² ˃ 0.5 and R² = 0.75, reflecting the predictive power of the model.

CONCLUSION

The QSAR model was applied for screening botanical components of herbal preparations used to treat BPH, and could predict the activity of some, among others, making reasonable attribution to the proposed effect of these preparations. Gamma tocopherol was found to be an active inhibitor, in consistence with many previous studies, anticipating the power of this model in the prediction of new candidate molecules and suggesting further investigations.

Tuning Properties for Blood-Brain Barrier Permeation: A Statistics-Based Analysis

In the effort to define a set of rules useful in tuning the properties for a successful blood-brain barrier (BBB) permeation, we statistically analyzed a set of 328 compounds and correlated their experimental in vivo logBB with a series of computed descriptors. Contingency tables were constructed, observed and expected distributions were calculated, and chi-square (χ²) distributions were evaluated. This allowed to point out a significant dependence of certain physicochemical properties in influencing the BBB permeation. Of over 15 computed descriptors, 9 resulted to be particularly important showing highly significant χ² distribution: polar surface area (χ² = 66.79; p = 1.08 × 10^-13), nitrogen and oxygen count (χ² = 51.17; p = 2.06 × 10^-10), logP (χ² = 47.38; p = 1.27 × 10^-9), nitrogen count (χ² = 38.29; p = 9.77 × 10^-8), logD (χ² = 36.80; p = 36.80), oxygen count (χ² = 35.83; p = 3.13 × 10^-7), ionization state (χ² = 33.02, p = 3.19 × 10^-7), hydrogen bond acceptors (χ² = 30.80; p = 3.36 × 10^-6), and hydrogen bond donors (χ² = 29.29; p = 6.81 × 10^-6). Other parameters describing the mass and size of the molecules (molecular weight: 11.18; p = 2.46 × 10^-2) resulted in being not significant since the population within the observed and expected distribution was similar. Depending on the combination of the significant descriptors, we set a three cases probabilistic scenario (BBB⁺, BBB^-, BBB⁺/BBB^-) that would prospectively be used to tune properties for BBB permeation.

The use of fast molecular descriptors and artificial neural networks approach in organochlorine compounds electron ionization mass spectra classification

Developing of theoretical tools can be very helpful for supporting new pollutant detection. Nowadays, a combination of mass spectrometry and chromatographic techniques are the most basic environmental monitoring methods. In this paper, two organochlorine compound mass spectra classification systems were proposed. The classification models were developed within the framework of artificial neural networks (ANNs) and fast 1D and 2D molecular descriptor calculations. Based on the intensities of two characteristic MS peaks, namely, [M] and [M-35], two classification criterions were proposed. According to criterion I, class 1 comprises [M] signals with the intensity higher than 800 NIST units, while class 2 consists of signals with the intensity lower or equal than 800. According to criterion II, class 1 consists of [M-35] signals with the intensity higher than 100, while signals with the intensity lower or equal than 100 belong to class 2. As a result of ANNs learning stage, five models for both classification criterions were generated. The external model validation showed that all ANNs are characterized by high predicting power; however, criterion I-based ANNs are much more accurate and therefore are more suitable for analytical purposes. In order to obtain another confirmation, selected ANNs were tested against additional dataset comprising popular sunscreen agents disinfection by-products reported in previous works.

The influence of drug properties and host factors on delayed onset of symptoms in drug-induced liver injury

BACKGROUND & AIMS

Most patients with drug-induced liver injury (DILI) manifest clinical symptoms while on therapy, while some patients manifest days or weeks after drug cessation (delayed onset). This challenges DILI causality assessment and diagnosis. Factors contributing to the delayed onset phenotype are currently unknown. We explored factors contributing to delayed onset of DILI by analysing culprit drug properties, host factors and their interactions in a large patient population from the Spanish DILI Registry.

METHODS

Clinical information from 388 patients (69 presented delayed onset) and drug properties of 43 causative drugs (45 active ingredients) were analysed. A two-tier regression-based model was used to assess host/drug interactions affecting the probability of delayed onset.

RESULTS

Antibacterial and anti-inflammatory drugs accounted for the delayed onset cases. Drug property of <50% hepatic metabolism (odds ratio [OR] 11.06, 95% confidence interval [95% CI]: 4.4-32.2, P = 0.0003), daily dose ≥1000 mg (OR: 2.77, 95% CI: 1.3-6.1, P = 0.0063) and the absence of pre-existing conditions in a patient (OR: 2.55, 95% CI: 1.3-4.9, P = 0.0043) were independently associated with delayed onset. The findings were consistent when externally validated using Latin American DILI Network cases (N = 131). Likewise, drug properties of mitochondrial liability and Pauling electronegativity were associated with delayed onset, but dependent on specific host factors such as age, sex and pre-existing cardiac diseases.

CONCLUSIONS

This study demonstrated that delayed onset, a specific DILI phenotype, is explained by complex interactions among drug properties and host factors and provided mechanistic hypotheses for future studies. These findings can help improve the diagnostic capability and causality assessment.

Application of Multivariate Adaptive Regression Splines (MARSplines) for Predicting Hansen Solubility Parameters Based on 1D and 2D Molecular Descriptors Computed from SMILES String

A new method of Hansen solubility parameters (HSPs) prediction was developed by combining the multivariate adaptive regression splines (MARSplines) methodology with a simple multivariable regression involving 1D and 2D PaDEL molecular descriptors. In order to adopt the MARSplines approach to QSPR/QSAR problems, several optimization procedures were proposed and tested. The effectiveness of the obtained models was checked via standard QSPR/QSAR internal validation procedures provided by the QSARINS software and by predicting the solubility classification of polymers and drug-like solid solutes in collections of solvents. By utilizing information derived only from SMILES strings, the obtained models allow for computing all of the three Hansen solubility parameters including dispersion, polarization, and hydrogen bonding. Although several descriptors are required for proper parameters estimation, the proposed procedure is simple and straightforward and does not require a molecular geometry optimization. The obtained HSP values are highly correlated with experimental data, and their application for solving solubility problems leads to essentially the same quality as for the original parameters. Based on provided models, it is possible to characterize any solvent and liquid solute for which HSP data are unavailable.

Discovery of Pyrazolo[4,3-c]quinolines Derivatives as Potential Anti-Inflammatory Agents through Inhibiting of NO Production

The synthesis and anti-inflammatory effects of certain pyrazolo[4,3-c]quinoline derivatives 2a–2r are described. The anti-inflammatory activities of these derivatives were evaluated by means of inhibiting nitric oxide (NO) production in lipopolysaccharide (LPS)-induced RAW 264.7 cells. Among them, 3-amino-4-(4-hydroxyphenylamino)-1H-pyrazolo[4,3-c]-quinoline (2i) and 4-(3-amino-1H-pyrazolo[4,3-c]quinolin-4-ylamino)benzoic acid (2m) exhibited significant inhibition of LPS-stimulated NO production with a potency approximately equal to that of the positive control, 1400 W. Important structure features were analyzed by quantitative structure–activity relationship (QSAR) analysis to give better insights into the structure determinants for predicting the inhibitory effects on the accumulation of nitric oxide for RAW 264.7 cells in response to LPS. In addition, our results indicated that their anti-inflammatory effects involve the inhibition of inducible nitric oxide synthase (iNOS) and cyclooxygenase 2 (COX-2) protein expression. Further studies on the structural optimization are ongoing.

Selection of effective cocrystals former for dissolution rate improvement of active pharmaceutical ingredients based on lipoaffinity index

New theoretical screening procedure was proposed for appropriate selection of potential cocrystal formers possessing the ability of enhancing dissolution rates of drugs. The procedure relies on the training set comprising 102 positive and 17 negative cases of cocrystals found in the literature. Despite the fact that the only available data were of qualitative character, performed statistical analysis using binary classification allowed to formulate quantitative criterions. Among considered 3679 molecular descriptors the relative value of lipoaffinity index, expressed as the difference between values calculated for active compound and excipient, has been found as the most appropriate measure suited for discrimination of positive and negative cases. Assuming 5% precision, the applied classification criterion led to inclusion of 70% positive cases in the final prediction. Since lipoaffinity index is a molecular descriptor computed using only 2D information about a chemical structure, its estimation is straightforward and computationally inexpensive. The inclusion of an additional criterion quantifying the cocrystallization probability leads to the following conjunction criterions H_mix<-0.18 and ΔLA>3.61, allowing for identification of dissolution rate enhancers. The screening procedure was applied for finding the most promising coformers of such drugs as Iloperidone, Ritonavir, Carbamazepine and Enthenzamide.

Discovery of Indeno[1,2-c]quinoline Derivatives as Potent Dual Antituberculosis and Anti-Inflammatory Agents

A series of indeno[1,2-c]quinoline derivatives were designed, synthesized and evaluated for their anti-tuberculosis (anti-TB) and anti-inflammatory activities. The minimum inhibitory concentration (MIC) of the newly synthesized compound was tested against Mycobacterium tuberculosis H₃₇R_V. Among the tested compounds, (E)-N′-[6-(4-hydroxypiperidin-1-yl)-11H-indeno[1,2-c]quinolin-11-ylidene]isonicotino-hydrazide (12), exhibited significant activities against the growth of M. tuberculosis (MIC values of 0.96 μg/mL) with a potency approximately equal to that of isoniazid (INH), an anti-TB drug. Important structure features were analyzed by quantitative structure–activity relationship (QSAR) analysis to give better insights into the structure determinants for predicting the anti-TB activity. The anti-inflammatory activity was induced by superoxide anion generation and neutrophil elastase (NE) release using the formyl-l-methionyl-l-leucyl-l-phenylalanine (fMLF)-activated human neutrophils method. Results indicated that compound 12 demonstrated a potent dual inhibitory effect on NE release and superoxide anion generation with IC₅₀ values of 1.76 and 1.72 μM, respectively. Our results indicated that compound 12 is a potential lead compound for the discovery of dual anti-TB and anti-inflammatory drug candidates. In addition, 6-[3-(hydroxymethyl)piperidin-1-yl]-9-methoxy-11H-indeno[1,2-c]quinolin-11-one (4g) showed a potent dual inhibitory effect on NE release and superoxide anion generation with IC₅₀ values of 0.46 and 0.68 μM, respectively, and is a potential lead compound for the discovery of anti-inflammatory drug candidates.

Predicting a Drug's Membrane Permeability: A Computational Model Validated With in Vitro Permeability Assay Data

Membrane permeability is a key property to consider during the drug design process, and particularly vital when dealing with small molecules that have intracellular targets as their efficacy highly depends on their ability to cross the membrane. In this work, we describe the use of umbrella sampling molecular dynamics (MD) computational modeling to comprehensively assess the passive permeability profile of a range of compounds through a lipid bilayer. The model was initially calibrated through in vitro validation studies employing a parallel artificial membrane permeability assay (PAMPA). The model was subsequently evaluated for its quantitative prediction of permeability profiles for a series of custom synthesized and closely related compounds. The results exhibited substantially improved agreement with the PAMPA data, relative to alternative existing methods. Our work introduces a computational model that underwent progressive molding and fine-tuning as a result of its synergistic collaboration with numerous in vitro PAMPA permeability assays. The presented computational model introduces itself as a useful, predictive tool for permeability prediction.

Comprehensive identification of sphingolipid species by in silico retention time and tandem mass spectral library

Background

Liquid chromatography coupled with electrospray ionization tandem mass spectrometry (LC–ESI–MS/MS) is used for comprehensive metabolome and lipidome analyses. Compound identification relies on similarity matching of the retention time (RT), precursor m/z, isotopic ratio, and MS/MS spectrum with reference compounds. For sphingolipids, however, little information on the RT and MS/MS references is available.

Results

Negative-ion ESI–MS/MS is a useful method for the structural characterization of sphingolipids. We created theoretical MS/MS spectra for 21 sphingolipid classes in human and mouse (109,448 molecules), with substructure-level annotation of unique fragment ions by MS-FINDER software. The existence of ceramides with β-hydroxy fatty acids was confirmed in mouse tissues based on cheminformatic- and quantum chemical evidences. The RT of sphingo- and glycerolipid species was also predicted for our LC condition. With this information, MS-DIAL software for untargeted metabolome profiling could identify 415 unique structures including 282 glycerolipids and 133 sphingolipids from human cells (HEK and HeLa) and mouse tissues (ear and liver).

Conclusions

MS-DIAL and MS-FINDER software programs can identify 42 lipid classes (21 sphingo- and 21 glycerolipids) with the in silico RT and MS/MS library. The library is freely available as Microsoft Excel files at the software section of our RIKEN PRIMe website (http://prime.psc.riken.jp/).

Electronic supplementary material

The online version of this article (doi:10.1186/s13321-017-0205-3) contains supplementary material, which is available to authorized users.

Molecular determinants of blood-brain barrier permeation

The blood-brain barrier (BBB) is a microvascular unit which selectively regulates the permeability of drugs to the brain. With the rise in CNS drug targets and diseases, there is a need to be able to accurately predict a priori which compounds in a company database should be pursued for favorable properties. In this review, we will explore the different computational tools available today, as well as underpin these to the experimental methods used to determine BBB permeability. These include in vitro models and the in vivo models that yield the dataset we use to generate predictive models. Understanding of how these models were experimentally derived determines our accurate and predicted use for determining a balance between activity and BBB distribution.

Estimating sensory irritation potency of volatile organic chemicals using QSARs based on decision tree methods for regulatory purpose

Volatile organic compounds (VOCs) are among the priority atmospheric pollutants that have high indoor and outdoor exposure potential. The toxicity assessment of VOCs to living ecosystems has received considerable attention in recent years. Development of computational methods for safety assessment of chemicals has been advocated by various regulatory agencies. The paper proposes robust and reliable quantitative structure-activity relationships (QSARs) for estimating the sensory irritation potency and screening of the VOCs. Here, decision tree (DT) based classification and regression QSARs models, such as single DT, decision tree forest (DTF), and decision tree boost (DTB) were developed using the sensory irritation data on VOCs in mice following the OECD principles. Structural diversity and nonlinearity in the data were evaluated through the Euclidean distance and Brock-Dechert-Scheinkman statistics. The constructed QSAR models were validated with external test data and the predictive performance of these models was established through a set of coefficients recommended in QSAR literature. The performance of all three classification and regression QSAR models was satisfactory, but DTF and DTB performed relatively better. The classification and regression QSAR models (DTF, DTB) rendered classification accuracies of 98.59 and 100 %, and yielded correlations (R(2)) of 0.950 and 0.971, respectively in complete data. The lipoaffinity index and SwHBa were identified as the most influential descriptors in proposed QSARs. The developed QSARs performed better than the previous studies. The developed models exhibited high statistical confidence and identified the structural properties of the VOCs responsible for their sensory irritation, and hence could be useful tools in screening of chemicals for regulatory purpose.

In silico prediction of chemical aquatic toxicity with chemical category approaches and substructural alerts

Aquatic toxicity is an important endpoint in the evaluation of chemically adverse effects on ecosystems.

Systems approaches and polypharmacology for drug discovery from herbal medicines: an example using licorice

ETHNOPHARMACOLOGICAL RELEVANCE

Licorice, one of the oldest and most popular herbal medicines in the world, has been widely used in traditional Chinese medicine as a cough reliever, anti-inflammatory, anti-anabrosis, immunomodulatory, anti-platelet, antiviral (hepatitis) and detoxifying agent. Licorice was used as an example to show drug discovery from herbal drugs using systems approaches and polypharmacology.

AIM OF THE STUDY

Herbal medicines are becoming more mainstream in clinical practice and show value in treating and preventing diseases. However, due to its extreme complexity both in chemical components and mechanisms of action, deep understanding of botanical drugs is still difficult. Thus, a comprehensive systems approach which could identify active ingredients and their targets in the crude drugs and more importantly, understand the biological basis for the pharmacological properties of herbal medicines is necessary.

MATERIALS AND METHODS

In this study, a novel systems pharmacology model that integrates oral bioavailability screening, drug-likeness evaluation, blood-brain barrier permeation, target identification and network analysis has been established to investigate the herbal medicines.

RESULTS

The comprehensive systems approach effectively identified 73 bioactive components from licorice and 91 potential targets for this medicinal herb. These 91 targets are closely associated with a series of diseases of respiratory system, cardiovascular system, and gastrointestinal system, etc. These targets are further mapped to drug-target and drug-target-disease networks to elucidate the mechanism of this herbal medicine.

CONCLUSION

This work provides a novel in silico strategy for investigation of the botanical drugs containing a huge number of components, which has been demonstrated by the well-studied licorice case. This attempt should be helpful for understanding definite mechanisms of action for herbal medicines and discovery of new drugs from plants.

Back to the future: can physical models of passive membrane permeability help reduce drug candidate attrition and move us beyond QSPR?

It is widely recognized that adsorption, distribution, metabolism, excretion, and toxicology liabilities kill the majority of drug candidates that progress to clinical trials. The development of computational models to predict small molecule membrane permeability is therefore of considerable scientific and public health interest. Empirical qualitative structure permeability relationship models of permeability have been a mainstay in industrial applications, but lack a deep understanding of the underlying biologic physics. Others and we have shown that implicit solvent models to predict passive permeability for small molecules exhibit mediocre predictive performance when validated across experimental test sets. Given the vast increase in computer power, more efficient parallelization schemes, and extension of current atomistic simulation codes to general use graphical processing units, the development and application of physical models based on all-atom simulations may now be feasible. Preliminary results from rigorous free energy calculations using all-atom simulations indicate that performance relative to implicit solvent models may be improved, but many outstanding questions remain. Here, we review the current state-of-the-art physical models for passive membrane permeability prediction and present a prospective look at promising new directions for all-atom approaches.

Testing physical models of passive membrane permeation

The biophysical basis of passive membrane permeability is well-understood, but most methods for predicting membrane permeability in the context of drug design are based on statistical relationships that indirectly capture the key physical aspects. Here, we investigate molecular mechanics-based models of passive membrane permeability and evaluate their performance against different types of experimental data, including parallel artificial membrane permeability assays (PAMPA), cell-based assays, in vivo measurements, and other in silico predictions. The experimental data sets we use in these tests are diverse, including peptidomimetics, congeneric series, and diverse FDA approved drugs. The physical models are not specifically trained for any of these data sets; rather, input parameters are based on standard molecular mechanics force fields, such as partial charges, and an implicit solvent model. A systematic approach is taken to analyze the contribution from each component in the physics-based permeability model. A primary factor in determining rates of passive membrane permeation is the conformation-dependent free energy of desolvating the molecule, and this measure alone provides good agreement with experimental permeability measurements in many cases. Other factors that improve agreement with experimental data include deionization and estimates of entropy losses of the ligand and the membrane, which lead to size-dependence of the permeation rate.

ADMET evaluation in drug discovery. 12. Development of binary classification models for prediction of hERG potassium channel blockage

Inhibition of the human ether-a-go-go related gene (hERG) potassium channel may result in QT interval prolongation, which causes severe cardiac side effects and is a major problem in clinical studies of drug candidates. The development of in silico tools to filter out potential hERG potassium channel blockers in early stages of the drug discovery process is of considerable interest. Here, a diverse set of 806 compounds with hERG inhibition data was assembled, and the binary hERG classification models using naive Bayesian classification and recursive partitioning (RP) techniques were established and evaluated. The naive Bayesian classifier based on molecular properties and the ECFP_8 fingerprints yielded 84.8% accuracy for the training set using the leave-one-out (LOO) cross-validation procedure and 85% accuracy for the test set of 120 molecules. For the two additional test sets, the model achieved 89.4% accuracy for the WOMBAT-PK test set, and 86.1% accuracy for the PubChem test set. The naive Bayesian classifiers gave better predictions than the RP classifiers. Moreover, the Bayesian classifier, employing molecular fingerprints, highlights the important structural fragments favorable or unfavorable for hERG potassium channel blockage, which offers extra valuable information for the design of compounds avoiding undesirable hERG activity.

Qualitative prediction of blood-brain barrier permeability on a large and refined dataset

The prediction of blood-brain barrier permeation is vitally important for the optimization of drugs targeting the central nervous system as well as for avoiding side effects of peripheral drugs. Following a previously proposed model on blood-brain barrier penetration, we calculated the cross-sectional area perpendicular to the amphiphilic axis. We obtained a high correlation between calculated and experimental cross-sectional area (r = 0.898, n = 32). Based on these results, we examined a correlation of the calculated cross-sectional area with blood-brain barrier penetration given by logBB values. We combined various literature data sets to form a large-scale logBB dataset with 362 experimental logBB values. Quantitative models were calculated using bootstrap validated multiple linear regression. Qualitative models were built by a bootstrapped random forest algorithm. Both methods found similar descriptors such as polar surface area, pKa, logP, charges and number of positive ionisable groups to be predictive for logBB. In contrast to our initial assumption, we were not able to obtain models with the cross-sectional area chosen as relevant parameter for both approaches. Comparing those two different techniques, qualitative random forest models are better suited for blood-brain barrier permeability prediction, especially when reducing the number of descriptors and using a large dataset. A random forest prediction system (n(trees) = 5) based on only four descriptors yields a validated accuracy of 88%.

Rank order entropy: why one metric is not enough

The use of Quantitative Structure-Activity Relationship models to address problems in drug discovery has a mixed history, generally resulting from the misapplication of QSAR models that were either poorly constructed or used outside of their domains of applicability. This situation has motivated the development of a variety of model performance metrics (r(2), PRESS r(2), F-tests, etc.) designed to increase user confidence in the validity of QSAR predictions. In a typical workflow scenario, QSAR models are created and validated on training sets of molecules using metrics such as Leave-One-Out or many-fold cross-validation methods that attempt to assess their internal consistency. However, few current validation methods are designed to directly address the stability of QSAR predictions in response to changes in the information content of the training set. Since the main purpose of QSAR is to quickly and accurately estimate a property of interest for an untested set of molecules, it makes sense to have a means at hand to correctly set user expectations of model performance. In fact, the numerical value of a molecular prediction is often less important to the end user than knowing the rank order of that set of molecules according to their predicted end point values. Consequently, a means for characterizing the stability of predicted rank order is an important component of predictive QSAR. Unfortunately, none of the many validation metrics currently available directly measure the stability of rank order prediction, making the development of an additional metric that can quantify model stability a high priority. To address this need, this work examines the stabilities of QSAR rank order models created from representative data sets, descriptor sets, and modeling methods that were then assessed using Kendall Tau as a rank order metric, upon which the Shannon entropy was evaluated as a means of quantifying rank-order stability. Random removal of data from the training set, also known as Data Truncation Analysis (DTA), was used as a means for systematically reducing the information content of each training set while examining both rank order performance and rank order stability in the face of training set data loss. The premise for DTA ROE model evaluation is that the response of a model to incremental loss of training information will be indicative of the quality and sufficiency of its training set, learning method, and descriptor types to cover a particular domain of applicability. This process is termed a "rank order entropy" evaluation or ROE. By analogy with information theory, an unstable rank order model displays a high level of implicit entropy, while a QSAR rank order model which remains nearly unchanged during training set reductions would show low entropy. In this work, the ROE metric was applied to 71 data sets of different sizes and was found to reveal more information about the behavior of the models than traditional metrics alone. Stable, or consistently performing models, did not necessarily predict rank order well. Models that performed well in rank order did not necessarily perform well in traditional metrics. In the end, it was shown that ROE metrics suggested that some QSAR models that are typically used should be discarded. ROE evaluation helps to discern which combinations of data set, descriptor set, and modeling methods lead to usable models in prioritization schemes and provides confidence in the use of a particular model within a specific domain of applicability.

Identification of novel functional inhibitors of acid sphingomyelinase

We describe a hitherto unknown feature for 27 small drug-like molecules, namely functional inhibition of acid sphingomyelinase (ASM). These entities named FIASMAs (Functional Inhibitors of Acid SphingoMyelinAse), therefore, can be potentially used to treat diseases associated with enhanced activity of ASM, such as Alzheimer's disease, major depression, radiation- and chemotherapy-induced apoptosis and endotoxic shock syndrome. Residual activity of ASM measured in the presence of 10 µM drug concentration shows a bimodal distribution; thus the tested drugs can be classified into two groups with lower and higher inhibitory activity. All FIASMAs share distinct physicochemical properties in showing lipophilic and weakly basic properties. Hierarchical clustering of Tanimoto coefficients revealed that FIASMAs occur among drugs of various chemical scaffolds. Moreover, FIASMAs more frequently violate Lipinski's Rule-of-Five than compounds without effect on ASM. Inhibition of ASM appears to be associated with good permeability across the blood-brain barrier. In the present investigation, we developed a novel structure-property-activity relationship by using a random forest-based binary classification learner. Virtual screening revealed that only six out of 768 (0.78%) compounds of natural products functionally inhibit ASM, whereas this inhibitory activity occurs in 135 out of 2028 (6.66%) drugs licensed for medical use in humans.

Prediction of blood-brain partitioning: a model based on molecular electronegativity distance vector descriptors

The objective of this paper is to build a reliable model based on the molecular electronegativity distance vector (MEDV) descriptors for predicting the blood-brain barrier (BBB) permeability and to reveal the effects of the molecular structural segments on the BBB permeability. Using 70 structurally diverse compounds, the partial least squares regression (PLSR) models between the BBB permeability and the MEDV descriptors were developed and validated by the variable selection and modeling based on prediction (VSMP) technique. The estimation ability, stability, and predictive power of a model are evaluated by the estimated correlation coefficient (r), leave-one-out (LOO) cross-validation correlation coefficient (q), and predictive correlation coefficient (R(p)). It has been found that PLSR model has good quality, r=0.9202, q=0.7956, and R(p)=0.6649 for M1 model based on the training set of 57 samples. To search the most important structural factors affecting the BBB permeability of compounds, we performed the values of the variable importance in projection (VIP) analysis for MEDV descriptors. It was found that some structural fragments in compounds, such as -CH(3), -CH(2)-, =CH-, =C, triple bond C-, -CH<, =C<, =N-, -NH-, =O, and -OH, are the most important factors affecting the BBB permeability.