Steered molecular dynamics simulation of the binding of the bovine auxilin J domain to the Hsc70 nucleotide-binding domain

The Hsp70 and Hsp40 chaperone machine plays critical roles in protein folding, membrane translocation, and protein degradation by binding and releasing protein substrates in a process that utilizes ATP. The activities of the Hsp70 family of chaperones are recruited and stimulated by the J domains of Hsp40 chaperones. However, structural information on the Hsp40–Hsp70 complex is lacking, and the molecular details of this interaction are yet to be elucidated. Here we used steered molecular dynamics (SMD) simulations to investigate the molecular interactions that occur during the dissociation of the auxilin J domain from the Hsc70 nucleotide-binding domain (NBD). The changes in energy observed during the SMD simulation suggest that electrostatic interactions are the dominant type of interaction. Additionally, we found that Hsp70 mainly interacts with auxilin through the surface residues Tyr866, Arg867, and Lys868 of helix II, His874, Asp876, Lys877, Thr879, and Gln881 of the HPD loop, and Phe891, Asn895, Asp896, and Asn903 of helix III. The conservative residues Tyr866, Arg867, Lys868, His874, Asp876, Lys877, and Phe891 were also found in a previous study to be indispensable to the catalytic activity of the DnaJ J domain and the binding of it with the NBD of DnaK. The in silico identification of the importance of auxilin residues Asn895, Asp896, and Asn903 agrees with previous mutagenesis and NMR data suggesting that helix III of the J domain of the T antigen interacts with Hsp70. Furthermore, our data indicate that Thr879 and Gln881 from the HPD loop are also important as they mediate the interaction between the bovine auxilin J domain and Hsc70.     

Steered molecular dynamics studies of titin I1 domain unfolding

The cardiac muscle protein titin, responsible for developing passive elasticity and extensibility of muscle, possesses about 40 immunoglobulin-like (Ig) domains in its I-band region. Atomic force microscopy (AFM) and steered molecular dynamics (SMD) have been successfully combined to investigate the reversible unfolding of individual Ig domains. However, previous SMD studies of titin I-band modules have been restricted to I27, the only structurally known Ig domain from the distal region of the titin I-band. In this paper we report SMD simulations unfolding I1, the first structurally available Ig domain from the proximal region of the titin I-band. The simulations are carried out with a view toward upcoming atomic force microscopy experiments. Both constant velocity and constant force stretching have been employed to model mechanical unfolding of oxidized I1, which has a disulfide bond bridging beta-strands C and E, as well as reduced I1, in which the disulfide bridge is absent. The simulations reveal that I1 is protected against external stress mainly through six interstrand hydrogen bonds between its A and B beta-strands. The disulfide bond enhances the mechanical stability of oxidized I1 domains by restricting the rupture of backbone hydrogen bonds between the A'- and G-strands. The disulfide bond also limits the maximum extension of I1 to approximately 220 A. Comparison of the unfolding pathways of I1 and I27 are provided and implications to AFM experiments are discussed.     

Steered molecular dynamics simulation of the Rieske subunit motion in the cytochrome bc(1) complex.

Crystallographic structures of the mitochondrial ubiquinol/cytochrome c oxidoreductase (cytochrome bc(1) complex) suggest that the mechanism of quinol oxidation by the bc(1) complex involves a substantial movement of the soluble head of the Rieske iron-sulfur protein (ISP) between reaction domains in cytochrome b and cytochrome c(1) subunits. In this paper we report the results of steered molecular dynamics simulations inducing, through an applied torque within 1 ns, a 56 degrees rotation of the soluble domain of ISP. For this purpose, a solvated structure of the bc(1) complex in a phospholipid bilayer (a total of 206,720 atoms) was constructed. A subset of 91,061 atoms was actually simulated with 45,131 moving atoms. Point charge distributions for the force field parametrization of heme groups and the Fe(2)S(2) cluster of the Rieske protein included in the simulated complex were determined. The simulations showed that rotation of the soluble domain of ISP is actually feasible. Several metastable conformations of the ISP during its rotation were identified and the interactions stabilizing the initial, final, and intermediate positions of the soluble head of the ISP domain were characterized. A pathway for proton conduction from the Q(o) site to the solvent via a water channel has been identified.     

Prediction of the binding mode between BMS-378806 and HIV-1 gp120 by docking and molecular dynamics simulation

BMS-378806 is a newly discovered small molecule that effectively blocks the binding of CD4 with gp120. The binding mode of this kind of inhibitor remains unknown. In this paper, AutoDock 3.0 in conjunction with molecular dynamics simulation, accommodating the receptor's flexibility, was used to explore the binding mode between BMS-378806 and gp120. Two structures, Mode I and Mode II, with the lowest docking energy were selected as different representative binding modes. The analysis of the results from the molecular dynamics simulation indicated that the binding of BMS-348806 in Mode II is more stable. The average structure of Mode II was analyzed and compared with the experimental data. The conclusion was that BMS-378806 inserts the azaindole ring deeply into the PHE43 cavity and makes contact with a number of residues in the cavity, on the cavity and near the cavity. This study benefits the understanding of the mechanism of this kind of inhibitor and may provide useful information for rational drug design.     

Improved prediction of HIV-1 protease-inhibitor binding energies by molecular dynamics simulations

Background  

The accurate prediction of enzyme-substrate interaction energies is one of the major challenges in computational biology. This study describes the improvement of protein-ligand binding energy prediction by incorporating protein flexibility through the use of molecular dynamics (MD) simulations.     

Steered molecular dynamics and mechanical functions of proteins

Atomic force microscopy of single molecules, steered molecular dynamics and the theory of stochastic processes have established a new field that investigates mechanical functions of proteins, such as ligand-receptor binding/unbinding and elasticity of muscle proteins during stretching. The combination of these methods yields information on the energy landscape that controls mechanical function and on the force-bearing components of proteins, as well as on the underlying physical mechanisms.     

Hsc70-Auxilin复合物的拉伸分子动力学模拟

Hsc70与auxilin蛋白组成的系统是Hsp70/Hsp40分子伴侣系统家族的一员,在热休克反应中发挥重要作用。本文为得出auxilin蛋白J结构域的关键氨基酸,首先采用由二硫键交联的Hsc70 R171C与auxilin D876C的复合物结晶结构作为初始模型,进行分子动力学模拟,通过比较平衡后的结合部位发现,将形成二硫键的氨基酸突变为原来的氨基酸结构在结合位点上与生化结果较为相近,之后利用此结构通过拉伸动力学模拟分析了auxilin蛋白J结构域与Hsc70的ATPase功能域的解离过程,并探讨了Hsc70与auxilin蛋白之间的相互作用力。结果表明位于HPD loop上的His874,Asp876,Thr879,螺旋Ⅲ上的Glu884,Asn895,Asp896,Ser899,Glu902,Asn903为关键氨基酸,这些数据符合之前核磁共振实验证实的T抗原J结构域的HPD基序和螺旋Ⅲ与Hsc70的ATPase功能域之间的相互作用。     

Conformational dynamics of RNA-peptide binding: a molecular dynamics simulation study

Molecular dynamics simulations of the binding of the heterochiral tripeptide KkN to the transactivation responsive (TAR) RNA of HIV-1 is presented, using an all-atom force field with explicit water. To obtain starting structures for the TAR-KkN complex, semirigid docking calculations were performed that employ an NMR structure of free TAR RNA. The molecular dynamics simulations show that the starting structures in which KkN binds to the major groove of TAR (as it is the case for the Tat-TAR complex of HIV-1) are unstable. On the other hand, the minor-groove starting structures are found to lead to several binding modes, which are stabilized by a complex interplay of stacking, hydrogen bonding, and electrostatic interactions. Although the ligand does not occupy the binding position of Tat protein, it is shown to hinder the interhelical motion of free TAR RNA. The latter is presumably necessary to achieve the conformational change of TAR RNA to bind Tat protein. Considering the time evolution of the trajectories, the binding process is found to be ligand-induced and cooperative. That is, the conformational rearrangement only occurs in the presence of the ligand and the concerted motion of the ligand and a large part of the RNA binding site is necessary to achieve the final low-energy binding state.     

HIV-1 protease substrate binding and product release pathways explored with coarse-grained molecular dynamics

We analyze the encounter of a peptide substrate with the native HIV-1 protease, the mechanism of substrate incorporation in the binding cleft, and the dissociation of products after substrate hydrolysis. To account for the substrate, we extend a coarse-grained model force field, which we previously developed to study the flap opening dynamics of HIV-1 protease on a microsecond timescale. Molecular and Langevin dynamics simulations show that the flaps need to open for the peptide to bind and that the protease interaction with the substrate influences the flap opening frequency and interval. On the other hand, release of the products does not require flap opening because they can slide out from the binding cleft to the sides of the enzyme. Our data show that in the protease-substrate complex the highest fluctuations correspond to the 17- and 39-turns and the substrate motion is anticorrelated with the 39-turn. Moreover, the active site residues and the flap tips move in phase with the peptide. We suggest some mechanistic principles for how the flexibility of the protein may be involved in ligand binding and release.     

Conformational stability of PCID2 upon DSS1 binding with molecular dynamics simulation

DSS1 is a small acidic intrinsically disordered protein (IDP) that can fold upon binding with PCID2 TREX-2. The resulting complex plays a key role in mRNA export. However, the binding mechanism between DSS1 and PCID2 is unsolved. Here, three independent 500-ns molecular dynamics (MD) simulations were performed to study the DSS1–PCID2 binding mechanism by comparing apo-PCID2 and bound PCID2. The results show that the conformational variation of bound PCID2 is smaller than that of apo-PCID2, especially in the binding domain of two helices (helix IV and VIII). The probability of coil formation between helix III and helix IV of bound PCID2 increases, and a short anti-parallel β-sheet forms upon DSS1 binding. The decomposition of binding free energy into protein and residue pairs suggests that electrostatic and hydrophobic interactions play key roles in the recognition between DSS1 and PCID2. There is a hydrophobic core of seven residues in DSS1 favorable to the binding of PCID2. These analytical methods can be used to reveal the recognition mechanisms of other IDPs and their partners.

Open image in new window

Figure Differences of second structure of PCID2 in bound and unbound states. The interaction surface between the helix VIII of PCID2 and helix of DSS1

     

Investigation of the effect of homocysteinylation of substance P on its binding to the NK1 receptor using molecular dynamics simulation

Substance P is a neurotransmitter or modulator in both the central and peripheral nervous systems. In this work, modifications of the lysine in SP by homocysteine and an acetyl group as well as the conformational dynamics of the native and modified SP peptides and their complexes with the NK1 receptor were studied via MD simulation. It was found that modifying SP stabilizes the peptide structure, but the modified SP peptides are less likely to bind to the NK1 receptor, so the resulting complexes are less stable. The RMSD of native SP (~0.33 nm) is about twice as large as that of the modified SP peptides (~0.18 nm), while the RMSD for the receptor complexed with native SP is ~0.3 nm, and that for the receptor complexed with either of the modified peptides is ~0.35 nm, which demonstrates the high stability of the modified SP peptides as well as the receptor complexed with native SP. Such behavior was also observed in other structural analyses. The binding free energies of the native and modified SP peptides with the NK1 receptor were also compared. The ΔG_bind values for the binding of homocysteinylated SP to the NK1 receptor and the binding of the acetylated SP and native SP to the NK1 receptor were ?38.89, ?64.46, and???264.52 kJ mol^?1, respectively. Modification of the lysine of SP decreases the binding affinity of the peptide to the NK1 receptor. In other words, homocysteinylation or acetylation of SP leads to weaker interactions of the peptide with the NK1 receptor compared to those between native SP and NK1. We propose that this phenomenon leads to increased levels of homocysteinylated SP in plasma in many diseases such as breast cancer.

Open image in new window

Graphical abstract Substance P (SP) is a neuropeptide which binds to the NK1 receptor. SP is of great pharmacological interest, as agonists and antagonists of SP can potentially be used to treat many chronic diseases. Therefore, in this work, the lysine (LYS) in SP was theoretically modified with a homocysteine or acetyl group to explore the effects of such a modification on the binding affinity of this peptide with the NK1 receptor and the structural dynamics of the resulting complex

     

Computational investigation of the HIV-1 Rev multimerization using molecular dynamics simulations and binding free energy calculations

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.     

Understanding the molecular mechanism of the broad and potent neutralization of HIV-1 by antibody VRC01 from the perspective of molecular dynamics simulation and binding free energy calculations

VRC01 is one of the most broadly and potently neutralizing HIV-1 antibodies known-it has been shown to neutralize 91?% of the tested primary isolate Env pseudoviruses by recognizing the viral envelope glycoprotein gp120. To explore the mechanism of HIV-1 neutralization by VRC01 and thus obtain valuable information for vaccine design, we performed molecular dynamics simulations and binding free energy calculations for apo-VRC01, apo-gp120, and the gp120-VRC01 complex. For gp120, residue energy decomposition analysis showed that the hotspot residues Asn280, Lys282, Asp368, Ile371, and Asp457 are located in three primary loops, including the CD4-binding loop, loop D, and loop V5. For VRC01, the hotspot residues Trp47, Trp50, Asn58, Arg61, Gln64, Trp100, and Tyr91 mainly come from CDR2 of the heavy chain. By decomposing the binding free energy into different components, intermolecular van der Waals interactions and nonpolar solvation were found to dominate the binding process. Principal component analysis of loops D and V5, which are related to neutralization resistance, indicated that these two areas have a larger conformational space in apo-gp120 compared to bound gp120. A comparison of three representative structures from the cluster analysis of loops D and V5 indicated that changes primarily occur at the tip of loop V5, and are caused by fluctuations in the terminal Glu1 residue of the antibody. This information can be used to guide the design of vaccines and small molecule inhibitors.     

Effects of ligand binding on the mechanical stability of protein GB1 studied by steered molecular dynamics simulation

Regulation of the mechanical properties of proteins plays an important role in many biological processes, and sheds light on the design of biomaterials comprised of protein. At present, strategies to regulate protein mechanical stability focus mainly on direct modulation of the force-bearing region of the protein. Interestingly, the mechanical stability of GB1 can be significantly enhanced by the binding of Fc fragments of human IgG antibody, where the binding site is distant from the force-bearing region of the protein. The mechanism of this long-range allosteric control of protein mechanics is still elusive. In this work, the impact of ligand binding on the mechanical stability of GB1 was investigated using steered molecular dynamics simulation, and a mechanism underlying the enhanced protein mechanical stability is proposed. We found that the external force causes deformation of both force-bearing region and ligand binding site. In other words, there is a long-range coupling between these two regions. The binding of ligand restricts the distortion of the binding site and reduces the deformation of the force-bearing region through a long-range allosteric communication, which thus improves the overall mechanical stability of the protein. The simulation results are very consistent with previous experimental observations. Our studies thus provide atomic-level insights into the mechanical unfolding process of GB1, and explain the impact of ligand binding on the mechanical properties of the protein through long-range allosteric regulation, which should facilitate effective modulation of protein mechanical properties.     

Insight into binding mechanisms of inhibitors MKP56, MKP73, MKP86, and MKP97 to HIV-1 protease by using molecular dynamics simulation

HIV-1 protease (PR) has been a significant target for design of potent inhibitors curing acquired immunodeficiency syndrome. Molecular dynamics simulations coupled with molecular mechanics Poisson–Boltzmann surface area method were performed to study interaction modes of four inhibitors MKP56, MKP73, MKP86, and MKP97 with PR. The results suggest that the main force controlling interactions of inhibitors with PR should be contributed by van der Waals interactions between inhibitors and PR. The cross-correlation analyses based on MD trajectories show that inhibitor binding produces significant effect on the flap dynamics of PR. Hydrogen bond analyses indicate that inhibitors can form stable hydrogen bonding interactions with the residues from the catalytic strands of PR. The contributions of separate residues to inhibitor bindings are evaluated by using residue-based free energy decomposition method and the results demonstrate that the CH–π and CH–CH interactions between the hydrophobic groups of inhibitors with residues drive the associations of inhibitors with PR. We expect that this study can provide a significant theoretical aid for design of potent inhibitors targeting PR.     

Pharmacophoric features of nucleosidic HIV-1RT inhibitors

A quantum pharmacological study on nucleosidic inhibitors of HIV-1RT has been performed. The main aim of this study is to discuss the pharmacophoric features (conformational and electrostatic) of nucleosidic inhibitors and compare them with normal substrate dNTP. Present study stresses on the need to refine nucleosidic drugs, as combination therapy to date is still one of the best remedies for AIDS. The results of ab initio HF calculations indicate very little effect of 3' substituent on ring puckering and suggest that potency regulation may be via very intricate phosphate-catalytic triad interactions. Our MESP maps also show charge complementarity between the drug and receptor.