Molecular docking and structure-based virtual screening studies of potential drug target,CAAX prenyl proteases,of Leishmania donovani

Targeting CAAX prenyl proteases of Leishmania donovani can be a good approach towards developing a drug molecule against Leishmaniasis. We have modeled the structure of CAAX prenyl protease I and II of L. donovani, using homology modeling approach. The structures were further validated using Ramachandran plot and ProSA. Active site prediction has shown difference in the amino acid residues present at the active site of CAAX prenyl protease I and CAAX prenyl protease II. The electrostatic potential surface of the CAAX prenyl protease I and II has revealed that CAAX prenyl protease I has more electropositive and electronegative potentials as compared CAAX prenyl protease II suggesting significant difference in their activity. Molecular docking with known bisubstrate analog inhibitors of protein farnesyl transferase and peptidyl (acyloxy) methyl ketones reveals significant binding of these molecules with CAAX prenyl protease I, but comparatively less binding with CAAX prenyl protease II. New and potent inhibitors were also found using structure-based virtual screening. The best docked compounds obtained from virtual screening were subjected to induced fit docking to get best docked configurations. Prediction of drug-like characteristics has revealed that the best docked compounds are in line with Lipinski’s rule. Moreover, best docked protein–ligand complexes of CAAX prenyl protease I and II are found to be stable throughout 20 ns simulation. Overall, the study has identified potent drug molecules targeting CAAX prenyl protease I and II of L. donovani whose drug candidature can be verified further using biochemical and cellular studies.     

Virtual screening,docking, and dynamics of potential new inhibitors of dihydrofolate reductase from Yersinia pestis

In the present work, we propose to design drugs that target the enzyme dihydrofolate redutase (DHFR) as a means of a novel drug therapy against plague. Potential inhibitors of DHFR from Yersinia pestis (YpDHFR) were selected by virtual screening and subjected to docking, molecular dynamics (MD) simulations, and Poisson–Boltzmann surface area method, in order to evaluate their interactions in the active sites of YpDHFR and human DHFR (HssDHFR). The results suggested selectivity for three compounds that were further used to propose the structures of six new potential selective inhibitors for YpDHFR.     

Inhibition of cysteine proteases by a natural biflavone: behavioral evaluation of fukugetin as papain and cruzain inhibitor

Cruzain is the major cysteine protease of Trypanosoma cruzi, the infectious agent responsible for Chagas disease, and cruzain inhibitors display considerable antitrypanosomal activity. In the present work we elucidated crystallographic data of fukugetin, a biflavone isolated from Garcinia brasiliensis, and investigated the role of this molecule as cysteine protease inhibitor. The kinetic analyses demonstrated that fukugetin inhibited cruzain and papain by a slow reversible type inhibition with K_I of 1.1 and 13.4 µM, respectively. However, cruzain inhibition was about 12 times faster than papain inhibition. Lineweaver–Burk plots demonstrated partial competitive inhibition for cruzain and hyperbolic mixed-type inhibition for papain. Furthermore, the docking results showed that the biflavone binds to ring C′ in the S2 pocket and to ring C in the S3 pocket through hydrophobic interactions and hydrogen bonds. Finally, fukugetin also presented inhibitory activity on proteases of the T. cruzi extract, with IC₅₀ of 7 µM.     

Mutations in exons 10 and 11 of human glucokinase result in conformational variations in the active site of the structure contributing to poor substrate binding – explains hyperglycemia in type 2 diabetic patients

Mutations in the glucokinase (GK) gene play a critical role in the establishment of type 2 diabetes. In our earlier study, R308K mutation in GK in a clinically proven type 2 diabetic patient showed, structural and functional variations that contributed immensely to the hyperglycemic condition. In the extension of this work, a cohort of 30 patients with established type 2 diabetic condition were chosen and the exons 10 and 11 of GK were PCR-amplified and sequenced. The sequence alignment showed A379S, D400Y, E300A, E395A, E395G, H380N, I348N, L301M, M298I, M381G, M402R, R308K, R394P, R397S, and S398R mutations in 12 different patients. The structural analysis of these mutated GKs, showed a variable number of β-α-β units, hairpins, β-bulges, strands, helices, helix–helix interactions, β-turns, and γ-turns along with the RMSD variations when compared to wild-type GK. Molecular modeling studies revealed that the substrate showed variable binding orientations and could not fit into the active site of these mutated structures; moreover, it was expelled out of the conformations. Therefore, these structural variations in GK due to mutations could be one of the strongest reasons for the hyperglycemic levels in these type 2 diabetic patients.     

Heterodimerization of the Entamoeba histolytica EhCPADH virulence complex through molecular dynamics and protein–protein docking

EhCPADH is a protein complex involved in the virulence of Entamoeba histolytica, the protozoan responsible for human amebiasis. It is formed by the EhCP112 cysteine protease and the EhADH adhesin. To explore the molecular basis of the complex formation, three-dimensional models were built for both proteins and molecular dynamics simulations (MDS) and docking calculations were performed. Results predicted that the pEhCP112 proenzyme and the mEhCP112 mature enzyme were globular and peripheral membrane proteins. Interestingly, in pEhCP112, the propeptide appeared hiding the catalytic site (C167, H329, N348); while in mEhCP112, this site was exposed and its residues were found structurally closer than in pEhCP112. EhADH emerged as an extended peripheral membrane protein with high fluctuation in Bro1 and V shape domains. 500 ns-long MDS and protein–protein docking predictions evidenced different heterodimeric complexes with the lowest free energy. pEhCP112 interacted with EhADH by the propeptide and C-terminal regions and mEhCP112 by the C-terminal through hydrogen bonds. In contrast, EhADH bound to mEhCP112 by 442–479 residues, adjacent to the target cell-adherence region (480–600 residues), and by the Bro1 domain (9–349 residues). Calculations of the effective binding free energy and per residue free energy decomposition showed that EhADH binds to mEhCP112 with a higher binding energy than to pEhCP112, mainly through van der Waals interactions and the nonpolar part of solvation energy. The EhADH and EhCP112 structural relationship was validated in trophozoites by immunofluorescence, TEM, and immunoprecipitation assays. Experimental findings fair agreed with in silico results.     

Structural basis of pesticide detection by enzymatic biosensing: a molecular docking and MD simulation study

Designing of rapid, facile, selective, and cost-effective biosensor technology is a growing area for the detection of various classes of pesticides. The biosensor with these features can be achieved only through the various bio-components using different transducers. This study, therefore, focuses on the usage of molecular docking, specificity tendencies, and capabilities of proteins for the detection of pesticides. Accordingly, the four transducers, acetylcholinesterase (ACH), cytochromes P450 (CYP), glutathione S-transferase (GST), and protein kinase C (PKC) were selected based on their applications including neurotransmitter, metabolism, detoxification enzyme, and protein phosphorylation. Then after molecular docking of the pesticides, fenobucarb, dichlorodiphenyltrichloroethane (DDT), and parathion onto each enzyme, the conformational behavior of the most stable complexes was further analyzed using 50 ns Molecular Dynamics (MD) simulations carried out under explicit water conditions. In the case of protein kinase C (PKC) and cytochrome P450 3A4 enzyme (CYP), the fenobucarb complex showed the most suitable combination of free energy of binding and inhibition constant ?4.42 kcal/mol (573.73 μM) and ?5.1 kcal/mol (183.49 μM), respectively. Parathion dominated for acetylcholinesterase (ACH) with ?4.57 kcal/mol (448.09 μM) and lastly dichlorodiphenyltrichloroethane for glutathione S-transferase (GST), ?5.43 kcal/mol (103.88 μM). The RMSD variations were critical for understanding the impact of pesticides as they distinctively influence the energetic attributes of the proteins. Overall, the outcomes from the extensive analysis provide an insight into the structural features of the proteins studied, thereby highlighting their potential use as a substrate in biorecognition sensing of pesticide compounds.     

Molecular dynamics and molecular docking studies on E166A point mutant,R274N/R276N double mutant,and E166A/R274N/R276N triple mutant forms of class A β-lactamases

Bacterial resistance to β-lactams antibiotics is a serious threat to human health. The most common cause of resistance to the β-lactams is the production of β-lactamase that inactivates β-lactams. Specifically, class A extended-spectrum β-lactamase produced by antibiotic resistant bacteria is capable of hydrolyzing extended-spectrum Cephalosporins and Monobactams. Mutations in class A β-lactamases play a crucial role in substrate and inhibitor specificity. In this present study, the E166A point mutant, R274N/R276N double mutant, and E166A/R274N/R276N triple mutant class A β-lactamases are analyzed. Molecular dynamics (MD) simulations are done to understand the consequences of mutations in class A β-lactamases. Root mean square deviation, root mean square fluctuation, radius of gyration, solvent accessibility surface area, hydrogen bond, and essential dynamics analysis results indicate notable loss in stability for mutant class A β-lactamases. MD simulations of native and mutant structures clearly confirm that the substitution of alanine at the position of 166, Asparagine at 274 and 276 causes more flexibility in 3D space. Molecular docking results indicate the mutation in class A β-lactamases which decrease the binding affinity of Cefpirome and Ceftobiprole which are third and fifth generation Cephalosporins, respectively. MD simulation of Ceftobiprole-native and mutant type Class A β-lactamases complexes reveal that E166A/R274N/R276N mutations alter the structure and notable loss in the stability for Ceftobirole-mutant type Class A β-lactamases complexes. Ceftobiprole is currently prescribed for patients with serious bacterial infections; this phenomenon is the probable cause for the effectiveness of Ceftobiprole in controlling bacterial infections.     

Molecular dynamics simulation and molecular docking studies of 1,4-Dihydropyridines as P-glycoprotein’s allosteric inhibitors

P-glycoprotein (P-gp) is a main factor contributing to multidrug resistance. The effect of this transporter protein on limiting the effectiveness of chemotherapy has been shown by various studies. In a previous report, we synthesized some 14-dihydropyridine (DHP) derivatives as inhibitors of human P-gp. In the present study, a computational approach has been exploited to reveal the main interactions between DHPs and P-gp. In order to do this, homology modeling was performed to obtain a model of the protein. Then, molecular dynamics simulation was used to refine the constructed model of P-gp in the presence of the lipids bilayer. Model validation was performed with several tools. Finally, molecular docking followed by MD simulation of ligand–protein complex was employed to elucidate the binding mode and the dynamical changes of protein with/without DHPs bound. The results emphasized that interaction of the residues Gln912, Ser909, Arg905, Ser474, Val472 with DHPs play a crucial role in the inhibitory of these ligands and this was in a relatively good accordance with the results reported in the experimental studies.     

Molecular dynamics simulation,binding free energy calculation and molecular docking of human D-amino acid oxidase (DAAO) with its inhibitors

Schizophrenia is a mental illness; most affected people live in developing countries, and neither appropriate treatment nor commercial drugs are currently available. One possibility is to inhibit human-d-amino acid oxidase (h-DAAO). In this study, molecular dynamic simulations of the monomer, dimer and tetramer forms of h-DAAO complexed with the inhibitor 3-hydroxyquinolin-2(1H)-one(2) were performed. Seven residues, Leu51, Gln53, Leu215, Tyr228, Ile230, Arg283 and Gly313, were identified as essential for interacting with the inhibitor. Molecular docking of h-DAAO with pyrrole, quinoline and kojic acid derivatives, representing 69 known or potential h-DAAO inhibitors, was also performed. The results indicated that the activity of the inhibitor can be improved by modifying the compounds to have a substituent group capable of interacting with the side chain of Tyr228. Van der Waals interactions of the inhibitor with the hydrophobic pocket of h-DAAO and electrostatic interactions or H-bonds with Arg283 and Gly313 were important elements in determining the efficiency of the inhibitor. These results provide information on the interaction between h-DAAO and its inhibitors at the molecular level and can aid in the design of novel inhibitors against h-DAAO for new drug development in the treatment of schizophrenia.     

分子伴侣HdeA与底物蛋白SurA作用机制的模拟研究

分子伴侣HdeA与底物蛋白间的相互作用可帮助底物蛋白复性,这是肠道致病菌得以在酸性环境中幸存的重要原因之一.为探究HdeA发挥伴侣活性的作用机制,本研究采用分子对接和分子动力学的方法,模拟了HdeA与底物蛋白SurA间的相互作用,计算了二者的结合自由能.通过分析HdeA-SurA复合物体系的作用模式、氢键作用以及能量分解的结果,确定了HdeA与底物蛋白SurA结合时发挥重要作用的关键氨基酸残基.该研究结果为以后采用实验手段探究HdeA与底物蛋白之间的作用提供了重要的理论参考,同时为今后设计与开发HdeA的抑制剂提供了理论指导依据.     

Structural and dynamic properties of thymopoietin mimetics

We propose a hypothesis that the T-cell receptor is a possible target of thymic hormones. We modelled the conformational dynamics of thymopentin and its structural variants in solution, as well as the interactions of these short peptides with the proposed molecular target. Thymopentin is a five-amino-acid fragment of the thymic hormone thymopoietin (residues 32 to 36) that reproduces the immunomodulatory activity of the complete hormone. Using molecular dynamics and flexible docking methods, we demonstrated high-affinity binding of thymopentin and its prospective mimetics with the T-cell receptor. The calculated biological activity spectra of thymopentin and its two promising modifications can be used in immunomodulatory activity screenings with live systems.     

Predicting substrate selectivity between UGT1A9 and UGT1A10 using molecular modelling and molecular dynamics approach

Uridine 5′-diphospho-glucuronosyltransferase-1A9 (UGT1A9) expressed in the liver, shows good sequence identity with UGT1A10, expressed in the intestine. Both uridine 5′-diphospho-glucuronosyltransferase (UGT) isoforms show comprehensive overlapping substrate selectivity but there are differences in stereoselectivity, regiospecificity and rate of glucuronidation of the substrates. Multiple sequence alignment analyses of UGT1A9 and UGT1A10 showed that 13% of the residues in N-terminal domain (NTD) are non-identical between them. Herein, authors attempted homology modelling of UGT1A9 and UGT1A10 and validation using software tools and reported mutagenic studies. A molecular docking study of the known substrates is performed on UGT1A9 and UGT1A10 homology models. The non-identical N-terminal residues ranging from 111 to 117 in UGT1A9 and UGT1A10 were identified to play a central role in the substrate selectivity. However, substrate binding is performed by Ser111, Gly115 and Leu117 in UGT1A10 and Gly111, Asp115 and Phe117 in UGT1A9. This study reports new residues in NTD, showing interaction with uridine 5′-diphospho-glucuronic acid which binds with C-terminal domain. Further, molecular dynamics simulations were carried out to study the role of non-identical residues in substrate identification. The study demonstrates the folding of the UGT enzyme, particularly, helix-loop-helix transition and movement of Nα3-2 helix, in response to substrate and co-substrate binding.     

Computational investigations of gram-negative bacteria phosphopantetheine adenylyltransferase inhibitors using 3D-QSAR,molecular docking and molecular dynamic simulations

Abstract

Phosphopantetheine adenylyltransferase (PPAT) has been recognized as a promising target to develop novel antimicrobial agents, which is a hexameric enzyme that catalyzes the penultimate step in coenzyme A biosynthesis. In this work, molecular modeling study was performed with a series of PPAT inhibitors using molecular docking, three-dimensional qualitative structure-activity relationship (3D-QSAR) and molecular dynamic (MD) simulations to reveal the structural determinants for their bioactivities. Molecular docking study was applied to understand the binding mode of PPAT with its inhibitors. Subsequently, 3D-QSAR model was constructed to find the features required for different substituents on the scaffolds. For the best comparative molecular field analysis (CoMFA) model, the Q² and R² values of which were calculated as 0.702 and 0.989, while they were calculated as 0.767 and 0.983 for the best comparative molecular similarity index analysis model. The statistical data verified the significance and accuracy of our 3D-QSAR models. Furthermore, MD simulations were carried out to evaluate the stability of the receptor–ligand contacts in physiological conditions, and the results were consistent with molecular docking studies and 3D-QSAR contour map analysis. Binding free energy was calculated with molecular mechanics generalized born surface area approach, the result of which coincided well with bioactivities and demonstrated that van der Waals accounted for the largest portion. Overall, our study provided a valuable insight for further research work on the recognition of potent PPAT inhibitors.

Communicated by Ramaswamy H. Sarma