Quantitative Assessment of Force Fields on Both Low-Energy Conformational Basins and Transition-State Regions of the (ϕ-ψ) Space

Correct folding of an α-helix and a β-hairpin using a polarized 2D torsional potential

A new modification to the AMBER force field that incorporates the coupled two-dimensional main chain torsion energy has been evaluated for the balanced representation of secondary structures. In this modified AMBER force field (AMBER03(2D)), the main chain torsion energy is represented by 2-dimensional Fourier expansions with parameters fitted to the potential energy surface generated by high-level quantum mechanical calculations of small peptides in solution. Molecular dynamics simulations are performed to study the folding of two model peptides adopting either α-helix or β-hairpin structures. Both peptides are successfully folded into their native structures using an AMBER03(2D) force field with the implementation of a polarization scheme (AMBER03(2D)p). For comparison, simulations using a standard AMBER03 force field with and without polarization, as well as AMBER03(2D) without polarization, fail to fold both peptides successfully. The correction to secondary structure propensity in the AMBER03 force field and the polarization effect are critical to folding Trpzip2; without these factors, a helical structure is obtained. This study strongly suggests that this new force field is capable of providing a more balanced preference for helical and extended conformations. The electrostatic polarization effect is shown to be indispensable to the growth of secondary structures.

Electrostatic frequency shifts in amide I vibrational spectra: direct parameterization against experiment

The interpretation of protein amide I infrared spectra has been greatly assisted by the observation that the vibrational frequency of a peptide unit reports on its local electrostatic environment. However, the interpretation of spectra remains largely qualitative due to a lack of direct quantitative connections between computational models and experimental data. Here, we present an empirical parameterization of an electrostatic amide I frequency map derived from the infrared absorption spectra of 28 dipeptides. The observed frequency shifts are analyzed in terms of the local electrostatic potential, field, and field gradient, evaluated at sites near the amide bond in molecular dynamics simulations. We find that the frequency shifts observed in experiment correlate very well with the electric field in the direction of the C=O bond evaluated at the position of the amide oxygen atom. A linear best-fit mapping between observed frequencies and electric field yield sample standard deviations of 2.8 and 3.7 cm(-1) for the CHARMM27 and OPLS-AA force fields, respectively, and maximum deviations (within our data set) of 9 cm(-1). These results are discussed in the broader context of amide I vibrational models and the effort to produce quantitative agreement between simulated and experimental absorption spectra.

A coupled two-dimensional main chain torsional potential for protein dynamics: generation and implementation

Main chain torsions of alanine dipeptide are parameterized into coupled 2-dimensional Fourier expansions based on quantum mechanical (QM) calculations at M06 2X/aug-cc-pvtz//HF/6-31G** level. Solvation effect is considered by employing polarizable continuum model. Utilization of the M06 2X functional leads to precise potential energy surface that is comparable to or even better than MP2 level, but with much less computational demand. Parameterization of the 2D expansions is against the full main chain torsion space instead of just a few low energy conformations. This procedure is similar to that for the development of AMBER03 force field, except unique weighting factor was assigned to all the grid points. To avoid inconsistency between quantum mechanical calculations and molecular modeling, the model peptide is further optimized at molecular mechanics level with main chain dihedral angles fixed before the calculation of the conformational energy on molecular mechanical level at each grid point, during which generalized Born model is employed. Difference in solvation models at quantum mechanics and molecular mechanics levels makes this parameterization procedure less straightforward. All force field parameters other than main chain torsions are taken from existing AMBER force field. With this new main chain torsion terms, we have studied the main chain dihedral distributions of ALA dipeptide and pentapeptide in aqueous solution. The results demonstrate that 2D main chain torsion is effective in delineating the energy variation associated with rotations along main chain dihedrals. This work is an implication for the necessity of more accurate description of main chain torsions in the future development of ab initio force field and it also raises a challenge to the development of quantum mechanical methods, especially the quantum mechanical solvation models.

Sequence-specific recognition of cancer drug-DNA adducts by HMGB1a repair protein

The efficacy of cancer drugs such as cisplatin (Cp) and oxaliplatin (Ox), which covalently bind to DNA to form drug-DNA adducts, is linked to their recognition by repair proteins such as HMGB1a. Previous experimental studies showed that HMGB1a's binding affinity for Cp- and Ox-DNA varies with the drug used and the local DNA sequence context of the adduct. We link this differential binding affinity to the free energy of deforming (bending and minor groove opening) the drug-DNA molecule during HMGB1a binding. Specifically, the minimal binding affinity of HMGB1a for Ox-DNA in the TGGA context is explained by its larger deformation free energy compared with Cp-DNA or Ox-DNA in other sequence contexts. Methyl groups on neighboring thymine bases in Ox-TGGA crowd the minor groove and sterically hinder the motion of the diaminocyclohexane ring of Ox, leading to this reduced deformability and resultant decrease in HMGB1a's binding affinity.

Assessment of the intrinsic conformational preferences of dipeptide amino acids in aqueous solution by combined umbrella sampling/MBAR statistics. A comparison with experimental results

The propensities of 19 amino acid dipeptides have been calculated by a distributed umbrella sampling molecular dynamics simulation procedure using the OPLS-AA force field. The potential of mean force maps was estimated with the multiple Bennett acceptance ratio statistics. The resulting propensities compare satisfactorily well with very recently published experimental data on equivalent systems. In particular, α conformation-probabilities for all of the dipeptides remain much lower than either β or P(II) propensities. This result is in agreement with most experimental data for dipeptides. However, it is also in contrast with most simulation studies performed so far with other force fields, where α conformations result even more probable than P(II) or β ones. We discuss the behavior of the OPLS-AA force field, which can be useful for the improvement of this model in reproducing the recent experimental observations on amino acid dipeptides.