MD simulation study of several peptides including a polyalanine,a helix(pdb:2I9M),and a leucine zipper were carried out to investigate hydrogen bond energetics using dynamic polarized protein-specific charge(DPPC)to account for the polarization effect in protein dynamics.Results show that the backbone hydrogen-bond strength is generally correlated with its specific local electrostatic environment,measured by the number of water molecules near the hydrogen bond in the first solvation shell.The correlation coefficient is found to be 0.89,0.78,and 0.80,respectively,for polyalanine,2I9M protein,and leucine zipper.In the polyalanine,the energies of the backbone hydrogen bonds are very similar to each other due to their similar local electrostatic environment.The current study helps demonstrate and support the understanding that hydrogen bonds are stronger in a hydrophobic surrounding than in a hydrophilic one.For comparison,the result from simulation using standard force field shows a much weaker correlation between hydrogen bond energy and local electrostatic environment due to the lack of polarization effect in the force field.
DUAN Li LGAO YaJI Chang GMEI YeZHANG Qing GTANG BoZHANG John Z.H
The stereodynamics of the reaction of Ca + HCl are calculated at three different collision energies based on the potential energy surface [Verbockhaven G et al. 2005 J. Chem. Phys. 122 204307] using quasi-classical trajectory theory. The polarization-dependent differential cross sections (PDDCSs) (2π/σ )(dσ 00 /dω t ), (2π/σ )(dσ 20 /dωt ), (2π/σ )(dσ 22+ /dωt ), (2π/σ )(dσ 21 /dω t ) and the distributions of P(θ r ), P(φr ), and P(θr ,φr ) are calculated. The results indicate that the rotational polarization of the CaCl product presents different characteristics for the different collision energies, and the effects of the collision energy on the vector potential, including the alignment, orientation, and PDDCSs, are not obvious.