Accurate sequence-dependent coarse-grained model for conformational and elastic properties of double-stranded DNA

We introduce MADna, a sequence-dependent coarse-grained model of double-stranded DNA (dsDNA), where each nucleotide is described by three beads localized at the sugar and base moieties, and at the phosphate group. The sequence dependence is included by considering a step-dependent parameterization of the bonded interactions, which are tuned in order to reproduce the values of key observables obtained from exhaustive atomistic simulations from literature. The predictions of the model are benchmarked against an independent set of all-atom simulations, showing that it captures with high fidelity the sequence dependence of conformational and elastic features beyond the single step considered in its formulation. A remarkably good agreement with experiments is found for both sequence-averaged and sequence-dependent conformational and elastic features, including the stretching and torsion moduli, the twist-stretch and twist-bend couplings, the persistence length and the helical pitch. Overall, for the inspected quantities, the model has a precision comparable to atomistic simulations, hence providing a reliable coarse-grained description for the rationalization of singlemolecule experiments and the study of cellular processes involving dsDNA. Owing to the simplicity of its formulation, MADna can be straightforwardly included in common simulation engines.

Observing the base-by-base search for native structure along transition paths during the folding of single nucleic acid hairpins

Biomolecular folding involves searching among myriad possibilities for the native conformation, but the elementary steps expected from theory for this search have never been detected directly. We probed the dynamics of folding at high resolution using optical tweezers, measuring individual trajectories as nucleic acid hairpins passed through the high-energy transition states that dominate kinetics and define folding mechanisms. We observed brief but ubiquitous pauses in the transition states, with a dwell time distribution that matched microscopic theories of folding quantitatively. The sequence dependence suggested that pauses were dominated by microbarriers from nonnative conformations during the search by each nucleotide residue for the native base-pairing conformation. Furthermore, the pauses were position dependent, revealing subtle local variations in energy–landscape roughness and allowing the diffusion coefficient describing the microscopic dynamics within the barrier to be found without reconstructing the shape of the energy landscape. These results show how high-resolution measurements can elucidate key microscopic events during folding to test fundamental theories of folding.

RNA Helix Thermodynamics: The End Game

Nearest neighbor parameters for estimating the folding stability of RNA secondary structures are in widespread use. For helices, current parameters penalize terminal AU base pairs relative to terminal GC base pairs. We curated an expanded database of helix stabilities determined by optical melting experiments. Analysis of the updated database shows that terminal penalties depend on the sequence identity of the adjacent penultimate base pair. New nearest neighbor parameters that include this additional sequence dependence accurately predict the measured values of 271 helices in an updated database with a correlation coefficient of 0.982. This refined understanding of helix ends facilitates fitting terms for base pair stacks with GU pairs. Prior parameter sets treated 5'GGUC/3'CUGG separately from other 5'GU/3'UG stacks. The improved understanding of helix end stability, however, makes the separate treatment unnecessary. Introduction of the additional terms was tested with three optical melting experiments. The average absolute difference between measured and predicted free energy changes at 37° C for these three duplexes containing terminal adjacent AU and GU pairs improved from 1.38 to 0.27 kcal/mol. This confirms the need for the additional sequence dependence in the model.

Progressive Domain Segregation in Early Embryonic Development and Underlying Correlation to Genetic and Epigenetic Changes

Chromatin undergoes drastic structural organization and epigenetic reprogramming during embryonic development. We present here a consistent view of the chromatin structural change, epigenetic reprogramming, and the corresponding sequence-dependence in both mouse and human embryo development. The two types of domains, identified earlier as forests (CGI-rich domains) and prairies (CGI-poor domains) based on the uneven distribution of CGI in the genome, become spatially segregated during embryonic development, with the exception of zygotic genome activation (ZGA) and implantation, at which point significant domain mixing occurs. Structural segregation largely coincides with DNA methylation and gene expression changes. Genes located in mixed prairie domains show proliferation and ectoderm differentiation-related function in ZGA and implantation, respectively. The chromatin of the ectoderm shows the weakest and the endoderm the strongest domain segregation in germ layers. This chromatin structure difference between different germ layers generally enlarges upon further differentiation. The systematic chromatin structure establishment and its sequence-based segregation strongly suggest the DNA sequence as a possible driving force for the establishment of chromatin 3D structures that profoundly affect the expression profile. Other possible factors correlated with or influencing chromatin structures, including transcription, the germ layers, and the cell cycle, are discussed for an understanding of concerted chromatin structure and epigenetic changes in development.

Understanding the Origin of Structural Diversity of DNA Double Helix

Deciphering the contribution of DNA subunits to the variability of its 3D structure represents an important step toward the elucidation of DNA functions at the atomic level. In the pursuit of that goal, our previous studies revealed that the essential conformational characteristics of the most populated “canonic” BI and AI conformational families of Watson–Crick duplexes, including the sequence dependence of their 3D structure, preexist in the local energy minima of the elemental single-chain fragments, deoxydinucleoside monophosphates (dDMPs). Those computations have uncovered important sequence-dependent regularity in the superposition of neighbor bases. The present work expands our studies to new minimal fragments of DNA with Watson–Crick nucleoside pairs that differ from canonic families in the torsion angles of the sugar-phosphate backbone (SPB). To address this objective, computations have been performed on dDMPs, cdDMPs (complementary dDMPs), and minimal fragments of SPBs of respective systems by using methods of molecular and quantum mechanics. These computations reveal that the conformations of dDMPs and cdDMPs having torsion angles of SPB corresponding to the local energy minima of separate minimal units of SPB exhibit sequence-dependent characteristics representative of canonic families. In contrast, conformations of dDMP and cdDMP with SPB torsions being far from the local minima of separate SPB units exhibit more complex sequence dependence.

The cause of on-target point mutations generated by CRISPR-Cas9 treatment in the yeast Xanthophyllomyces dendrorhous

ABSTRACTRecognizing outcomes of DNA repair induced by CRISPR-Cas9 cutting is vital for precise genome editing. Reported DNA repair outcomes after Cas9 cutting include deletions/insertions and low frequency of genomic rearrangements and nucleotide substitutions. Thus far, substitution mutations caused by CRISPR-Cas9 has not attracted much attention. Here, we identified on-target point mutations induced by CRISPR-Cas9 treatment in the yeast Xanthophyllomyces dendrorhous by Sanger and Illumina sequencing. Different from previous studies, our findings suggested that the on-target mutations are not random and they cannot render the gRNA effective. Moreover, these point mutations showed strong sequence dependence that is not consistent with the observations in Hela cells, in which CRISPR-mediated substitutions were considered lacking sequence dependence and conversion preferences. Furthermore, this study demonstrated that the NHEJ components Ku70, Ku80, Mre11, or RAD50, and the overlapping roles of non-essential DNA polymerases were necessary for the emergence of point mutations, increasing the knowledge on CRISPR-Cas9 mediated DNA repair.

Sequence dependence of transient Hoogsteen base-pairing in DNA

Hoogsteen (HG) base-pairing is characterized by a 180° rotation of the purine base with respect to the Watson-Crick-Franklin (WCF) motif. Recently, it has been found that both conformations coexist in a dynamical equilibrium and that several biological functions require HG pairs. This relevance has motivated experimental and computational investigations of the base-pairing transition. However, a systematic simulation of sequence variations has remained out of reach. Here, we employ advanced path-based methods to perform unprecedented free-energy calculations. Our methodology enables us to study the different mechanisms of purine rotation, either remaining inside or after flipping outside of the double helix. We study seven different sequences, which are neighbor variations of a well-studied A·T pair in A6-DNA. We observe the known effect of A·T steps favoring HG stability, and find evidence of triple-hydrogen-bonded neighbors hindering the inside transition. More importantly, we identify a dominant factor: the direction of the A rotation, with the 6-ring pointing either towards the longer or shorter segment of the chain, respectively relating to a lower or higher barrier. This highlights the role of DNA's relative flexibility as a modulator of the WCF/HG dynamic equilibrium. Additionally, we provide a robust methodology for future HG proclivity studies.

Electrical characteristics of amyloid beta peptides in vertical junctions

Assembled amyloid beta (Aβ) peptides have been considered pathological assemblies involved in human brain diseases, and the electron transfer or electron transport characteristics of Aβ are important for the formation of structured assemblies. Here, we report the electrical characteristics of surface-assembled Aβ peptides similar to those observed in Alzheimer’s patients. These characteristics correlate to their electron transfer characteristics. Electrical current–voltage plots of Aβ vertical junction devices show the Aβ sequence dependence of the current densities at both Aβ monomers (mono-Aβs) and Aβ oligomers (oli-Aβs), while Aβ sequence dependence is not clearly observed in the electrical characteristics of Aβ planar field effect transistors (FETs). In particular, surface oligomerization of Aβ peptides drastically decreases the activity of electron transfer, which presents a change in the electron transport pathway in the Aβ vertical junctions. Electron transport at oli-Aβ junctions is symmetric (tunneling/tunneling) due to the weak and voltage-independent coupling of the less redox-reactive oli-Aβ to the contacts, while that at mono-Aβ junctions is asymmetric (hopping/tunneling) due to redox levels of mono-Aβ voltage-dependently coupled with contact electrodes. Consequently, through vertical junctions, the sequence- and conformation-dependent electrical characteristics of Aβs can reveal their electron transfer activities.