Insight into the structural requirements of aminopyrimidine derivatives for good potency against both purified enzyme and whole cells of M. tuberculosis: combination of HQSAR,CoMSIA, and MD simulation studies

The Mycobacterium tuberculosis protein kinase B (PknB) is critical for growth and survival of M. tuberculosis within the host. The series of aminopyrimidine derivatives show impressive activity against PknB (IC₅₀ < .5 μM). However, most of them show weak or no cellular activity against M. tuberculosis (MIC > 63 μM). Consequently, the key structural features related to activity against of both PknB and M. tuberculosis need to be investigated. Here, two- and three-dimensional quantitative structure–activity relationship (2D and 3D QSAR) analyses combined with molecular dynamics (MD) simulations were employed with the aim to evaluate these key structural features of aminopyrimidine derivatives. Hologram quantitative structure–activity relationship (HQSAR) and CoMSIA models constructed from IC₅₀ and MIC values of aminopyrimidine compounds could establish the structural requirements for better activity against of both PknB and M. tuberculosis. The NH linker and the R₁ substituent of the template compound are not only crucial for the biological activity against PknB but also for the biological activity against M. tuberculosis. Moreover, the results obtained from MD simulations show that these moieties are the key fragments for binding of aminopyrimidine compounds in PknB. The combination of QSAR analysis and MD simulations helps us to provide a structural concept that could guide future design of PknB inhibitors with improved potency against both the purified enzyme and whole M. tuberculosis cells.     

Exploration of potential RSK2 inhibitors by pharmacophore modelling,structure-based 3D-QSAR,molecular docking study and molecular dynamics simulation

Abstract

The p90 ribosomal s6 kinase 2 (RSK2) is a promising target because of its over expression and activation in human cancer cells and tissues. Over the last few years, significant efforts have been made in order to develop RSK2 inhibitors to treat myeloma, prostatic cancer, skin cancer and etc., but with limited success so far. In this paper, pharmacophore modelling, molecular docking study and molecular dynamics (MD) simulation have been performed to explore the novel inhibitors of RSK2. Pharmacophore models were developed by 95 molecules having pIC₅₀ ranging from 4.577 to 9.000. The pharmacophore model includes one hydrogen bond acceptor (A), one hydrogen bond donor (D), one hydrophobic feature (H) and one aromatic ring (R). It is the best pharmacophore hypothesis that has the highest correlation coefficient (R² = 0.91) and cross validation coefficient (Q² = 0.71) at 5 component PLS factor. It was evaluated using enrichment analysis and the best model was used for virtual screening. The constraints used in this study were docking score, ADME properties, binding free energy estimates and IFD Score to screen the database. Ultimately, 12 hits were identified as potent and novel RSK2 inhibitors. A 15 ns molecular dynamics (MD) simulation was further employed to validate the reliability of the docking results.     

Analysis of dissociation process for gas hydrates by molecular dynamics simulation

The dissociation processes of methane and carbon dioxide hydrates were investigated by molecular dynamics simulation. The simulations were performed with 368 water molecules and 64 gas molecules using NPT ensembles. The TraPPE (single-site) and 5-site models were adopted for methane molecules. The EPM2 (3-site) and SPC/E models were used for carbon dioxide and water molecules, respectively. The simulations were carried out at 270 K and 5.0 MPa for hydrate stabilisation. Then, temperature was increased up to 370 K. The temperature increasing rates were 0.1–20 TK/s. The gas hydrates dissociated during increasing temperature or at 370 K. The potential models of methane molecule did not much influence the dissociation process of methane hydrate. The mechanisms of dissociation process were analysed with the coordination numbers and mean square displacements. It was found that the water cages break down first, then the gas molecules escape from the water cages. The methane hydrate was more stable than the carbon dioxide hydrate at the calculated conditions.     

In silico studies on potential MCF-7 inhibitors: a combination of pharmacophore and 3D-QSAR modeling,virtual screening,molecular docking,and pharmacokinetic analysis

Gallic acid and its derivatives exhibit a diverse range of biological applications, including anti-cancer activity. In this work, a data-set of forty-six molecules containing the galloyl moiety, and known to show anticarcinogenic activity against the MCF-7 human cancer cell line, have been chosen for pharmacophore modeling and 3D-Quantitative Structure Activity Relationship (3D-QSAR) studies. A tree-based partitioning algorithm has been used to find common pharmacophore hypotheses. The QSAR model was generated for three, four, and five featured hypotheses with increasing PLS factors and analyzed. Results for five featured hypotheses with three acceptors and two aromatic rings were the best out of all the possible combinations. On analyzing the results, the most robust (R^2?=?.8990) hypothesis with a good predictive power (Q^2?=?.7049) was found to be AAARR.35. A good external validation (R² = .6109) was also obtained. In order to design new MCF-7 inhibitors, the QSAR model was further utilized in pharmacophore-based virtual screening of a large database. The predicted IC₅₀ values of the identified potential MCF-7 inhibitors were found to lie in the micromolar range. Molecular docking into the colchicine domain of tubulin was performed in order to examine one of the probable mechanisms. This revealed various interactions between the ligand and the active site protein residues. The present study is expected to provide an effective guide for methodical development of potent MCF-7 inhibitors.     

Structural insights of Staphylococcus aureus FtsZ inhibitors through molecular docking, 3D-QSAR and molecular dynamics simulations

Filamentous temperature-sensitive protein Z (FtsZ) is a protein encoded by the FtsZ gene that assembles into a Z-ring at the future site of the septum of bacterial cell division. Structurally, FtsZ is a homolog of eukaryotic tubulin but has low sequence similarity; this makes it possible to obtain FtsZ inhibitors without affecting the eukaryotic cell division. Computational studies were performed on a series of substituted 3-arylalkoxybenzamide derivatives reported as inhibitors of FtsZ activity in Staphylococcus aureus. Quantitative structure-activity relationship models (QSAR) models generated showed good statistical reliability, which is evident from r²_ncv and r²_loo values. The predictive ability of these models was determined and an acceptable predictive correlation (r²_Pred) values were obtained. Finally, we performed molecular dynamics simulations in order to examine the stability of protein-ligand interactions. This facilitated us to compare free binding energies of cocrystal ligand and newly designed molecule B1. The good concordance between the docking results and comparative molecular field analysis (CoMFA)/comparative molecular similarity indices analysis (CoMSIA) contour maps afforded obliging clues for the rational modification of molecules to design more potent FtsZ inhibitors.     

Pharmacophore mapping,molecular docking and QSAR studies of structurally diverse compounds as CYP2B6 inhibitors

Pharmacophore mapping, molecular docking and quantitative structure–activity relationship (QSAR) studies were carried out for a structurally diverse set of 48 compounds as CYP2B6 inhibitors. The generated best pharmacophore hypotheses from the three methods of conformer generation (FAST, BEST and conformer algorithm based on energy screening and recursive buildup) indicate the importance of two features, namely, hydrogen bond acceptor [electron-rich centre] and ring aromaticity. The distance between the two centres of the important features for ideal inhibitors varied from 5.82 to 6.03 Å. The chemometric tools used for the QSAR analysis were genetic function approximation (GFA) and genetic partial least squares. The developed QSAR models indicate the importance of an electron-rich centre, size of molecule, impact of branching and ring system and distribution of charges in the molecular surface. The docking study confirms the importance of an electron-rich centre for binding with the iron atom of the cytochrome enzyme. A GFA model with spline option was found to be the best model based on internal validation as well as the r ² _{m (overall)} criterion (Q ² = 0.772, r ² _{m (overall)} = 0.774). According to the external prediction statistics (R ² _pred = 0.876), another GFA-derived model with spline option outperforms the remaining models.     

Potential inhibitors of the enzyme acetylcholinesterase and juvenile hormone with insecticidal activity: study of the binding mode via docking and molecular dynamics simulations

Abstract

Models validation in QSAR, pharmacophore, docking and others can ensure the accuracy and reliability of future predictions in design and selection of molecules with biological activity. In this study, pyriproxyfen was used as a pivot/template to search the database of the Maybridge Database for potential inhibitors of the enzymes acetylcholinesterase and juvenile hormone as well. The initial virtual screening based on the 3D shape resulted in 2000 molecules with Tanimoto index ranging from 0.58 to 0.88. A new reclassification was performed on the overlapping of positive and negative charges, which resulted in 100 molecules with Tanimoto's electrostatic score ranging from 0.627 to 0.87. Using parameters related to absorption, distribution, metabolism and excretion and the pivot molecule, the molecules selected in the previous stage were evaluated regarding these criteria, and 21 were then selected. The pharmacokinetic and toxicological properties were considered and for 12 molecules, the DEREK software not fired any alert of toxicity, which were thus considered satisfactory for prediction of biological activity using the Web server PASS. In the molecular docking with insect acetylcholinesterase, the Maybridge3_002654 molecule had binding affinity of ?11.1?kcal/mol, whereas in human acetylcholinesterase, the Maybridge4_001571molecule show in silico affinity of ?10.2?kcal/mol, and in the juvenile hormone, the molecule MCULE-8839595892 show in silico affinity value of ?11.6?kcal/mol. Subsequent long-trajectory molecular dynamics studies indicated considerable stability of the novel molecules compared to the controls.

Abbreviations QSAR quantitative structure–activity relationships

PASS prediction of activity spectra for substances

Communicated by Ramaswamy H. Sarma     

Identification of InhA inhibitors: A combination of virtual screening,molecular dynamics simulations and quantum chemical studies

In the present work, multiple pharmacophore-based virtual screening of the SPECS natural product database was carried out to identify novel inhibitors of the validated biological target, InhA. The pharmacophore models were built from the five different groups of the co-crystallized ligands present within the active site. The generated models with the same features from each group were pooled and subjected to the test set validation, receiver–operator characteristic analysis and Güner–Henry studies. A set of five hypotheses with sensitivity > 0.5, specificity > 0.5, area under curve (AUC) > 0.7, and goodness of hit score > 0.7 were retrieved and exploited for the virtual screening. The common hits (87 molecules) obtained from these hypotheses were processed via drug-likeness filters. The filtered molecules (27 molecules) were compared for the binding modes and the top scored molecules (12 molecules) along with the reference (triclosan (TCL), docking score = ?11.65 kcal/mol) were rescored and reprioritized via molecular mechanics-generalized Born surface area approach. Eventually, the stability of reprioritized (10 molecules) docked complexes was scrutinized via molecular dynamics simulations. Moreover, the quantum chemical studies of the dynamically stable compounds (9 molecules) were performed to understand structural features essential for the activity. Overall, the protocol resulted in the recognition of nine lead compounds that can be targeted against InhA.     

Molecular dynamics-based investigation of InhA substrate binding loop for diverse biological activity of direct InhA inhibitors

The closed conformation of substrate binding loop (SBL) is considered significant for biological activity of direct InhA inhibitors (DIIs). However, exact interactions of SBL with inhibitors are not characterized yet to emphasize over SBL conformations. The seven InhA-DII complexes are analyzed using molecular dynamics simulation to deduce the mechanism for closed and open conformation of SBL. MMGBSA binding energy calculations and decompositions help to identify Ala198, Met199, Ile202, Val203, Ile215, and Leu218 in SBL region as the key residues. The interactions of DIIs with SBL residues particularly Ile202, Val203, Ile215, and Leu218 are found considerable for closed SBL conformation. This difference is accounted for closed state of SBL in 2X23, and open/moderately open state in other complexes. This study substantiates the loop ordering property of DIIs as the basis for high-affinity InhA inhibitors under the molecular recognition phenomena. This property can be used as a parameter to identify potential DIIs using virtual screening approaches.     

Inhibition of the Mycobacterium tuberculosis enoyl acyl carrier protein reductase InhA by arylamides

InhA, the enoyl acyl carrier protein reductase (ENR) from Mycobacterium tuberculosis, is one of the key enzymes involved in the type II fatty acid biosynthesis pathway of M. tuberculosis. We report here the discovery, through high-throughput screening, of a series of arylamides as a novel class of potent InhA inhibitors. These direct InhA inhibitors require no mycobacterial enzymatic activation and thus circumvent the resistance mechanism to antitubercular prodrugs such as INH and ETA that is most commonly observed in drug-resistant clinical isolates. The crystal structure of InhA complexed with one representative inhibitor reveals the binding mode of the inhibitor within the InhA active site. Further optimization through a microtiter synthesis strategy followed by in situ activity screening led to the discovery of a potent InhA inhibitor with in vitro IC(50)=90 nM, representing a 34-fold potency improvement over the lead compound.     

Identification of natural compound inhibitors against PfDXR: A hybrid structure-based molecular modeling approach and molecular dynamics simulation studies

In the present contribution, multicomplex-based pharmacophore studies were carried out on the structural proteome of Plasmodium falciparum 1-deoxy-D -xylulose-5-phosphate reductoisomerase. Among the constructed models, a representative model with complementary features, accountable for the inhibition was used as a primary filter for the screening of database molecules. Auxiliary evaluations of the screened molecules were performed via drug-likeness and molecular docking studies. Subsequently, the stability of the docked inhibitors was envisioned by molecular dynamics simulations, principle component analysis, and molecular mechanics-Poisson-Boltzmann surface area-based free binding energy calculations. The stability assessment of the hits was done by comparing with the reference (beta-substituted fosmidomycin analog, LC5) to prioritize more potent candidates. All the complexes showed stable dynamic behavior while three of them displayed higher binding free energy compared with the reference. The work resulted in the identification of the compounds with diverse scaffolds, which could be used as initial leads for the design of novel PfDXR inhibitors.     

3D-QSAR pharmacophore-based virtual screening,molecular docking and molecular dynamics simulation toward identifying lead compounds for NS2B–NS3 protease inhibitors

NS2B–NS3 protease has been identified to serve as lead drug design target due to its significant role in West Nile viral (WNV) and dengue virus (DENV) reproduction and replication. There are currently no approved chemotherapeutic drugs and effective vaccines to inhibit DENV and WNV infections. In this work, 3D-QSAR pharmacophore model has been developed to discover potential inhibitory candidates. Validation through Fischer’s model and decoy test indicate that the developed 3D pharmacophore model is highly predictive for DENV inhibitors, which was then employed to screen ZINC chemical library to obtain reasonable hits. Following ADMET filtering, 15 hits were subjected to further filter through molecular docking and CoMFA modeling. Finally, top three hits were identified as lead compounds or potential inhibitory candidates with IC₅₀ values of ～0.4637?µM and fitness of ～57.73. It is implied from CoMFA modeling that substituents at the side site of benzotriazole such as a p-nitro group (e.g. biphenyl head) and a carbonyl (e.g. carboxylate function) at the side site of furan or amino group may improve bioactivity of ZINC85645245, respectively. Molecular dynamics simulations (MDS) were performed to discover new interactions and reinforce the binding modes from docking for the hits also. The QSAR and MDS results obtained from this work should be useful in determining structural requirements for inhibitor development as well as in designing more potential inhibitors for NS2B–NS3 protease.     

Computational analyses of interactions between ALK-5 and bioactive ligands: insights for the design of potential anticancer agents

Activin Receptor-Like Kinase 5 (ALK-5) is related to some types of cancer, such as breast, lung, and pancreas. In this study, we have used molecular docking, molecular dynamics simulations, and free energy calculations in order to explore key interactions between ALK-5 and six bioactive ligands with different ranges of biological activity. The motivation of this work is the lack of crystal structure for inhibitor–protein complexes for this set of ligands. The understanding of the molecular structure and the protein–ligand interaction could give support for the development of new drugs against cancer. The results show that the calculated binding free energy using MM-GBSA, MM-PBSA, and SIE is correlated with experimental data with r² = 0.88, 0.80, and 0.94, respectively, which indicates that the calculated binding free energy is in excellent agreement with experimental data. In addition, the results demonstrate that H bonds with Lys232, Glu245, Tyr249, His283, Asp351, and one structural water molecule play an important role for the inhibition of ALK-5. Overall, we discussed the main interactions between ALK-5 and six inhibitors that may be used as starting points for designing new molecules to the treatment of cancer.     

Discovery of novel N-methyl carbazole tethered rhodanine derivatives as direct inhibitors of Mycobacterium tuberculosis InhA

InhA (Enoyl-ACP reductase) plays a crucial role in the biosynthetic pathway of cell wall synthesis in Mycobacterium tuberculosis (Mtb). Isoniazid (INH) is an important first-line drug, which inhibits InhA. The rapid increase in resistance to INH and currently marketed drugs as well as emergence of MDR-TB and XDR-TB has complicated the diagnosis and treatment of Mtb with ever increasing threat to human kind. Herein, we report novel N-methyl carbazole derivatives as potential anti-TB compounds acting directly via InhA inhibition. All the synthesized final compounds were screened against Mtb virulent cell line H₃₇Rv and investigated the InhA enzyme inhibition. Interestingly, compound 9e displayed promising inhibition (91%) at 50 µM concentration and IC₅₀ of 2.82 µM against InhA. To understand the ligand receptor interaction between compound 9e and InhA, molecular docking and molecular dynamics experiments were performed. The computational results were in agreement with the observed experimental data. Further, the cytotoxicity studies on mammalian cells revealed that all the compounds were safe.     

2D-QSAR study,molecular docking,and molecular dynamics simulation studies of interaction mechanism between inhibitors and transforming growth factor-beta receptor I (ALK5)

Transforming growth factor type 1 receptor (ALK5) is kinase associated with a wide variety of pathological processes, and inhibition of ALK5 is a good strategy to treat many kinds of cancer and fibrotic diseases. Recently, a series of compounds have been synthesized as ALK5 inhibitors. However, the study of their selectivity against other potential targets remains elusive. In this research, a data-set of ALK5 inhibitors were collected and studied based on the combination of 2D-QSAR, molecular docking and molecular dynamics simulation. The quality of QSAR models were assessed statistically by F, R², and R²_ADJ, proved to be credible. The cross-validations for the models (q²_LOO = 0.571 and 0.629, respectively) showed their robustness, while the external validations (r²_test = 0.703 and 0.764, respectively) showed their predictive power. Besides, the predicted binding free energy results calculated by MM/GBSA method were in accordance with the experimental data, and the van der Waals energy term was the factor that had the most significant impact on ligand binding. What is more, several important residues were found to significantly affect the binding affinity. Finally, based on our analyses above, a proposed series of molecules were designed.