Local states of chromatin compaction at transcription start sites control transcription levels.

The 'open' and 'compact' regions of chromatin are considered to be regions of active and silent transcription, respectively. However, individual genes produce transcripts at different levels, suggesting that transcription output does not depend on the simple open-compact conversion of chromatin, but on structural variations in chromatin itself, which so far have remained elusive. In this study, weakly crosslinked chromatin was subjected to sedimentation velocity centrifugation, which fractionated the chromatin according to its degree of compaction. Open chromatin remained in upper fractions, while compact chromatin sedimented to lower fractions depending on the level of nucleosome assembly. Although nucleosomes were evenly detected in all fractions, histone H1 was more highly enriched in the lower fractions. H1 was found to self-associate and crosslinked to histone H3, suggesting that H1 bound to H3 interacts with another H1 in an adjacent nucleosome to form compact chromatin. Genome-wide analyses revealed that nearly the entire genome consists of compact chromatin without differences in compaction between repeat and non-repeat sequences; however, active transcription start sites (TSSs) were rarely found in compact chromatin. Considering the inverse correlation between chromatin compaction and RNA polymerase binding at TSSs, it appears that local states of chromatin compaction determine transcription levels.

Free energy landscape of RNA binding dynamics in start codon recognition by eukaryotic ribosomal pre-initiation complex.

Specific interaction between the start codon, 5'-AUG-3', and the anticodon, 5'-CAU-3', ensures accurate initiation of translation. Recent studies show that several near-cognate start codons (e.g. GUG and CUG) can play a role in initiating translation in eukaryotes. However, the mechanism allowing initiation through mismatched base-pairs at the ribosomal decoding site is still unclear at an atomic level. In this work, we propose an extended simulation-based method to evaluate free energy profiles, through computing the distance between each base-pair of the triplet interactions involved in recognition of start codons in eukaryotic translation pre-initiation complex. Our method provides not only the free energy penalty for mismatched start codons relative to the AUG start codon, but also the preferred pathways of transitions between bound and unbound states, which has not been described by previous studies. To verify the method, the binding dynamics of cognate (AUG) and near-cognate start codons (CUG and GUG) were simulated. Evaluated free energy profiles agree with experimentally observed changes in initiation frequencies from respective codons. This work proposes for the first time how a G:U mismatch at the first position of codon (GUG)-anticodon base-pairs destabilizes the accommodation in the initiating eukaryotic ribosome and how initiation at a CUG codon is nearly as strong as, or sometimes stronger than, that at a GUG codon. Our method is expected to be applied to study the affinity changes for various mismatched base-pairs.

Protocol for calculating binding free energy of RNA:RNA interactions through molecular dynamics simulations using adaptive biasing force technique.

The adaptive biasing force (ABF) technique allows sampling to proceed in a flat free energy surface when performing molecular dynamics (MD) simulations. Here, we present a protocol to perform MD simulations using the ABF technique and apply it to calculate the binding free energy of an RNA:RNA interaction. We describe steps for server setup, test running software, and building molecular models. We then detail procedures for running and configuring ABF-MD simulations and analyzing binding free energy and structural change. For complete details on the use and execution of this protocol, please refer to Fujita et al.1 and Kameda et al.2.

Molecular dynamics analysis of biomolecular systems including nucleic acids.

With the recent progress in structural biology and genome biology, structural dynamics of molecular systems that include nucleic acids has attracted attention in the context of gene regulation. The structure-function relationship is an important topic that highlights the importance of the physicochemical properties of nucleotides, as well as that of amino acids in proteins. Simulations are a useful tool for the detailed analysis of molecular dynamics that complement experiments in molecular biology; however, molecular simulation of nucleic acids is less well developed than that of proteins partly due to the physical nature of nucleic acids. In this review, we briefly describe the current status and future directions of the field as a guide to promote collaboration between experimentalists and computational biologists.

Translational recoding by chemical modification of non-AUG start codon ribonucleotide bases.

In contrast to prokaryotes wherein GUG and UUG are permissive start codons, initiation frequencies from non-AUG codons are generally low in eukaryotes, with CUG being considered as strongest. Here, we report that combined 5-cytosine methylation (5mC) and pseudouridylation (Ψ) of near-cognate non-AUG start codons convert GUG and UUG initiation strongly favored over CUG initiation in eukaryotic translation under a certain context. This prokaryotic-like preference is attributed to enhanced NUG initiation by Ψ in the second base and reduced CUG initiation by 5mC in the first base. Molecular dynamics simulation analysis of tRNAiMet anticodon base pairing to the modified codons demonstrates that Ψ universally raises the affinity of codon:anticodon pairing within the ribosomal preinitiation complex through partially mitigating discrimination against non-AUG codons imposed by eukaryotic initiation factor 1. We propose that translational control by chemical modifications of start codon bases can offer a new layer of proteome diversity regulation and therapeutic mRNA technology.

Structural dynamics of DNA depending on methylation pattern.

DNA methylation is associated with a number of biological phenomena, and plays crucial roles in epigenetic regulation of eukaryotic gene expression. It is also suggested that DNA methylation alters the mechanical properties of DNA molecules, which is likely to affect epigenetic regulation. However, it has not been systematically investigated how methylation changes the structural and dynamic features of DNA. In this research, to elucidate the effects of methylation on DNA mechanics, a fully atomic molecular dynamics simulation of double-stranded DNA with several methylation patterns was performed. Through the analysis of the relative positioning of the nucleotides (base-step variables), characteristic changes in terms of local flexibility were observed, which further affected the overall DNA geometry and stiffness. These findings may serve as a basis for a discussion on methylation-dependent DNA dynamics in physiological conditions.

Histone Tail Dynamics in Partially Disassembled Nucleosomes During Chromatin Remodeling.

Nucleosomes are structural units of the chromosome consisting of DNA wrapped around histone proteins, and play important roles in compaction and regulation of the chromatin structure. While the structure and dynamics of canonical nucleosomes have been studied extensively, those of nucleosomes in intermediate states, that occur when their structure or positioning is modulated, have been less understood. In particular, the dynamic features of partially disassembled nucleosomes have not been discussed in previous studies. Using all-atom molecular dynamics simulations, in this study, we investigated the dynamics and stability of nucleosome structures lacking a histone-dimer. DNA in nucleosomes lacking a histone H2A/H2B dimer was drastically deformed due to loss of local interactions between DNA and histones. In contrast, conformation of DNA in nucleosomes lacking H3/H4 was similar to the canonical nucleosome, as the H2A C-terminal domain infiltrated the space originally occupied by the dissociated H3/H4 histones and restricted DNA dynamics in close proximity. Our results suggest that, besides histone chaperones, the intrinsic dynamics of nucleosomes support the exchange of H2A/H2B, which is significantly more frequent than that of H3/H4.

The 1-Particle-per-k-Nucleotides (1PkN) Elastic Network Model of DNA Dynamics with Sequence-Dependent Geometry.

Coarse-grained models of DNA have made important contributions to the determination of the physical properties of genomic DNA, working as a molecular machine for gene regulation. In this study, to analyze the global dynamics of long DNA sequences with consideration of sequence-dependent geometry, we propose elastic network models of DNA where each particle represents k nucleotides (1-particle-per-k-nucleotides, 1PkN). The models were adjusted according to profiles of the anisotropic fluctuations obtained from our previous 1-particle-per-1-nucleotide (1P1N) model, which was proven to reproduce such profiles of all-atom models. We confirmed that the 1P3N and 1P4N models are suitable for the analysis of detailed dynamics such as local twisting motion. The models are intended for the analysis of large structures, e.g., 10-nm fibers in the nucleus, and nucleoids of mitochondrial or phage DNA at low computational costs. As an example, we surveyed the physical characteristics of the whole mitochondrial human and Plasmodium falciparum genomes.