Revealing innovative JAK1 and JAK3 inhibitors: a comprehensive study utilizing QSAR, 3D-Pharmacophore screening, molecular docking, molecular dynamics, and MM/GBSA analyses.

The heterocycle compounds, with their diverse functionalities, are particularly effective in inhibiting Janus kinases (JAKs). Therefore, it is crucial to identify the correlation between their complex structures and biological activities for the development of new drugs for the treatment of rheumatoid arthritis (RA) and cancer. In this study, a diverse set of 28 heterocyclic compounds selective for JAK1 and JAK3 was employed to construct quantitative structure-activity relationship (QSAR) models using multiple linear regression (MLR). Artificial neural network (ANN) models were employed in the development of QSAR models. The robustness and stability of the models were assessed through internal and external methodologies, including the domain of applicability (DoA). The molecular descriptors incorporated into the model exhibited a satisfactory correlation with the receptor-ligand complex structures of JAKs observed in X-ray crystallography, making the model interpretable and predictive. Furthermore, pharmacophore models ADRRR and ADHRR were designed for each JAK1 and JAK3, proving effective in discriminating between active compounds and decoys. Both models demonstrated good performance in identifying new compounds, with an ROC of 0.83 for the ADRRR model and an ROC of 0.75 for the ADHRR model. Using a pharmacophore model, the most promising compounds were selected based on their strong affinity compared to the most active compounds in the studied series each JAK1 and JAK3. Notably, the pharmacokinetic, physicochemical properties, and biological activities of the selected compounds (As compounds ZINC79189223 and ZINC66252348) were found to be consistent with their therapeutic effects in RA, owing to their non-toxic, cholinergic nature, absence of P-glycoprotein, high gastrointestinal absorption, and ability to penetrate the blood-brain barrier. Furthermore, ADMET properties were assessed, and molecular dynamics and MM/GBSA analysis revealed stability in these molecules.

Computational 3D Modeling-Based Identification of Inhibitors Targeting Cysteine Covalent Bond Catalysts for JAK3 and CYP3A4 Enzymes in the Treatment of Rheumatoid Arthritis.

This work aimed to find new inhibitors of the CYP3A4 and JAK3 enzymes, which are significant players in autoimmune diseases such as rheumatoid arthritis. Advanced computer-aided drug design techniques, such as pharmacophore and 3D-QSAR modeling, were used. Two strong 3D-QSAR models were created, and their predictive power was validated by the strong correlation (R2 values > 80%) between the predicted and experimental activity. With an ROC value of 0.9, a pharmacophore model grounded in the DHRRR hypothesis likewise demonstrated strong predictive ability. Eight possible inhibitors were found, and six new inhibitors were designed in silico using these computational models. The pharmacokinetic and safety characteristics of these candidates were thoroughly assessed. The possible interactions between the inhibitors and the target enzymes were made clear via molecular docking. Furthermore, MM/GBSA computations and molecular dynamics simulations offered insightful information about the stability of the binding between inhibitors and CYP3A4 or JAK3. Through the integration of various computational approaches, this study successfully identified potential inhibitor candidates for additional investigation and efficiently screened compounds. The findings contribute to our knowledge of enzyme-inhibitor interactions and may help us create more effective treatments for autoimmune conditions like rheumatoid arthritis.

In silico molecular investigations of pyridine N-Oxide compounds as potential inhibitors of SARS-CoV-2: 3D QSAR, molecular docking modeling, and ADMET screening.

The new coronavirus SARS-CoV-2 virus is causing a severe pneumonia in human, provoking the serious outbreak epidemic CoV-2. Since its appearance in Wuhan, China on December 2019, CoV-2 becomes the biggest challenge the world is facing today, including the discovery of antiviral drug for SARS-CoV-2. In this study, the potential inhibitory of a class of human SARS inhibitors, namely pyridine N-oxide derivatives, against CoV-2 was addressed by quantitative structure-activity relationship 3 D-QSAR. The reliable CoMSIA developed model of 110 pyridine N-oxide based-antiviral compounds, showed Q2= 0.54 and rext2=0.71. The molecular surflex-docking was applied to identify the crystal structure of CoV-2 main protease 3CLpro (PDB: 6LU7) and two potentially and largely used antiviral molecules, namely chloroquine, hydroxychloroquine. The obtained free energy affinity and ADMET properties indicate that among the series of model antiviral compounds examined, the new antiviral compound A5 could be an excellent antiviral drug inhibitor against COVID-19. The inhibition activity of pyridine N-oxyde compounds against CoV-2 was compared with the activity of two common antiviral drug, namely chloroquine (CQ) and hydroxychloroquine (HCQ). DFT method was also used to define the sites of reactivity of pyridine N-oxyde derivatives as well as CQ and HCQ.Communicated by Ramaswamy H. Sarma.

QSAR study of unsymmetrical aromatic disulfides as potent avian SARS-CoV main protease inhibitors using quantum chemical descriptors and statistical methods.

In silico research was executed on forty unsymmetrical aromatic disulfide derivatives as inhibitors of the SARS Coronavirus (SARS-CoV-1). Density functional theory (DFT) calculation with B3LYP functional employing 6-311 â€‹+ â€‹G(d,p) basis set was used to calculate quantum chemical descriptors. Topological, physicochemical and thermodynamic parameters were calculated using ChemOffice software. The dataset was divided randomly into training and test sets consisting of 32 and 8 compounds, respectively. In attempt to explore the structural requirements for bioactives molecules with significant anti-SARS-CoV activity, we have built valid and robust statistics models using QSAR approach. Hundred linear pentavariate and quadrivariate models were established by changing training set compounds and further applied in test set to calculate predicted IC50 values of compounds. Both built models were individually validated internally as well as externally along with Y-Randomization according to the OECD principles for the validation of QSAR model and the model acceptance criteria of Golbraikh and Tropsha's. Model 34 is chosen with higher values of R2, R2 test and Q2cv (R2 â€‹= â€‹0.838, R2 test â€‹= â€‹0.735, Q2 cv â€‹= â€‹0.757). It is very important to notice that anti-SARS-CoV main protease of these compounds appear to be mainly governed by five descriptors, i.e. highest occupied molecular orbital energy (EHOMO), energy of molecular orbital below HOMO energy (EHOMO-1), Balaban index (BI), bond length between the two sulfur atoms (S1S2) and bond length between sulfur atom and benzene ring (S2Bnz). Here the possible action mechanism of these compounds was analyzed and discussed, in particular, important structural requirements for great SARS-CoV main protease inhibitor will be by substituting disulfides with smaller size electron withdrawing groups. Based on the best proposed QSAR model, some new compounds with higher SARS-CoV inhibitors activities have been designed. Further, in silico prediction studies on ADMET pharmacokinetics properties were conducted.

Identification of a novel dual-target scaffold for 3CLpro and RdRp proteins of SARS-CoV-2 using 3D-similarity search, molecular docking, molecular dynamics and ADMET evaluation.

The new SARS-CoV-2 coronavirus is the causative agent of the COVID-19 pandemic outbreak that affected whole the world with more than 6 million confirmed cases and over 370,000 deaths. At present, there are no effective treatments or vaccine for this disease, which constitutes a serious global health crisis. As the pandemic still spreading around the globe, it is of interest to use computational methods to identify potential inhibitors for the virus. The crystallographic structures of 3CLpro (PDB: 6LU7) and RdRp (PDB 6ML7) were used in virtual screening of 50000 chemical compounds obtained from the CAS Antiviral COVID19 database using 3D-similarity search and standard molecular docking followed by ranking and selection of compounds based on their binding affinity, computational techniques for the sake of details on the binding interactions, absorption, distribution, metabolism, excretion, and toxicity prediction; we report three 4-(morpholin-4-yl)-1,3,5-triazin-2-amine derivatives; two compounds (2001083-68-5 and 2001083-69-6) with optimal binding features to the active site of the main protease and one compound (833463-19-7) with optimal binding features to the active site of the polymerase for further consideration to fight COVID-19. The structural stability and dynamics of lead compounds at the active site of 3CLpro and RdRp were examined using molecular dynamics (MD) simulation. Essential dynamics demonstrated that the three complexes remain stable during simulation of 20 ns, which may be suitable candidates for further experimental analysis. As the identified leads share the same scaffold, they may serve as promising leads in the development of dual 3CLpro and RdRp inhibitors against SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Discovery of Potent SARS-CoV-2 Inhibitors from Approved Antiviral Drugs via Docking and Virtual Screening.

BACKGROUND: Coronavirus Disease 2019 (COVID-19) pandemic continues to threaten patients, societies and healthcare systems around the world. There is an urgent need to search for possible medications. OBJECTIVE: This article intends to use virtual screening and molecular docking methods to find potential inhibitors from existing drugs that can respond to COVID-19. METHODS: To take part in the current research investigation and to define a potential target drug that may protect the world from the pandemic of corona disease, a virtual screening study of 129 approved drugs was carried out which showed that their metabolic characteristics, dosages used, potential efficacy and side effects are clear as they have been approved for treating existing infections. Especially 12 drugs against chronic hepatitis B virus, 37 against chronic hepatitis C virus, 37 against human immunodeficiency virus, 14 anti-herpesvirus, 11 anti-influenza, and 18 other drugs currently on the market were considered for this study. These drugs were then evaluated using virtual screening and molecular docking studies on the active site of the (SARS-CoV-2) main protease (6lu7). Once the efficacy of the drug is determined, it can be approved for its in vitro and in vivo activity against the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), which can be beneficial for the rapid clinical treatment of patients. These drugs were considered potentially effective against SARS-CoV-2 and those with high molecular docking scores were proposed as novel candidates for repurposing. The N3 inhibitor cocrystallized with protease (6lu7) and the anti-HIV protease inhibitor Lopinavir were used as standards for comparison. RESULTS: The results suggest the effectiveness of Beclabuvir, Nilotinib, Tirilazad, Trametinib and Glecaprevir as potent drugs against SARS-CoV-2 since they tightly bind to its main protease. CONCLUSION: These promising drugs can inhibit the replication of the virus; hence, the repurposing of these compounds is suggested for the treatment of COVID-19. No toxicity measurements are required for these drugs since they were previously tested prior to their approval by the FDA. However, the assessment of these potential inhibitors as clinical drugs requires further in vivo tests of these drugs.

Discovering Anti-Cancer Drugs via Computational Methods.

New drug discovery has been acknowledged as a complicated, expensive, time-consuming, and challenging project. It has been estimated that around 12 years and 2.7 billion USD, on average, are demanded for a new drug discovery via traditional drug development pipeline. How to reduce the research cost and speed up the development process of new drug discovery has become a challenging, urgent question for the pharmaceutical industry. Computer-aided drug discovery (CADD) has emerged as a powerful, and promising technology for faster, cheaper, and more effective drug design. Recently, the rapid growth of computational tools for drug discovery, including anticancer therapies, has exhibited a significant and outstanding impact on anticancer drug design, and has also provided fruitful insights into the area of cancer therapy. In this work, we discussed the different subareas of the computer-aided drug discovery process with a focus on anticancer drugs.

In Silico Exploration of Aryl Halides Analogues as Checkpoint Kinase 1 Inhibitors by Using 3D QSAR, Molecular Docking Study, and ADMET Screening.

Purpose: In this review, a set of aryl halides analogs were identified as potent checkpoint kinase 1 (Chk1) inhibitors through a series of computer-aided drug design processes, to develop models with good predictive ability, highlight the important interactions between the ligand and the Chk1 receptor protein and determine properties of the new proposed drugs as Chk1 inhibitors agents. Methods: Three-dimensional quantitative structure-activity relationship (3D-QSAR) modeling, molecular docking and absorption, distribution, metabolism, excretion and toxicity (ADMET) approaches are used to determine structure activity relationship and confirm the stable conformation on the receptor pocket. Results: The statistical analysis results of comparative -molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) models that employed for a training set of 24 compounds gives reliable values of Q2 (0.70 and 0.94, respectively) and R2 (0.68 and 0.96, respectively). Conclusion: Computer-aided drug design tools used to develop models that possess good predictive ability, and to determine the stability of the observed and predicted molecules in the receptor pocket, also in silico of pharmacokinetic (ADMET) results shows good properties and bioavailability for these new proposed Chk1 inhibitors agents.

QSAR study of anti-Human African Trypanosomiasis activity for 2-phenylimidazopyridines derivatives using DFT and Lipinski's descriptors.

The quantitative structure-activity relationship (QSAR) of sixty 2-phenylimidazopyridines derivatives with anti-Human African Trypanosomiasis (anti-HAT) activity has been studied by using the density functional theory (DFT) and statistical methods. Becke's three-parameter hybrid method and the Lee-Yang-Parr B3LYP functional employing 6-31G(d) basis set are used to calculate quantum chemical descriptors using Gaussian 03W software, and the five Lipinski's parameters were calculated using ChemOffice software. In order to obtain robust and reliable QSAR model, the original dataset was randomly divided into training and prediction sets comprising 48 and 12 compounds, respectively. An optimal model for the training set with significant statistical quality was established. The same model was further applied to predict pEC50 values of the 12 compounds in the test set, further showing that this QSAR model has high predictive ability. It is very interesting to find that the anti-HAT of these compounds appear to be mainly governed by four factors, i.e., the number of H-bond donors, the lowest unoccupied molecular orbital energy, the molecular weight and the octanol/water partition coefficient. Here the possible action mechanism of these compounds was analysed and discussed, in particular, important structural requirements for great anti-HAT activity will be by increasing molecular size and substitute the 2-phenylimidazopyridines derivatives with polar, ionic, stronger accepting electron ability group and heteroatoms attached to one or more hydrogen atoms. Based on this proposed QSAR model, some new compounds with higher anti-HAT activities have been theoretically designed. Such results can offer useful theoretical references for future experimental works.

3D QSAR studies, molecular docking and ADMET evaluation, using thiazolidine derivatives as template to obtain new inhibitors of PIM1 kinase.

Proviral Integration site for Moloney murine leukemia virus-1 (PIM1) belongs to the serine/threonine kinase family of Ca2+-calmodulin-dependent protein kinase (CAMK) group, which is involved in cell survival and proliferation as well as a number of other signal transduction pathways. Thus, PIM1 is regarded as a promising target for treatment of cancers. In the present paper, a three-dimensional quantitative structure activity relationship (3D-QSAR) and molecular docking were performed to investigate the binding between PIM1 and thiazolidine inhibitors in order to design potent inhibitors. The comparative molecular similarity indices analysis (CoMSIA) was developed using twenty-six molecules having pIC50 ranging from 8.854 to 6.011 (IC50 in nM). The best CoMSIA model gave significant statistical quality. The determination coefficient (R2) and Leave-One-Out cross-validation coefficient (Q2) are 0.85 and 0.58, respectively. Furthermore, the predictive ability of this model was evaluated by external validation((n = 11, R2test = 0.72, and MAE = 0.170 log units). The graphical contour maps could provide structural features to improve inhibitory activity. Furthermore, a good consistency between contour maps and molecular docking strongly demonstrates that the molecular modeling is reliable. Based on these satisfactory results, we designed several new potent PIM1 inhibitors and their inhibitory activities were predicted by the molecular models. Additionally, those newly designed inhibitors, showed promising results in the preliminary in silico ADMET evaluations, compared to the best inhibitor from the studied dataset. The results expand our understanding of thiazolidines as inhibitors of PIM1 and could be of great help in lead optimization for early drug discovery of highly potent inhibitors.

QSAR study and rustic ligand-based virtual screening in a search for aminooxadiazole derivatives as PIM1 inhibitors.

BACKGROUND: Quantitative structure-activity relationship (QSAR) was carried out to study a series of aminooxadiazoles as PIM1 inhibitors having pki ranging from 5.59 to 9.62 (k i in nM). The present study was performed using Genetic Algorithm method of variable selection (GFA), multiple linear regression analysis (MLR) and non-linear multiple regression analysis (MNLR) to build unambiguous QSAR models of 34 substituted aminooxadiazoles toward PIM1 inhibitory activity based on topological descriptors. RESULTS: Results showed that the MLR and MNLR predict activity in a satisfactory manner. We concluded that both models provide a high agreement between the predicted and observed values of PIM1 inhibitory activity. Also, they exhibit good stability towards data variations for the validation methods. Furthermore, based on the similarity principle we performed a database screening to identify putative PIM1 candidates inhibitors, and predict their inhibitory activities using the proposed MLR model. CONCLUSIONS: This approach can be easily handled by chemists, to distinguish, which ones among the future designed aminooxadiazoles structures could be lead-like and those that couldn't be, thus, they can be eliminated in the early stages of drug discovery process.

Computer aided drug design based on 3D-QSAR and molecular docking studies of 5-(1H-indol-5-yl)-1,3,4-thiadiazol-2-amine derivatives as PIM2 inhibitors: a proposal to chemists.

PIM2 kinase plays a crucial role in the cell cycle events including survival, proliferation, and differentiation in normal and neoplastic neuronal cells. Thus, it is regarded as an essential target for cancer pharmaceutical. Design of novel 5-(1H-indol-5-yl)-1,3,4-thiadiazol-2-amine derivatives with enhanced PIM2 inhibitory activity. A series of twenty-five PIM2 inhibitors reported in the literature containing 5-(1H-indol-5-yl)-1,3,4-thiadiazol-2-amines scaffold was studied by using two computational techniques, namely, three-dimensional quantitative structure activity relationship (3D-QSAR) and molecular docking. The comparative molecular field analysis (CoMFA) and comparative molecular similarity indexes analysis (CoMSIA) studies were developed using nineteen molecules having pIC50 ranging from 8.222 to 4.157. The best generated CoMFA and CoMSIA models exhibit conventional determination coefficients R2 of 0.91 and 0.90 as well as the Leave One Out cross-validation determination coefficients Q2 of 0.68 and 0.62, respectively. Moreover, the predictive ability of those models was evaluated by the external validation using a test set of six compounds with predicted determination coefficients Rtest 2 of 0.96 and 0.96, respectively. Besides, y-randomization test was also performed to validate our 3D-QSAR models. The most and the least active compounds were docked into the active site of the protein (PDB ID: 4 × 7q) to confirm those obtained results from 3D-QSAR models and elucidate the binding mode between this kind of compounds and the PIM2 enzyme. These satisfactory results are not offered help only to understand the binding mode of 5-(1H-indol-5-yl)-1,3,4-thiadiazol series compounds into this kind of targets, but provide information to design new potent PIM2 inhibitors.

QSAR studies on PIM1 and PIM2 inhibitors using statistical methods: a rustic strategy to screen for 5-(1H-indol-5-yl)-1,3,4-thiadiazol analogues and predict their PIM inhibitory activity.

BACKGROUND: Quantitative structure activity relationship was carried out to study a series of PIM1 and PIM2 inhibitors. The present study was performed on twenty-five substituted 5-(1H-indol-5-yl)-1,3,4-thiadiazols as PIM1 and PIM2 inhibitors having pIC50 ranging from 5.55 to 9 µM and from 4.66 to 8.22 µM, respectively, using genetic function algorithm for variable selection and multiple linear regression analysis (MLR) to establish unambiguous and simple QSAR models based on topological molecular descriptors. RESULTS: Results showed that the MLR predict activity in a satisfactory manner for both activities. Consequently, the aim of the current study is twofold, first, a simple linear QSAR model was developed, which could be easily handled by chemist to screen chemical databases, or design for new potent PIM1 and PIM2 inhibitors. Second, the outcomes extracted from the current study were exploited to predict the PIM inhibitory activity of some studied compound analogues. CONCLUSIONS: The goal of this study is to develop easy and convenient QSAR model could be handled by everyone to screen chemical databases or to design newly PIM1 and PIM2 inhibitors derived from 5-(1H-indol-5-yl)-1,3,4-thiadiazol. Graphical abstract Flow chart of the methodology used in this work.