Optimization of the CHARMM Additive Force Field for DNA: Improved Treatment of the BI/BII Conformational Equilibrium

Citations Over TimeTop 10% of 2011 papers

Abstract

The B-form of DNA can populate two different backbone conformations: BI and BII, defined by the difference between the torsion angles ε and ζ (BI = ε-ζ 0). BI is the most populated state, but the population of the BII state, which is sequence dependent, is significant and accumulating evidence shows that BII affects the overall structure of DNA, and thus influences protein-DNA recognition. This work presents a reparametrization of the CHARMM27 additive nucleic acid force field to increase the sampling of the BII form in MD simulations of DNA. In addition, minor modifications of sugar puckering were introduced to facilitate sampling of the A form of DNA under the appropriate environmental conditions. Parameter optimization was guided by quantum mechanical data on model compounds, followed by calculations on several DNA duplexes in the condensed phase. The selected optimized parameters were then validated against a number of DNA duplexes, with the most extensive tests performed on the EcoRI dodecamer, including comparative calculations using the Amber Parm99bsc0 force field. The new CHARMM model better reproduces experimentally observed sampling of the BII conformation, including sampling as a function of sequence. In addition, the model reproduces the A form of the 1ZF1 duplex in 75 % ethanol, and yields a stable Z-DNA conformation of duplex (GTACGTAC) in its crystal environment. The resulting model, in combination with a recent reoptimization of the CHARMM27 force field for RNA, will be referred to as CHARMM36.

Related Papers

• → An adaptive umbrella sampling procedure in conformational analysis using molecular dynamics and its application to glycol(1992)149 cited
• → Molecular dynamics simulations of a DMSO/water mixture using the AMBER force field(2018)20 cited
• → Development and Optimization of a New Force Field for Flexible Aluminosilicates, Enabling Fast Molecular Dynamics Simulations on Parallel Architectures(2012)28 cited
• → Secondary structure dependence of amyloid‐β(1–40) on simulation techniques and force field parameters(2021)24 cited
• → Molecular dynamics simulations of an apoliprotein A–I derived peptide in explicit water(2008)7 cited