Discover potential inhibitors for PFKFB3 using 3D-QSAR,virtual screening,molecular docking and molecular dynamics simulation

Abstract

The 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase-3 (PFKFB3) is a master regulator of glycolysis in cancer cells by synthesizing fructose-2,6-bisphosphate (F-2,6-BP), a potent allosteric activator of phosphofructokinase-1 (PFK-1), which is a rate-limiting enzyme of glycolysis. PFKFB3 is an attractive target for cancer treatment. It is valuable to discover promising inhibitors by using 3D-QSAR pharmacophore modeling, virtual screening, molecular docking and molecular dynamics simulation. Twenty molecules with known activity were used to build 3D-QSAR pharmacophore models. The best pharmacophore model was ADHR called Hypo1, which had the highest correlation value of 0.98 and the lowest RMSD of 0.82. Then, the Hypo1 was validated by cost value method, test set method and decoy set validation method. Next, the Hypo1 combined with Lipinski's rule of five and ADMET properties were employed to screen databases including Asinex and Specs, total of 1,048,159 molecules. The hits retrieved from screening were docked into protein by different procedures including HTVS, SP and XP. Finally, nine molecules were picked out as potential PFKFB3 inhibitors. The stability of PFKFB3-lead complexes was verified by 40?ns molecular dynamics simulation. The binding free energy and the energy contribution of per residue to the binding energy were calculated by MM-PBSA based on molecular dynamics simulation.     

Pharmacophore generation,atom-based 3D-QSAR,molecular docking and molecular dynamics simulation studies on benzamide analogues as FtsZ inhibitors

FtsZ is an appealing target for the design of antimicrobial agent that can be used to defeat the multidrug-resistant bacterial pathogens. Pharmacophore modelling, molecular docking and molecular dynamics (MD) simulation studies were performed on a series of three-substituted benzamide derivatives. In the present study a five-featured pharmacophore model with one hydrogen bond acceptors, one hydrogen bond donors, one hydrophobic and two aromatic rings was developed using 97 molecules having MIC values ranging from .07 to 957 μM. A statistically significant 3D-QSAR model was obtained using this pharmacophore hypothesis with a good correlation coefficient (R² = .8319), cross validated coefficient (Q² = .6213) and a high Fisher ratio (F = 103.9) with three component PLS factor. A good correlation between experimental and predicted activity of the training (R² = .83) and test set (R² = .67) molecules were displayed by ADHRR.1682 model. The generated model was further validated by enrichment studies using the decoy test and MAE-based criteria to measure the efficiency of the model. The docking studies of all selected inhibitors in the active site of FtsZ protein showed crucial hydrogen bond interactions with Val 207, Asn 263, Leu 209, Gly 205 and Asn-299 residues. The binding free energies of these inhibitors were calculated by the molecular mechanics/generalized born surface area VSGB 2.0 method. Finally, a 15 ns MD simulation was done to confirm the stability of the 4DXD–ligand complex. On a wider scope, the prospect of present work provides insight in designing molecules with better selective FtsZ inhibitory potential.     

Cloning and Molecular Modeling of Duodenase with Respect to Evolution of Substrate Specificity within Mammalian Serine Proteases That Have Lost a Conserved Active-Site Disulfide Bond

Mammalian serine proteases such as the chromosome 14 (Homo sapiens, Mus musculus) located granzymes, chymases, cathepsin G, and related enzymes including duodenase collectively represent a special group within the chymotrypsin family which we refer to here as "granases". Enzymes of this group have lost the ancient active-site disulfide bond Cys191-Cys220 (bovine chymotrypsinogen A numbering) which is strongly conserved in classic serine proteases such as pancreatic, blood coagulation, and fibrinolysis proteases and others (granzymes A, M, K and leukocyte elastases). We sequenced the cDNA encoding bovine (Bos taurus) duodenase, a granase with unusual dual trypsin-like and chymotrypsin-like specificity. The sequence revealed a 17-residue signal peptide and two-residue (GlyLys) activation peptide typical for granases. Production of the mature enzyme is apparently accompanied by further proteolytic processing of the C-terminal pentapeptide extension of duodenase. Similar C-terminal processing is known for another dual-specific granase, human cathepsin G. Using phylogenetic analysis based on 39 granases we retraced the evolution of residues 189 and 226 crucial for serine protease primary specificity. The analysis revealed that while there is no obvious link between mutability of residue 189 and the appearance of novel catalytic properties in granases, the mutability of residue 226 evidently gives rise to different specificity subgroups within this enzyme group. The architecture of the extended substrate-binding site of granases and structural basis of duodenase dual specificity based on molecular dynamic method are discussed. We conclude that the marked selectivity of granases that is crucial to their role as regulatory proteases has evolved through the fine-tuning of specificity at three levels--primary, secondary, and conformational.     

3D-QSAR,molecular docking,and molecular dynamic simulations for prediction of new Hsp90 inhibitors based on isoxazole scaffold

Heat shock protein 90(Hsp90), as a molecular chaperone, play a crucial role in folding and proper function of many proteins. Hsp90 inhibitors containing isoxazole scaffold are currently being used in the treatment of cancer as tumor suppressers. Here in the present studies, new compounds based on isoxazole scaffold were predicted using a combination of molecular modeling techniques including three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamic (MD) simulations. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were also done. The steric and electrostatic contour map of CoMFA and CoMSIA were created. Hydrophobic, hydrogen bond donor and acceptor of CoMSIA model also were generated, and new compounds were predicted by CoMFA and CoMSIA contour maps. To investigate the binding modes of the predicted compounds in the active site of Hsp90, a molecular docking simulation was carried out. MD simulations were also conducted to evaluate the obtained results on the best predicted compound and the best reported Hsp90 inhibitors in the 3D-QSAR model. Findings indicate that the predicted ligands were stable in the active site of Hsp90.     

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     

3D-QSAR and molecular recognition of Klebsiella pneumoniae NDM-1 inhibitors

New Delhi metallo-β-lactamase-1 (NDM-1) as a target for the development of anti-superbug agents, plays an important role in the resistance of β-lactam antibiotics and has received worldwide attention. Sulfhydryl propionic acid derivatives can effectively inhibit the catalytic activity of NDM-1, but the quantitative structure–activity relationship (QSAR) and inhibitor-target recognition mechanism both remain unclear. In this work, CoMFA and CoMSIA models of sulfhydryl propionic acid inhibitors with high predictive ability were obtained, from which the effect of different substituents on the inhibitory activity against NDM-1 were revealed at the molecular level. Then, two 120-nanosecond comparative molecular dynamics (MD) simulations for NDM-1 enzyme and NDM-1-inhibitor complex systems were performed to study the recognition and inhibition mechanism of sulfhydryl propionic acid derivatives. The binding of inhibitors alters the entire H-bond network of the NDM-1 system accompanied by the formation of strong interactions with I35, W93, H120, H122, D124, H189 and H250, further weakens the recognition of NDM-1 by its endogenic substrates. Finally, the results of free energy landscape and conformation cluster analyses show that NDM-1 underwent a significant conformational change after the association with sulfhydryl propionic acid inhibitors. Our findings can provide theoretical support and help for anti-superbug agents design based on the structures of NDM-1 and sulfhydryl propionic acid derivatives.     

Molecular Dynamics Analysis of Chymosin in Solution and in a Crystalline Environment

The method of molecular dynamics in explicit solvent was applied to test the hypothesis of the existence of a self-inhibited form of chymosin in solution. The paths and energies were calculated for chymosin in solution and in a crystalline environment. The modeling revealed that the intermolecular contacts of chymosin in crystal have negligible influence on the energy stabilization of its self-inhibited conformation. On the other hand, upon molecular dynamics simulation of the active and self-inhibited forms in solution their conformational energies proved to be quite close and the potential barrier between them relatively low. All this supports the possibility of chymosin to adopt spontaneously the self-inhibited conformation in solution, and indicates that it is one of the really existing enzyme forms rather than a crystal packing artifact. The results obtained open novel approaches to studying the specificity of chymosin as well as other aspartic proteinases.     

Molecular Modeling of Cytochrome P450 2B1: Mode of Membrane Insertion and Substrate Specificity

A molecular model of a mammalian membrane-bound cytochrome P450, rat P450 2B1, was constructed in order to elucidate its mode of attachment to the endoplasmic reticulum and the structural basis of substrate specificity. The model was primarily derived from the structure of P450BM-3, which as a class II P450 is the most functionally similar P450 of known structure. However, model development was also guided by the conserved core regions of P450cam and P450terp. To optimally align the P450 2B1 and P450BM-3 sequences, multiple alignment was performed using sequences of five P450s in the II family, followed by minor adjustments on the basis of secondary structure predictions. The resulting P450 2B1 homology model structure was refined by molecular dynamics heating, equilibration, simulation, and energy minimization. The model suggests that the F–G loop serves as both a hydrophobic membrane anchor and entrance channel for hydrophobic substrates from the membrane to the P450 active site. To assess the mode of substrate binding, benzphetamine, testosterone, and benzo[a]pyrene were docked into the active site. The hydrophobic substrate-binding pocket is consistent with the preferences of this P450 toward hydrophobic substrates, while the presence of an acidic Glu-105 in this pocket is consistent with the preference of this P450 for the cationic substrate benzphetamine. This model is thus consistent with several known experimental properties of this P450, such as membrane attachment and substrate selectivity.     

3D-QSAR with CoMFA Model of Prolylendopeptidase Substrates

Abstract

The prolylendopeptidase (PEP) is the proteolytic enzyme, which plays an essential role in the regulation of some processes in central nervous system, such as memory, learning and behavior. It was shown that PEP activity changes at different diseases, like Parkinsons or Alzheimer's diseases, and some PEP inhibitors are used in therapy. At present time the discovery of new types of PEP inhibitors are the actual task.

In this study the structure of PEP active site was analyzed by 3D-QSAR with CoMFA methods using of 12 PEP substrates. The designed pharmacophore model assumes that substrates interact with PEP active site by pyrrolidol ring of proline residue and by hydrogen bonding.

The 3-D-QSAR + CoMFA model of PEP substrates propose that the hydrophobic bonds play the essential role in substrate interaction with enzyme. This model reveals the important steric and electrostatic areas around the molecules and the presence of substituents controls the PEP activity for substrates. Analysis of obtained data allows to assume, that substrate binding in PEP active site causes essential perturbations of substrate structure. This effect mainly depends on chemical nature of the amino acid side chain, located near to proline.     

w-转氨酶分子改造研究进展

ω-转氨酶能催化羰基化合物发生不对称还原胺化反应,在制备手性胺类化合物方面具有较好的应用前景。由于底物结合区域特殊的空间结构,野生型ω-转氨酶在合成大位阻手性胺方面的应用受到了限制。此外,在立体选择性和稳定性方面这一类酶也存在一些不足,目前满足工业应用需求的ω-转氨酶仍较为有限。文中首先介绍了ω-转氨酶的结构特征和催化机制,并探讨S型和R型酶在结构特征方面的主要差异。然后对ω-转氨酶的分子改造研究进行了综述,重点阐述了基于结构特征和催化机制进行的分子改造研究,包括底物特异性改造、立体选择性改造和稳定性改造三方面。最后,对ω-转氨酶分子改造研究进展进行总结和展望。