Quantitative Structure-activity Relationship Models Research Articles

First report on the QSAR modelling and multistep virtual screening of the inhibitors of nonstructural protein Nsp14 of SARS-CoV-2: Reducing unnecessary chemical synthesis and experimental tests

Corona Virus Disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), poses a serious threat to human health and life safety. How to effectively prevent and treat COVID-19 is crucial. In this study, we used the inhibitors of nonstructural protein Nsp14 of SARS-CoV-2 to perform the quantitative structure activity relationship (QSAR) modelling for the first time. Based on different dataset division strategies, we selected partial least square (PLS) and multiple linear regression (MLR) methods to develop easily interpretable and reproducible QSAR models with 2D molecular descriptors. All models complied with the strict QSAR validation principles of OECD and internationally recognized validation metrics. The best model contained two molecular descriptors with the following statistical parameters: R2 = 0.7796, QLOO2= 0.7373, Rtest2 = 0.8539 and CCCtest = 0.9073. Obviously, the model exhibited good prediction performance and can be used for quickly predicting the inhibitory activity of unknown compounds against Nsp14. Mechanistic interpretation identified the detailed relationship between molecular structure information and inhibitory activity. The best QSAR model was used to predict the inhibitory activity of 263 true external compounds without experimental values against Nsp14, and the prediction reliability was analyzed and discussed. Molecular docking and ADMET analyses were conducted for compounds with higher similarity to the modelling compounds. Finally, two compounds were identified as potential candidate drugs of targeting Nsp14. The current work lays a solid theoretical foundation for the discovery of inhibitors targeting Nsp14, and has an important reference significance for the development of anti-COVID-19 drugs.

A machine learning-based QSAR model reveals important molecular features for understanding the potential inhibition mechanism of ionic liquids to acetylcholinesterase

The broad application of ionic liquids (ILs) has been hindered by uncertainties surrounding their ecotoxicity. In this work, a Quantitative Structure-Activity Relationship (QSAR) model was devised to predict the inhibition of ILs towards the activity of AChE, employing both Random Forest (RF) and eXtreme Gradient Boosting (XGBoost) machine learning approaches. Fourteen kings of essential molecular feature descriptors were screened from an initial roster of 244 descriptors through the application of a feature importance index and they showed a significant impact on the activity of AChE activity. The two models based solely on the 14 most critical molecular descriptors could maintain model's robustness and reliability. The correlation analysis between these 14 descriptors and the inhibition of AChE activity revealed the potential impact of the molecular characteristics on ILs toxicity. The results underscored the main influence of cations in ILs on the inhibitory activity towards the AChE enzyme. Specifically, cations exhibiting hydrophobicity properties were found to exert more potent inhibitory effects on the AChE enzyme. In addition, some other properties of the cations, such as the degree of branching, atomic weight and partial charge also modulated their inhibition potential. This study enhances the comprehension of the structure-activity relationship between ILs and AChE inhibition, providing a reference for designing safer and greener ILs.

Aminopyrimidine Derivatives as Multiflavivirus Antiviral Compounds Identified from a Consensus Virtual Screening Approach.

Around three billion people are at risk of infection by the dengue virus (DENV) and potentially other flaviviruses. Worldwide outbreaks of DENV, Zika virus (ZIKV), and yellow fever virus (YFV), the lack of antiviral drugs, and limitations on vaccine usage emphasize the need for novel antiviral research. Here, we propose a consensus virtual screening approach to discover potential protease inhibitors (NS3pro) against different flavivirus. We employed an in silico combination of a hologram quantitative structure-activity relationship (HQSAR) model and molecular docking on characterized binding sites followed by molecular dynamics (MD) simulations, which filtered a data set of 7.6 million compounds to 2,775 hits. Lastly, docking and MD simulations selected six final potential NS3pro inhibitors with stable interactions along the simulations. Five compounds had their antiviral activity confirmed against ZIKV, YFV, DENV-2, and DENV-3 (ranging from 4.21 ± 0.14 to 37.51 ± 0.8 μM), displaying aggregator characteristics for enzymatic inhibition against ZIKV NS3pro (ranging from 28 ± 7 to 70 ± 7 μM). Taken together, the compounds identified in this approach may contribute to the design of promising candidates to treat different flavivirus infections.

Quantitative Predictive Studies of Multiple Biological Activities of TRPV1 Modulators.

TRPV1 channel agonists and antagonists, which have powerful analgesic effects without the addictive qualities associated with traditional analgesics, have become a focus area for the development of novel analgesics. In this study, quantitative structure-activity relationship (QSAR) models for three bioactive endpoints (Ki, IC50, and EC50) were successfully constructed using four machine learning algorithms: SVM, Bagging, GBDT, and XGBoost. These models were based on 2922 TRPV1 modulators and incorporated four types of molecular descriptors: Daylight, E-state, ECFP4, and MACCS. After the rigorous five-fold cross-validation and external test set validation, the optimal models for the three endpoints were obtained. For the Ki endpoint, the Bagging-ECFP4 model had a Q2 value of 0.778 and an R2 value of 0.780. For the IC50 endpoint, the XGBoost-ECFP4 model had a Q2 value of 0.806 and an R2 value of 0.784. For the EC50 endpoint, the SVM-Daylight model had a Q2 value of 0.784 and an R2 value of 0.809. These results demonstrate that the constructed models exhibit good predictive performance. In addition, based on the model feature importance analysis, the influence between substructure and biological activity was also explored, which can provide important theoretical guidance for the efficient virtual screening and structural optimization of novel TRPV1 analgesics. And subsequent studies on novel TRPV1 modulators will be based on the feature substructures of the three endpoints.

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.

Control list of high-priority chemicals based on 5-HT-RI functionality and the human health interference effects selective CNN-GRU deep learning model

The antidepressant drug known as 5-HT reuptake inhibitor (5-HT-RI) was commonly detected in biological tissues and result in significant adverse health effects. Homology modeling was used to characterize the functionalities (efficacy and resistance), and the adverse outcome pathway was used to characterize its human health interferences (olfactory toxicity, neurotoxicity, and gut microbial interference). The convolutional neural network coupled with the gated recurrent unit (CNN-GRU) deep learning method was used to construct a comprehensive model of 5-HT-RI functionality and human health interference effects selectivity with small sample data. The architecture with 2 SE, 320 neuronal nodes and 6-folds cross-validation showed the best applicability. The results showed that the confidence interval of the constructed model reached 90 % indicating that the model had reliable prediction ability and generalization ability. Based on the CNN-GRU deep learning model, seven high-priority chemicals with a weak comprehensive effect, including D-VEN, (1R,4S)-SER, S-FLX, CTP, S-CTP, NEF, and VEN, were screened. Based on the molecular three-dimensional structure information, a comprehensive-effect three-dimensional quantitative structure-activity relationship (3D-QSAR) model was constructed to confirm the reliability of the constructed control list of 5-HT-RI high-priority chemicals. Analysis with the ranking of calculated values based on the molecular dynamics method and predicted values based on the CNN-GRU deep learning model, we found that the consistency of the three methods was above 85 %. Additionally, by analyzing the sensitivity, molecular electrostatic potential, polar surface area of the comprehensive-effect CNN-GRU deep learning model, and the electrostatic field of the 3D-QSAR models, we found that the significant effects of five key characteristics (DM, Qyy, Qxz, I, and BP), molecular electronegativity, and polarity significantly affected the high-priority degree of 5-HT-RI. In this study, we provided reasonable and reliable prediction tools and discussed theoretical methods for the risk assessment of functionality and human health interference of emerging pollutants such as 5-HT-RI.

QSAR studies of quinoline alkaloids camptothecin derivatives for prediction anticancer activity using linear and nonlinear methods

AbstractQuinoline alkaloid camptothecin (CPT) derivatives are compounds with a wide range of inhibitory activities and they are the basis of several groups of drugs. CPT is one of the prominent lead compounds in anticancer drug development. It inhibits DNA topoisomerase I (topo I) enzyme and has shown remarkable anticancer activity against lung, colon, rectum, ovarian, breast, and stomach cancers. Quantitative structure–activity relationship (QSAR) is basically a statistical approach correlating the response activity data with descriptors encoding chemical information. In the current research, a QSAR study was performed for modeling and predicting the anticancer activity (pIC50) of 76 CPT derivatives as an inhibitor of DNA topo I using the genetic algorithm multiple linear regression (GA‐MLR) and back propagation artificial neural network (BP‐ANN) methods. The results of this study indicate that the use of constitutional and geometrical descriptors provide good estimate for pIC50. The obtained results of statistical criteria, internal and external validation show the superiority of BP‐ANN model than GA‐MLR to predict the pIC50 values of the investigated compounds. QSAR model developed in this study can be used for design and development of novel potent CPT derivatives and for predicting their anticancer activity.

In silico study for discovering novel thiazole derivatives as anti-breast cancer agents (MCF-7)

The work validated the pharmacokinetic characteristics and docking approach, and discovered novel thiazole compounds with MCF-7 breast anti-cancer efficacy using quantitative structure-activity relationships (QSAR) models. Using an experimental dataset of 53 thiazole derivatives with IC50 values, QSAR models were developed using multivariate linear regression (QSARMLR) and artificial neural networks (QSARANN). Eight descriptions, with statistical values of R2train = 0.875 and Q2LOO = 0.834, comprised the successfully developed QSARMLR model. Based on the descriptors of the QSARMLR model, which has statistical values of R2 = 0.918, Q2test = 0.934, and Q2CV = 0.916, the QSARANN neural network model with the architectural network of I(8)-HL(9)-O(1) has also been constructed. Using the AD and outliers analysis, the models were utilized to forecast 120 new design derivatives based on the thiazole framework, and the IC50 values of 10 new thiazole derivatives were determined. Additionally, the derivatives were assessed for resistance to estrogen-positive breast cancer by docking them onto the Polo-like kinases (Plk1) receptor and screening them for pharmacokinetic features in accordance with Lipinski and Ghose guidelines. Consequently, it was discovered that the thiazole TAZ5 had promising derivative activity against the primary MCF-7 breast cancer cell line.

Exploring Citrus sinensis Phytochemicals as Potential Inhibitors for Breast Cancer Genes BRCA1 and BRCA2 Using Pharmacophore Modeling, Molecular Docking, MD Simulations, and DFT Analysis.

Structure-activity relationship (SAR) is considered to be an effective in silico approach when discovering potential antagonists for breast cancer due to gene mutation. Major challenges are faced by conventional SAR in predicting novel antagonists due to the discovery of diverse antagonistic compounds. Methodologyand Results: In predicting breast cancer antagonists, a multistep screening of phytochemicals isolated from the seeds of the Citrus sinensis plant was applied using feasible complementary methodologies. A three-dimensional quantitative structure-activity relationship (3D-QSAR) model was developed through the Flare project, in which conformational analysis, pharmacophore generation, and compound alignment were done. Ten hit compounds were obtained through the development of the 3D-QSAR model. For exploring the mechanism of action of active compounds against cocrystal inhibitors, molecular docking analysis was done through Molegro software (MVD) to identify lead compounds. Three new proteins, namely, 1T15, 3EU7, and 1T29, displayed the best Moldock scores. The quality of the docking study was assessed by a molecular dynamics simulation. Based on binding affinities to the receptor in the docking studies, three lead compounds (stigmasterol P8, epoxybergamottin P28, and nobiletin P29) were obtained, and they passed through absorption, distribution, metabolism, and excretion (ADME) studies via the SwissADME online service, which proved that P28 and P29 were the most active allosteric inhibitors with the lowest toxicity level against breast cancer. Then, density functional theory (DFT) studies were performed to measure the active compound's reactivity, hardness, and softness with the help of Gaussian 09 software. This multistep screening of phytochemicals revealed high-reliability antagonists of breast cancer by 3D-QSAR using flare, docking analysis, and DFT studies. The present study helps in providing a proper guideline for the development of novel inhibitors of BRCA1 and BRCA2.

Establishing short-term occupational exposure limits (STELs) for sensory irritants using predictive and in silico respiratory rate depression (RD50) models

Sensory irritation is a health endpoint that serves as the critical effect basis for many occupational exposure limits (OELs). Schaper 1993 described a significant relationship with high correlation between the measured exposure concentration producing a 50% respiratory rate decrease (RD50) in a standard rodent assay and the American Conference of Governmental Industrial Hygienists (ACGIH®) Threshold Limit Values (TLVs®) as time-weighted averages (TWAs) for airborne chemical irritants. The results demonstrated the potential use of the RD50 values for deriving full-shift TWA OELs protective of irritant responses. However, there remains a need to develop a similar predictive model for deriving workplace short-term exposure limits (STELs) for sensory irritants. The aim of our study was to establish a model capable of correlating the relationship between RD50 values and published STELs to prospectively derive short-term exposure OELs for sensory irritants. A National Toxicology Program (NTP) database that included chemicals with both an RD50 and established STELs was used to fit several linear regression models. A strong correlation between RD50s and STELs was identified, with a predictive equation of ln (STEL) (ppm) = 0.86 * ln (RD50) (ppm) − 2.42 and an R2 value of 0.75. This model supports the use of RD50s to derive STELs for chemicals without existing exposure recommendations. Further, for data-poor sensory irritants, predicted RD50 values from in silico quantitative structure activity relationship (QSAR) models can be used to derive STELs. Hence, in silico methods and statistical modeling can present a path forward for establishing reliable OELs and improving worker safety and health.

QSAR models for the ozonation of diverse volatile organic compounds at different temperatures.

In order to assess the fate and persistence of volatile organic compounds (VOCs) in the atmosphere, it is necessary to determine their oxidation rate constants for their reaction with ozone (kO3). However, given that experimental values of kO3 are only available for a few hundred compounds and their determination is expensive and time-consuming, developing predictive models for kO3 is of great importance. Thus, this study aimed to develop reliable quantitative structure-activity relationship (QSAR) models for 302 values of 149 VOCs across a broad temperature range (178-409 K). The model was constructed based on the combination of a simplified molecular-input line-entry system (SMILES) and temperature as an experimental condition, namely quasi-SMILES. In this study, temperature was incorporated in the models as an independent feature. The hybrid optimal descriptor generated from the combination of quasi-SMILES and HFG (hydrogen-filled graph) was used to develop reliable, accurate, and predictive QSAR models employing the CORAL software. The balance between the correlation method and four different target functions (target function without considering IIC or CII, target function using each IIC or CII, and target function based on the combination of IIC and CII) was used to improve the predictability of the QSAR models. The performance of the developed models based on different target functions was compared. The correlation intensity index (CII) significantly enhanced the predictability of the model. The best model was selected based on the numerical value of Rm2 of the calibration set (split #1, Rtrain2 = 0.9834, Rcalibration2 = 0.9276, Rvalidation2 = 0.9136, and calibration = 0.8770). The promoters of increase/decrease for log kO3 were also computed based on the best model. The presence of a double bond (BOND10000000 and $10 000 000 000), absence of halogen (HALO00000000), and the nearest neighbor codes for carbon equal to 321 (NNC-C⋯321) are some significant promoters of endpoint increase.

Low-cost quantum mechanical descriptors for data efficient skin sensitization QSAR models

Quantitative Structure Activity Relationship modelling methodologies need to incorporate relevant mechanistic information to have high predictive performance and validity. Electrophilic reactivity is a common mechanistic feature of skin sensitization endpoints which could be concisely characterized with electronic descriptors which is key to enabling the modelling of small datasets in this domain. However, quantum mechanical methodologies have previously featured high computational costs which would exclude the use of large datasets. Consequently, we investigate the use of electronic descriptors calculated using the Hartree Fock with 3 corrections (Hf-3c) method, a low-cost ab initio methodology that has higher chemical accuracy than previous semiempirical methodologies for modelling in vitro skin sensitization assay outcomes. We also model the Ames assay as a surrogate for determining skin sensitization outcomes. The quantum chemical descriptors calculated using the Hf-3c method with conductor-like polarizable continuum model (CPCM) implicit solvation found improved QSAR model performance for the in vitro Ames (n = 6049, 0.770 AUC), KeratinoSens (n = 164, 0.763 AUC), and Direct Peptide Reactivity Assay (n = 122, 0.750 AUC) datasets, with their combination producing high predictive performance for unseen in vivo Local Lymph Node Assay (n = 86, 0.789 AUC) and Human Repeated Insult Patch Test (n = 86, 0.791 AUC) assay toxicant outcomes.

Computational Investigation of Corrosion Inhibition Properties of Hydrazine and Its Derived Molecules for Mild Steel

Corrosion of mild steel has gained a lot of attention by many industries and academia. This study was aimed to find the corrosion inhibition properties of hydrazine and its 11 derived molecules for mild steel using computational method, namely, the density functional theory with Beck 3‐parameter Lee Yang and Parr functional along with the 6‐311++G (d, p) basis set using simple quantum chemical parameters. The main quantum chemical parameters were energies of highest occupied and lowest unoccupied molecular orbitals and their energy gap, dipole moment, hardness, softness, electrophilicity index, the fractions of electrons transmitted, the change energy back‐donation, electron‐donating power, electron‐accepting power, molecular electrostatic potential, and Fukui function as well as polarizability. The results indicated that hydrazine derivatives have good anticorrosion nature which is consistent with the experimental findings. Out of the 12 target inhibitor molecules, N′‐(3‐nitrobenzylidene) hydrazine carbodithioic acid (NBHCA) and N′‐(4‐dimethylaminobenzylidene) hydrazine carbodithioic acid (DABHCA) are the most effective inhibitors. According to the value of the quantum chemical parameters, EHOMO, ΔEgap, Χ, η, σ, ε, ΔN, α, and ΔEb-d DABHCA showed the top effective anticorrosion inhibitor, and by the value of the quantum chemical parameters, ELUMO, TE, ω−, and ω+ NBHCA showed an effective anticorrosion inhibitor. The quantum chemical parameter data showed that NBHCA and DABHCA are nearly equivalent in their values, which is consistent with experimental values. Furthermore, a multiple‐linear regression equation was proposed as a quantitative structural activity relationship model to relate the calculated molecular descriptors of hydrazine derivatives to their experimental inhibition efficiencies obtained from literature. As a result, this model showed an excellent accuracy with an R‐squared of 0.924 and this uses to predict inhibition efficiencies of another molecules.

A review on the structural characterization of nanomaterials for nano-QSAR models.

Quantitative structure-activity relationship (QSAR) models are routinely used to predict the properties and biological activity of chemicals to direct synthetic advances, perform massive screenings, and even to register new substances according to international regulations. Currently, nanoscale QSAR (nano-QSAR) models, adapting this methodology to predict the intrinsic features of nanomaterials (NMs) and quantitatively assess their risks, are blooming. One of the challenges is the characterization of the NMs. This cannot be done with a simple SMILES representation, as for organic molecules, because their chemical structure is complex, including several layers and many inorganic materials, and their size and geometry are key features. In this review, we survey the literature for existing predictive models for NMs and discuss the variety of calculated and experimental features used to define and describe NMs. In the light of this research, we propose a classification of the descriptors including those that directly describe a component of the nanoform (core, surface, or structure) and also experimental features (related to the nanomaterial's behavior, preparation, or test conditions) that indirectly reflect its structure.

ARKA: a framework of dimensionality reduction for machine-learning classification modeling, risk assessment, and data gap-filling of sparse environmental toxicity data.

Due to the lack of experimental toxicity data for environmental chemicals, there arises a need to fill data gaps by in silico approaches. One of the most commonly used in silico approaches for toxicity assessment of small datasets is the Quantitative Structure-Activity Relationship (QSAR), which generates predictive models for the efficient prediction of query compounds. However, the reliability of the predictions from QSARs derived from small datasets is often questionable from a statistical point of view. This is due to the presence of a larger number of descriptors as compared to the number of training compounds, which reduces the degree of freedom of the developed model. To reduce the overall prediction error for a particular QSAR model, we have proposed here the computation of the novel Arithmetic Residuals in K-groups Analysis (ARKA) descriptors. We have reduced the number of modeling descriptors in a supervised manner by partitioning them into K classes (K = 2 here) depending on the higher mean normalized values of the descriptors to a particular response class, thus preventing the loss of chemical information. A scatter plot of the data points using the values of two ARKA descriptors (ARKA_2 vs. ARKA_1) can potentially identify activity cliffs, less confident data points, and less modelable data points. We have used here five representative environmentally relevant endpoints (skin sensitization, earthworm toxicity, milk/plasma partitioning, algal toxicity, and rodent carcinogenicity of hazardous chemicals) with graded responses to which the ARKA framework was applied for classification modeling. On comparing the performance of the models generated using conventional QSAR descriptors and the ARKA descriptors, the prediction quality of the models derived from ARKA descriptors was found, based on multiple graded-data validation metrics-derived decision criteria, much better than the models derived from QSAR descriptors signifying the potential of ARKA descriptors in ecotoxicological classification modeling of small data sets. Additionally, this holds true for the Read-Across approach as well, since the Read-Across predictions using ARKA descriptors supersede the predictions generated from QSAR descriptors. For the ease of users, a Java-based expert system has been developed that computes the ARKA descriptors from the input of QSAR descriptors.

Identification of apigenin-4'-glucoside as bacterial DNA gyrase inhibitor by QSAR modeling, molecular docking, DFT, molecular dynamics, and in vitro confirmation studies.

It is well known that antibiotic resistance is a major health hazard. To eradicate antibiotic-resistant bacterial infections, it is essential to find a novel antibacterial agent. Hence, in this study, a quantitative structure-activity relationship (QSAR) model was developed using 43 DNA gyrase inhibitors, and 700 natural compounds were screened for their antibacterial properties. Based on molecular docking and absorption, distribution, metabolism, excretion, and toxicity (ADMET) studies, the top three leads viz., apigenin-4'-glucoside, 8-deoxygartanin, and cryptodorine were selected and structurally optimized using density functional theory (DFT) studies. The optimized structures were redocked, and molecular dynamic (MD) simulations were performed. Binding energies were calculated by molecular mechanics/Poisson-Boltzmann surface area solvation (MM-PBSA). Based on the above studies, apigenin-4'-glucoside was identified as a potent antibacterial lead. Further in vitro confirmation studies were performed using the plant Lawsonia inermis containing apigenin-4'-glucoside to confirm the antibacterial activity. For QSAR modeling, 2D descriptors were calculated by PaDEL-Descriptors v2.21 software, and the model was developed using the DTClab QSAR tool. Docking was performed using PyRx v0.8 software. ORCA v5.0.1 computational package was used to optimize the structures. The job type used in optimization was equilibrium structure search using the DFT hybrid functional ORCA method B3LYP. The basis set was 6-311G (3df, 3pd) plus four polarization functions for all atoms. Accurate docking was performed for optimized leads using the iGEMDOCK v2.1 tool with a genetic algorithm by 10 solutions each of 80 generations. Molecular dynamic simulations were performed using GROMACS 2020.04 software with CHARMM36 all-atom force field.

Analysing the effect caused by increasing the molecular volume in M1-AChR receptor agonists and antagonists: a structural and computational study.

M1 muscarinic acetylcholine receptor (M1-AChR), a member of the G protein-coupled receptors (GPCR) family, plays a crucial role in learning and memory, making it an important drug target for Alzheimer's disease (AD) and schizophrenia. M1-AChR activation and deactivation have shown modifying effects in AD and PD preclinical models, respectively. However, understanding the pharmacology associated with M1-AChR activation or deactivation is complex, because of the low selectivity among muscarinic subtypes, hampering their therapeutic applications. In this regard, we constructed two quantitative structure-activity relationship (QSAR) models, one for M1-AChR agonists (total and partial), and the other for the antagonists. The binding mode of 59 structurally different compounds, including agonists and antagonists with experimental binding affinity values (pKi), were analyzed employing computational molecular docking over different structures of M1-AChR. Furthermore, we considered the interaction energy (Einter), the number of rotatable bonds (NRB), and lipophilicity (ilogP) for the construction of the QSAR model for agonists (R2 = 89.64, QLMO2 = 78, and Qext2 = 79.1). For the QSAR model of antagonists (R2 = 88.44, QLMO2 = 82, and Qext2 = 78.1) we considered the Einter, the fraction of sp3 carbons fCsp3, and lipophilicity (MlogP). Our results suggest that the ligand volume is a determinant to establish its biological activity (agonist or antagonist), causing changes in binding energy, and determining the affinity for M1-AChR.

QSAR of SARS-CoV-2 Main Protease Inhibitors Utilizing Theoretical Molecular Descriptors

Background: COVID-19 is caused by a novel strain of severe acute respiratory syndrome coronaviruses (SARS-CoV-2). It has claimed casualties around the world since the end of 2019 due to its high virulence and quick multiplicity in the human body. Hence, there has been a requirement to develop effective remedial measures to mitigate the mortality. Scientists have been able to develop corona vaccines to provide immunity, but there are no specific small-molecule chemotherapeutics to combat the novel coronavirus which has spread to the whole world due to its contagiousness. In the viral genome exploration, it has been found that the main protease, also known as chymotrypsin-like cysteine protease ([Mpro] or 3C-like protease [3CLpro]) is responsible for the novel coronavirus replication, transcription, and host immunity destruction. Objectives: Therefore, the main protease has been selected as one of the major targets for the design of new inhibitors. The protein crystallographic and molecular docking studies on SARS-CoV-2 Mpro inhibitors and some quantitative structure-activity relationship (QSAR) studies have been carried out on SARSCoV main protease inhibitors to get some lead molecules for SARS-CoV-2 inhibition. However, there is hardly any QSAR done on the diverse data of SARS-CoV-2 main protease inhibitors. In view of it, QSAR studies have been attempted on SARS-CoV-2 Mpro inhibitors utilizing theoretical molecular descriptors solely computed from the structures of novel corona viral main protease inhibitors. Methods: Methods: As the number of structural descriptors is more than the observations, a genetic algorithm coupled with multiple linear methods has been applied for the development of QSAR models taking diverse SARS-CoV-2 Mpro inhibitors. Results: The developed best QSAR model showing R2 , Q2 Loo, and R2 pred values of 0.7389, 0.6666, and 0.6453 respectively has been further validated on an external data set where a good correlation (r = 0.787) has been found. Conclusion: Therefore, this model may be useful for the design of new SARS-CoV-2 main protease inhibitors.

Evaluation of polymer-preservative interactions for preservation efficacy: molecular dynamics simulation and QSAR approaches.

Preservatives are critical ingredients in various pharmaceutical and consumer products. In particular, a high efficacy preservative system is essential in enhancing the shelf-life and safety of these products. However, the development of such a preservative system heavily relies on experimental approaches. In this study, molecular dynamics (MD) simulation was complemented with quantitative structure-activity relationship (QSAR) modelling to comprehensively evaluate polymer-preservative interactions between three different polymers (polyethylene terephthalate, PET; polypropylene, PP; and cellulose) and a series of preservatives from the classes of aliphatic, aromatic, and organic acids. First, adsorption of preservatives onto polymer surfaces was simulated in an aqueous environment. The preservatives did not adhere to hydrophilic cellulose, but most preservatives were adsorbed by PET and PP in distinct configurations. Interaction energies (IEs) between the preservatives and the polymers generally increase from cellulose to PP and PET. The diffusion coefficients of preservatives are dependent on polymer nature, preservative structure, and their resulting molecular interactions. Linear and low molecular weight preservatives exhibit higher diffusion coefficients in polymers. For a particular preservative, diffusion coefficients increased in the order of cellulose < PET < PP. Finally, using MD properties and molecular descriptors of preservatives, QSAR models were developed to identify key descriptors of preservatives and predict their IEs and diffusion coefficients in polymers. This study demonstrates a computational approach for identifying critical materials properties, and predicting polymer-preservative molecular interactions in water. Such an approach streamlines the rational selection and design of high efficacy preservative systems for various pharmaceutical, food and cosmetic products. Furthermore, the integrated computational strategy also reduces trial-and-error experimental efforts, thereby accelerating the development of high efficacy preservative systems.

2D-QSAR and molecular docking study on nitrofuran analogues as antitubercular agents

&lt;abstract&gt;&lt;sec&gt; &lt;title&gt;Background&lt;/title&gt; &lt;p&gt;Resistance to most of the antitubercular drugs has been on rising trends due to the misuse of existing drugs. This has encouraged us to explore a novel scaffold that has the potential for quick antimicrobial action with minimum side effects. Nitrofurans have attracted us due to their extensive biological activities, such as antibacterial and antifungal activities.&lt;/p&gt; &lt;/sec&gt;&lt;sec&gt; &lt;title&gt;Objective&lt;/title&gt; &lt;p&gt;The antitubercular activities of 126 nitrofuran derivatives have been investigated by using indicator parameters and topological and structural fragment descriptors.&lt;/p&gt; &lt;/sec&gt;&lt;sec&gt; &lt;title&gt;Methods&lt;/title&gt; &lt;p&gt;The different quantitative structure activity relationship (QSAR) models have been created and validated by using two different methodologies: combinatorial protocol in multiple linear regression (CP-MLR) and partial least-squares (PLS) analysis.&lt;/p&gt; &lt;/sec&gt;&lt;sec&gt; &lt;title&gt;Results&lt;/title&gt; &lt;p&gt;The 16 descriptors identified in CP-MLR are from six different classes: Constitutional, Functional, Atom Centered Fragments, Topological, Galvez, and 2D autocorrelation. Indicator parameters and Dragon descriptors suggested that the presence of a furan ring substituted by nitro group is essential for antitubercular activity. Further descriptors from constitutional, and functional classes suggest that the number of double bonds, number of sulphur atoms and number of fragments like thiazole, morpholine and thiophene should be minimum, along with the positive influence of Kier-Hall electrotopological states (Ss) for improved activity. The ACF class descriptors, GALVEZ class descriptors, and 2D-AUTO descriptor GATS4p have also shown positive influence on the antitubercular activity. The TOPO class descriptor T(O…S) suggests that the minimum gap between sulphur and oxygen is favorable for activity.&lt;/p&gt; &lt;/sec&gt;&lt;sec&gt; &lt;title&gt;Conclusions&lt;/title&gt; &lt;p&gt;The models acknowledged in the study have explained the variance between 72 to 76% in the training set and in the prediction of the test set compounds. Also, compounds &lt;bold&gt;122&lt;/bold&gt;, &lt;bold&gt;123&lt;/bold&gt; and &lt;bold&gt;82&lt;/bold&gt; were found to possess good binding affinity towards nitroreductase.&lt;/p&gt; &lt;/sec&gt;&lt;/abstract&gt;