Register-Shifted Structures in Base-Flipped Uracil-Damaged DNA.

We report the occurrence of register-shifted structures in simulations of uracil-containing dsDNA. These occur when the 3' base vicinal to uracil is thymine in U:A base-paired DNA. Upon base flipping of uracil, this 3' thymine hydrogen bonds with the adenine across the uracil instead of its complementary base. The register-shifted structure is persistent and sterically blocks re-entry of uracil into the helix stack. Register shifting might be important for DNA repair since the longer exposure of the lesion in register-shifted structures could facilitate enzymatic recognition and repair.

Orthogonal Halogen-Bonding-Driven 3D Supramolecular Assembly of Right-Handed Synthetic Helical Peptides.

Peptide-mediated self-assembly is a prevalent method for creating highly ordered supramolecular architectures. Herein, we report the first example of orthogonal C-X⋅⋅⋅X-C/C-X⋅⋅⋅π halogen bonding and hydrogen bonding driven crystalline architectures based on synthetic helical peptides bearing hybrids of l-sulfono-γ-AApeptides and natural amino acids. The combination of halogen bonding, intra-/intermolecular hydrogen bonding, and intermolecular hydrophobic interactions enabled novel 3D supramolecular assembly. The orthogonal halogen bonding in the supramolecular architecture exerts a novel mechanism for the self-assembly of synthetic peptide foldamers and gives new insights into molecular recognition, supramolecular design, and rational design of biomimetic structures.

Hydrogen-Bonding-Driven 3D Supramolecular Assembly of Peptidomimetic Zipper.

Hydrogen-bonding-driven three-dimensional (3D) assembly of a peptidomimetic zipper has been established for the first time by using an α/AApeptide zipper that assembles into a de novo lattice arrangement through two layers of hydrogen-bonded linker-directed interactions. Via a covalently bridged 1D 413-helix, drastic enhancement in stability has been achieved in the formed 3D crystalline supramolecular architecture as evidenced by gas-sorption studies. As the first example of an unnatural peptidic zipper, the dimensional augmentation of the zipper differs from metal-coordinated strategies, and may have general implications for the preparation of peptidic functional materials for a variety of future applications.

De Novo Left-Handed Synthetic Peptidomimetic Foldamers.

The development of peptidomimetic helical foldamers with a wide repertoire of functions is of significant interest. Herein, we report the X-ray crystal structures of a series of homogeneous l-sulfono-γ-AA foldamers and elucidate their folding conformation at the atomic level. Single-crystal X-ray crystallography revealed that this class of oligomers fold into unprecedented dragon-boat-shaped and unexpected left-handed helices, which are stabilized by the 14-hydrogen-bonding pattern present in all sequences. These l-sulfono-γ-AApeptides have a helical pitch of 5.1 Šand exactly four side chains per turn, and the side chains lie perfectly on top of each other along the helical axis. 2D NMR spectroscopy, computational simulations, and CD studies support the folding conformation in solution. Our results provide a structural basis at the atomic level for the design of novel biomimetics with a precise arrangement of functional groups in three dimensions.

Right-Handed Helical Foldamers Consisting of De Novo d-AApeptides.

New types of foldamer scaffolds are formidably challenging to design and synthesize, yet highly desirable as structural mimics of peptides/proteins with a wide repertoire of functions. In particular, the development of peptidomimetic helical foldamers holds promise for new biomaterials, catalysts, and drug molecules. Unnatural l-sulfono-γ-AApeptides were recently developed and shown to have potential applications in both biomedical and material sciences. However, d-sulfono-γ-AApeptides, the enantiomers of l-sulfono-γ-AApeptides, have never been studied due to the lack of high-resolution three-dimensional structures to guide structure-based design. Herein, we report the first synthesis and X-ray crystal structures of a series of 2:1 l-amino acid/d-sulfono-γ-AApeptide hybrid foldamers, and elucidate their folded conformation at the atomic level. Single-crystal X-ray crystallography indicates that this class of oligomers folds into well-defined right-handed helices with unique helical parameters. The helical structures were consistent with data obtained from solution 2D NMR, CD studies, and molecular dynamics simulations. Our findings are expected to inspire the structure-based design of this type of unique folding biopolymers for biomaterials and biomedical applications.

Free Energy Coupling between DNA Bending and Base Flipping.

Free energy simulations are presented to probe the energetic coupling between DNA bending and the flipping of a central thymine in double stranded DNA 13mers. The energetics are shown to depend on the neighboring base pairs, and upstream C or T or downstream C tended to make flipping more costly. Flipping to the major groove side was generally preferred. Bending aids flipping, by pushing the system up in free energy, but for small and intermediate bending angles the two were uncorrelated. At higher bending angles, bending and flipping became correlated, and bending primed the system for base flipping toward the major groove. Flipping of the 6-4 pyrimidine-pyrimidone and pyrimidine dimer photoproducts is shown to be more facile than for undamaged DNA. For the damages, major groove flipping was preferred, and DNA bending was much facilitated in the 6-4 pyrimidine-pyrimidone damaged system. Aspects of the calculations were verified by structural analyses of protein-DNA complexes with flipped bases.

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.

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.

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.

Secondary Structure Adopted by the Gly-Gly-X Repetitive Regions of Dragline Spider Silk.

Solid-state NMR and molecular dynamics (MD) simulations are presented to help elucidate the molecular secondary structure of poly(Gly-Gly-X), which is one of the most common structural repetitive motifs found in orb-weaving dragline spider silk proteins. The combination of NMR and computational experiments provides insight into the molecular secondary structure of poly(Gly-Gly-X) segments and provides further support that these regions are disordered and primarily non-ß-sheet. Furthermore, the combination of NMR and MD simulations illustrate the possibility for several secondary structural elements in the poly(Gly-Gly-X) regions of dragline silks, including ß-turns, 310-helicies, and coil structures with a negligible population of α-helix observed.

Importance of ß2-ß3 Loop Motion for the Increased Binding and Decreased Selectivity of the ΔLL Mutant of the Human Papillomavirus Type 6 E2 Protein.

The binding affinity of the human papillomavirus type 6 E2 protein is strongly mediated by the sequence of the DNA linker region, with high affinity for the AATT linker and low affinity for the CCGG linker. When two terminal leucine residues are removed from the protein, the level of binding to both strands increases, but unequally, resulting in a significant decrease in selectivity for the AATT linker strand. To rationalize this behavior, we performed molecular dynamics simulations of the wild-type and mutant protein in the apo state and bound to DNA with high-affinity AATT and low-affinity CCGG linker strands. While no stable contacts were made between the ß2-ß3 loop and DNA in the wild type, this loop was repositioned in the mutant complexes and formed electrostatic contacts with the DNA backbone. More contacts were formed when the mutant was bound to the CCGG linker strand than to the AATT linker strand, resulting in a more favorable change in interaction energy for the CCGG strand. In addition, significant differences in correlated motions were found, which further explained the differences in binding. The simulations suggest that ß2-ß3 loop motions are responsible for the increased affinity and decreased selectivity of the mutant protein.

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.

RevErbα preferentially deforms DNA by induced fit.

Using free-energy simulations, we have shown that RevErbα-induced DNA deformation preferentially occurs by induced fit rather than by conformational selection, even though the DNA is only slightly distorted in the complex. Our study shows that information on the sequence of binding events is needed to establish whether conformational selection or induced fit is operative, and that the presence of multiple apo-state structures might not be enough to distinguish between these binding models

Free energy simulation of helical transitions.

An umbrella sampling method for the calculation of free energies for helical transitions is presented. The method biases structures toward helices of a desired radius and pitch. Although computationally complex, the method has negligible overhead in actual applications. To illustrate the method, calculations of the helical free energy landscape of several peptides are presented for both the CHARMM and the AMBER force fields.

DNA binding and bending by Sac7d is stepwise.

Which came first? Using thermodynamic cycles constructed from experimental measures of the binding free energy and free energy simulations, we have shown that binding of Sulfolobus acidocaldarius protein Sac7d to DNA occurs before DNA bending, thus indicating that a conformational selection mechanism is unlikely to be operative in Sac7d binding.

Systematic assessment of the flexibility of uracil damaged DNA.

Uracil is a common DNA lesion which is recognized and removed by uracil DNA-glycosylase (UDG) as a part of the base excision repair pathway. Excision proceeds by base flipping, and UDG efficiency is thought to depend on the ease of deformability of the bases neighboring the lesion. We used molecular dynamics simulations to assess the flexibility of a large library of dsDNA strands, containing all tetranucleotide motifs with U:A, U:G, T:A or C:G base pairs. Our study demonstrates that uracil damaged DNA largely follows trends in flexibility of undamaged DNA. Measured bending persistence lengths, groove widths, step parameters and base flipping propensities demonstrate that uracil increases the flexibility of DNA, and that U:G base paired strands are more flexible than U:A strands. Certain sequence contexts are more deformable than others, with a key role for the 3' base next to uracil. Flexibilities are large when this base is an A or G, and repressed for a C or T. A 5' T adjacent to the uracil strongly promotes flexibility, but other 5' bases are less influential. DNA bending is correlated to step deformations and base flipping, and bending aids flipping. Our study implies that the link between substrate flexibility and UDG efficiency is widely valid, helps explain why UDG prefers to bind U:G base paired strands, and suggests that the DNA bending angle of the UDG-substrate complex is optimal for base flipping.Communicated by Ramaswamy H. Sarma.

Uracil-DNA glycosylase efficiency is modulated by substrate rigidity.

Uracil DNA-glycosylase (UNG) is a DNA repair enzyme that removes the highly mutagenic uracil lesion from DNA using a base flipping mechanism. Although this enzyme has evolved to remove uracil from diverse sequence contexts, UNG excision efficiency depends on DNA sequence. To provide the molecular basis for rationalizing UNG substrate preferences, we used time-resolved fluorescence spectroscopy, NMR imino proton exchange measurements, and molecular dynamics simulations to measure UNG specificity constants (kcat/KM) and DNA flexibilities for DNA substrates containing central AUT, TUA, AUA, and TUT motifs. Our study shows that UNG efficiency is dictated by the intrinsic deformability around the lesion, establishes a direct relationship between substrate flexibility modes and UNG efficiency, and shows that bases immediately adjacent to the uracil are allosterically coupled and have the greatest impact on substrate flexibility and UNG activity. The finding that substrate flexibility controls UNG efficiency is likely significant for other repair enzymes and has major implications for the understanding of mutation hotspot genesis, molecular evolution, and base editing.

The natural DNA bending angle in the lac repressor headpiece-O1 operator complex is determined by protein-DNA contacts and water release.

We performed molecular dynamics simulations of the lac repressor headpiece-O1 operator complex for natural, over and underbent DNA to assess the factors that determine the natural DNA bending angle. At the natural angle, the specific and nonspecific contacts between the protein and DNA are optimized. Protein-DNA contacts show different angle dependences in the right and left sites, with the left site generally getting weaker and the right site getting stronger as the bending angle increases. Two entropic factors were identified as well: at the natural bending angle, water release and the quasiharmonic protein configurational entropy are maximized. The gain in protein configurational entropy might stem from an entropy-entropy compensation mechanism, in which a reduction in protein fluctuations is offset by a loss in correlations between the right and left sites.

Secondary structure of peptides mimicking the Gly-rich regions of major ampullate spidroin protein 1 and 2.

Spider dragline silk has highly desirable material properties, possessing high extensibility, strength, and biocompatibility. Before it is spun, the constituent proteins are stored in a concentrated dope that is void of fibrils. To investigate the structural properties of the amorphous fiber regions in the dope, computer simulations were performed on model peptides representing the N. clavipes Gly-rich regions. Analysis of the secondary structure found predominantly turns, bends and coils; a small 31-helical population decreased with increasing concentration. Interestingly, the population of 31-helices saw a large increase in octanol. These results indicate that the unusual 31-helical secondary structure of the Gly-rich region of the fiber is a consequence of the spinning process, and that the low dielectric environment of the fiber may assist in favoring this structure.

Cy3-DNA stacking interactions strongly depend on the identity of the terminal basepair.

We characterized the effect of the first basepair on the conformational dynamics of the fluorescent dye Cy3 attached to the 5' end of double-stranded DNA using gaussian-mixture adaptive umbrella sampling simulations. In the simulations, the sampling of all five dihedral angles along the linker was enhanced, so that both stacked and unstacked states were sampled. The affinity of Cy3 for a T·A basepair (with the dye attached to T) was found to be significantly less than for the other basepairs. This was verified experimentally by measuring the activation energies for cis-trans isomerization of the dye. The simulation and experimental results indicate the existence of partially unstacked conformations amenable to photoisomerization. The simulations also showed that stacking of Cy3 straightens the DNA while stabilizing the first basepair. Our findings indicate that fluorescence is modulated by Cy3-DNA interactions in a sequence-dependent manner.