DNA Binding and Bending by Sac7d is Stepwise

C-B-A Test of DNA Force Fields

The DNA duplex may be locally strongly bent in complexes with proteins, for example, with polymerases or in a nucleosome. At such bends, the DNA helix is locally in the noncanonical forms A (with a narrow major groove and a large amount of north sugars) or C (with a narrow minor groove and a large share of BII phosphates). To model the formation of such complexes by molecular dynamics methods, the force field is required to reproduce these conformational transitions for a naked DNA. We analyzed the available experimental data on the B-C and B-A transitions under the conditions easily implemented in modeling: in an aqueous NaCl solution. We selected six DNA duplexes which conformations at different salt concentrations are known reliably enough. At low salt concentrations, poly(GC) and poly(A) are in the B-form, classical and slightly shifted to the A-form, respectively. The duplexes ATAT and GGTATACC have a strong and salt concentration dependent bias toward the A-form. The polymers poly(AC) and poly(G) take the C- and A-forms, respectively, at high salt concentrations. The reproduction of the behavior of these oligomers can serve as a test for the balance of interactions between the base stacking and the conformational flexibility of the sugar-phosphate backbone in a DNA force field. We tested the AMBER bsc1 and CHARMM36 force fields and their hybrids, and we failed to reproduce the experiment. In all the force fields, the salt concentration dependence is very weak. The known B-philicity of the AMBER force field proved to result from the B-philicity of its excessively strong base stacking. In the CHARMM force field, the B-form is a result of a fragile balance between the A-philic base stacking (especially for G:C pairs) and the C-philic backbone. Finally, we analyzed some recent simulations of the LacI-, SOX-4-, and Sac7d-DNA complex formation in the framework of the AMBER force field.

Evidence for a bind-then-bend mechanism for architectural DNA binding protein yNhp6A

The yeast Nhp6A protein (yNhp6A) is a member of the eukaryotic HMGB family of chromatin factors that enhance apparent DNA flexibility. yNhp6A binds DNA nonspecifically with nM affinity, sharply bending DNA by >60°. It is not known whether the protein binds to unbent DNA and then deforms it, or if bent DNA conformations are ‘captured’ by protein binding. The former mechanism would be supported by discovery of conditions where unbent DNA is bound by yNhp6A. Here, we employed an array of conformational probes (FRET, fluorescence anisotropy, and circular dichroism) to reveal solution conditions in which an 18-base-pair DNA oligomer indeed remains bound to yNhp6A while unbent. In 100 mM NaCl, yNhp6A-bound DNA unbends as the temperature is raised, with no significant dissociation of the complex detected up to ∼45°C. In 200 mM NaCl, DNA unbending in the intact yNhp6A complex is again detected up to ∼35°C. Microseconds-resolved laser temperature-jump perturbation of the yNhp6a–DNA complex revealed relaxation kinetics that yielded unimolecular DNA bending/unbending rates on timescales of 500 μs−1 ms. These data provide the first direct observation of bending/unbending dynamics of DNA in complex with yNhp6A, suggesting a bind-then-bend mechanism for this protein.

Biasing Simulations of DNA Base Pair Parameters with Application to Propellor Twisting in AT/AT, AA/TT, and AC/GT Steps and Their Uracil Analogs

An accurate and efficient implementation of the six DNA base pair parameters as order parameters for enhanced sampling simulations is presented. The parameter definitions are defined by vector algebra operations on a reduced atomic set of the base pair, and correlate very well with standard definitions. Application of the model is illustrated by umbrella sampling simulations of propeller twisting within AT/AT, AA/TT, and AC/GT steps and their uracil analogs. Strong correlations are found between propeller twisting and a number of conformational parameters, including buckle, opening, BI/BII backbone configuration, and sugar puckering. The thymine methyl group is observed to notably alter the local conformational free energy landscape, with effects within and directly upstream of the thymine containing base pair.

Anisotropy of B-DNA Groove Bending

DNA bending is critical for DNA packaging, recognition, and repair, and occurs toward either the major or the minor groove. The anisotropy of B-DNA groove bending was quantified for eight DNA sequences by free energy simulations employing a novel reaction coordinate. The simulations show that bending toward the major groove is preferred for non-A-tracts while the A-tract has a high tendency of bending toward the minor groove. Persistence lengths were generally larger for bending toward the minor groove, which is thought to originate from differences in groove hydration. While this difference in stiffness is one of the factors determining the overall preference of bending direction, the dominant contribution is shown to be a free energy offset between major and minor groove bending. The data suggests that, for the A-tract, this offset is largely determined by inherent structural properties, while differences in groove hydration play a large role for non-A-tracts. By quantifying the energetics of DNA groove bending and rationalizing the origins of the anisotropy, the calculations provide important new insights into a key biological process.

Molecular Dynamics Simulations of 5-Hydroxycytosine Damaged DNA

Oxidation of cytosine is a leading cause of mutations and can lead to cancer. Here we report molecular dynamics simulations that characterized the structure and flexibility of 5-hydroxycytosine damaged DNA. A total of four systems were studied: undamaged DNA, damaged DNA base paired to a matching guanine, damaged DNA base paired to a mismatching adenine, and the corresponding undamaged mismatched strand. The simulations showed high spatial similarity between undamaged and damaged DNA; however, the matched damaged strand had greater overtwisting flexibility, and for both the matched and unmatched strands sugar puckering was much more flexible at the damaged site. The mismatch introduced larger changes, notably a loss in hydrogen bonding and a gain in stacking interactions, as well as effects on base pair and step geometry and solvation. Implications for damage recognition are discussed.

Enhanced sampling simulations of DNA step parameters

A novel approach for the selection of step parameters as reaction coordinates in enhanced sampling simulations of DNA is presented. The method uses three atoms per base and does not require coordinate overlays or idealized base pairs. This allowed for a highly efficient implementation of the calculation of all step parameters and their Cartesian derivatives in molecular dynamics simulations. Good correlation between the calculated and actual twist, roll, tilt, shift, and slide parameters is obtained, while the correlation with rise is modest. The method is illustrated by its application to the methylated and unmethylated 5'-CATGTGACGTCACATG-3' double stranded DNA sequence. One-dimensional umbrella simulations indicate that the flexibility of the central CG step is only marginally affected by methylation.

Coupled binding-bending-folding: The complex conformational dynamics of protein-DNA binding studied by atomistic molecular dynamics simulations

BACKGROUND

Protein-DNA binding often involves dramatic conformational changes such as protein folding and DNA bending. While thermodynamic aspects of this behavior are understood, and its biological function is often known, the mechanism by which the conformational changes occur is generally unclear. By providing detailed structural and energetic data, molecular dynamics simulations have been helpful in elucidating and rationalizing protein-DNA binding.

SCOPE OF REVIEW

This review will summarize recent atomistic molecular dynamics simulations of the conformational dynamics of DNA and protein-DNA binding. A brief overview of recent developments in DNA force fields is given as well.

MAJOR CONCLUSIONS

Simulations have been crucial in rationalizing the intrinsic flexibility of DNA, and have been instrumental in identifying the sequence of binding events, the triggers for the conformational motion, and the mechanism of binding for a number of important DNA-binding proteins.

GENERAL SIGNIFICANCE

Molecular dynamics simulations are an important tool for understanding the complex binding behavior of DNA-binding proteins. With recent advances in force fields and rapid increases in simulation time scales, simulations will become even more important for future studies. This article is part of a Special Issue entitled Recent developments of molecular dynamics.