Large-scale classification of P-glycoprotein inhibitors using SMILES-based descriptors

SMILES-based QSAR virtual screening to identify potential therapeutics for COVID-19 by targeting 3CLpro and RdRp viral proteins

The COVID-19 pandemic has prompted the medical systems of many countries to develop effective treatments to combat the high rate of infection and death caused by the disease. Within the array of proteins found in SARS-CoV-2, the 3 chymotrypsin-like protease (3CLpro) holds significance as it plays a crucial role in cleaving polyprotein peptides into distinct functional nonstructural proteins. Meanwhile, RNA-dependent RNA polymerase (RdRp) takes center stage as the key enzyme tasked with replicating the viral genomic RNA within host cells. These proteins, 3CLpro and RdRp, are deemed optimal subjects for QSAR modeling due to their pivotal functions in the viral lifecycle. In this study, SMILES-based QSAR classification models were developed for a dataset of 2377 compounds that were defined as either active or inactive against 3CLpro and RdRp. Pharmacophore (PH4) and QSAR modeling were used for the virtual screening on 60.2 million compounds including ZINC, ChEMBL, Molport, and MCULE databases to identify new potent inhibitors against 3CLpro and RdRp. Then, a filter was established based on typical molecular characteristics to identify drug-like molecules. The molecular docking was also performed to evaluate the binding affinity of 156 AND 51 potential inhibitors to 3CLpro and RdRp, respectively. Among the 15 hits identified based on molecular docking scores, M3, N2, and N4 were identified as promising inhibitors due to their good synthetic accessibility scores (3.07, 3.11, and 3.29 out of 10 for M3, N2, and N4 respectively). These compounds contain amine functional groups, which are known for their crucial role in the binding interactions between drugs and their targets. Consequently, these hits have been chosen for further biological assay studies to validate their activity. They may represent novel 3CLpro and RdRp inhibitors possessing drug-like properties suitable for COVID-19 therapy.

Correlation intensity index-index of ideality of correlation: A hyphenated target function for furtherance of MAO-B inhibitory activity assessment

Monoamine oxidases are the enzymes involved in the management of brain homeostasis through oxidative deamination of monoamines such as neurotransmitters, tyramine etc. The excessive production of monoamine oxidase-B specifically results in numerous neurodegenerative disorders like Alzheimer's and Parkinson's diseases. Inhibitors of monoamine oxidase-B are applied in the management of these disorders. Here in this article we have developed robust hybrid descriptor based QSAR models related to 123 monoamine oxidase-B inhibitors through CORAL software by means of Monte Carlo optimization method. Three target functions were applied to prepare QSAR models and three splits were made for each target function. The most reliable, robust and better predictive QSAR models were developed with TF3 (correlation intensity index -index of ideality of correlation). Correlation intensity index showed positive effect on QSAR models. The structural features obtained from the QSAR modeling were incorporated in newly designed molecules and exhibited positive effect on their endpoint. Significant binding interactions were represented by these molecules in docking studies. Molecule B5 displayed prominent pIC50 (8.3) and binding affinity (-11.5 kcal mol-1) towards monoamine oxidase-B.

Recent Advances on P-Glycoprotein (ABCB1) Transporter Modelling with In Silico Methods

ABC transporters play a critical role in both drug bioavailability and toxicity, and with the discovery of the P-glycoprotein (P-gp), this became even more evident, as it plays an important role in preventing intracellular accumulation of toxic compounds. Over the past 30 years, intensive studies have been conducted to find new therapeutic molecules to reverse the phenomenon of multidrug resistance (MDR) ), that research has found is often associated with overexpression of P-gp, the most extensively studied drug efflux transporter; in MDR, therapeutic drugs are prevented from reaching their targets due to active efflux from the cell. The development of P-gp inhibitors is recognized as a good way to reverse this type of MDR, which has been the subject of extensive studies over the past few decades. Despite the progress made, no effective P-gp inhibitors to reverse multidrug resistance are yet on the market, mainly because of their toxic effects. Computational studies can accelerate this process, and in silico models such as QSAR models that predict the activity of compounds associated with P-gp (or analogous transporters) are of great value in the early stages of drug development, along with molecular modelling methods, which provide a way to explain how these molecules interact with the ABC transporter. This review highlights recent advances in computational P-gp research, spanning the last five years to 2022. Particular attention is given to the use of machine-learning approaches, drug-transporter interactions, and recent discoveries of potential P-gp inhibitors that could act as modulators of multidrug resistance.

Machine learning models for classification tasks related to drug safety

In this review, we outline the current trends in the field of machine learning-driven classification studies related to ADME (absorption, distribution, metabolism and excretion) and toxicity endpoints from the past six years (2015-2021). The study focuses only on classification models with large datasets (i.e. more than a thousand compounds). A comprehensive literature search and meta-analysis was carried out for nine different targets: hERG-mediated cardiotoxicity, blood-brain barrier penetration, permeability glycoprotein (P-gp) substrate/inhibitor, cytochrome P450 enzyme family, acute oral toxicity, mutagenicity, carcinogenicity, respiratory toxicity and irritation/corrosion. The comparison of the best classification models was targeted to reveal the differences between machine learning algorithms and modeling types, endpoint-specific performances, dataset sizes and the different validation protocols. Based on the evaluation of the data, we can say that tree-based algorithms are (still) dominating the field, with consensus modeling being an increasing trend in drug safety predictions. Although one can already find classification models with great performances to hERG-mediated cardiotoxicity and the isoenzymes of the cytochrome P450 enzyme family, these targets are still central to ADMET-related research efforts.

The Monte Carlo method to build up models of the hydrolysis half-lives of organic compounds

The hydrolysis of organic chemicals such as pesticides, pollutants, or drugs can affect the fate and behaviour of environmental contaminants, so it is of interest to evaluate the stability of substances in water for various purposes. For the registration of organic compounds in Europe, information on hydrolysis must be presented. However, the experimental measurements of all chemicals would require enormous resources, and computational models may become attractive. Applying the CORAL software (http://www.insilico.eu/coral) quantitative structure-property relationships (QSPRs) were built up to model hydrolysis. The 2D-optimal descriptor is calculated with so-called correlation weights for attributes of simplified molecular input-line entry systems (SMILES). The correlation weights are obtained as results of the special Monte Carlo optimization. The nature of (five- and six-member) rings is an important component of this approach. Another important component is the atom pair proportions for nitrogen, oxygen, and sulphur. The statistical quality of the best model is: n = 44, r² = 0.74 (training set); n = 14, r² = 0.75 (calibration set); and n = 12, r² = 0.80 (validation set).

Artificial Intelligence in Drug Discovery: A Comprehensive Review of Data-driven and Machine Learning Approaches

As expenditure on drug development increases exponentially, the overall drug discovery process requires a sustainable revolution. Since artificial intelligence (AI) is leading the fourth industrial revolution, AI can be considered as a viable solution for unstable drug research and development. Generally, AI is applied to fields with sufficient data such as computer vision and natural language processing, but there are many efforts to revolutionize the existing drug discovery process by applying AI. This review provides a comprehensive, organized summary of the recent research trends in AI-guided drug discovery process including target identification, hit identification, ADMET prediction, lead optimization, and drug repositioning. The main data sources in each field are also summarized in this review. In addition, an in-depth analysis of the remaining challenges and limitations will be provided, and proposals for promising future directions in each of the aforementioned areas.

The Monte Carlo Method as a Tool to Build up Predictive QSPR/QSAR

Background:

The Monte Carlo method has a wide application in various scientific researches. For the development of predictive models in a form of the quantitative structure-property / activity relationships (QSPRs/QSARs), the Monte Carlo approach also can be useful. The CORAL software provides the Monte Carlo calculations aimed to build up QSPR/QSAR models for different endpoints.

Methods:

Molecular descriptors are a mathematical function of so-called correlation weights of various molecular features. The numerical values of the correlation weights give the maximal value of a target function. The target function leads to a correlation between endpoint and optimal descriptor for the visible training set. The predictive potential of the model is estimated with the validation set, i.e. compounds that are not involved in the process of building up the model.

Results:

The approach gave quite good models for a large number of various physicochemical, biochemical, ecological, and medicinal endpoints. Bibliography and basic statistical characteristics of several CORAL models are collected in the present review. In addition, the extended version of the approach for more complex systems (nanomaterials and peptides), where behaviour of systems is defined by a group of conditions besides the molecular structure is demonstrated.

Conclusion:

The Monte Carlo technique available via the CORAL software can be a useful and convenient tool for the QSPR/QSAR analysis.

In silico modelling, identification of crucial molecular fingerprints, and prediction of new possible substrates of human organic cationic transporters 1 and 2

The cation membrane transporters are crucial to regulate movement of foreign molecules within the body. The present study found out structural fingerprints within molecules to be recognized as substrate/non-substrate against these transporters.

Idealization of correlations between optimal simplified molecular input-line entry system-based descriptors and skin sensitization

The Index of Ideality of Correlation (IIC) is a new criterion of the predictive potential for quantitative structure-property/activity relationships. The value of the IIC is a mathematical function sensitive to the value of the correlation coefficient and dispersion (expressed via mean absolute error). The IIC has been applied to develop QSAR models for skin sensitization achieving good predictive potential. The 'ideal correlation' is based on elementary fragments of simplified molecular input-line entry system (SMILES) and on the taking into account of the total numbers of nitrogen, oxygen, sulphur and phosphorus in the molecule.

Use of Simplified Molecular Input Line Entry System and molecular graph based descriptors in prediction and design of pancreatic lipase inhibitors

Aim: The inhibition of pancreatic lipase (PL) enzyme is the most explored strategy for the treatment of obesity. The present study describes the development of quantitative structure–activity relationship (QSAR) models for a diverse set of 293 PL inhibitors by means of the Monte Carlo optimization technique. Methodology & results: The hybrid optimal descriptors were used to build QSAR models with three subsets of three splits. The developed QSAR models were further validated with corresponding external sets. The best QSAR model has the following statistical particulars: R2 = 0.752, Q LOO 2 = 0 . 736 for the test set and R2 = 0.768, Q F 1 2 = 0 . 628 , Q F 2 2 = 0 . 621 for the validation set. Conclusion: The developed QSAR models were robust, stable and predictive and led to the design of novel PL inhibitors.

Hydroxyethylamine derivatives as HIV-1 protease inhibitors: a predictive QSAR modelling study based on Monte Carlo optimization

Application of HIV-1 protease inhibitors (as an anti-HIV regimen) may serve as an attractive strategy for anti-HIV drug development. Several investigations suggest that there is a crucial need to develop a novel protease inhibitor with higher potency and reduced toxicity. Monte Carlo optimized QSAR study was performed on 200 hydroxyethylamine derivatives with antiprotease activity. Twenty-one QSAR models with good statistical qualities were developed from three different splits with various combinations of SMILES and GRAPH based descriptors. The best models from different splits were selected on the basis of statistically validated characteristics of the test set and have the following statistical parameters: r² = 0.806, Q² = 0.788 (split 1); r² = 0.842, Q² = 0.826 (split 2); r² = 0.774, Q² = 0.755 (split 3). The structural attributes obtained from the best models were analysed to understand the structural requirements of the selected series for HIV-1 protease inhibitory activity. On the basis of obtained structural attributes, 11 new compounds were designed, out of which five compounds were found to have better activity than the best active compound in the series.

Monte Carlo method based QSAR modelling of natural lipase inhibitors using hybrid optimal descriptors

Obesity is one of the most provoking health burdens in the developed countries. One of the strategies to prevent obesity is the inhibition of pancreatic lipase enzyme. The aim of this study was to build QSAR models for natural lipase inhibitors by using the Monte Carlo method. The molecular structures were represented by the simplified molecular input line entry system (SMILES) notation and molecular graphs. Three sets - training, calibration and test set of three splits - were examined and validated. Statistical quality of all the described models was very good. The best QSAR model showed the following statistical parameters: r² = 0.864 and Q² = 0.836 for the test set and r² = 0.824 and Q² = 0.819 for the validation set. Structural attributes for increasing and decreasing the activity (expressed as pIC₅₀) were also defined. Using defined structural attributes, the design of new potential lipase inhibitors is also presented. Additionally, a molecular docking study was performed for the determination of binding modes of designed molecules.