共查询到20条相似文献，搜索用时 0 毫秒
1.
Leukocyte function associated antigen-1 (LFA-1) plays a critical role in T cell migration and has been recognized as a therapeutic target for immune disorders. Several classes of small molecule antagonists have been developed to block LFA-1 interaction with intercellular adhesion molecule-1 (ICAM-1). Recent structural studies show that the antagonists bind to an allosteric site in the I-domain of LFA-1. However, it is not yet clear how these small molecules work as antagonists since no significant conformational change is observed in the I-domain-antagonist complex structures. Here we present a computational study suggesting how these allosteric antagonists affect the dynamics of the I-domain. The lowest frequency vibrational mode calculated from an LFA-1 I-domain structure shows large scale \"coil-down\" motion of the C-terminal alpha7 helix, which may lead to the open form of the I-domain. The presence of an allosteric antagonist greatly reduces this motion of the alpha7 helix as well as other parts of the I-domain. Thus, our study suggests that allosteric antagonists work by eliminating breathing motion that leads to the open conformation of the I-domain. 相似文献
2.
Antiestradiol antibody 57-2 binds 17beta-estradiol (E2) with moderately high affinity (K(a) = 5 x 10(8) M(-1)). The structurally related natural estrogens estrone and estriol as well synthetic 17-deoxy-estradiol and 17alpha-estradiol are bound to the antibody with 3.7-4.9 kcal mol(-1) lower binding free energies than E2. Free energy perturbation (FEP) simulations and the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method were applied to investigate the factors responsible for the relatively low cross-reactivity of the antibody with these four steroids, differing from E2 by the substituents of the steroid D-ring. In addition, computational alanine scanning of the binding site residues was carried out with the MM-PBSA method. Both the FEP and MM-PBSA methods reproduced the experimental relative affinities of the five steroids in good agreement with experiment. On the basis of FEP simulations, the number of hydrogen bonds formed between the antibody and steroids, which varied from 0 to 3 in the steroids studied, determined directly the magnitude of the steroid-antibody interaction free energies. One hydrogen bond was calculated to contribute about 3 kcal mol(-1) to the interaction energy. Because the relative binding free energies of estrone (two antibody-steroid hydrogen bonds), estriol (three hydrogen bonds), 17-deoxy-estradiol (no hydrogen bonds), and 17alpha-estradiol (two hydrogen bonds) are close to each other and clearly lower than that of E2 (three hydrogen bonds), the water-steroid interactions lost upon binding to the antibody make an important contribution to the binding free energies. The MM-PBSA calculations showed that the binding of steroids to the antiestradiol antibody is driven by van der Waals interactions, whereas specificity is solely due to electrostatic interactions. In addition, binding of steroids to the antiestradiol antibody 57-2 was compared to the binding to the antiprogesterone antibody DB3 and antitestosterone antibody 3-C4F5, studied earlier with the MM-PBSA method. 相似文献
3.
p53 is a homotetrameric tumor suppressor protein that is found to be mutated in most human cancers. Some of these mutations, particularly mutations to R337, fall in the tetramerization domain and cause defects in tetramer formation leading to loss of function. Mutation to His at this site has been found to destabilize the tetramer in a pH-dependent fashion. In structures of the tetramerization domain determined by crystallography, R337 from one monomer makes a salt bridge with D352 from another monomer, apparently helping to stabilize the tetramer. Here we present molecular dynamics simulations of wild-type p53 and the R337His mutant at several different pH and salt conditions. We find that the 337-352 salt bridge is joined by two other charged side chains, R333 and E349. These four residues do not settle into a fixed pattern of salt bridging, but continue to exchange salt-bridging partners on the nanosecond time scale throughout the simulation. This unusual system of fluid salt bridging may explain the previous finding from alanine scanning experiments that R333 contributes significantly to protein stability, even though in the crystal structure it is extended outward into solvent. This extended conformation of R333 appears to be the result of a specific crystal contact and, this contact being absent in the simulation, R333 turns inward to join its interaction partners. When R337 is mutated to His but remains positively charged, it maintains the original interaction with D352, but the newly observed interaction with E349 is weakened, accounting for the reduced stability of R337H even under mildly acidic conditions. When this His is deprotonated, the interaction with D352 is also lost, accounting for the further destabilization observed under mildly alkaline conditions. Simulations were carried out using both explicit and implicit solvent models, and both displayed similar behavior of the fluid salt-bridging cluster, suggesting that implicit solvent models can capture at least the qualitative features of this phenomenon as well as explicit solvent. Simulations under strongly acidic conditions in implicit solvent displayed the beginnings of the unfolding process, a destabilization of the hydrophobic dimer-dimer interface. Computational alanine scanning using the molecular mechanics Poisson-Boltzmann surface area method showed significant correlation to experimental unfolding data for charged and polar residues, but much weaker correlation for hydrophobic residues. 相似文献
4.
We have applied the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method (J. Srinivasan, T. E. Cheatham, P. Cieplak, P. A. Kollman, and D. A. Case, Journal of the American Chemical Society, 1998, Vol. 120, pp. 9401-9409) to study the interaction of an RNA aptamer with theophylline and its analogs. The MM-PBSA free energy analysis provides a reasonable absolute binding free energy for the RNA aptamer-theophylline complex formation. Energetic analysis reveals that the van der Waals interaction and the nonpolar contribution to solvation provide the basis for the favorable absolute free energy of complex. This trend is similar to other protein-ligand interactions studied previously. The MM-PBSA method also ranks the relative binding energies of five theophylline analogs approximately correctly, but not as well as the more conventional thermodynamic integration calculations, which were carried out to convert theophylline into its analogs. The comparison of MM-PBSA with TI suggests that the MM-PBSA method has some difficulties with the first-solvation-shell energetics. 相似文献
5.
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. 相似文献
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.
The cytolethal distending toxin (CDT) is a widespread bacterial toxin that consists of an active subunit CdtB with nuclease activity and two ricin-like lectin domains, CdtA and CdtC, that are involved in the delivery of CdtB into the host cell. The three subunits form a tripartite complex that is required to achieve the fully active holotoxin. In the present study we investigate the assembly and dynamic properties of the CDT holotoxin using molecular dynamics simulations and binding free energy calculations. The results have revealed that CdtB likely adopts a different conformation in the unbound state with a closed DNA binding site. The two characterized structural elements of the aromatic patch and groove on the CdtA and CdtC protein surfaces exhibit high mobility, and free energy calculations show that the heterodimeric complex CdtA-CdtC, as well as the CdtA-CdtB and CdtB-CdtC sub-complexes are less energetically stable as compared to the binding in the tripartite complex. Analysis of the dynamical cross-correlation map reveals information on the correlated motions and long-range interplay among the CDT subunits associated with complex formation. Finally, the estimated binding free energies of subunit interactions are presented, together with the free energy decomposition to determine the contributions of residues for both binding partners, providing insight into the protein-protein interactions in the CDT holotoxin. 相似文献
8.
Ming Liu Xiao Jing Cong Ping Li Jian Jun Tan Wei Zu Chen Cun Xin Wang 《Biopolymers》2009,91(9):700-709
Human immunodeficiency virus type 1 integrase (IN) is an essential enzyme in the life cycle of this virus and also an important target for the study of anti‐HIV drugs. In this work, the binding modes of the wild type IN core domain and the two mutants, that is, W132G and C130S, with the 4‐hydroxycoumarin compound NSC158393 were evaluated by using the “relaxed complex” molecular docking approach combined with molecular dynamics (MD) simulations. Based on the monomer MD simulations, both of the two substitutions affect not only the stability of the 128–136 peptides, but also the flexibility of the functional 140s loop. In principle, NSC158393 binds the 128–136 peptides of IN; however, the specific binding modes for the three systems are various. According to the binding mode of NSC158393 with WT, NSC158393 can effectively interfere with the stability of the IN dimer by causing a steric hindrance around the monomer interface. Additionally, through the comparative analysis of the MD trajectories of the wild type IN and the IN‐NSC158393 complex, we found that NSC15893 may also exert its inhibitory function by diminishing the mobility of the function loop of IN. Three key binding residues, that is, W131, K136, and G134, were discovered by energy decomposition calculated with the Molecular Mechanics Generalized Born Surface Area method. Characterized by the largest binding affinity, W131 is likely to be indispensable for the ligand binding. All the above results are consistent with experiment data, providing us some helpful information for understanding the mechanism of the coumarin‐based inhibitors. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 700–709, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com 相似文献
9.
10.
The pH-driven opening and closure of beta-lactoglobulin EF loop, acting as a lid and closing the internal cavity of the protein, has been studied by molecular dynamics (MD) simulations and free energy calculations based on molecular mechanics/Poisson-Boltzmann (PB) solvent-accessible surface area (MM/PBSA) methodology. The forms above and below the transition pH differ presumably only in the protonation state of residue Glu89. MM/PBSA calculations are able to reproduce qualitatively the thermodynamics of the transition. The analysis of MD simulations using a combination of MM/PBSA methodology and the colony energy approach is able to highlight the driving forces implied in the transition. The analysis suggests that global rearrangements take place before the equilibrium local conformation is reached. This conclusion may bear general relevance to conformational transitions in all lipocalins and proteins in general. 相似文献
11.
Histone acetylation is a very important regulatory mechanism in gene expression in the chromatin context. A new protein family‐YEATS domains have been found as a novel histone acetylation reader, which could specific recognize the histone lysine acetylation. AF9 is an important one in the YEATS family. Focused on the AF9‐H3K9ac (K9 acetylation) complex (ALY) (PDB code: 4TMP) and a serials of mutants, MUT (the acetyllsine of H3K9ac was mutated to lysine), F59A, G77A, and D103A, we applied molecular dynamics simulation and molecular mechanics Poisson?Boltzmann (MM‐PBSA) free energy calculations to examine the role of AF9 protein in recognition interaction. The simulation results and analysis indicate that some residues of the protein have significant influence on recognition and binding to H3K9ac peptides and hydrophobic surface show the hydrophobic interactions play an important role in the binding. Our work can give important information to understand how the protein AF9 recognizes the peptides H3K9ac. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 779–786, 2016. 相似文献
12.
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. 相似文献
13.
14.
Elena Moreno‐Castillo Yoanna M. lvarez‐Ginarte Mario E. Valds‐Tresanco Luis A. Montero‐Cabrera Ernesto Moreno Pedro A. Valiente 《Journal of molecular recognition : JMR》2020,33(9)
Alzheimer's disease is a progressive neurodegenerative disorder characterized by the abnormal processing of the Tau and the amyloid precursor proteins. The unusual aggregation of Tau is based on the formation of intermolecular β‐sheets through two motifs: 275VQIINK280 and 306VQIVYK311. Phenylthiazolyl‐hydrazides (PTHs) are capable of inhibiting/disassembling Tau aggregates. However, the disaggregation mechanism of Tau oligomers by PTHs is still unknown. In this work, we studied the disruption of the oligomeric form of the Tau motif 306VQIVYK311 by PTHs through molecular docking, molecular dynamics, and free energy calculations. We predicted hydrophobic interactions as the major driving forces for the stabilization of Tau oligomer, with V306 and I308 being the major contributors. Nonpolar component of the binding free energy is essential to stabilize Tau‐PTH complexes. PTHs disrupted mainly the van der Waals interactions between the monomers, leading to oligomer destabilization. Destabilization of full Tau filament by PTHs and emodin was not observed in the sampled 20 ns; however, in all cases, the nonpolar component of the binding free energy is essential for the formation of Tau filament‐PTH and Tau filament‐emodin. These results provide useful clues for the design of more effective Tau‐aggregation inhibitors. 相似文献
15.
16.
17.
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. 相似文献
18.
Meng-Nan Jiang Xiao-Ping Zhou Dong-Ru Sun Huan Gao Qing-Chuan Zheng Hong-Xing Zhang 《Journal of biomolecular structure & dynamics》2013,31(14):3705-3717
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, R2, and R2ADJ, proved to be credible. The cross-validations for the models (q2LOO = 0.571 and 0.629, respectively) showed their robustness, while the external validations (r2test = 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. 相似文献
19.
We describe the application of the molecular dynamics (MD) and molecular mechanics-generalized Born/surface area (MM-GB/SA) approaches to the simulation of the different biological activity of diethylstilbestrol (DES) on two highly homologous nuclear receptors-estrogen receptor alpha (ER-alpha) and estrogen-related receptor gamma (ERR-gamma). DES exerts an agonistic effect against ER-alpha and an antagonistic effect against ERR-gamma. Using the x-ray crystal structures of ER-alpha in the canonical agonist bound form (PDB code: 3ERD) and antagonist bound form (PDB code: 3ERT), ERR-gamma homology models have been constructed for the receptor in two different conformations. MM-GB/SA binding free energy calculations of DES in the ER-alpha and ERR-gamma structures suggest that DES exhibits a greater free energy of binding in the agonist bound conformation of ER-alpha, while the antagonist bound conformation is preferred for ERR-gamma. Further dissection of the free energy contributions coupled with calculation of the ligand binding pocket volume suggests that the van der Waals interactions for DES within the smaller binding pocket volume of ERR-gamma are less favorable and this is the main factor for DES antagonism in ERR-gamma. This approach has potential general applicability to the prediction of the biological activity of nuclear receptor ligands. 相似文献
20.