Base-pair Steps Research Articles

Exploration of Configurational and Topological Properties of Minichromosomes using Elastic Energy Optimizations at the DNA Base-Pair Level

Advances in experimental biophysics continue to broaden our understanding of protein-DNA interactions, highlighting structural characteristics that may account for modes of genetic expression and regulation. These interactions influence the configuration of DNA that enters and exits the protein and may lead to DNA supercoiling. Previous studies of minichromosomes showcased discrepancies in linking number and topology across short circular chains containing nucleosomes with ∼145 base-pairs of DNA wrapped less than two superhelical turns around a core of eight histone proteins. The increasing collection of high-resolution nucleosome structures has been valuable in developing computational models to survey differences in minichromosomes with different histone binding surfaces or DNA sequences. The models offer an approximation of the energetically favorable configurations of protein-free DNA, which depend on the spatial arrangements of the entry and exit base-pairs and/or the nucleotide sequence. Elastic energy minimization calculations performed at the level of base-pair steps starting from approximately closed circular DNA models generated by Monte Carlo sampling have revealed the sensitivity of minichromosomal configurations to free-DNA chain length, entry-exit spatial arrangements, and imposed base-pair twist. Producing optimized circular structures on a mesoscale level may shed light on certain topological behaviors, such as strand crossings at the entry and exit regions of linker DNA, and may reconcile the linking number paradox that surround minichromosomes.

Computational and NMR studies of RNA duplexes with an internal pseudouridine-adenosine base pair

Pseudouridine (Ψ) is the most common chemical modification present in RNA. In general, Ψ increases the thermodynamic stability of RNA. However, the degree of stabilization depends on the sequence and structural context. To explain experimentally observed sequence dependence of the effect of Ψ on the thermodynamic stability of RNA duplexes, we investigated the structure, dynamics and hydration of RNA duplexes with an internal Ψ-A base pair in different nearest-neighbor sequence contexts. The structures of two RNA duplexes containing 5′-GΨC/3′-CAG and 5′-CΨG/3′-GAC motifs were determined using NMR spectroscopy. To gain insight into the effect of Ψ on duplex dynamics and hydration, we performed molecular dynamics (MD) simulations of RNA duplexes with 5′-GΨC/3′-CAG, 5′-CΨG/3′-GAC, 5′-AΨU/3′-UAA and 5′-UΨA/3′-AAU motifs and their unmodified counterparts. Our results showed a subtle impact from Ψ modification on the structure and dynamics of the RNA duplexes studied. The MD simulations confirmed the change in hydration pattern when U is replaced with Ψ. Quantum chemical calculations showed that the replacement of U with Ψ affected the intrinsic stacking energies at the base pair steps depending on the sequence context. The calculated intrinsic stacking energies help to explain the experimentally observed sequence dependent changes in the duplex stability from Ψ modification.

Inconsistencies in ab initio evaluations of non-additive contributions of DNA stacking energies

We evaluated the non-additive contributions of the inter-molecular interactions in B-DNA stacking by using diffusion Monte Carlo methods with fixed node approximations (FNDMC). For some base-pair steps, we found that their non-additive contributions evaluated by FNDMC significantly differ from those by any other ab initio methods, while there are no remarkable findings on their stacking energies themselves. The apparently unexpected results of non-additivity raise issues in both FNDMC and correlated wavefunction methods. For the latter, it can be partly attributed to the imperfect complete basis set (CBS) correction scheme due to the limitation of the computational costs. On the other hand, the striking contrast between the stacking and non-additivity behaviors was found in FNDMC. This might imply that the error cancellations of the fixed node biases in FNDMC work well for the stacking energies, while not for the non-additivity contributions involving charge transfers caused by hydrogen bonds bridging Watson-Crick base pairs.

Sequence-dependent RNA helix conformational preferences predictably impact tertiary structure formation

Structured RNAs and RNA complexes underlie biological processes ranging from control of gene expression to protein translation. Approximately 50% of nucleotides within known structured RNAs are folded into Watson-Crick (WC) base pairs, and sequence changes that preserve these pairs are typically assumed to preserve higher-order RNA structure and binding of macromolecule partners. Here, we report that indirect effects of the helix sequence on RNA tertiary stability are, in fact, significant but are nevertheless predictable from a simple computational model called RNAMake-∆∆G. When tested through the RNA on a massively parallel array (RNA-MaP) experimental platform, blind predictions for >1500 variants of the tectoRNA heterodimer model system achieve high accuracy (rmsd 0.34 and 0.77 kcal/mol for sequence and length changes, respectively). Detailed comparison of predictions to experiments support a microscopic picture of how helix sequence changes subtly modulate conformational fluctuations at each base-pair step, which accumulate to impact RNA tertiary structure stability. Our study reveals a previously overlooked phenomenon in RNA structure formation and provides a framework of computation and experiment for understanding helix conformational preferences and their impact across biological RNA and RNA-protein assemblies.

Cooperative DNA binding by proteins through DNA shape complementarity.

Localized arrays of proteins cooperatively assemble onto chromosomes to control DNA activity in many contexts. Binding cooperativity is often mediated by specific protein–protein interactions, but cooperativity through DNA structure is becoming increasingly recognized as an additional mechanism. During the site-specific DNA recombination reaction that excises phage λ from the chromosome, the bacterial DNA architectural protein Fis recruits multiple λ-encoded Xis proteins to the attR recombination site. Here, we report X-ray crystal structures of DNA complexes containing Fis + Xis, which show little, if any, contacts between the two proteins. Comparisons with structures of DNA complexes containing only Fis or Xis, together with mutant protein and DNA binding studies, support a mechanism for cooperative protein binding solely by DNA allostery. Fis binding both molds the minor groove to potentiate insertion of the Xis β-hairpin wing motif and bends the DNA to facilitate Xis-DNA contacts within the major groove. The Fis-structured minor groove shape that is optimized for Xis binding requires a precisely positioned pyrimidine-purine base-pair step, whose location has been shown to modulate minor groove widths in Fis-bound complexes to different DNA targets.

How B-DNA Dynamics Decipher Sequence-Selective Protein Recognition

The rules governing sequence-specific DNA–protein recognition are under a long-standing debate regarding the prevalence of base versus shape readout mechanisms to explain sequence specificity and of the conformational selection versus induced fit binding paradigms to explain binding-related conformational changes in DNA. Using a combination of atomistic simulations on a subset of representative sequences and mesoscopic simulations at the protein–DNA interactome level, we demonstrate the prevalence of the shape readout model in determining sequence-specificity and of the conformational selection paradigm in defining the general mechanism for binding-related conformational changes in DNA. Our results suggest that the DNA uses a double mechanism to adapt its structure to the protein: it moves along the easiest deformation modes to approach the bioactive conformation, while final adjustments require localized rearrangements at the base-pair step and backbone level. Our study highlights the large impact of B-DNA dynamics in modulating DNA–protein binding.

Insight into the recognition mechanism of EcoRV

EcoRV, a restriction enzyme in Escherichia coli, destroys invading foreign DNA by cleaving it at the center step of a GATATC sequence. In the EcoRV–cognate DNA crystallographic complex, a sharp kink of 50° has been found at the center base‐pair step (TA). In a previous research, we have examined the interplay between the intrinsic propensity of the cognate sequence to kink and the induction by the enzyme, and have understood the mechanism of action of the enzyme and the intrinsically behavior of the cognate sequence to undergo a kink allowing the proper placement of the atoms for cleavage. We have also demonstrated that the enzyme requires 2 magnesium ions at the catalytic site to perform the cleavage.In the present study, we perform multiple independent molecular dynamics simulations to investigate the behavior of a mutant EcoRV (E45A), in the presence of one or two magnesium ions and complexed to three DNA sequences: the cognate sequence, GATATC (TA); the non‐cognate sequence, GAATTC (AT); and with the cognate sequence methylated on the first adenine GACH3TATC (TA‐CH3). The results give further insight into the recognition and cleavage mechanism of EcoRV.This abstract is from the Experimental Biology 2019 Meeting. There is no full text article associated with this abstract published in The FASEB Journal.

Effect of Tetrameric Base-Pair Context on the Sequence-Dependent Configurations of DNA Minicircles

The biophysical characteristics of various DNA configurations, ranging from changes in localized base-pair twisting to global writhe, are known to be influenced by nucleotide sequence and composition. Decades worth of experimental findings have shown how regions of high adenine content exhibit high degrees of bending or how pyrimidine-purine dimers exhibit a high twist/low twist polymorphism, examples critical in protein interactions and genetic expression. Calculations and subsequent modeling in the past have been limited to dimeric base pair steps but the recent influx of high-resolution structural data has shown these steps to be influenced by both their 3′ and 5′ nearest neighbors, such as conformational differences of the internal AA dimers in the context of A-AA-A and A-AA-C tetramers. The growing database of high-resolution nucleic acid-containing structures can be harnessed to generate intrinsic rigid-body parameters and elastic constants of successive base pairs in all tetrameric contexts. Optimizing the elastic energy of DNA minicircles at the dinucleotide level has been used to determine the effects of sequence on global shape and elasticity. These optimizations make use of an array of intrinsic dimer models ranging from a simple linear model, a curved three-state model, or one of higher complexity with rigid-body parameters and elastic constants unique to each dimer, as well as impose spatial constraints to mimic behavior observed in protein-DNA interactions. These models have been further refined to include nearest neighbor interactions. Elastic energy minimization calculations were used on multiple circular configurations with sequences of high A⋅T content. Comparative analysis between optimization results are useful in exploring how purine and pyrimidine nucleotides before and after a dinucleotide sequence may affect polymer shape and elasticity.

Investigating Structure and Topology of Protein-Mediated DNA Loops via Computational Modeling of Elastic Energy

Many biological processes involve deformations of DNA 3D structure. DNA looping generated by a protein bound to two distant sites on the double helix is a familiar mechanism employed in regulation of cellular processes. We employ a computational method to obtain minimum energy structures of DNA loops mediated by the Lac repressor protein, whereby the potential energy of elastic deformation is optimized. By evaluating the energy required to deform DNA, we can predict the likely spatial pathway of the DNA strand. By adjusting the optimization parameters, we can evaluate the effects of torsional stress on loop configuration and global topology. For instance, altering the position or orientation of the end constraints provides insight into mechanisms that proteins may use to control DNA folding. Alternatively, changing the chain length of the constrained DNA molecule alters the total twist, which is reflected in the supercoiling of the molecule. Adjusting the intrinsic helical repeat in the optimization simulates changes in the cellular environment such as found with increased salt content. Using knowledge of the cost of deforming base-pair steps of different types we simulate the effect of nucleotide sequence on the overall structure of constrained DNA chains.

Hexameric helicase G40P unwinds DNA in single base pair steps.

Most replicative helicases are hexameric, ring-shaped motor proteins that translocate on and unwind DNA. Despite extensive biochemical and structural investigations, how their translocation activity is utilized chemo-mechanically in DNA unwinding is poorly understood. We examined DNA unwinding by G40P, a DnaB-family helicase, using a single-molecule fluorescence assay with a single base pair resolution. The high-resolution assay revealed that G40P by itself is a very weak helicase that stalls at barriers as small as a single GC base pair and unwinds DNA with the step size of a single base pair. Binding of a single ATPγS could stall unwinding, demonstrating highly coordinated ATP hydrolysis between six identical subunits. We observed frequent slippage of the helicase, which is fully suppressed by the primase DnaG. We anticipate that these findings allow a better understanding on the fine balance of thermal fluctuation activation and energy derived from hydrolysis.

Interchromophoric Interactions Determine the Maximum Brightness Density in DNA Origami Structures.

An ideal point light source is as small and as bright as possible. For fluorescent point light sources, homogeneity of the light sources is important as well as that the fluorescent units inside the light source maintain their photophysical properties, which is compromised by dye aggregation. Here we propose DNA origami as a rigid scaffold to arrange dye molecules in a dense pixel array with high control of stoichiometry and dye-dye interactions. In order to find the highest labeling density in a DNA origami structure without influencing dye photophysics, we alter the distance of two ATTO647N dyes in single base pair steps and probe the dye-dye interactions on the single-molecule level. For small distances strong quenching in terms of intensity and fluorescence lifetime is observed. With increasing distance, we observe reduced quenching and molecular dynamics. However, energy transfer processes in the weak coupling regime still have a significant impact and can lead to quenching by singlet-dark-state-annihilation. Our study fills a gap of studying the interactions of dyes relevant for superresolution microscopy with dense labeling and for single-molecule biophysics. Incorporating these findings in a 3D DNA origami object will pave the way to bright and homogeneous DNA origami nanobeads.

An In-Depth Look at DNA Crystals through the Prism of Molecular Dynamics Simulations

Summary X-ray crystallography is the primary tool for biomolecular structural determination. However, contacts formed through the crystal lattice are known to affect structures, especially for small and flexible molecules such as DNA oligomers, by introducing significant structural changes in comparison to solution. Furthermore, why molecules crystallize in certain symmetry groups, which role crystallization additives play, and whether they are just innocuous and unspecific crystallization catalysts remain unclear. By using one of the currently best-performing DNA force fields and applying significant computational effort, we described the nature of intermolecular forces that stabilize B-DNA crystals in various symmetry groups and solvent environments with an unprecedented level of detail. We showed a tight coupling between the lattice stability and the type of crystallization additives and that certain symmetry groups are stable only in the presence of a specific additive. Additives and crystal contacts induce small but non-negligible changes in the physical properties of DNA.

Investigations of Stacked DNA Base-Pair Steps: Highly Accurate Stacking Interaction Energies, Energy Decomposition, and Many-Body Stacking Effects.

The stacking energies of 10 unique B-DNA base-pair steps were calculated with highly accurate quantum chemistry and used as reference values in a thorough benchmark of (dispersion-corrected) DFT, wave function methods, tight-binding methods, and different force fields, including charge variants thereof. The reference values were computed using a focal-point energy function based on extrapolated explicitly correlated MP2-F12 and conventional CCSD(T) data at the triple-ζ level. A collection of 29 different density functionals, sometimes with multiple dispersion corrections (D3(BJ), D3M(BJ), and VV10) were evaluated, including recent functionals like B97M-V, ωB97M-V, and SCAN-D3(BJ), which perform excellently. The double-hybrid DSD-BLYP-NL was found to be the best DFT method. Common wave function methods (MP2, SCS-MP2, and MP2.5) and the SNS-MP2 protocol were tested as well, where only the latter and DLPNO-CCSD(T)/CBS were competitive with DFT. The tight-binding methods DFTB3-D3 and GFN-xTB revealed a comparatively low accuracy. The AMBER force field outperformed CHARMM and GROMOS but still showed systematic gas-phase overbinding, which could be traced back to the electrostatic term, as revealed by comparison of different sets of point charges. High-order SAPT, e.g., SAPT2 + 3δ(MP2), was not only benchmarked but also used to study the nature of the stacking interactions to high accuracy. The δ(MP2) term turned out to be crucially important to reach high accuracy. Finally, we investigated four-body stacking effects with DLPNO-CCSD(T) and DFT, which were found to be significant and strongest for the CpC base-pair step where they reached almost 30% of the total stacking energy.

Investigating the sequence-dependent mechanical properties of DNA nicks for applications in twisted DNA nanostructure design.

DNA nick can be used as a design motif in programming the shape and reconfigurable deformation of synthetic DNA nanostructures, but its mechanical properties have rarely been systematically characterized at the level of base sequences. Here, we investigated sequence-dependent mechanical properties of DNA nicks through molecular dynamics simulation for a comprehensive set of distinct DNA oligomers constructed using all possible base-pair steps with and without a nick. We found that torsional rigidity was reduced by 28–82% at the nick depending on its sequence and location although bending and stretching rigidities remained similar to those of regular base-pair steps. No significant effect of a nick on mechanically coupled deformation such as the twist-stretch coupling was observed. These results suggest that the primary structural role of nick is the relaxation of torsional constraint by backbones known to be responsible for relatively high torsional rigidity of DNA. Moreover, we experimentally demonstrated the usefulness of quantified nick properties in self-assembling DNA nanostructure design by constructing twisted DNA origami structures to show that sequence design of nicks successfully controls the twist angle of structures. Our study illustrates the importance as well as the opportunities of considering sequence-dependent properties in structural DNA nanotechnology.

Role of Inosine⁻Uracil Base Pairs in the Canonical RNA Duplexes.

Adenosine to inosine (A–I) editing is the most common modification of double-stranded RNA (dsRNA). This change is mediated by adenosine deaminases acting on RNA (ADARs) enzymes with a preference of U>A>C>G for 5′ neighbor and G>C=A>U or G>C>U=A for 3′ neighbor. A–I editing occurs most frequently in the non-coding regions containing repetitive elements such as ALUs. It leads to disruption of RNA duplex structure, which prevents induction of innate immune response. We employed standard and biased molecular dynamics (MD) simulations to analyze the behavior of RNA duplexes with single and tandem inosine–uracil (I–U) base pairs in different sequence context. Our analysis showed that the I–U pairs induce changes in base pair and base pair step parameters and have different dynamics when compared with standard canonical base pairs. In particular, the first I–U pair from tandem I–U/I–U systems exhibited increased dynamics depending on its neighboring 5′ base. We discovered that UII sequence, which is frequently edited, has lower flexibility compared with other sequences (AII, GII, CII), hence it only modestly disrupts dsRNA. This might indicate that the UAA motifs in ALUs do not have to be sufficiently effective in preventing immune signaling.

Predictive Model for the Sequence-Dependent Fluorogenic Response of Forced-Intercalation Peptide Nucleic Acid.

The forced-intercalation peptide nucleic acid (FIT-PNA) concept, introduced by Seitz and co-workers, is based on replacing a nucleobase of the PNA sequence with a cyanine dye (such as thiazole orange). The cyanine dye is thus a surrogate base that is forced to intercalate in the duplex (e.g., PNA:DNA). This allows single-mismatch sensitivity as the introduction of a mismatch in the vicinity of the dye increases freedom of motion and leads to a significant depletion of its fluorescence because of the free rotation of the monomethine bond separating the two π-systems of the cyanine dye. Herein, we designed and synthesized six FIT-PNA probes, featuring bisquinoline (BisQ), a red-emitting cyanine dye recently developed in our laboratory for FIT-PNAs. By following PNA–DNA duplex fluorescence, we found new sequence-based factors governing the fluorescence response to the mismatched FIT-PNA:DNA duplex. Fluorogenic properties are correlated with the π-stacking energy of three distinctive base pair steps (BPSs) in the PNA:DNA duplex. The first two are the two BPSs opposite BisQ, whereas the third is the BPS of the mismatch position, which presumably becomes unstacked due to the mismatch. We suggest a predictive model for FIT-PNA single-mismatch detection mechanism, a model that can be used in future research to improve FIT-PNA design.

DEPENDENCE OF HOMO–LUMO GAP OF DNA BASE PAIR STEPS ON TWIST ANGLE: A DENSITY FUNCTIONAL APPROACH

Electronic structure calculations of all 10 unique base pair (bp) steps have been calculated to study the interaction energies of the bp steps, density of states (DOS), projected density of states (pDOS) using the density functional theory (DFT). Plane wave basis with ultrasoft pseudo-potential method has been used within the local density approximation (LDA) for the exchange correlation functional. Electron densities of the bp steps corresponding to HOMO and LUMO level have been calculated to understand the difference in stacking energies of the bp steps. The variation of HOMO–LUMO gap (g) of all possible bp steps on twist angle has been studied. We have observed that out of the 10 bp steps, 4 purine–purine bp steps (d(AA), d(GG), d(AG) and d(GA)), show significant variation of [Formula: see text] on twist angle. The observed variation on twist angle of d(AA) bp step has been explained by the calculated DOS and electron densities.

Kinetic and mechanistic insight into the formation of amphetamine using the Leuckart–Wallach reaction and interaction of the drug with GpC·CpG base-pair step of DNA: a DFT study

Zahedi, E (reprint author), Islamic Azad Univ, Chem Dept, Shahrood Branch, Shahrood, Iran. e_zahedi@iau-shahrood.ac.ir

Sequence-dependent response of DNA to torsional stress: a potential biological regulation mechanism.

Torsional restraints on DNA change in time and space during the life of the cell and are an integral part of processes such as gene expression, DNA repair and packaging. The mechanical behavior of DNA under torsional stress has been studied on a mesoscopic scale, but little is known concerning its response at the level of individual base pairs and the effects of base pair composition. To answer this question, we have developed a geometrical restraint that can accurately control the total twist of a DNA segment during all-atom molecular dynamics simulations. By applying this restraint to four different DNA oligomers, we are able to show that DNA responds to both under- and overtwisting in a very heterogeneous manner. Certain base pair steps, in specific sequence environments, are able to absorb most of the torsional stress, leaving other steps close to their relaxed conformation. This heterogeneity also affects the local torsional modulus of DNA. These findings suggest that modifying torsional stress on DNA could act as a modulator for protein binding via the heterogeneous changes in local DNA structure.

Sequence-dependent nucleosome sliding in rotation-coupled and uncoupled modes revealed by molecular simulations.

While nucleosome positioning on eukaryotic genome play important roles for genetic regulation, molecular mechanisms of nucleosome positioning and sliding along DNA are not well understood. Here we investigated thermally-activated spontaneous nucleosome sliding mechanisms developing and applying a coarse-grained molecular simulation method that incorporates both long-range electrostatic and short-range hydrogen-bond interactions between histone octamer and DNA. The simulations revealed two distinct sliding modes depending on the nucleosomal DNA sequence. A uniform DNA sequence showed frequent sliding with one base pair step in a rotation-coupled manner, akin to screw-like motions. On the contrary, a strong positioning sequence, the so-called 601 sequence, exhibits rare, abrupt transitions of five and ten base pair steps without rotation. Moreover, we evaluated the importance of hydrogen bond interactions on the sliding mode, finding that strong and weak bonds favor respectively the rotation-coupled and -uncoupled sliding movements.