Hydration structure and dynamics of B- and Z-DNA in the presence of counterions via molecular dynamics simulations

Following our previous attempts at understanding the structural and dynamical properties of water and counterions hydrating nucleic acids, we have performed molecular dynamics simulations for B- and Z-DNA. In these simulations, the nucleic acids were held rigid. In the case of B-DNA, one turn of B-DNA double helix was considered in the presence of 1500 water molecules and 20 counterions (K⁺). The simulations were performed for 4.0 ps after equilibrating the system. For Z-DNA, we considered one turn of the double helix in the presence of 1851 water molecules and 24 counterions (K⁺). The simulations were carried out for 3.5 ps after equilibration. The average temperature of these simulations was ～ 360 K for Z-DNA and ～ 345 K for B-DNA. In these simulations the hydrogen atoms were explicitly taken into account. For both simulations, a fifth-order predictor-corrector was used for solving the translational equations of motion. The rotational motion of the water molecules was represented in terms of quaternion algebra and the rotational equations of motion were solved with a second-order quaternion method using a sixth-order predictor-corrector method. A time step of 0.5 · 10^?15 s was used in these simulations. The structural and the dynamical properties of water solvating the counterions, and the phosphate groups of the DNA, were computed to understand the hydration structure. Diffusion coefficients and velocity correlation functions were calculated for both ions and the water molecules. The velocity correlation functions for the ions exhibit a caged behavior. The dipole correlation functions for the water molecules indicate that the water molecules close to the helix retain the memory of their initial orientations for longer periods of time than those away from the helix. During the time period of our simulation (3–4 ps) the ion probability distributions show a well-defined pattern and suggest limited mobility for the ions, being close to the helix.     

Molecular dynamics simulations investigation of neocarzinostatin chromophore-releasing pathways from the holo-NCS protein

The enediyne ring chromophore with strong DNA cleavage activity of neocarzinostatin is labile and therefore stabilization by forming the complex (carrying protein + chromophore: holo-NCS). Holo-NCS has gained much attention in clinical use as well as for drug delivery systems, but the chromophore-releasing mechanism to trigger binding to the target DNA with high affinity and producing DNA damage remain unclear. Three possible pathways were initially determined by conventional MD, essential dynamics and essential dynamics sampling. One of the paths runs along the naphthoate moiety; another runs along the amino sugar moiety; the third along the enediyne ring. Further, calculated forces and time by FPMD (force-probe molecular dynamics) suggest that the opening of the naphthoate moiety is most favorable pathway and Leu45, Phe76 and Phe78 all are key residues for chromophore release. In addition, conformational analyses indicate that the chromophore release is only local motions for the protein.     

Molecular dynamics simulations of the whey protein beta-lactoglobulin.

Molecular dynamics simulations have been used to model the motions and conformational behavior of the whey protein bovine beta-lactoglobulin. Simulations were performed for the protein by itself and complexed to a single retinol ligand located in a putative interior binding pocket. In the absence of the retinol ligand, the backbone loops around the opening of this interior pocket shifted inward to partially close off this cavity, similar to the shifts observed in previously reported molecular dynamics simulations of the uncomplexed form of the homologous retinol binding protein. The protein complexed with retinol does not exhibit the same conformational shifts. Conformational changes of this type could serve as a recognition signal allowing in vivo discrimination between the free and retinol complexed forms of the beta-lactoglobulin molecule. The unusual bending of the single alpha-helix observed in the simulations of retinol binding protein were not observed in the present calculations.     

Weighted-ensemble Brownian dynamics simulations for protein association reactions.

A new method, weighted-ensemble Brownian dynamics, is proposed for the simulation of protein-association reactions and other events whose frequencies of outcomes are constricted by free energy barriers. The method features a weighted ensemble of trajectories in configuration space with energy levels dictating the proper correspondence between "particles" and probability. Instead of waiting a very long time for an unlikely event to occur, the probability packets are split, and small packets of probability are allowed to diffuse almost immediately into regions of configuration space that are less likely to be sampled. The method has been applied to the Northrup and Erickson (1992) model of docking-type diffusion-limited reactions and yields reaction rate constants in agreement with those obtained by direct Brownian simulation, but at a fraction of the CPU time (10(-4) to 10(-3), depending on the model). Because the method is essentially a variant of standard Brownian dynamics algorithms, it is anticipated that weighted-ensemble Brownian dynamics, in conjunction with biophysical force models, can be applied to a large class of association reactions of interest to the biophysics community.     

Molecular dynamics simulations of the TrkH membrane protein

TrkH is a transmembrane protein that mediates uptake of K(+) through the cell membrane. Despite the recent determination of its crystallographic structure, the nature of the permeation mechanism is still unknown, that is, whether K(+) ions move across TrkH by active transport or passive diffusion. Here, molecular dynamics simulations and the umbrella sampling technique have been employed to shed light on this question. The existence of binding site S3 and two alternative binding sites have been characterized. Analysis of the coordination number renders values that are almost constant, with a full contribution from the carbonyls of the protein only at S3. This observation contrasts with observations of K(+) channels, where the contribution of the protein to the coordination number is roughly constant in all four binding sites. An intramembrane loop is found immediately after the selectivity filter at the intracellular side of the protein, which obstructs the permeation pathway, and this is reflected in the magnitude of the energy barriers.     

Molecular dynamics simulations of the photoactive protein nitrile hydratase

Nitrile hydratase (NHase) is an enzyme used in the industrial biotechnological production of acrylamide. The active site, which contains nonheme iron or noncorrin cobalt, is buried in the protein core at the interface of two domains, α and β. Hydrogen bonds between βArg-56 and αCys-114 sulfenic acid (αCEA114) are important to maintain the enzymatic activity. The enzyme may be inactivated by endogenous nitric oxide (NO) and activated by absorption of photons of wavelength λ < 630 nm. To explain the photosensitivity and to propose structural determinants of catalytic activity, differences in the dynamics of light-active and dark-inactive forms of NHase were investigated using molecular dynamics (MD) modeling. To this end, a new set of force field parameters for nonstandard NHase active sites have been developed. The dynamics of the photodissociated NO ligand in the enzyme channel was analyzed using the locally enhanced sampling method, as implemented in the MOIL MD package. A series of 1 ns trajectories of NHases shows that the protonation state of the active site affects the dynamics of the catalytic water and NO ligand close to the metal center. MD simulations support the catalytic mechanism in which a water molecule bound to the metal ion directly attacks the nitrile carbon.     

Electrophoresis in protein crystal: nonequilibrium molecular dynamics simulations

Electrophoresis of a mixture of NaCl and CaCl₂ in a lysozyme crystal is investigated using nonequilibrium molecular dynamics (MD) simulations. Upon exposure to an electric field, the stability of lysozyme is found to decrease slightly. This finding is demonstrated by increases in the root mean-square deviations of the heavy atoms of lysozyme, in the solvent-accessible surface area of hydrophobic residues, and in the number of hydrogen bonds between lysozyme and water. The solvent-accessible surface area of hydrophilic residues changes marginally, and the number of hydrogen bonds between lysozyme molecules decreases. Water molecules tend to align preferentially parallel to the electric field, and the dipole moment along the pore axis increases linearly with increasing field strength. Two pronounced layered structures are observed for Na⁺ and Ca²⁺ in the vicinity of protein surface, but only one enriched layer is observed for Cl⁻. The number distributions of all ions are nearly independent of the electric field. The water coordination numbers of all ions are smaller in the crystal than in aqueous bulk solution; however, the reverse is found for the Cl⁻ coordination numbers of cations. Both the water and the Cl⁻ coordination numbers are insensitive to the electric field. Ion diffusivities in the crystal are ∼2 orders of magnitude smaller than those in aqueous bulk solution. The drift velocities of ions increase proportionally to the electric field, particularly at high strengths, and depend on ionic charge and coordination with oppositely charged ions. Electrical current exhibits a linear relationship with the field strength. The zero-field electrical conductivity is estimated to be 0.56 S/m, which is very close to 0.61 S/m as predicted by the Nernst-Einstein equation.     

Molecular dynamics simulations of a protein on hydrophobic and hydrophilic surfaces.

Molecular dynamics simulations have been used to investigate the behavior of the peripheral membrane protein, cytochrome c, covalently tethered to hydrophobic (methyl-terminated) and hydrophilic (thiol-terminated) self-assembled monolayers (SAMs). The simulations predict that the protein will undergo minor structural changes when it is tethered to either surface, and the structures differ qualitatively on the two surfaces: the protein is less spherical on the hydrophilic SAM where the polar surface residues reach out to interact with the SAM surface. The protein is completely excluded from the hydrophobic SAM but partially dissolves in the hydrophilic SAM. Consequently, the surface of the thiol-terminated SAM is considerably less ordered than that of the methyl-terminated SAM, although a comparable, high degree of order is maintained in the bulk of both SAMs: the chains exhibit collective tilts in the nearest-neighbor direction at angles of 20 degrees and 17 degrees with respect to the surface normal in the hydrophobic and the hydrophilic SAMs, respectively. On the hydrophobic SAM the protein is oriented so that the heme plane is more nearly parallel to the surface, whereas on the hydrophilic surface it is more nearly perpendicular. The secondary structure of the protein, dominated by alpha helices, is not significantly affected, but the structure of the loops as well as the helix packing is slightly modified by the surfaces.     

Molecular dynamics simulations of protein-tyrosine phosphatase 1B. I. ligand-induced changes in the protein motions.

Activity of enzymes, such as protein tyrosine phosphatases (PTPs), is often associated with structural changes in the enzyme, resulting in selective and stereospecific reactions with the substrate. To investigate the effect of a substrate on the motions occurring in PTPs, we have performed molecular dynamics simulations of PTP1B and PTP1B complexed with a high-affinity peptide DADEpYL, where pY stands for phosphorylated tyrosine. The peptide sequence is derived from the epidermal growth factor receptor (EGFR988-993). Simulations were performed in water for 1 ns, and the concerted motions in the protein were analyzed using the essential dynamics technique. Our results indicate that the predominately internal motions in PTP1B occur in a subspace of only a few degrees of freedom. Upon substrate binding, the flexibility of the protein is reduced by approximately 10%. The largest effect is found in the protein region, where the N-terminal of the substrate is located, and in the loop region Val198-Gly209. Displacements in the latter loop are associated with the motions in the WPD loop, which contains a catalytically important aspartic acid. Estimation of the pKa of the active-site cysteine along the trajectory indicates that structural inhomogeneity causes the pKa to vary by approximately +/-1 pKa unit. In agreement with experimental observations, the active-site cysteine is negatively charged at physiological pH.     

Tryptophan dynamics of the FK506 binding protein: time-resolved fluorescence and simulations.

The FK506-binding protein (FKBP12) is important in the immunosuppressant action of FK506 and rapamycin. We have investigated Trp side chain dynamics in FKBP12, with and without a bound immunosuppressant, by measuring the Trp time-resolved fluorescence anisotropy decay r(t). The r(t) for W59 in aqueous uncomplexed FKBP12 at 20 degrees C is well described by a single exponential with a recovered initial anisotropy, r(eff)o, of 0.192 and an overall rotational correlation time for the protein, phi p, of 4.7 ns; r(eff)o = 0.214 and phi p = 4.2 ns for the FKBP12/FK506 complex. Using an expression for the order parameter squared, namely S2 = r(eff)o/rTo, where rTo is the vitrified steady-state excitation anisotropy, we recovered an S2 of 0.75 for W59 fluorescence in uncomplexed FKBP12 and S2 approximately equal to 1 in the FKBP12/FK506 complex. Results obtained for the FKBP12/rapamycin complex are similar to those found for the FKBP12/FK506 complex. Minimum perturbation mapping simulations were performed on the free and complexed forms of FKBP12 and the results were generally in agreement with the experimental data.     

Finite size effects in simulations of protein aggregation

It is becoming increasingly clear that the soluble protofibrillar species that proceed amyloid fibril formation are associated with a range of neurodegenerative disorders such as Alzheimer's and Parkinson diseases. Computer simulations of the processes that lead to the formation of these oligomeric species are starting to make significant contributions to our understanding of the determinants of protein aggregation. We simulate different systems at constant concentration but with a different number of peptides and we study the how the finite number of proteins affects the underlying free energy of the system and therefore the relative stability of the species involved in the process. If not taken into account, this finite size effect can undermine the validity of theoretical predictions regarding the relative stability of the species involved and the rates of conversion from one to the other. We discuss the reasons that give rise to this finite size effect form both a probabilistic and energy fluctuations point of view and also how this problem can be dealt by a finite size scaling analysis.     

Molecular dynamics simulations and membrane protein structure quality

Despite a growing repertoire of membrane protein structures (currently ∼120 unique structures), considerations of low resolution and crystallization in the absence of a lipid bilayer require the development of techniques to assess the global quality of membrane protein folds. This is also the case for assessment of, e.g. homology models of human membrane proteins based on structures of (distant) bacterial homologues. Molecular dynamics (MD) simulations may be used to help evaluate the quality of a membrane protein structure or model. We have used a structure of the bacterial ABC transporter MsbA which has the correct transmembrane helices but an incorrect handedness and topology of their packing to test simulation methods of quality assessment. An MD simulation of the MsbA model in a lipid bilayer is compared to a simulation of another bacterial ABC transporter, BtuCD. The latter structure has demonstrated good conformational stability in the same bilayer environment and over the same timescale (20 ns) as for the MsbA model simulation. A number of comparative analyses of the two simulations were performed to assess changes in the structural integrity of each protein. The results show a significant difference between the two simulations, chiefly due to the dramatic structural deformations of MsbA. We therefore propose that MD could become a useful quality control tool for membrane protein structural biology. In particular, it provides a way in which to explore the global conformational stability of a model membrane protein fold.     

Dominant-negative effects in prion diseases: insights from molecular dynamics simulations on mouse prion protein chimeras

Mutations in the prion protein (PrP) can cause spontaneous prion diseases in humans (Hu) and animals. In transgenic mice, mutations can determine the susceptibility to the infection of different prion strains. Some of these mutations also show a dominant-negative effect, thus halting the replication process by which wild type mouse (Mo) PrP is converted into Mo scrapie. Using all-atom molecular dynamics (MD) simulations, here we studied the structure of HuPrP, MoPrP, 10?Hu/MoPrP chimeras, and 1 Mo/sheepPrP chimera in explicit solvent. Overall, ～2?μs of MD were collected. Our findings suggest that the interactions between α1 helix and N-terminal of α3 helix are critical in prion propagation, whereas the β2–α2 loop conformation plays a role in the dominant-negative effect.

An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:4.     

Brownian dynamics simulations of probe and self-diffusion in concentrated protein and DNA solutions.

We have developed a Brownian dynamics algorithm for simulating probe and self-diffusion in concentrated solutions of DNA and protein. In these simulations, proteins are represented as spheres with radii given by their hydrodynamic radii, while DNA is modeled as a wormlike chain of hydrodynamically equivalent spherical frictional elements. The molecular interaction potentials employed by the program allow for intramolecular stretching and bending motions of the DNA chains, short-range Lennard-Jones interactions, and long-range electrostatic interactions. To test the program, we have carried out simulations of bovine serum albumin (BSA) probe diffusion and DNA self-diffusion in solutions of short-chain DNA as a function of both DNA concentration and solution ionic strength. In addition, we report on simulations of BSA self-diffusion as a function of BSA concentration and ionic strength. Based on a comparison to available experimental data, we find that our simulations accurately predict these transport properties under conditions of physiological salt concentration and predict the stronger concentration dependence observed at lower salt concentrations. These results are discussed in light of the nature of the intermolecular interactions in such systems and the approximations and limitations of the simulation algorithm.     

Time-resolved fluorescence spectroscopy and molecular dynamics simulations point out the effects of pressure on the stability and dynamics of the porcine odorant-binding protein

The effects of hydrostatic pressure on the structure and stability of porcine odorant-binding protein (pOBP) in the presence and absence of the odorant molecule 2-isobutyl-3-methoxypyrazine (IBMP) were studied by steady-state and time-resolved fluorescence spectroscopy as well as by molecular dynamics simulation. The authors found that the application of moderate values of hydrostatic pressure to pOBP solutions perturbed the microenvironment of Trp(16) and disrupted its highly quenched complex with Met(39). In addition, compared with the protein in the absence of IBMP, the MD simulations experiments carried out at different pressures highlighted the role of this ligand in stabilizing the Trp(16)/Met(39) interaction even at 2000 bar. The obtained results will assist for the tailoring of this protein as specific sensing element in a new class of fluorescence-based biosensors for the detection of explosives.     

分子动力学模拟尿素对水溶液中蛋白质构象转变的影响

采用分子动力学方法和全原子模型研究尿素和水分子对模型蛋白S-肽链结构转化的影响。模拟结果显示S-肽链的变性速率常数k值随着尿素浓度的增加而先降低后升高,在尿素浓度为2.9 mol/L时达到最低值。模拟了不同尿素浓度下尿素-肽链、水-肽链以及肽链分子氢键的形成状况。结果表明:尿素浓度较低时,尿素分子与S-肽链的极性氨基酸侧链形成氢键,但不破坏其分子内的骨架氢键,尿素在S-肽链水化层外形成限制性空间,增强了S-肽链的稳定性。随着尿素的升高,尿素分子进入S-肽链内部并与其内部氨基酸残基形成氢键,导致S-肽链的骨架氢键丧失,S-肽链发生去折叠。上述模拟结果与文献报道的实验结果一致,从分子水平上揭示了尿素对蛋白质分子结构变化的影响机制,对于研究和发展蛋白质折叠及稳定化技术具有指导意义。     

Computer simulations of protein folding by targeted molecular dynamics

We have performed 128 folding and 45 unfolding molecular dynamics runs of chymotrypsin inhibitor 2 (CI2) with an implicit solvation model for a total simulation time of 0.4 microseconds. Folding requires that the three-dimensional structure of the native state is known. It was simulated at 300 K by supplementing the force field with a harmonic restraint which acts on the root-mean-square deviation and allows to decrease the distance to the target conformation. High temperature and/or the harmonic restraint were used to induce unfolding. Of the 62 folding simulations started from random conformations, 31 reached the native structure, while the success rate was 83% for the 66 trajectories which began from conformations unfolded by high-temperature dynamics. A funnel-like energy landscape is observed for unfolding at 475 K, while the unfolding runs at 300 K and 375 K as well as most of the folding trajectories have an almost flat energy landscape for conformations with less than about 50% of native contacts formed. The sequence of events, i.e., secondary and tertiary structure formation, is similar in all folding and unfolding simulations, despite the diversity of the pathways. Previous unfolding simulations of CI2 performed with different force fields showed a similar sequence of events. These results suggest that the topology of the native state plays an important role in the folding process.