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

The searching for agents for Alzheimer's disease treatment via the system of self-consistent models

Robust quantitative structure-activity relationships (QSARs) for hBACE-1 inhibitors (pIC50) for a large database (n = 1706) are established. New statistical criteria of the predictive potential of models are suggested and tested. These criteria are the index of ideality of correlation (IIC) and the correlation intensity index (CII). The system of self-consistent models is a new approach to validate the predictive potential of QSAR-models. The statistical quality of models obtained using the CORAL software (http://www.insilico.eu/coral) for the validation sets is characterized by the average determination coefficient R²_v= 0.923, and RMSE =0.345. Three new promising molecular structures which can become inhibitors hBACE-1 are suggested.

The QSAR-search of effective agents towards coronaviruses applying the Monte Carlo method

Perhaps there is some similarity between the coronavirus of 2017 and the COVID-19. Consequently, a predictive model for the antiviral activity for the Middle East respiratory syndrome coronavirus (MERS-CoV, 2017) could be useful for designing the strategy and tactics in the struggle with coronaviruses in general and with COVID 19 in particular. Quantitative structure-activity relationships (QSARs) of inhibitory activity to MERS-CoV were developed. The index of ideality of correlation was applied to build up these models for the antiviral activity. The statistical quality of the best model is quite good (r² = 0.84). A mechanistic interpretation of these models based on the molecular features with strong positive (i.e. promoters for endpoint increase) and strong negative (i.e. promoters for endpoint decrease) influence on the inhibitory activity is suggested. A collection of possible biologically active compounds, constructed using data on the above molecular features which are statistically reliable promoters of increase or decrease of the activity, is presented.

Can the Monte Carlo method predict the toxicity of binary mixtures?

Risk assessment of toxicants mainly is a result of experiments with single substances. However, toxicity in natural ecosystems typically does not result from single toxicant exposure but is rather a result of exposure to mixtures of toxicants. It is not surprising a mixture of toxicity is a subject of eco-toxicological interest for several decades. A quantitative structure-activity relationships (QSAR)-based approach is an attractive approach to assessing the joint effects in the binary mixtures. The validity of the proposed approach was demonstrated by comparing the predicted values against the experimentally determined values. Simplified molecular input-line entry system (SMILES) is used for the representation of the molecular structures of components of two-component mixtures to build up QSAR. The SMILES-based models are improving if the Monte Carlo optimization aimed to define 2D-optimal descriptors apply the so-called index of ideality of correlation (IIC), which is a mathematical function of both the correlation coefficient and mean absolute error calculated for the positive and negative difference between observed and calculated values of toxicity. The average statistical quality of these models (for the validation set) is n=25, R²=0.95, and RMSE=0.375.

Discovery of New Hydroxyethylamine Analogs against 3CLpro Protein Target of SARS-CoV-2: Molecular Docking, Molecular Dynamics Simulation, and Structure-Activity Relationship Studies

The novel coronavirus, SARS-CoV-2, has caused a recent pandemic called COVID-19 and a severe health threat around the world. In the current situation, the virus is rapidly spreading worldwide, and the discovery of a vaccine and potential therapeutics are critically essential. The crystal structure for the main protease (Mpro) of SARS-CoV-2, 3-chymotrypsin-like cysteine protease (3CLpro), was recently made available and is considerably similar to the previously reported SARS-CoV. Due to its essentiality in viral replication, it represents a potential drug target. Herein, a computer-aided drug design (CADD) approach was implemented for the initial screening of 13 approved antiviral drugs. Molecular docking of 13 antivirals against the 3-chymotrypsin-like cysteine protease (3CLpro) enzyme was accomplished, and indinavir was described as a lead drug with a docking score of -8.824 and a XP Gscore of -9.466 kcal/mol. Indinavir possesses an important pharmacophore, hydroxyethylamine (HEA), and thus, a new library of HEA compounds (>2500) was subjected to virtual screening that led to 25 hits with a docking score more than indinavir. Exclusively, compound 16 with a docking score of -8.955 adhered to drug-like parameters, and the structure-activity relationship (SAR) analysis was demonstrated to highlight the importance of chemical scaffolds therein. Molecular dynamics (MD) simulation analysis performed at 100 ns supported the stability of 16 within the binding pocket. Largely, our results supported that this novel compound 16 binds with domains I and II, and the domain II-III linker of the 3CLpro protein, suggesting its suitability as a strong candidate for therapeutic discovery against COVID-19.

SARS-CoV Mpro inhibitory activity of aromatic disulfide compounds: QSAR model

The main protease (M^pro) of SARS-associated coronavirus (SARS-CoV) had caused a high rate of mortality in 2003. Current events (2019-2020) substantiate important challenges for society due to coronaviruses. Consequently, advancing models for the antiviral activity of therapeutic agents is a necessary component of the fast development of treatment for the virus. An analogy between anti-SARS agents suggested in 2017 and anti-coronavirus COVID-19 agents are quite probable. Quantitative structure-activity relationships for SARS-CoV are developed and proposed in this study. The statistical quality of these models is quite good. Mechanistic interpretation of developed models is based on the statistical and probability quality of molecular alerts extracted from SMILES. The novel, designed structures of molecules able to possess anti-SARS activities are suggested. For the final assessment of the designed molecules inhibitory potential, developed from the obtained QSAR model, molecular docking studies were applied. Results obtained from molecular docking studies were in a good correlation with the results obtained from QSAR modeling.

In silico enhancement of azo dye adsorption affinity for cellulose fibre through mechanistic interpretation under guidance of QSPR models using Monte Carlo method with index of ideality correlation

Azo dyes are a group of chemical moieties joined by azo (-N=N-) group with potential usefulness in different industrial applications. But these dyes are not devoid of hazardous consequence because of poor affinity for the fibre and discharge into the water stream. The chemical aspects of 72 azo dyes towards cellulose fibre in terms of their affinity by QSPR have been explored in the present work. We have employed two approaches, namely balance of correlation without IIC (TF₁) and balance of correlation with IIC (TF₂), to generate 16 QSAR models from 8 splits. The determination coefficient of calibration and validation set was found higher when the QSPR models were developed using the index of ideality correlation (IIC) parameter (TF₂). The model developed with TF₂ for split 3 was considered as a prominent model because the determination coefficient of the validation set was maximum (r ² = 0.9468). The applicability domain (AD) was also analysed based on 'statistical defect', d(A) for a SMILES attribute. The mechanistic interpretation was done by identifying the SMILES attributes responsible for the promoter of endpoint increase and promoter of endpoint decrease. These SMILES attributes were applied to design 15 new dyes with higher affinity for cellulose fibre.

An Overview of Antiretroviral Agents for Treating HIV Infection in Paediatric Population

Paediatric Acquired ImmunoDeficiency Syndrome (AIDS) is a life-threatening and infectious disease in which the Human Immunodeficiency Virus (HIV) is mainly transmitted through Mother-To- Child Transmission (MTCT) during pregnancy, labour and delivery, or breastfeeding. This review provides an overview of the distinct therapeutic alternatives to abolish the systemic viral replication in paediatric HIV-1 infection. Numerous classes of antiretroviral agents have emerged as therapeutic tools for downregulation of different steps in the HIV replication process. These classes encompass Non- Nucleoside Analogue Reverse Transcriptase Inhibitors (NNRTIs), Nucleoside/Nucleotide Analogue Reverse Transcriptase Inhibitors (NRTIs/NtRTIs), INtegrase Inhibitors (INIs), Protease Inhibitors (PIs), and Entry Inhibitors (EIs). Co-administration of certain antiretroviral drugs with Pharmacokinetic Enhancers (PEs) may boost the effectiveness of the primary therapeutic agent. The combination of multiple antiretroviral drug regimens (Highly Active AntiRetroviral Therapy - HAART) is currently the standard therapeutic approach for HIV infection. So far, the use of HAART offers the best opportunity for prolonged and maximal viral suppression, and preservation of the immune system upon HIV infection. Still, the frequent administration of high doses of multiple drugs, their inefficient ability to reach the viral reservoirs in adequate doses, the development of drug resistance, and the lack of patient compliance compromise the complete HIV elimination. The development of nanotechnology-based drug delivery systems may enable targeted delivery of antiretroviral agents to inaccessible viral reservoir sites at therapeutic concentrations. In addition, the application of Computer-Aided Drug Design (CADD) approaches has provided valuable tools for the development of anti-HIV drug candidates with favourable pharmacodynamics and pharmacokinetic properties.

QSPR/QSAR: State-of-Art, Weirdness, the Future

Ability of quantitative structure-property/activity relationships (QSPRs/QSARs) to serve for epistemological processes in natural sciences is discussed. Some weirdness of QSPR/QSAR state-of-art is listed. There are some contradictions in the research results in this area. Sometimes, these should be classified as paradoxes or weirdness. These points are often ignored. Here, these are listed and briefly commented. In addition, hypotheses on the future evolution of the QSPR/QSAR theory and practice are suggested. In particular, the possibility of extending of the QSPR/QSAR problematic by searching for the "statistical similarity" of different endpoints is suggested and illustrated by an example for relatively "distanced each from other" endpoints, namely (i) mutagenicity, (ii) anticancer activity, and (iii) blood-brain barrier.

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.

Good and bad molecular fingerprints for human rhinovirus 3C protease inhibition: identification, validation, and application in designing of new inhibitors through Monte Carlo-based QSAR study

HRV 3 C protease (HRV 3Cpro) is an important target for common cold and upper respiratory tract infection. Keeping in view of the non-availability of drug for the treatment, newer computer-based modelling strategies should be applied to rationalize the process of antiviral drug discovery in order to decrease the valuable time and huge expenditure of the process. The present work demonstrates a structure wise optimization using Monte Carlo-based QSAR method that decomposes ligand compounds (in SMILES format) into several molecular fingerprints/descriptors. The current state-of-the-art in QSAR study involves the balance of correlation approach using four different sets: training, invisible training, calibration, and validation. The final models were also validated through mean absolute error, index of ideality of correlation, Y-randomization and applicability domain analysis. R2 and Q2 values for the best model were 0.8602, 0.8507 (training); 0.8435, 0.8331 (invisible training); 0.7424, 0.7020 (calibration); 0.5993, 0.5216 (validation), respectively. The process identified some molecular substructures as good and bad fingerprints depending on their effect to increase or decrease the HRV 3Cpro inhibition. Finally, new inhibitors were designed based on the fundamental concept to replace the bad fragments with the good fragments as well as including more good fragments into the structure. The study points out the importance of the fingerprint based drug design strategy through Monte Carlo optimization method in the modelling of HRV 3Cpro inhibitors.

Insight into structural features of phenyltetrazole derivatives as ABCG2 inhibitors for the treatment of multidrug resistance in cancer

ABCG2 is the principal ABC transporter involved in the multidrug resistance of breast cancer. Looking at the current demand in the development of ABCG2 inhibitors for the treatment of multidrug-resistant cancer, we have explored structural requirements of phenyltetrazole derivatives for ABCG2 inhibition by combining classical QSAR, Bayesian classification modelling and molecular docking studies. For classical QSAR, structural descriptors were calculated from the free software tool PaDEL-descriptor. Stepwise multiple linear regression (SMLR) was used for model generation. A statistically significant model was generated and validated with different parameters (For training set: r = 0.825; Q² = 0.570 and for test set: r = 0.894, r²_pred = 0.783). The predicted model was found to satisfy the Golbraikh and Trospha criteria for model acceptability. Bayesian classification modelling was also performed (ROC scores were 0.722 and 0.767 for the training and test sets, respectively). Finally, the binding interactions of phenyltetrazole type inhibitor with the ABCG2 receptor were mapped with the help of molecular docking study. The result of the docking analysis is aligned with the classical QSAR and Bayesian classification studies. The combined modelling study will guide the medicinal chemists to act faster in the drug discovery of ABCG2 inhibitors for the management of resistant breast cancer.

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.

Exploration of good and bad structural fingerprints for inhibition of indoleamine-2,3-dioxygenase enzyme in cancer immunotherapy using Monte Carlo optimization and Bayesian classification QSAR modeling

Indoleamine-2,3-dioxygenase 1 (IDO1) is an extrahepatic, heme-containing and tryptophan-catalyzing enzyme responsible for causing blockade of T-cell proliferation and differentiation by depleting tryptophan level in cancerous cells. Therefore, inhibition of IDO1 may be a useful strategy for immunotherapy against cancer. In this study, 448 structurally diverse IDO1 inhibitors with a wide range of activity has been taken into consideration for classification QSAR analysis through Monte Carlo Optimization by using different splits as well as different combinations of SMILES-based, graph-based and hybrid descriptors. The best model from Monte Carlo optimization was interpreted to find out the good and bad structural fingerprints for IDO1 and further justified by using Bayesian classification QSAR modeling. Among the three splits in Monte Carlo optimization, the statistics of the best model was obtained from Split 3: sensitivity = 0.87, specificity = 0.91, accuracy = 0.89 and MCC = 0.78. In Bayesian classification modeling, the ROC scores for training and test set were found to be 0.91 and 0.86, respectively. The combined modeling analysis revealed that the presence of aryl hydrazyl sulphonyl moiety, furazan ring, halogen substitution, nitro group and hetero atoms in aromatic system can be very useful in designing IDO1 inhibitors. All the good and bad structural fingerprints for IDO1 were identified and are justified by correlating these fragments to the inhibition of IDO1 enzyme. These structural fingerprints will guide the researchers in this field to design better inhibitors against IDO1 enzyme for cancer immunotherapy.Communicated by Ramaswamy H. Sarma.

Structural exploration of hydroxyethylamines as HIV-1 protease inhibitors: new features identified

The current study deals with chemometric modelling strategies (Naïve Bayes classification, hologram-based quantitative structure-activity relationship (HQSAR), comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA)) to explore the important features of hydroxylamine derivatives for exerting potent human immunodeficiency virus-1 (HIV-1) protease inhibition. Depending on the statistically validated reliable and robust quantitative structure-activity relationship (QSAR) models, important and crucial structural features have been identified that may be responsible for enhancing the activity profile of these hydroxylamine compounds. Arylsulfonamide function along with methoxy or fluoro substitution is important for enhancing activity. Bulky steric substitution at the sulfonamide nitrogen disfavours activity whereas smaller hydrophobic substitution at the same position is found to be favourable. Apart from the crucial oxazolidinone moiety, pyrrolidine, cyclic urea and methyl ester functions are also responsible for increasing the HIV-1 protease inhibitory profile. Observations derived from these modelling studies may be utilized further in designing promising HIV-1 protease inhibitors of this class.

Identification of structural requirements and prediction of inhibitory activity of natural flavonoids against Zika virus through molecular docking and Monte Carlo based QSAR Simulation

There has been growing interest in the research of flavonoids due to their potential antiviral activities. Recently, some natural flavonoids have shown potential inhibitory activity against zika virus NS3-NS2B protease. In order to accelerate the drug discovery efforts using flavonoids, a Monte Carlo simulation-based QSAR method has been applied to find out the structural requirements for the inhibitory activity. The best QSAR model was obtained using SMILES descriptors and HSG descriptors with ¹EC connectivity with the following statistical parameters: R ² = 0.9569 and Q ² = 0.9050 for the test set. The best model was further utilised for the prediction of inhibitory activity of some other natural flavonoids. Four flavonoids (amentoflavone, fisetin, isorhamnetin and theaflavin-3-gallate) have shown higher predicted inhibitory activity and further validated by performing docking analysis. This study may help in understanding and performing natural flavonoids-based drug discovery against zika virus.