Radial Distribution Function descriptors for predicting affinity for vitamin D receptor

Heterogeneous Photodegradation Behavior of Liquid Crystal Monomers in Dust: Quantitative Structure-Activity Relationship and Product Identification

The heterogeneous photodegradation behavior of liquid crystal monomers (LCMs) in standard dust (standard reference material, SRM 2583) and environmental dust was investigated. The measured photodegradation ratios for 23 LCMs in SRM and environmental dust in 12 h were 11.1 ± 1.8 to 23.2 ± 1.1% and 8.7 ± 0.5 to 24.0 ± 2.8%, respectively. The degradation behavior of different LCM compounds varied depending on their structural properties. A quantitative structure-activity relationship model for predicting the degradation ratio of LCMs in SRM dust was established, which revealed that the molecular descriptors related to molecular polarizability, electronegativity, and molecular mass were closely associated with LCMs' photodegradation. The photodegradation products of the LCM compound 4'-propoxy-4-biphenylcarbonitrile (PBIPHCN) in dust, including •OH oxidation, C-O bond cleavage, and ring-opening products, were identified by nontarget analysis, and the corresponding degradation pathways were suggested. Some of the identified products, such as 4'-hydroxyethoxy-4-biphenylcarbonitrile, showed predicted toxicity (with an oral rat lethal dose of 50%) comparable to that of PBIPHCN. The half-lives of the studied LCMs in SRM dust were estimated at 32.2-82.5 h by fitting an exponential decay curve to the observed photodegradation data. The photodegradation mechanisms of LCMs in dust were revealed for the first time, enhancing the understanding of LCMs' environmental behavior and risks.

Machine learning methods applied for the prediction of biological activities of triple reuptake inhibitors

Major depressive disorder (MDD) is characterized by a series of disabling symptoms like anhedonia, depressed mood, lack of motivation for daily tasks and self-extermination thoughts. The monoamine deficiency hypothesis states that depression is mainly caused by a deficiency of monoamine at the synaptic cleft. Thus, major efforts have been made to develop drugs that inhibit serotonin (SERT), norepinephrine (NET) and dopamine (DAT) transporters and increase the availability of these monoamines. Current gold standard treatment of MDD uses drugs that target one or more monoamine transporters. Triple reuptake inhibitors (TRIs) can target SERT, NET, and DAT simultaneously, and are believed to have the potential to be early onset antidepressants. Quantitative structure-activity relationship models were developed using machine learning algorithms in order to predict biological activities of a series of triple reuptake inhibitor compounds that showed in vitro inhibitory activity against multiple targets. The results, using mostly interpretable descriptors, showed that the internal and external predictive ability of the models are adequate, particularly of the DAT and NET by Random Forest and Support Vector Machine models. The current work shows that models developed from relatively simple, chemically interpretable descriptors can predict the activity of TRIs with similar structure in the applicability domain using ML methods.Communicated by Ramaswamy H. Sarma.

Absorption matters: A closer look at popular oral bioavailability rules for drug approvals

This study examines how two popular drug-likeness concepts used in early development, Lipinski Rule of Five (Ro5) and Veber's Rules, possibly affected drug profiles of FDA approved drugs since 1997. Our findings suggest that when all criteria are applied, relevant compounds may be excluded, addressing the harmfulness of blindly employing these rules. Of all oral drugs in the period used for this analysis, around 66 % conform to the RO5 and 85 % to Veber's Rules. Molecular Weight and calculated LogP showed low consistent values over time, apart from being the two least followed rules, challenging their relevance. On the other hand, hydrogen bond related rules and the number of rotatable bonds are amongst the most followed criteria and show exceptional consistency over time. Furthermore, our analysis indicates that topological polar surface area and total count of hydrogen bonds cannot be used as interchangeable parameters, contrary to the original proposal. This research enhances the comprehension of drug profiles that were FDA approved in the post-Lipinski period. Medicinal chemists could utilize these heuristics as a limited guide to direct their exploration of the oral bioavailability chemical space, but they must also steer the wheel to break these rules and explore different regions when necessary.

Strategic sampling with stochastic surface walking for machine learning force fields in iron's bcc-hcp phase transitions

This study developed a machine learning-based force field for simulating the bcc-hcp phase transitions of iron. By employing traditional molecular dynamics sampling methods and stochastic surface walking sampling methods, combined with Bayesian inference, we construct an efficient machine learning potential for iron. By using SOAP descriptors to map structural data, we find that the machine learning force field exhibits good coverage in the phase transition space. Accuracy evaluation shows that the machine learning force field has small errors compared to DFT calculations in terms of energy, force, and stress evaluations, indicating excellent reproducibility. Additionally, the machine learning force field accurately predicts the stable crystal structure parameters, elastic constants, and bulk modulus of bcc and hcp phases of iron, and demonstrates good performance in predicting higher-order derivatives and phase transition processes, as evidenced by comparisons with DFT calculations and existing experimental data. Therefore, our study provides an effective tool for investigating the phase transitions of iron using machine learning methods, offering new insights and approaches for materials science and solid-state physics research.

In silico characterization of aryl benzoyl hydrazide derivatives as potential inhibitors of RdRp enzyme of H5N1 influenza virus

RNA-dependent RNA polymerase (RdRp) is a potential therapeutic target for the discovery of novel antiviral agents for the treatment of life-threatening infections caused by newly emerged strains of the influenza virus. Being one of the most conserved enzymes among RNA viruses, RdRp and its inhibitors require further investigations to design novel antiviral agents. In this work, we systematically investigated the structural requirements for antiviral properties of some recently reported aryl benzoyl hydrazide derivatives through a range of in silico tools such as 2D-quantitative structure-activity relationship (2D-QSAR), 3D-QSAR, structure-based pharmacophore modeling, molecular docking and molecular dynamics simulations. The 2D-QSAR models developed in the current work achieved high statistical reliability and simultaneously afforded in-depth mechanistic interpretability towards structural requirements. The structure-based pharmacophore model developed with the docked conformation of one of the most potent compounds with the RdRp protein of H5N1 influenza strain was utilized for developing a 3D-QSAR model with satisfactory statistical quality validating both the docking and the pharmacophore modeling methodologies performed in this work. However, it is the atom-based alignment of the compounds that afforded the most statistically reliable 3D-QSAR model, the results of which provided mechanistic interpretations consistent with the 2D-QSAR results. Additionally, molecular dynamics simulations performed with the apoprotein as well as the docked complex of RdRp revealed the dynamic stability of the ligand at the proposed binding site of the receptor. At the same time, it also supported the mechanistic interpretations drawn from 2D-, 3D-QSAR and pharmacophore modeling. The present study, performed mostly with open-source tools and webservers, returns important guidelines for research aimed at the future design and development of novel anti-viral agents against various RNA viruses like influenza virus, human immunodeficiency virus-1, hepatitis C virus, corona virus, and so forth.

In silico prediction of polyethylene-aqueous and air partition coefficients of organic contaminants using linear and nonlinear approaches

Low-density polyethylene (LDPE) passive sampling is very attractive for use in determining chemicals concentrations. Crucial to the measurement is the coefficient (K_PE) describing partitioning between LDPE and environmental matrices. 255, 117 and 190 compounds were collected for the development of datasets in three different matrices, i.e., water, air and seawater, respectively. Further, 3 pp-LFER models and 9 QSPR models based on classical multiple linear regression (MLR) coupled with prevalent nonlinear algorithms (artificial neural network, ANN and support vector machine, SVM) were performed to predict LDPE-water (K_PE-W), LDPE-air (K_PE-A) and LDPE-seawater (K_PE-SW) partition coefficients. These developed models have satisfying predictability (R²_adj: 0.805-0.966, 0.963-0.991 and 0.817-0.941; RMSE_tra: 0.233-0.565, 0.200-0.406 and 0.260-0.459) and robustness (Q²_ext: 0.840-0.943, 0.968-0.984 and 0.797-0.842; RMSE_ext: 0.308-0.514, 0.299-0.426 and 0.407-0.462) in three datasets (water, air and seawater), respectively. In particular, the reasonable mechanism interpretations revealed that the molecular size, hydrophobicity, polarizability, ionization potential, and molecular stability were the most relevant properties, for governing chemicals partitioning between LDPE and environmental matrices. The application domains (ADs) assessed here exhibited the satisfactory applicability. As such, the derived models can act as intelligent tools to predict unknown K_PE values and fill the experimental gaps, which was further beneficial for the construction of enormous and reliable database to facilitate a distinct understanding of the distribution for organic contaminants in total environment.

Performance of radial distribution function-based descriptors in the chemoinformatic studies of HIV-1 protease

Aim: This letter investigates the role of radial distribution function-based descriptors for in silico design of new drugs. Methodology: The multiple linear regression models for HIV-1 protease and its complexes with a series of inhibitors were constructed. A detailed analysis of major atomic contributions to the radial distribution function descriptor weighted by the number of valence shell electrons identified residues Arg8, Asp29 and residues of the catalytic triad as crucial for the correlation with the inhibition constant, together with residues Asp30 and Ile50, whose mutations are known to cause an emergence of drug resistant variants. Conclusion: This study demonstrates an easy and fast assessment of the activity of potential drugs and the derivation of structural information of their complexes with the receptor or enzyme.

Structural refinement and prediction of potential CCR2 antagonists through validated multi-QSAR modeling studies

Chemokines trigger numerous inflammatory responses and modulate the immune system. The interaction between monocyte chemoattractant protein-1 and chemokine receptor 2 (CCR2) may be the cause of atherosclerosis, obesity, and insulin resistance. However, CCR2 is also implicated in other inflammatory diseases such as rheumatoid arthritis, multiple sclerosis, asthma, and neuropathic pain. Therefore, there is a paramount importance of designing potent and selective CCR2 antagonists despite a number of drug candidates failed in clinical trials. In this article, 83 CCR2 antagonists by Jhonson and Jhonson Pharmaceuticals have been considered for robust validated multi-QSAR modeling studies to get an idea about the structural and pharmacophoric requirements for designing more potent CCR2 antagonists. All these QSAR models were validated and statistically reliable. Observations resulted from different modeling studies correlated and validated results of other ones. Finally, depending on these QSAR observations, some new molecules were proposed that may exhibit higher activity against CCR2.

Benchmarking ligand-based virtual High-Throughput Screening with the PubChem database

With the rapidly increasing availability of High-Throughput Screening (HTS) data in the public domain, such as the PubChem database, methods for ligand-based computer-aided drug discovery (LB-CADD) have the potential to accelerate and reduce the cost of probe development and drug discovery efforts in academia. We assemble nine data sets from realistic HTS campaigns representing major families of drug target proteins for benchmarking LB-CADD methods. Each data set is public domain through PubChem and carefully collated through confirmation screens validating active compounds. These data sets provide the foundation for benchmarking a new cheminformatics framework BCL::ChemInfo, which is freely available for non-commercial use. Quantitative structure activity relationship (QSAR) models are built using Artificial Neural Networks (ANNs), Support Vector Machines (SVMs), Decision Trees (DTs), and Kohonen networks (KNs). Problem-specific descriptor optimization protocols are assessed including Sequential Feature Forward Selection (SFFS) and various information content measures. Measures of predictive power and confidence are evaluated through cross-validation, and a consensus prediction scheme is tested that combines orthogonal machine learning algorithms into a single predictor. Enrichments ranging from 15 to 101 for a TPR cutoff of 25% are observed.

A quantitative structure-retention relationship for the prediction of retention indices of the essential oils of Ammoides atlantica

A simple, descriptive and interpretable model, based on a quantitative structure-retention relationship (QSRR), was developed using the genetic algorithm-multiple linear regression (GA-MLR) approach for the prediction of the retention indices (RI) of essential oil components. By molecular modeling, three significant descriptors related to the RI values of the essential oils were identified. A data set was selected consisting of the retention indices for 32 essential oil molecules with a range of more than 931 compounds. Then, a suitable set of the molecular descriptors was calculated and the important descriptors were selected with the aid of the genetic algorithm and multiple regression method. A model with a low prediction error and a good correlation coefficient was obtained. This model was used for the prediction of the RI values of some essential oil components which were not used in the modeling procedure.