The InhA inhibitors play key role in mycolic acid synthesis by preventing the fatty acid biosynthesis pathway. In this present article, Pharmacophore modelling and molecular docking study followed by ...in silico virtual screening could be considered as effective strategy to identify newer enoyl-ACP reductase inhibitors. Pyrrolidine carboxamide derivatives were opted to generate pharmacophore models using HypoGen algorithm in Discovery studio 2.1. Further it was employed to screen Zinc and Minimaybridge databases to identify and design newer potent hit molecules. The retrieved newer hits were further evaluated for their drug likeliness and docked against enoyl acyl carrier protein reductase. Here, novel pyrazolo1,5-apyrimidine analogues were designed and synthesized with good yields. Structural elucidation of synthesized final molecules was perform through IR, MASS,
1
H-NMR,
13
C-NMR spectroscopy and further tested for its in vitro anti-tubercular activity against H37Rv strain using Microplate Alamar blue assay (MABA) method. Most of the synthesized compounds displayed strong anti-tubercular activities. Further, these potent compounds were gauged for MDR-TB, XDR-TB and cytotoxic study.
Alzheimer's disease (AD) is one of the most significant neurodegenerative disorders and its symptoms mostly appear in aged people. Catechol-o-methyltransferase (COMT) is one of the known target ...enzymes responsible for AD. With the use of 23 known inhibitors of COMT, a query has been generated and validated by screening against the database of 1500 decoys to obtain the GH score and enrichment value. The crucial features of the known inhibitors were evaluated by the online ZINC Pharmer to identify new leads from a ZINC database. Five hundred hits were retrieved from ZINC Pharmer and by ADMET (absorption, distribution, metabolism, excretion, and toxicity) filtering by using FAF-Drug-3 and 36 molecules were considered for molecular docking. From the COMT inhibitors, opicapone, fenoldopam, and quercetin were selected, while ZINC63625100_413 ZINC39411941_412, ZINC63234426_254, ZINC63637968_451, and ZINC64019452_303 were chosen for the molecular dynamics simulation analysis having high binding affinity and structural recognition. This study identified the potential COMT inhibitors through pharmacophore-based inhibitor screening leading to a more complete understanding of molecular-level interactions.
In this work, the binding mechanism of new Polyketide Synthase 13 (Pks13) inhibitors has been studied through molecular dynamics simulation and free energy calculations. The drug Tam1 and its ...analogs, belonging to the benzofuran class, were submitted to 100 ns simulations, and according to the results obtained for root mean square deviation, all the simulations converged from approximately 30 ns. For the analysis of backbone flotation, the root mean square fluctuations were plotted for the Cα atoms; analysis revealed that the greatest fluctuation occurred in the residues that are part of the protein lid domain. The binding free energy value (ΔG
bind
) obtained for the Tam16 lead molecule was of −51.43 kcal/mol. When comparing this result with the ΔG
bind
values for the remaining analogs, the drug Tam16 was found to be the highest ranked: this result is in agreement with the experimental results obtained by Aggarwal and collaborators, where it was verified that the IC
50
for Tam16 is the smallest necessary to inhibit the Pks13 (IC
50
= 0.19 μM). The energy decomposition analysis suggested that the residues which most interact with inhibitors are: Ser1636, Tyr1637, Asn1640, Ala1667, Phe1670, and Tyr1674, from which the greatest energy contribution to Phe1670 was particularly notable. For the lead molecule Tam16, a hydrogen bond with the hydroxyl of the phenol not observed in the other analogs induced a more stable molecular structure. Aggarwal and colleagues reported this hydrogen bonding as being responsible for the stability of the molecule, optimizing its physic-chemical, toxicological, and pharmacokinetic properties.
Thrombin is a key component for chemotherapeutic and antithrombotic therapy development. As the physiologic and pathologic roles of the light chain still remain vague, here, we continue previous ...efforts to understand the impacts of the disease-associated single deletion of LYS9 in the light chain. By combining supervised and unsupervised machine learning methodologies and more traditional structural analyses on data from 10 μs molecular dynamics simulations, we show that the conformational ensemble of the ΔK9 mutant is significantly perturbed. Our analyses consistently indicate that LYS9 deletion destabilizes both the catalytic cleft and regulatory functional regions and result in some conformational changes that occur in tens to hundreds of nanosecond scaled motions. We also reveal that the two forms of thrombin each prefer a distinct binding mode of a Na
+
ion. We expand our understanding of previous experimental observations and shed light on the mechanisms of the LYS9 deletion associated bleeding disorder by providing consistent but more quantitative and detailed structural analyses than early studies in literature. With a novel application of supervised learning, i.e. the decision tree learning on the hydrogen bonding features in the wild-type and ΔK9 mutant forms of thrombin, we predict that seven pairs of critical hydrogen bonding interactions are significant for establishing distinct behaviors of wild-type thrombin and its ΔK9 mutant form. Our calculations indicate the LYS9 in the light chain has both localized and long-range allosteric effects on thrombin, supporting the opinion that light chain has an important role as an allosteric effector.
Kinesin spindle protein (KSP) belongs to the kinesin superfamily of microtubule-based motor proteins. KSP is responsible for the establishment of the bipolar mitotic spindle which mediates cell ...division. Inhibition of KSP expedites the blockade of the normal cell cycle during mitosis through the generation of monoastral MT arrays that finally cause apoptotic cell death. As KSP is highly expressed in proliferating/cancer cells, it has gained considerable attention as a potential drug target for cancer chemotherapy. Therefore, this study envisaged to design novel KSP inhibitors by employing computational techniques/tools such as pharmacophore modelling, virtual database screening, molecular docking and molecular dynamics. Initially, the pharmacophore models were generated from the data-set of highly potent KSP inhibitors and the pharmacophore models were validated against in house test set ligands. The validated pharmacophore model was then taken for database screening (Maybridge and ChemBridge) to yield hits, which were further filtered for their drug-likeliness. The potential hits retrieved from virtual database screening were docked using CDOCKER to identify the ligand binding landscape. The top-ranked hits obtained from molecular docking were progressed to molecular dynamics (AMBER) simulations to deduce the ligand binding affinity. This study identified MB-41570 and CB-10358 as potential hits and evaluated these experimentally using in vitro KSP ATPase inhibition assays.
Brown plant hopper (BPH) is one of the major destructive insect pests of rice, causing severe yield loss. Thirty-two BPH resistance genes have been identified in cultivated and wild species of rice ...Although, molecular mechanism of rice plant resistance against BPH studied through map-based cloning, due to non-existence of NMR/crystal structures of Bph14 protein, recognition of leucine-rich repeat (LRR) domain and its interaction with different ligands are poorly understood. Thus, in the present study, in silico approach was adopted to predict three-dimensional structure of LRR domain of Bph14 using comparative modelling approach followed by interaction study with jasmonic and salicylic acids. LRR domain along with LRR-jasmonic and salicylic acid complexes were subjected to dynamic simulation using GROMACS, individually, for energy minimisation and refinement of the structure. Final binding energy of jasmonic and salicylic acid with LRR domain was calculated using MM/PBSA. Free-energy landscape analysis revealed that overall stability of LRR domain of Bph14 is not much affected after forming complex with jasmonic and salicylic acid. MM/PBSA analysis revealed that binding affinities of LRR domain towards salicylic acid is higher as compared to jasmonic acid. Interaction study of LRR domain with salicylic acid and jasmonic acid reveals that THR987 of LRR form hydrogen bond with both complexes. Thus, THR987 plays active role in the Bph14 and phytochemical interaction for inducing resistance in rice plant against BPH. In future, Bph14 gene and phytochemicals could be used in BPH management and development of novel resistant varieties for increasing rice yield.
5-HydroxyTriptamine 2A antagonists are potential targets for treatment of various cerebrovascular and cardiovascular disorders. In this study, we have developed and performed a unique screening ...pipeline for filtering ZINC database compounds on the basis of similarities to known antagonists to determine novel small molecule antagonists of 5-HydroxyTriptamine 2A. The screening pipeline is based on 2D similarity, 3D dissimilarity and a combination of 2D/3D similarity. The shortlisted compounds were docked to a 5-HydroxyTriptamine 2A homology-based model, and complexes with low binding energies (287 complexes) were selected for molecular dynamics (MD) simulations in a lipid bilayer. The MD simulations of the shortlisted compounds in complex with 5-HydroxyTriptamine 2A confirmed the stability of the complexes and revealed novel interaction insights. The receptor residues S239, N343, S242, S159, Y370 and D155 predominantly participate in hydrogen bonding. π-π stacking is observed in F339, F340, F234, W151 and W336, whereas hydrophobic interactions are observed amongst V156, F339, F234, V362, V366, F340, V235, I152 and W151. The known and potential antagonists shortlisted by us have similar overlapping molecular interaction patterns. The 287 potential 5-HydroxyTriptamine 2A antagonists may be experimentally verified.
The Escherichia coli uracil/H + symporter UraA, known as the representative nucleobase/cation symporter 2(NCS2) protein, gets involved in several crucial physiological processes for most living ...organisms on Earth, such as the uptake of nucleobases and transport of vitamin C. Some experiments proposed a working model to explain proton-coupling and uracil transporting process of UraA on the basis of the crystal structure of NCS2 protein, but the details of conformational changes remained unknown. Thus, in order to make clear conformational changes caused by the protonation and deprotonation process of some conserved proton-coupled residues, the molecular dynamics simulation was used to study the conformation of UraA complexes in different protonation states. The results demonstrated that the protonation of residue Glu241 and Glu290 resulted in the whole conformational transition from the inward-open to the outward-open state. It can be concluded that Glu290 was crucial in a network of hydrogen-bonds in the middle of the core domain involving another essential residue, mainly including tyr288 in TM8, Tyr342, Ser338 in TM12, and the network of hydrogen-bonds was the key to maintain the stability of conformation. Protonation of Glu290 affects the stability of network of H-bond and changed the domains TM3 TM10 TM12. Thus, Glu290 may play a vital role as a 'proton trigger' that affects spatial structural of amino and residues near substrate binding side leading to an outward-open conformation transition.
Protein kinase monopolar spindle 1 plays an important role in spindle assembly checkpoint at the onset of mitosis. Over expression of MPS1 correlated with a wide range of human tumors makes it an ...attractive target for finding an effective and specific inhibitor. In this work, we performed molecular dynamics simulations of protein MPS1 itself as well as protein bound systems with the inhibitor and natural substrate based on crystal structures. The reported orally bioavailable 1 h-pyrrolo 3,2-c pyridine inhibitors of MPS1 maintained stable binding in the catalytic site, while natural substrate ATP could not stay. Comparative study of stability and flexibility of three systems reveals position shifting of β-sheet region within the catalytic site, which indicates inhibition mechanism was through stabilizing the β-sheet region. Binding free energies calculated with MM-GB/PBSA method shows different binding affinity for inhibitor and ATP. Finally, interactions between protein and inhibitor during molecular dynamic simulations were measured and counted. Residue Gly605 and Leu654 were suggested as important hot spots for stable binding of inhibitor by molecular dynamic simulation. Our results reveal an important position shifting within catalytic site for non-inhibited proteins. Together with hot spots found by molecular dynamic simulation, the results provide important information of inhibition mechanism and will be referenced for designing novel inhibitors.