Multiple QSAR and molecular modelling for identification of potent human adenovirus inhibitors

This study has investigated docking-based 2D- and 3D-quantitative structure-activity relationships (QSARs) for a range of 53 hydroxybenzamide analogues as anti- Human adenoviruses (HAdVs). The best 3D-QSAR (Schrodinger, LLC, NY, 2020) and 2D-QSAR models were obtained for the training set and were found to be statistically significant, with cross-validated coefficients (q²) of 0.6775 and 0.7875, and coefficients of determination (r²) of 0.8106 and 0.8122, respectively. Our in-silico docking and virtual screening studies revealed significant higher binding affinity of dataset molecule 34 (-141.444 ​kcal/mol) and hit ZINC01088642 (-114.357 ​kcal/mol) with 4PIE protein than the standard drugs. In in-silico ADME/toxicity studies, molecule 34 and proposed hit ZINC01088642 were found safe with good intestinal absorption, aqueous solubility, medium blood–brain barrier (BBB), no eye corrosion, no skin irritancy, and non-mutagenic profiles. Molecular dynamics analysis showed good stability of complex, hit ZINC01088642 with protein, 4PIE over the simulation period of 20 ns. We believe that further experimental, as well as in-vitro investigation, will shed more lights on the identification of ZINC01088642 as a potential human adenovirus agent.     

Computational Approaches for the Design of Novel Anticancer Compounds Based on Pyrazolo[3,4-d]pyrimidine Derivatives as TRAP1 Inhibitor

In the present in-silico study, various computational techniques were applied to determine potent compounds against TRAP1 kinase. The pharmacophore hypothesis DHHRR_1 consists of important features required for activity. The 3D QSAR study showed a statistically significant model with R² = 0.96 and Q² = 0.57. Leave one out (LOO) cross-validation (R² CV = 0.58) was used to validate the QSAR model. The molecular docking study showed maximum XP docking scores (−11.265, −10.532, −10.422, −10.827, −10.753 kcal/mol) for potent pyrazole analogs (42, 46, 49, 56, 43), respectively, with significant interactions with amino acid residues (ASP 594, CYS 532, PHE 583, SER 536) against TRAP1 kinase receptors (PDB ID: 5Y3N). Furthermore, the docking results were validated using the 100 ns MD simulations performed for the selected five docked complexes. The selected inhibitors showed relatively higher binding affinities than the TRAP1 inhibitor molecules present in the literature. The ZINC database was used for a virtual screening study that screened ZINC05297837, ZINC05434822, and ZINC72286418, which showed similar binding interactions to those shown by potent ligands. Absorption, distribution, metabolism, and excretion (ADME) analysis showed noticeable results. The results of the study may be helpful for the further development of potent TRAP1 inhibitors     

Homology modeling and molecular dynamic simulation of UDP-N-acetylmuramoyl-l-alanine-d-glutamate ligase (MurD) from Mycobacterium tuberculosis H37Rv using in silico approach

The present study aimed to identify the prospective inhibitors of MurD, a cytoplasmic enzyme that catalyzes the addition of d-glutamate to the UDP-N-acetylmuramoyl-l-alanine nucleotide precursor in Mycobacterium tuberculosis (MTB), using virtual screening, docking studies, pharmacokinetic analysis, Molecular Dynamic (MD) simulation, and Molecular Mechanics Generalized Born and Surface Area (MM-GBSA) analyses. The three dimensional (3D) structure was determined based on the homology technique using a template from Streptococcus agalactiae. The modeled structure had three binding sites, namely; substrate binding site (Val18, Thr19, Asp39, Asp40, Gly75, Asn147, Gln171 and His192), the ATP binding site (Gly123, Lys124, Thr125, Thr126, Glu166, Asp283, and Arg314) and the glutamic acid binding site (Arg382, Ser463, and Tyr470). These residues mentioned above play a critical role in the catalytic activity of the enzyme, and their inhibition could serve as a stumbling block to the normal function of the enzyme. A total of 10,344 obtained from virtual screened of Zinc and PubChem databases. These compounds further screened for Lipinski rule of five, docking studies and pharmacokinetic analysis. Four compounds with good binding energies (ZINC11881196 = −10.33 kcal/mol, ZINC12247644 = −8.90 kcal/mol, ZINC14995379 =−8.42 kcal/mol, and PubChem6185 = −8.20 kcal/mol), better than the binding energies of the ATP (−2.31 kcal/mol) and the ligand with known IC₅₀, Aminothiazole (−7.11 kcal/mol) were selected for the MD simulation and MM-GBSA analyses. The result of the analyses showed that all the four ligands formed a stable complex and had the binding free energies better than the binding energy of ATP. Therefore, these ligands considered as suitable prospective inhibitors of the MurD after experimental validation.     

In vitro and in silico interaction of faba bean (Vicia faba L.) seed extract with xanthine oxidase and evaluation of antioxidant activity as a therapeutic potential

Acetone extract of faba bean (Vicia faba L.) was found to be highest total phenol and flavonoid content among all extracts. Antioxidant activity for inhibition percentage (free radical scavenging activity) had 86.47% for acetone extract, and 97.36% for ascorbic acid respectively. IC₅₀ value of ascorbic acid and acetone extact were found to be 9 μg/mL ± 0.20 and 30 μg/mL ± 0.21. Faba bean seeds had catechin, epicatechin, gallic acid and ellagic acid which on molecular docking study revealed that it binds effectively with xanthine oxidase by binding energy of –7.78, –6.11, –6.39, –5.78 kcal/mol respectively compared to allopurinol drug having binding energy of –4.94 kcal/mol. Gallic acid, ellagic acid, catechin, epicatechin (polyphenols) and allopurinol bind other than catalytic residues (Glu-1261) of xanthine oxidase. In vitro and in silico analysis recommended that mode of enzyme inhibition was mixed type.     

Identification of novel Plasmodium falciparum PI4KB inhibitors as potential anti-malarial drugs: Homology modeling,molecular docking and molecular dynamics simulations

The current study was set to discover selective Plasmodium falciparum phosphatidylinositol-4-OH kinase type III beta (pfPI4KB) inhibitors as potential antimalarial agents using combined structure-based and ligand-based drug discovery approach. A comparative model of pfPI4KB was first constructed and validated using molecular docking techniques. Performance of Autodock4.2 and Vina4 software in predicting the inhibitor-PI4KB binding mode and energy was assessed based on two Test Sets: Test Set I contained five ligands with resolved crystal structures with PI4KB, while Test Set II considered eleven compounds with known IC50 value towards PI4KB. The outperformance of Autodock as compared to Vina was reported, giving a correlation coefficient (R²) value of 0.87 and 0.90 for Test Set I and Test Set II, respectively. Pharmacophore-based screening was then conducted to identify drug-like molecules from ZINC database with physicochemical similarity to two potent pfPI4KB inhibitors –namely cpa and cpb. For each query inhibitor, the best 1000 hits in terms of TanimotoCombo scores were selected and subjected to molecular docking and molecular dynamics (MD) calculations. Binding energy was then estimated using molecular mechanics–generalized Born surface area (MM-GBSA) approach over 50 ns MD simulations of the inhibitor-pfPI4KB complexes. According to the calculated MM-GBSA binding energies, ZINC78988474 and ZINC20564116 were identified as potent pfPI4KB inhibitors with binding energies better than those of cpa and cpb, with ΔG_binding ≥ −34.56 kcal/mol. The inhibitor-pfPI4KB interaction and stability were examined over 50 ns MD simulation; as well the selectivity of the identified inhibitors towards pfPI4KB over PI4KB was reported.     

Theoretical study of the hydrogen chemisorption on a ZnO surface

Reactions of a hydrogen molecule with a ZnO surface are studied by an ab initio method. For simulating the ZnO (10 1 0) surface, one ZnO molecule both with and without a Madelung potential is used. Since the electrostatic potential due to the ionic layer decreases exponentially, the effect of the layers deeper than the second one can be neglected. The Madelung potential is, therefore, expressed by the 32 point charges of ±0.5 situated on the first and second layers. Several low-lying states of ZnO and the ZnO + H₂ system have been calculated by the symmetry-adapted cluster (SAC ) and SAC –CI methods. It is found that the ¹Σ⁺ state of ZnO is the ground state and catalytic active and the other states are inactive. ZnO (¹Σ⁺) reacts with H₂ and dissociatively adsorbs it with making Zn? H and O? H bonds. This occurs both with and without the Madelung potential. Without the Madelung potential, the heat of reaction is 81.3 kcal/mol and the reaction barrier is 14.0 kcal/mol. With the Madelung potential, the heat of reaction decreases to 73.5 kcal/mol and the barrier decreases to 11.5 kcal/mol. The mechanism of this reaction is the electron donation from the 2pπ orbital of O to the antibonding σ_u MO of H₂ and the back-donation from the bonding σ_g MO of H₂ to the LUMO of ZnO. In the intermediate stage of the reaction, the dipole of ZnO works to increase the overlap of the active MOS to make the reaction easier. Throughout the reaction, the in-plane 2pπ orbital of O and the HOMO of ZnO are inactive and work to keep the ZnO bond stable during the catalytic process.     

Combined molecular dynamics and continuum solvent approaches (MM-PBSA/GBSA) to predict noscapinoid binding to γ-tubulin dimer

γ-tubulin plays crucial role in the nucleation and organization of microtubules during cell division. Recent studies have also indicated its role in the regulation of microtubule dynamics at the plus end of the microtubules. Moreover, γ-tubulin has been found to be over-expressed in many cancer types, such as carcinomas of the breast and glioblastoma multiforme. These studies have led to immense interest in the identification of chemical leads that might interact with γ-tubulin and disrupt its function in order to explore γ-tubulin as potential chemotherapeutic target. Recently a colchicine-interacting cavity was identified at the interface of γ-tubulin dimer that might also interact with other similar compounds. In the same direction we theoretically investigated binding of a class of compounds, noscapinoids (noscapine and its derivatives) at the interface of the γ-tubulin dimer. Molecular interaction of noscapine and two of its derivatives, amino-noscapine and bromo-noscapine, was investigated by molecular docking, molecular dynamics simulation and binding free energy calculation. All noscapinoids displayed stable interaction throughout simulation of 25 ns. The predictive binding free energy (ΔG_bind) indicates that noscapinoids bind strongly with the γ-tubulin dimer. However, bromo-noscapine showed the best binding affinity (ΔG_bind = –37.6 kcal/mol) followed by noscapine (ΔG_bind = –29.85 kcal/mol) and amino-noscapine (ΔG_bind = –23.99 kcal/mol) using the MM-PBSA method. Similarly using the MM-GBSA method, bromo-noscapine showed highest binding affinity (ΔG_bind = –43.64 kcal/mol) followed by amino-noscapine (ΔG_bind = –37.56 kcal/mol) and noscapine (ΔG_bind = –34.57 kcal/mol). The results thus generate compelling evidence that these noscapinoids may hold great potential for preclinical and clinical evaluation.     

Structure based virtual screening of the Ebola virus trimeric glycoprotein using consensus scoring

Ebola virus (EBOV) causes zoonotic viral infection with a potential risk of global spread and a highly fatal effect on humans. Till date, no drug has gotten market approval for the treatment of Ebola virus disease (EVD), and this perhaps allows the use of both experimental and computational approaches in the antiviral drug discovery process. The main target of potential vaccines that are recently undergoing clinical trials is trimeric glycoprotein (GP) of the EBOV and its exact crystal structure was used in this structure based virtual screening study, with the aid of consensus scoring to select three possible hit compounds from about 36 million compounds in MCULE’s database. Amongst these three compounds, (5R)-5-[[5-(4-chlorophenyl)-1,2,4-oxadiazol-3-yl]methyl]-N-[(4-methoxyphenyl)methyl]-4,5-dihydroisoxazole-3-carboxamide (SC-2, C₂₁H₁₉ClN₄O₄) showed good features with respect to drug likeness, ligand efficiency metrics, solubility, absorption and distribution properties and non-carcinogenicity to emerge as the most promising compound that can be optimized to lead compound against the GP EBOV. The binding mode showed that SC-2 is well embedded within the trimeric chains of the GP EBOV with molecular interactions with some amino acids. The SC-2 hit compound, upon its optimization to lead, might be a good potential candidate with efficacy against the EBOV pathogen and subsequently receive necessary approval to be used as antiviral drug for the treatment of EVD.     

Computational insights into factor affecting the potency of diaryl sulfone analogs as Escherichia coli dihydropteroate synthase inhibitors

Dihydropteroate synthase (DHPS) is an alluring target for designing novel drug candidates to prevent infections caused by pathogenic Escherichia coli strains. Diaryl Sulfone (SO) compounds are found to inhibit DHPS competitively with respect to the substrate pABA (p-aminobenzoate). The extra aromatic ring of diaryl sulfone compounds found to stabilize them in highly flexible pABA binding loops. In this present study, a statistically significant 3D-QSAR model was developed using a data set of diaryl sulfone compounds. The favourable and unfavourable contributions of substitutions in sulfone compounds were illustrated by contour plot obtained from the developed 3D-QSAR model. Molecular docking calculations were performed to investigate the putative binding mode of diaryl sulfone compounds at the catalytic pocket. DFT calculations were carried out using SCF approach, B3LYP- 6-31 G (d) basis set to compute the HOMO, LUMO energies and their respective location at pABA binding pocket. Further, the developed model was validated by FEP (Free Energy Perturbation) calculations. The calculated relative free energy of binding between the highly potent and less potent sulfone compound was found to be −3.78 kcal/ mol which is comparable to the experimental value of −5.85 kcal/mol. A 10 ns molecular dynamics simulation of inhibitor and DHPS confirmed its stability at pABA catalytic site. Outcomes of the present work provide deeper insight in designing novel drug candidates for pathogenic Escherichia coli strains.     

A selectivity study of benzenesulfonamide derivatives on human carbonic anhydrase II/IX by 3D-QSAR,Molecular Docking and Molecular Dynamics Simulation

Nowadays, different approaches have been pursued with the intent to develop sulfonamide-like carbonic anhydrase inhibitors that possess better selectivity profiles toward the different human isoforms of the enzyme. Here, we used conventional 3D-QSAR methods, including comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), and Topomer CoMFA, to construct three-dimensional quantitative structure-activity relationship (3D-QSAR) models for benzenesulfonamide derivatives as human carbonic anhydrase (hCA) II/IX inhibitors. The theoretical models had good reliability (R²>0.75) and predictability (Q²>0.55), and the contour maps could graphically present the contributions of the force fields for activity and identify the structural divergence between human carbonic anhydrase II inhibitors and human carbonic anhydrase IX inhibitors. Consequently, we explored the selectivity of inhibitor for human carbonic anhydrase II and IX through molecular docking, and the difference of activity coincides with the potential binding mode well. According to the results of the predicted values and the molecule docking, we found that the inhibitors published in the literature had stronger inhibition on the hCA IX; based on the theoretical models, we designed seven new compounds with good potential activity and reasonably good ADMET profile, which could selectively inhibit hCA IX. Molecular Dynamics Simulation showed that newly-designed compound D7 had good selectivity on hCA IX. The findings from 3D-QSAR and docking studies maybe helpful in the rational drug design of isoform-selective inhibitors.     

Exploring facile synthesis and cholinesterase inhibiting potential of heteroaryl substituted imidazole derivatives for the treatment of Alzheimer’s disease

Alzheimer’s disease (AD) is a neurodegenerative disorder and cholinesterase (ChE) enzymes are considered as crucial targets for the treatment of AD. Herein, a series of heteroaryl substituted imidazole derivatives (5a-5x) was prepared using amino acid catalyzed, one-pot facile synthetic approach. In this context, the catalytic potentials of different amino acids were investigated and 15 mol% of glutamic acid was identified as the most suitable catalyst to obtain the target products in good yields up to 90 %. These structurally exciting heterocyclic hybrids were screened against acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) enzymes. This series displayed moderate to excellent inhibitory potential against AChE with IC₅₀ values > 25 µM and the most active compound was 3-(4-(1-(3,5-dimethylphenyl)-4,5-diphenyl-1H-imidazol-2-yl)-1-phenyl-1H-pyrazol-3-yl)–2H-chromen-2-one (5x) with IC₅₀ value of 25.83 ± 0.25 µM.This inhibitory potential was attributed to hydrophobicity as the major contributory factor. The most potent compound against BChE was 1,3-diphenyl-4-(1,4,5-triphenyl-1H-imidazol-2-yl)-1H-pyrazole (5a) with IC₅₀ value of 0.35 ± 0.02 µM followed by other potent compounds 5p, 5 m, 5x, 5b, 5c, 5e and 5f with IC₅₀ values < 10 µM. SAR studies further revealed that coumarinyl moiety at R¹ position in the imidazolylpyrazole skeleton significantly improved the overall cholinesterase inhibitory potential. However, a simple phenyl ring attached at this R¹ site was highly effective and selective for BChE inhibition (5a) over AChE. Docking data also demonstrated the interaction of 5x and AChE with a docking score of 7564 and atomic contact energy (ACE) value of –291.90 kcal/mol whereas docking score for 5a against BChE was 7096 with ACE value of –332.95 kcal/mol. The results altogether suggest further investigations of the heteroaryl substituted imidazole core skeleton in search of potential leads towards designing of new anti-cholinesterase drugs for the treatment of AD.     

Pharmacophore generation,atom-based 3D-QSAR,HQSAR and activity cliff analyses of benzothiazine and deazaxanthine derivatives as dual A2A antagonists/MAO‑B inhibitors

Dual inhibition of A_2A and MAO-B is an emerging strategy in neurodegenerative diseases, such as Alzheimer’s disease (AD) and Parkinson’s disease (PD). In this study, atom-based three-dimensional quantitative structure–activity relationship (3D-QSAR) and hologram quantitative structure–activity relationship (HQSAR) models were generated with benzothiazine and deazaxanthine derivatives. Based on activity against A_2A and MAO-B, two statistically signi?cant 3D-QSAR models (r² = 0.96, q² = 0.76 and r² = 0.91, q² = 0.63) and HQSAR models (r² = 0.93, q² = 0.68 and r² = 0.97, q² = 0.58) were developed. In an activity cliff analysis, structural outliers were identified by calculating the Mahalanobis distance for a pair of compounds with A_2A and MAO-B inhibitory activities. The generated 3D-QSAR and HQSAR models, activity cliff analysis, molecular docking and dynamic studies for dual target protein inhibitors provide key structural scaffolds that serve as building blocks in designing drug-like molecules for neurodegenerative diseases.     

Molecular modeling studies of human immunodeficiency virus type 1 protease inhibitors using three‐dimensional quantitative structure‐activity relationship,virtual screening,and docking simulations

In this paper, two 3‐dimensional quantitative structure‐activity relationship models for 60 human immunodeficiency virus (HIV)‐1 protease inhibitors were established using random sampling analysis on molecular surface and translocation comparative molecular field vector analysis (Topomer CoMFA). The non–cross‐validation (r²), cross‐validation (q²), correlation coefficient of external validation (Q²_ext), and F of 2 models were 0.94, 0.80, 0.79, and 198.84 and 0.94, 0.72, 0.75, and 208.53, respectively. The results indicated that 2 models were reasonable and had good prediction ability. Topomer Search was used to search R groups in the ZINC database, 20 new compounds were designed, and the Topomer CoMFA model was used to predicate the biological activity. The results showed that 18 new compounds were more active than the template molecule. So the Topomer Search is effective in screening and can guide the design of new HIV/AIDS drugs. The mechanism of action was studied by molecular docking, and it showed that the protease inhibitors and Ile50, Asp25, and Arg8 sites of HIV‐1 protease have interactions. These results have provided an insight for the design of new potent inhibitors of HIV‐1 protease.     

Computer-aided molecular design of (E)-N-Aryl-2-ethene-sulfonamide analogues as microtubule targeted agents in prostate cancer

Microtubules are tube-shaped, filamentous and cytoskeletal proteins that are essential in all eukaryotic cells. Microtubule is an attractive and promising target for anticancer agents. In this study, three-dimensional quantitative structure activity relationships (3D-QSAR) including comparative molecular field analysis, CoMFA, and comparative molecular similarity indices analysis, CoMSIA, were performed on a set of 45 (E)-N-Aryl-2-ethene-sulfonamide analogues as microtubule-targeted anti-prostate cancer agents. Automated grid potential analysis, AutoGPA module in Molecular Operating Environment 2009.10 (MOE) as a new 3D-QSAR approach with the pharmacophore-based alignment was carried out on the same dataset. AutoGPA-based 3D-QSAR model yielded better prediction parameters than CoMFA and CoMSIA. Based on the contour maps generated from the models, some key features were identified in (E)-N-Aryl-2-arylethene-sulfonamide analogues that were responsible for the anti-cancer activity. Virtual screening was performed based on pharmacophore modeling and molecular docking to identify the new inhibitors from ZINC database. Seven top ranked compounds were found based on Gold score fitness function. In silico ADMET studies were performed on compounds retrieved from virtual screening in compliance with the standard ranges.     

The 3D Structure of the Binding Pocket of the Human Oxytocin Receptor for Benzoxazine Antagonists,Determined by Molecular Docking,Scoring Functions and 3D-QSAR Methods

Summary Molecular docking and 3D-QSAR studies were performed to determine the binding mode for a series of benzoxazine oxytocin antagonists taken from the literature. Structural hypotheses were generated by docking the most active molecule to the rigid receptor by means of AutoDock 3.05. The cluster analysis yielded seven possible binding conformations. These structures were refined by using constrained simulated annealing, and the further ligands were aligned in the refined receptor by molecular docking. A good correlation was found between the estimated ΔG_bind and the pK_i values for complex F. The Connolly-surface analysis, CoMFA and CoMSIA models q_CoMFA² = 0.653, q_CoMSA² = 0.630 and r_pred,CoMFA² = 0.852 , r_pred,CoMSIA² = 0.815) confirmed the scoring function results. The structural features of the receptor–ligand complex and the CoMFA and CoMSIA fields are in closely connected. These results suggest that receptor–ligand complex F is the most likely binding hypothesis for the studied benzoxazine analogs.     

Interaction of the prototypical α-ketoamide inhibitor with the SARS-CoV-2 main protease active site in silico: Molecular dynamic simulations highlight the stability of the ligand-protein complex

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes an illness known as COVID-19, which has been declared a global pandemic with over 2 million confirmed cases and 137,000 deaths in 185 countries and regions at the time of writing (16 April 2020), over a quarter of these cases being in the United States. In the absence of a vaccine, or an approved effective therapeutic, there is an intense interest in repositioning available drugs or designing small molecule antivirals. In this context, in silico modelling has proven to be an invaluable tool. An important target is the SARS-CoV-2 main protease (M^pro), involved in processing translated viral proteins. Peptidomimetic α-ketoamides represent prototypical inhibitors of M^pro. A recent attempt at designing a compound with enhanced pharmacokinetic properties has resulted in the synthesis and evaluation of the α-ketoamide 13b analogue. Here, we performed molecular docking and molecular dynamics simulations to further characterize the interaction of α-ketoamide 13b with the active site of the SARS-CoV-2 M^pro. We included the widely used antibiotic, amoxicillin, for comparison. Our findings indicate that α-ketoamide 13b binds more tightly (predicted GlideScore = -8.7 and -9.2 kcal/mol for protomers A and B, respectively), to the protease active site compared to amoxicillin (-5.0 and -4.8 kcal/mol). Further, molecular dynamics simulations highlight the stability of the interaction of the α-ketoamide 13b ligand with the SARS-CoV-2 M^pro (ΔG = -25.2 and -22.3 kcal/mol for protomers A and B). In contrast, amoxicillin interacts unfavourably with the protease (ΔG = +32.8 kcal/mol for protomer A), with unbinding events observed in several independent simulations. Overall, our findings are consistent with those previously observed, and highlight the need to further explore the α-ketoamides as potential antivirals for this ongoing COVID-19 pandemic.     

Imidazo[1,2-a]pyrazine inhibitors of phosphoinositide 3-kinase alpha (PI3Kα): 3D-QSAR analysis utilizing the Hybrid Monte Carlo algorithm to refine receptor-ligand complexes for molecular alignment

Phosphoinositide 3-kinase alpha (PI3Kα) is a lipid kinase involved in several cellular functions such as cell growth, proliferation, differentiation and survival, and its anomalous regulation leads to cancerous conditions. PI3Kα inhibition completely blocks the cancer signalling pathway, hence it can be explored as an important therapeutic target for cancer treatment. In the present study, docking analysis of 49 selective imidazo[1,2-a]pyrazine inhibitors of PI3Kα was carried out using the QM-Polarized ligand docking (QPLD) program of the Schrödinger software, followed by the refinement of receptor–ligand conformations using the Hybrid Monte Carlo algorithm in the Liaison program, and alignment of refined conformations of inhibitors was utilized for the development of an atom-based 3D-QSAR model in the PHASE program. Among the five generated models, the best model was selected corresponding to PLS factor 2, displaying the highest value of Q²_test (0.650). The selected model also displayed high values of r²_train (0.917), F-value (166.5) and Pearson-r (0.877) and a low value of SD (0.265). The contour plots generated for the selected 3D-QSAR model were correlated with the results of docking simulations. Finally, this combined information generated from 3D-QSAR and docking analysis was used to design new congeners.     

Drug promiscuity: Exploring the polypharmacology potential of 1, 3, 6-trisubstituted 1, 4-diazepane-7-ones as an inhibitor of the ‘god father’ of immune checkpoint

High production cost, instability, low tumor penetration are some of the shortcomings that have characterized and undermined the use of antibodies as a target for Cytotoxic T-lymphocytes associated protein 4 (CTLA-4). Design and discovery of small molecule inhibitors have therefore become a sine qua non in targeting immune proteins implicated in immune disorders. In this study, we utilized a drug repositioning approach to explore the characteristic feature of unrelated proteins to have similar binding sites and the promiscuity of drugs to repurpose an existing drug to target CTLA-4. CTLA-4 and Kallikrein-7 were found to have similar binding sites, we therefore used 1, 3, 6-trisubstituted 1, 4-diazepane-7-ones (TDSO) which is an inhibitor of Kallikrein-7 as our lead compound. High throughput screening using TDSO as a lead compound resulted in 9 hits with ZINC04515726 and ZINC08985213 having the highest binding score. We went ahead to investigate the interaction of these compounds with CTLA-4 by conducting a molecular dynamic simulation. Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) estimations revealed that TDSO had the highest binding energy value of -28.51Kcal/mol, with ZINC04515726 and ZINC08985213 having -23.76Kcal/mol and -21.03Kcal/mol respectively. The per-residue decomposition highlighted Tyr24, Ala25, Gly28, Ala30, Tyr53 and Asn72 as having significantly high electrostatic energy contributions and the main contributing residues to the binding of TDSO, ZINC04515726 and ZINC08985213 to Cytotoxic T lymphocytes CTLA-4. Summarily, from the results gathered, we proposed that TDSO can be an effective immune check point small molecule inhibitor against the suppression of T-cell activation, proliferation, and tumor cell eradication.     

Strategies to design pyrazolyl urea derivatives for p38 kinase inhibition: a molecular modeling study

The p38 protein kinase is a serine–threonine mitogen activated protein kinase, which plays an important role in inflammation and arthritis. A combined study of 3D-QSAR and molecular docking has been undertaken to explore the structural insights of pyrazolyl urea p38 kinase inhibitors. The 3D-QSAR studies involved comparative molecular field analysis (CoMFA) and comparative molecular similarity indices (CoMSIA). The best CoMFA model was derived from the atom fit alignment with a cross-validated r ² (q ²) value of 0.516 and conventional r ² of 0.950, while the best CoMSIA model yielded a q ² of 0.455 and r ² of 0.979 (39 molecules in training set, 9 molecules in test set). The CoMFA and CoMSIA contour maps generated from these models provided inklings about the influence of interactive molecular fields in the space on the activity. GOLD, Sybyl (FlexX) and AutoDock docking protocols were exercised to explore the protein–inhibitor interactions. The integration of 3D-QSAR and molecular docking has proffered essential structural features of pyrazolyl urea inhibitors and also strategies to design new potent analogues with enhanced activity.