共查询到20条相似文献，搜索用时 78 毫秒
1.
Ramin Ekhteiari Salmas Serdar Durdagi Mehmet Fuat Gulhan Merve Duruyurek Huda I. Abdullah 《Journal of biomolecular structure & dynamics》2018,36(3):609-620
The objective of the present study was to evaluate the effects of propolis, pollen, and caffeic acid phenethyl ester (CAPE) on tyrosine hydroxylase (TH) activity and total RNA levels of Nω-nitro-L-arginine methyl ester (L-NAME) inhibition of nitric oxide synthase in the heart, adrenal medulla, and hypothalamus of hypertensive male Sprague dawley rats. The TH activity in the adrenal medulla, heart, and hypothalamus of the rats was significantly increased in the L-NAME group vs. control (p < 0.05). Treatment with L-NAME led to a significant increase in blood pressure (BP) in the L-NAME group compared to control (p < 0.05). These data suggest that propolis, pollen, and CAPE may mediate diminished TH activity in the heart, adrenal medulla, and hypothalamus in hypertensive rats. The decreased TH activity may be due to the modulation and synthesis of catecholamines and BP effects. In addition, the binding mechanism of CAPE within the catalytic domain of TH was investigated by means of molecular modeling approaches. These data suggest that the amino acid residues, Glu429 and Ser354 of TH may play a pivotal role in the stabilization of CAPE within the active site as evaluated by molecular dynamics (MD) simulations. Gibbs binding free energy (ΔGbinding) of CAPE in complex with TH was also determined by post-processing MD analysis approaches (i.e. Poisson-Boltzmann Surface Area (MM-PBSA) method). 相似文献
2.
Khaled M. Tumbi Prajwal P. Nandekar Naeem Shaikh Siddharth S. Kesharwani Abhay T. Sangamwar 《Journal of molecular recognition : JMR》2014,27(3):138-150
The discovery of novel anticancer molecules 5F‐203 (NSC703786) and 5‐aminoflavone (5‐AMF, NSC686288) has addressed the issues of toxicity and reduced efficacy by targeting over expressed Cytochrome P450 1A1 (CYP1A1) in cancer cells. CYP1A1 metabolizes these compounds into their reactive metabolites, which are proven to mediate their anticancer effect through DNA adduct formation. However, the drug metabolite–DNA binding has not been explored so far. Hence, understanding the binding characteristics and molecular recognition for drug metabolites with DNA is of practical and fundamental interest. The present study is aimed to model binding preference shown by reactive metabolites of 5F‐203 and 5‐AMF with DNA in forming DNA adducts. To perform this, three different DNA crystal structures covering sequence diversity were selected, and 12 DNA‐reactive metabolite complexes were generated. Molecular dynamics simulations for all complexes were performed using AMBER 11 software after development of protocol for DNA‐reactive metabolite system. Furthermore, the MM‐PBSA/GBSA energy calculation, per‐nucleotide energy decomposition, and Molecular Electrostatic Surface Potential analysis were performed. The results obtained from present study clearly indicate that minor groove in DNA is preferable for binding of reactive metabolites of anticancer compounds. The binding preferences shown by reactive metabolites were also governed by specific nucleotide sequence and distribution of electrostatic charges in major and minor groove of DNA structure. Overall, our study provides useful insights into the initial step of mechanism of reactive metabolite binding to the DNA and the guidelines for designing of sequence specific DNA interacting anticancer agents. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
3.
4.
The ecdysone receptor is a nuclear hormone receptor that plays a pivotal role in the insect metamorphosis and development. To address the molecular mechanisms of binding and selectivity, the interactions of two typical agonists Ponasterone A and 20-Hydroxyecdysone with Drosophila melanogaster (DME) and Leptinotarsa decemlineata ecdysone (LDE) receptors were investigated by homology modeling, molecular docking, molecular dynamic simulation, and thermodynamic analysis. We discover that 1) the L5-loop, L11-loop, and H12 helix for DME, L7-loop, and L11-loop for LDE are more flexible, which affect the global dynamics of the ligand-binding pocket, thus facilitating the ligand recognition of ecdysone receptor; 2) several key residues (Thr55/Thr37, Phe109/Phe91, Arg95/Arg77, Arg99/Arg81, Phe108/Leu90, and Ala110/Val92) are responsible for the binding of the proteins; 3) the binding-free energy is mainly contributed by the van der Waals forces as well as the electrostatic interactions of ligand and receptor; 4) the computed binding-free energy difference between DME-C1 and LDE-C1 is –4.65 kcal/mol, explains that C1 can form many more interactions with the DME; 5) residues Phe108/Leu90 and Ala110/Val92 have relatively position and orientation difference in the two receptors, accounting most likely for the ligand selectivity of ecdysone receptor from different orders of insects. This study underscores the expectation that different insect pests should be able to discriminate among compounds from different as yet undiscovered compounds, and the results firstly show a structural and functional relay between the agonists and receptors (DME and LDE), which can provide an avenue for the development of target-specific insecticides.
Communicated by Ramaswamy H. Sarma 相似文献
5.
The 2013 Nobel Prize in Chemistry has convinced the world that how important the role that computational sciences play in chemical and materials sciences. In this review, computational methods and rational molecule design, including quantum mechanics and molecular mechanics methods, have been applied to study electronic structures and the interactions in a number of important applications at molecular level. The applications which include bioactive compounds, drug candidates and photoactive molecules at Swinburne University in the past several years are discussed. The research is in close collaboration with world class experimental groups from spectroscopy, organic and medicinal synthesis laboratories and most recently to γ-ray spectroscopy as well as other theory groups in the world. Ionisation spectra of biomolecules and bioactive compounds including amino acids, DNA bases, cyclic dipeptides, drug candidates, complexes and photoactive molecules are discussed. Most recent projects such as infrared spectral studies of ferrocene, rational design of organic dyes in solar cell applications, and recent development in γ-ray spectra of positron annihilation in molecules are highlighted. 相似文献
6.
The HIV Rev protein mediates the nuclear export of viral mRNA, and is thereby essential for the production of late viral proteins in the replication cycle. Rev forms a large organized multimeric protein-protein complex for proper functioning. Recently, the three-dimensional structures of a Rev dimer and tetramer have been resolved and provide the basis for a thorough structural analysis of the binding interaction. Here, molecular dynamics (MD) and binding free energy calculations were performed to elucidate the forces thriving dimerization and higher order multimerization of the Rev protein. It is found that despite the structural differences between each crystal structure, both display a similar behavior according to our calculations. Our analysis based on a molecular mechanics-generalized Born surface area (MM/GBSA) and a configurational entropy approach demonstrates that the higher order multimerization site is much weaker than the dimerization site. In addition, a quantitative hot spot analysis combined with a mutational analysis reveals the most contributing amino acid residues for protein interactions in agreement with experimental results. Additional residues were found in each interface, which are important for the protein interaction. The investigation of the thermodynamics of the Rev multimerization interactions performed here could be a further step in the development of novel antiretrovirals using structure based drug design. Moreover, the variability of the angle between each Rev monomer as measured during the MD simulations suggests a role of the Rev protein in allowing flexibility of the arginine rich domain (ARM) to accommodate RNA binding. 相似文献
7.
Qiao Xue Qing‐Chuan Zheng Ji‐Long Zhang Ying‐Lu Cui Hong‐Xing Zhang 《Biopolymers》2014,101(8):849-860
Filoviruses often cause terrible infectious disease which has not been successfully dealt with pharmacologically. All filoviruses encode a unique protein termed VP35 which can mask doubled‐stranded RNA to deactivate interferon. The interface of VP35–dsRNA would be a feasible target for structure‐based antiviral agent design. To explore the essence of VP35–dsRNA interaction, molecular dynamics simulation combined with MM‐GBSA calculations were performed on Marburg virus VP35–dsRNA complex and several mutational complexes. The energetic analysis indicates that nonpolar interactions provide the main driving force for the binding process. Although the intermolecular electrostatic interactions play important roles in VP35–dsRNA interaction, the whole polar interactions are unfavorable for binding which result in a low binding affinity. Compared with wild type VP35, the studied mutants F228A, R271A, and K298A have obviously reduced binding free energies with dsRNA reflecting in the reduction of polar or nonpolar interactions. The results also indicate that the loss of binding affinity for one dsRNA strand would abolish the total binding affinity. Three important residues Arg271, Arg294, and Lys298 which makes the largest contribution for binding in VP35 lose their binding affinity significantly in mutants. The uncovering of VP35–dsRNA recognition mechanism will provide some insights for development of antiviral drug. © 2014 Wiley Periodicals, Inc. Biopolymers 101: 849–860, 2014. 相似文献
8.
9.
Fangfang Wang Wei Yang Yonghui Shi 《Journal of biomolecular structure & dynamics》2013,31(9):1929-1940
The discovery of clinically relevant inhibitors of retinoic acid receptor-related orphan receptor-gamma-t (RORγt) for autoimmune diseases therapy has proven to be a challenging task. In the present work, to find out the structural features required for the inhibitory activity, we show for the first time a three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulations for a series of novel thiazole/thiophene ketone amides with inhibitory activity at the RORγt receptor. The optimum CoMFA and CoMSIA models, derived from ligand-based superimposition I, exhibit leave-one-out cross-validated correlation coefficient (R2cv) of .859 and .805, respectively. Furthermore, the external predictive abilities of the models were evaluated by a test set, producing the predicted correlation coefficient (R2pred) of .7317 and .7097, respectively. In addition, molecular docking analysis was applied to explore the binding modes between the inhibitors and the receptor. MD simulation and MM/PBSA method were also employed to study the stability and rationality of the derived conformations, and the binding free energies in detail. The QSAR models and the results of molecular docking, MD simulation, binding free energies corroborate well with each other and further provide insights regarding the development of novel RORγt inhibitors with better activity. 相似文献
10.
The initiation of microtubule assembly within cells is guided by a cone shaped multi‐protein complex, γ‐tubulin ring complex (γTuRC) containing γ‐tubulin and atleast five other γ‐tubulin‐complex proteins (GCPs), i.e., GCP2, GCP3, GCP4, GCP5, and GCP6. The rim of γTuRC is a ring of γ‐tubulin molecules that interacts, via one of its longitudinal interfaces, with GCP2, GCP3, or GCP4 and, via other interface, with α/β?tubulin dimers recruited for the microtubule lattice formation. These interactions however, are not well understood in the absence of crystal structure of functional reconstitution of γTuRC subunits. In this study, we elucidate the atomic interactions between γ‐tubulin and GCP4 through computational techniques. We simulated two complexes of γ‐tubulin‐GCP4 complex (we called dimer1 and dimer2) for 25 ns to obtain a stable complex and calculated the ensemble average of binding free energies of ?158.82 and ?170.19 kcal/mol for dimer1 and ?79.53 and ?101.50 kcal/mol for dimer2 using MM‐PBSA and MM‐GBSA methods, respectively. These highly favourable binding free energy values points to very robust interactions between GCP4 and γ‐tubulin. From the results of the free‐energy decomposition and the computational alanine scanning calculation, we identified the amino acids crucial for the interaction of γ‐tubulin with GCP4, called hotspots. Furthermore, in the endeavour to identify chemical leads that might interact at the interface of γ‐tubulin‐GCP4 complex; we found a class of compounds based on the plant alkaloid, noscapine that binds with high affinity in a cavity close to γ‐tubulin‐GCP4 interface compared with previously reported compounds. All noscapinoids displayed stable interaction throughout the simulation, however, most robust interaction was observed for bromo‐noscapine followed by noscapine and amino‐noscapine. This offers a novel chemical scaffold for γ‐tubulin binding drugs near γ‐tubulin‐GCP4 interface. Proteins 2015; 83:827–843. © 2015 Wiley Periodicals, Inc. 相似文献
11.
Single chain fragment variables (ScFvs) have been extensively employed in studying the protein-protein interactions. ScFvs derived from phage display libraries have an additional advantage of being generated against a native antigen, circumventing loss of information on conformational epitopes. In the present study, an attempt has been made to elucidate human chorionic gonadotropin (hCG)-luteinizing hormone (LH) receptor interactions by using a neutral and two inhibitory ScFvs against hCG. The objective was to dock a computationally derived model of these ScFvs onto the crystal structure of hCG and understand the differential roles of the mapped epitopes in hCG-LH receptor interactions. An anti-hCG ScFv, whose epitope was mapped previously using biochemical tools, served as the positive control for assessing the quality of docking analysis. To evaluate the role of specific side chains at the hCG-ScFv interface, binding free energy as well as residue interaction energies of complexes in solution were calculated using molecular mechanics Poisson-Boltzmann/surface area method after performing the molecular dynamic simulations on the selected hCG-ScFv models and validated using biochemical and SPR analysis. The robustness of these calculations was demonstrated by comparing the theoretically determined binding energies with the experimentally obtained kinetic parameters for hCG-ScFv complexes. Superimposition of hCG-ScFv model onto a model of hCG complexed with the 51-266 residues of LH receptor revealed importance of the residues previously thought to be unimportant for hormone binding and response. This analysis provides an alternate tool for understanding the structure-function analysis of ligand-receptor interactions. 相似文献
12.
Effective virtual screening relies on our ability to make accurate prediction of protein-ligand binding, which remains a great challenge. In this work, utilizing the molecular-mechanics Poisson-Boltzmann (or Generalized Born) surface area approach, we have evaluated the binding affinity of a set of 156 ligands to seven families of proteins, trypsin β, thrombin α, cyclin-dependent kinase (CDK), cAMP-dependent kinase (PKA), urokinase-type plasminogen activator, β-glucosidase A, and coagulation factor Xa. The effect of protein dielectric constant in the implicit-solvent model on the binding free energy calculation is shown to be important. The statistical correlations between the binding energy calculated from the implicit-solvent approach and experimental free energy are in the range of 0.56-0.79 across all the families. This performance is better than that of typical docking programs especially given that the latter is directly trained using known binding data whereas the molecular mechanics is based on general physical parameters. Estimation of entropic contribution remains the barrier to accurate free energy calculation. We show that the traditional rigid rotor harmonic oscillator approximation is unable to improve the binding free energy prediction. Inclusion of conformational restriction seems to be promising but requires further investigation. On the other hand, our preliminary study suggests that implicit-solvent based alchemical perturbation, which offers explicit sampling of configuration entropy, can be a viable approach to significantly improve the prediction of binding free energy. Overall, the molecular mechanics approach has the potential for medium to high-throughput computational drug discovery. 相似文献
13.
Tomić A Abramić M Spoljarić J Agić D Smith DM Tomić S 《Journal of molecular recognition : JMR》2011,24(5):804-814
Human dipeptidyl peptidase III (DPP III) is a zinc-exopeptidase with implied roles in protein catabolism, pain modulation, and defense against oxidative stress. To understand the mode of ligand binding into its active site, we performed molecular modeling, site-directed mutagenesis, and biochemical analyses. Using the recently determined crystal structure of the human DPP III we built complexes between both, the wild-type (WT) protein and its mutant H568N with the preferred substrate Arg-Arg-2-naphthylamide (RRNA) and a competitive inhibitor Tyr-Phe-hydroxamate (Tyr-Phe-NHOH). The mutation of the conserved His568, structurally equivalent to catalytically important His231 in thermolysin, to Asn, resulted in a 1300-fold decrease of k(cat) for RRNA hydrolysis and in significantly lowered affinity for the inhibitor. Molecular dynamics simulations revealed the key protein-ligand interactions as well as the ligand-induced reorganization of the binding site and its partial closure. Simultaneously, the non-catalytic domain was observed to stretch and the opening at the wide side of the inter-domain cleft became enhanced. The driving force for these changes was the formation of the hydrogen bond between Asp372 and the bound ligand. The structural and dynamical differences, found for the ligand binding to the WT enzyme and the H568N mutant, and the calculated binding free energies, agree well with the measured affinities. On the basis of the obtained results we suggest a possible reaction mechanism. In addition, this work provides a foundation for further site-directed mutagenesis experiments, as well as for modeling the reaction itself. 相似文献
14.
David Rinaldo Claudio Vita Martin J. Field 《Journal of biomolecular structure & dynamics》2013,31(3):281-297
Abstract Proteins with the ability to specifically bind strontium would potentially be of great use in the field of nuclear waste management. Unfortunately, no such peptides or proteins are known—indeed, it is uncertain whether they exist under natural conditions due to low environmental concentrations of strontium. To investigate the possibility of devising such molecules, one of us (CV), in a previous experimental study [J. Biol. Inorg. Chem. 8, 33440 (2003)], proposed starting from an EF-hand motif of the protein calmodulin and mutating some residues to change the motif's specificity for calcium into one for strontium. In this paper, which represents a theoretical complement to the experimental work, we analyzed small-molecule crystallographic structures and performed quantum chemical calculations to identify possible mutations. We then constructed seven mutant sequences of the EF-hand motif and analyzed their dynamical and binding behaviors using molecular dynamics simulations and free-energy calculations (using the MM/PBSA method). As a result of these analyzes we were able to isolate some characteristics that could lead to mutant peptides with enhanced strontium affinity. 相似文献
15.
HIV-1 protease has been an important drug target for the antiretroviral treatment of HIV infection. The efficacy of protease drugs is impaired by the rapid emergence of resistant virus strains. Understanding the molecular basis and evaluating the potency of an inhibitor to combat resistance are no doubt important in AIDS therapy. In this study, we first identified residues that have significant contributions to binding with six substrates using molecular dynamics simulations and Molecular Mechanics Generalized Born Surface Area calculations. Among the critical residues, Asp25, Gly27, Ala28, Asp29, and Gly49 are well conserved, with which the potent drugs should form strong interactions. We then calculated the contribution of each residue to binding with eight FDA approved drugs. We analyzed the conservation of each protease residue and also compared the interaction between the HIV protease and individual residues of the drugs and substrates. Our analyses showed that resistant mutations usually occur at less conserved residues forming more favorable interactions with drugs than with substrates. To quantitatively integrate the binding free energy and conservation information, we defined an empirical parameter called free energy/variability (FV) value, which is the product of the contribution of a single residue to the binding free energy and the sequence variability at that position. As a validation, the FV value was shown to identify single resistant mutations with an accuracy of 88%. Finally, we evaluated the potency of a newly approved drug, darunavir, to combat resistance and predicted that darunavir is more potent than amprenavir but may be susceptible to mutations on Val32 and Ile84. 相似文献
16.
Venken T Krnavek D Münch J Kirchhoff F Henklein P De Maeyer M Voet A 《Proteins》2011,79(11):3221-3235
VIRus Inhibitory Peptide (VIRIP), a 20 amino acid peptide, binds to the fusion peptide (FP) of human immunodeficiency virus type 1 (HIV-1) gp41 and blocks viral entry. VIRIP derivatives with improved antiviral activity have been developed, and one of those derivatives has recently proven effective and safe in a phase 1/2 clinical trial. Here, molecular dynamics were executed in combination with molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) free energy calculations to explore the binding interaction between VIRIP derivatives and gp41 FP. A promising correlation between antiviral activity and simulated binding free energy was established thanks to restriction of the flexibility of the peptides, inclusion of configurational entropy calculations, and the use of multiple internal dielectric constants for the MM/PBSA calculations depending on the amino acid sequence. Based on these results, a virtual screening experiment was carried out to design VIRIP analogs with further improved antiretroviral activity. A selection of peptides was tested for inhibitory activity and several VIRIP derivatives were identified with significantly enhanced activity compared to the reference peptides. The results demonstrate that computational modeling strategies using an adapted MM/PBSA methodology improve the accuracy of binding free energy calculations of peptide complexes compared to the classic MM/PBSA protocol. As such, this virtual screening approach generated HIV-1 gp41 FP inhibitors with improved antiviral activity that could be useful for future clinical applications. 相似文献
17.
Human heparanase is an endo-beta-D-glycosidase that cleaves heparan sulphate (HS) chains in the extracellular matrix and basement membrane. It is known that the cleavage of HS by heparanase results in cell invasion and metastasis of cancer. Therefore, heparanase is considered an important target for cancer drug development. The three-dimensional structure of heparanase would be useful in the rational design of inhibitors targeted to the enzyme; however, the three-dimensional structure has not yet been determined. In our effort to design inhibitors, we developed a three-dimensional structure of heparanase using a homology-modeling approach. The homology-built structure is consistent to previous bioinformatics and site-mutation experimental results. The heparanase features a (alpha/beta)(8) TIM-barrel fold with two glutamate residues (Glu225 and Glu343) located in the active-site cleft. This feature supports the putative mechanism of proton donor and nucleophilic sites. Docking simulations yielded 41 complex structures, which indicate that the bound inhibitor could block ligand binding into the catalytic site. A free energy of binding model was established for 25 heparanase inhibitors with a training set of 25 heparanase inhibitors using the linear response MM-PBSA approach (LR-MM-PBSA). The correlation between calculated and experimental activity was 0.79 and the reliability of the model was validated with leave-one-out cross-validation method. Its predictive capability was further validated using a test set of 16 inhibitors similar to the training set of inhibitors. The correlation between the predicted and observed activities is significantly improved by the protein \"induced-fit\" that accounts for the flexibility of the receptor. These interaction and pharmacophore elements provide a unique insight to the rational design of new ligands targeted to the enzyme. 相似文献
18.
19.
Yudibeth Sixto-López Rolando Alberto Rodríguez-Fonseca Martha Cecilia Rosales-Hernández Marlet Martínez-Archundia José Antonio Gómez-Vidal 《Journal of biomolecular structure & dynamics》2017,35(13):2794-2814
Histone deacetylases (HDACs) are a family of proteins involved in the deacetylation of histones and other non-histones substrates. HDAC6 belongs to class II and shares similar biological functions with others of its class. Nevertheless, its three-dimensional structure that involves the catalytic site remains unknown for exploring the ligand recognition properties. Therefore, in this contribution, homology modeling, 100-ns-long Molecular Dynamics (MD) simulation and docking calculations were combined to explore the conformational complexity and binding properties of the catalytic domain 2 from HDAC6 (DD2-HDAC6), for which activity and affinity toward five different ligands have been reported. Clustering analysis allowed identifying the most populated conformers present during the MD simulation, which were used as starting models to perform docking calculations with five DD2-HDAC6 inhibitors: Cay10603 (CAY), Rocilinostat (RCT), Tubastatin A (TBA), Tubacin (TBC), and Nexturastat (NXT), and then were also submitted to 100-ns-long MD simulations. Docking calculations revealed that the five inhibitors bind at the DD2-HDAC6 binding site with the lowest binding free energy, the same binding mode is maintained along the 100-ns-long MD simulations. Overall, our results provide structural information about the molecular flexibility of apo and holo DD2-HDAC6 states as well as insight of the map of interactions between DD2-HDAC6 and five well-known DD2-HDAC6 inhibitors allowing structural details to guide the drug design. Finally, we highlight the importance of combining different theoretical approaches to provide suitable structural models for structure-based drug design. 相似文献
20.
Quantum mechanical/molecular mechanical free-energy simulations were performed to understand the deacylation reaction catalysed by sedolisin (a serine-carboxyl peptidase) and to elucidate the catalytic mechanism and the role of the active-site residues during the process. The results given here demonstrate that Asp170 may act as a general acid/base catalyst for the deacylation reaction. It is also shown that the electrostatic oxyanion hole interactions involving Asp170 may be less effective in transition state stabilisation for the deacylation step in the sedolisin-catalysed reaction compared to the general acid/base mechanism. The proton transfer processes during the enzyme-catalysed process were examined, and their role in the catalysis was discussed. 相似文献