Transcriptome-scale analysis uncovers conserved residues in the hydrophobic core of the bacterial RNA chaperone Hfq required for small regulatory RNA stability.

The RNA chaperone Hfq plays crucial roles in bacterial gene expression and is a major facilitator of small regulatory RNA (sRNA) action. The toroidal architecture of the Hfq hexamer presents three well-characterized surfaces that allow it to bind sRNAs to stabilize them and engage target transcripts. Hfq-interacting sRNAs are categorized into two classes based on the surfaces they use to bind Hfq.

Can We Ever Develop an Ideal RNA Force Field? Lessons Learned from Simulations of the UUCG RNA Tetraloop and Other Systems.

Molecular dynamics (MD) simulations are an important and well-established tool for investigating RNA structural dynamics, but their accuracy relies heavily on the quality of the employed force field (). In this work, we present a comprehensive evaluation of widely used pair-additive and polarizable RNA s using the challenging UUCG tetraloop (TL) benchmark system. Extensive standard MD simulations, initiated from the NMR structure of the 14-mer UUCG TL, revealed that most s did not maintain the native state, instead favoring alternative loop conformations.

How Binding Site Flexibility Promotes RNA Scanning by TbRGG2 RRM: A Molecular Dynamics Simulation Study.

RNA recognition motifs (RRMs) are a key class of proteins that primarily bind single-stranded RNAs. In this study, we applied standard atomistic molecular dynamics simulations to obtain insights into the intricate binding dynamics between uridine-rich RNAs and TbRGG2 RRM using the recently developed OL3-Stafix AMBER force field, which improves the description of single-stranded RNA molecules. Complementing structural experiments that unveil a primary binding mode with a single uridine bound, our simulations uncover two supplementary binding modes in which adjacent nucleotides encroach upon the binding pocket.

Refinement of the Sugar Puckering Torsion Potential in the AMBER DNA Force Field.

The transition from B-DNA to A-DNA occurs in many protein-DNA interactions or in DNA/RNA hybrid duplexes, and thus plays a role in many important biomolecular processes that convey the biological function of DNA. However, the stability of A-DNA is severely underestimated in current AMBER force fields such as OL15, OL21 or bsc1, potentially leading to unstable or deformed protein-DNA complexes. In this study, we refine the deoxyribose dihedral potential to increase the stability of the north (N) puckering present in A-DNA.

Computer Folding of Parallel DNA G-Quadruplex: Hitchhiker's Guide to the Conformational Space.

Guanine quadruplexes (GQs) play crucial roles in various biological processes, and understanding their folding pathways provides insight into their stability, dynamics, and functions. This knowledge aids in designing therapeutic strategies, as GQs are potential targets for anticancer drugs and other therapeutics. Although experimental and theoretical techniques have provided valuable insights into different stages of the GQ folding, the structural complexity of GQs poses significant challenges, and our understanding remains incomplete.

The Drosophila RNA binding protein Hrp48 binds a specific RNA sequence of the msl-2 mRNA 3' UTR to regulate translation.

Repression of msl-2 mRNA translation is essential for viability of Drosophila melanogaster females to prevent hypertranscription of both X chromosomes. This translational control event is coordinated by the female-specific protein Sex-lethal (Sxl) which recruits the RNA binding proteins Unr and Hrp48 to the 3' untranslated region (UTR) of the msl-2 transcript and represses translation initiation. The mechanism exerted by Hrp48 during translation repression and its interaction with msl-2 are not well understood.

Molecular Simulations to Investigate the Impact of N6-Methylation in RNA Recognition: Improving Accuracy and Precision of Binding Free Energy Prediction.

N6-Methyladenosine (mA) is a prevalent RNA post-transcriptional modification that plays crucial roles in RNA stability, structural dynamics, and interactions with proteins. The YT521-B (YTH) family of proteins, which are notable mA readers, functions through its highly conserved YTH domain. Recent structural investigations and molecular dynamics (MD) simulations have shed light on the mechanism of recognition of mA by the YTHDC1 protein.

N-terminal domain of polypyrimidine-tract binding protein is a dynamic folding platform for adaptive RNA recognition.

The N-terminal RNA recognition motif domain (RRM1) of polypyrimidine tract binding protein (PTB) forms an additional C-terminal helix α3, which docks to one edge of the β-sheet upon binding to a stem-loop RNA containing a UCUUU pentaloop. Importantly, α3 does not contact the RNA. The α3 helix therefore represents an allosteric means to regulate the conformation of adjacent domains in PTB upon binding structured RNAs.

Phosphoric acid salts of amino acids as a source of oligopeptides on the early Earth.

Because of their unique proton-conductivity, chains of phosphoric acid molecules are excellent proton-transfer catalysts. Here we demonstrate that this property could have been exploited for the prebiotic synthesis of the first oligopeptide sequences on our planet. Our results suggest that drying highly diluted solutions containing amino acids (like glycine, histidine and arginine) and phosphates in comparable concentrations at elevated temperatures (ca.

Comprehensive Assessment of Force-Field Performance in Molecular Dynamics Simulations of DNA/RNA Hybrid Duplexes.

Mixed double helices formed by RNA and DNA strands, commonly referred to as hybrid duplexes or hybrids, are essential in biological processes like transcription and reverse transcription. They are also important for their applications in CRISPR gene editing and nanotechnology. Yet, despite their significance, the hybrid duplexes have been seldom modeled by atomistic molecular dynamics methodology, and there is no benchmark study systematically assessing the force-field performance.

Structure of an internal loop motif with three consecutive U•U mismatches from stem-loop 1 in the 3'-UTR of the SARS-CoV-2 genomic RNA.

The single-stranded RNA genome of SARS-CoV-2 is highly structured. Numerous helical stem-loop structures interrupted by mismatch motifs are present in the functionally important 5'- and 3'-UTRs. These mismatches modulate local helical geometries and feature unusual arrays of hydrogen bonding donor and acceptor groups.

Mechanical Stability and Unfolding Pathways of Parallel Tetrameric G-Quadruplexes Probed by Pulling Simulations.

Guanine quadruplex (GQ) is a noncanonical nucleic acid structure formed by guanine-rich DNA and RNA sequences. Folding of GQs is a complex process, where several aspects remain elusive, despite being important for understanding structure formation and biological functions of GQs. Pulling experiments are a common tool for acquiring insights into the folding landscape of GQs.

Photoinduced charge separation and DNA self-repair depend on sequence directionality and stacking pattern.

Charge separation is one of the most common consequences of the absorption of UV light by DNA. Recently, it has been shown that this process can enable efficient self-repair of cyclobutane pyrimidine dimers (CPDs) in specific short DNA oligomers such as the GAT[double bond, length as m-dash]T sequence. The mechanism was characterized as sequential electron transfer through the nucleobase stack which is controlled by the redox potentials of nucleobases and their sequence.

Molecular dynamics simulations reveal the parallel stranded d(GGGA)GGG DNA quadruplex folds via multiple paths from a coil-like ensemble.

G-quadruplexes (G4s) are non-canonical nucleic acid structures that fold through complex processes. Characterization of the G4 folding landscape may help to elucidate biological roles of G4s but is challenging both experimentally and computationally. Here, we achieved complete folding of a three-quartet parallel DNA G4 with (GGGA)GGG sequence using all-atom explicit-solvent enhanced-sampling molecular dynamics (MD) simulations.

DNA Quadruplex Structure with a Unique Cation Dependency.

DNA quadruplex structures provide an additional layer of regulatory control in genome maintenance and gene expression and are widely used in nanotechnology. We report the discovery of an unprecedented tetrastranded structure formed from a native G-rich DNA sequence originating from the telomeric region of Caenorhabditis elegans. The structure is defined by multiple properties that distinguish it from all other known DNA quadruplexes.

The tautomer-specific excited state dynamics of 2,6-diaminopurine using resonance-enhanced multiphoton ionization and quantum chemical calculations.

2,6-Diaminopurine (2,6-dAP) is an alternative nucleobase that potentially played a role in prebiotic chemistry. We studied its excited state dynamics in the gas phase by REMPI, IR-UV hole burning, and ps pump-probe spectroscopy and performed quantum chemical calculations at the SCS-ADC(2) level of theory to interpret the experimental results. We found the 9H tautomer to have a small barrier to ultrafast relaxation via puckering of its 6-membered ring.

Simple Adjustment of Intranucleotide Base-Phosphate Interaction in the OL3 AMBER Force Field Improves RNA Simulations.

Molecular dynamics (MD) simulations represent an established tool to study RNA molecules. The outcome of MD studies depends, however, on the quality of the force field (). Here we suggest a correction for the widely used AMBER OL3 by adding a simple adjustment of the nonbonded parameters.

Prebiotic Synthesis of 3',5'-Cyclic Adenosine and Guanosine Monophosphates through Carbodiimide-Assisted Cyclization.

3',5'-Cyclic nucleotides play a fundamental role in modern biochemical processes and have been suggested to have played a central role at the origin of terrestrial life. In this work, we suggest that a formamide-based systems chemistry might account for their availability on the early Earth. In particular, we demonstrate that in a liquid formamide environment at elevated temperatures 3',5'-cyclic nucleotides are obtained in good yield and selectivity upon intramolecular cyclization of 5'-phosphorylated nucleosides in the presence of carbodiimides.

Crystallization as a selection force at the polymerization of nucleotides in a prebiotic context.

Accumulation and selection of nucleotides is one of the most challenging problems surrounding the origin of the first RNA molecules on our planet. In the current work we propose that guanosine 3',5' cyclic monophosphate could selectively crystallize upon evaporation of an acidic prebiotic pool containing various other nucleotides. The conditions of the evaporative crystallization are fully compatible with the subsequent acid catalyzed polymerization of this cyclic nucleotide reported in earlier studies and may be relevant in a broad range of possible prebiotic environments.

The Reactivity-Enhancing Role of Water Clusters in Ammonia Aqueous Solutions.

Among the many prototypical acid-base systems, ammonia aqueous solutions hold a privileged place, owing to their omnipresence in various planets and their universal solvent character. Although the theoretical optimal water-ammonia molar ratio to form NH and OH ion pairs is 50:50, our molecular dynamics simulations show that the tendency of forming these ionic species is inversely (directly) proportional to the amount of ammonia (water) in ammonia aqueous solutions, up to a water-ammonia molar ratio of ∼75:25. Here we prove that the reactivity of these liquid mixtures is rooted in peculiar microscopic patterns emerging at the H-bonding scale, where the highly orchestrated motion of 5 solvating molecules modulates proton transfer events through local electric fields.

Complexity of Guanine Quadruplex Unfolding Pathways Revealed by Atomistic Pulling Simulations.

Guanine quadruplexes (GQs) are non-canonical nucleic acid structures involved in many biological processes. GQs formed in single-stranded regions often need to be unwound by cellular machinery, so their mechanochemical properties are important. Here, we performed steered molecular dynamics simulations of human telomeric GQs to study their unfolding.

Assessing the Current State of Amber Force Field Modifications for DNA─2023 Edition.

Advances in molecular dynamics (MD) software alongside enhanced computational power and hardware have allowed for MD simulations to significantly expand our knowledge of biomolecular structure, dynamics, and interactions. Furthermore, it has allowed for the extension of conformational sampling times from nanoseconds to the microsecond level and beyond. This has not only made convergence of conformational ensembles through comprehensive sampling possible but consequently exposed deficiencies and allowed the community to overcome limitations in the available force fields.

Atomistic Picture of Opening-Closing Dynamics of DNA Holliday Junction Obtained by Molecular Simulations.

Holliday junction (HJ) is a noncanonical four-way DNA structure with a prominent role in DNA repair, recombination, and DNA nanotechnology. By rearranging its four arms, HJ can adopt either closed or open state. With enzymes typically recognizing only a single state, acquiring detailed knowledge of the rearrangement process is an important step toward fully understanding the biological function of HJs.

Sensitivity of the RNA Structure to Ion Conditions as Probed by Molecular Dynamics Simulations of Common Canonical RNA Duplexes.

RNA molecules play a key role in countless biochemical processes. RNA interactions, which are of highly diverse nature, are determined by the fact that RNA is a highly negatively charged polyelectrolyte, which leads to intimate interactions with an ion atmosphere. Although RNA molecules are formally single-stranded, canonical (Watson-Crick) duplexes are key components of folded RNAs.