Structural, Dynamical, and Entropic Differences between SARS-CoV and SARS-CoV-2 s2m Elements Using Molecular Dynamics Simulations

The functional role of the highly conserved stem-loop II motif (s2m) in SARS-CoV and SARS-CoV-2 in the viral lifecycle remains enigmatic and an intense area of research. Structure and dynamics of the s2m are key to establishing a structure-function connection, yet a full set of atomistic resolution coordinates is not available for SARS-CoV-2 s2m. Our work constructs three-dimensional coordinates consistent with NMR solution phase data for SARS-CoV-2 s2m and provides a comparative analysis with its counterpart SARS-CoV s2m. We employed initial coordinates based on PDB ID 1XJR for SARS-CoV s2m and two models for SARS-CoV-2 s2m: one based on 1XJR in which we introduced the mutations present in SARS-CoV-2 s2m and the second based on the available SARS-CoV-2 NMR NOE data supplemented with knowledge-based methods. For each of the three systems, 3.5 μs molecular dynamics simulations were used to sample the structure and dynamics, and principal component analysis (PCA) reduced the ensembles to hierarchal conformational substates for detailed analysis. Dilute solution simulations of SARS-CoV s2m demonstrate that the GNRA-like terminal pentaloop is rigidly defined by base stacking uniquely positioned for possible kissing dimer formation. However, the SARS-CoV-2 s2m simulation did not retain the reported crystallographic SARS-CoV motifs and the terminal loop expands to a highly dynamic "nonaloop." Increased flexibility and structural disorganization are observed for the larger terminal loop, where an entropic penalty is computed to explain the experimentally observed reduction in kissing complex formation. Overall, both SARS-CoV and SARS-CoV-2 s2m elements have a similarly pronounced L-shape due to different motif interactions. Our study establishes the atomistic three-dimensional structure and uncovers dynamic differences that arise from s2m sequence changes, which sets the stage for the interrogation of different mechanistic pathways of suspected biological function.

Deciphering nucleotide modification-induced structure and stability changes

Nucleotide modification in RNA controls a bevy of biological processes, including RNA degradation, gene expression, and gene editing. In turn, misregulation of modified nucleotides is associated with a host of chronic diseases and disorders. However, the molecular mechanisms driving these processes remain poorly understood. To partially address this knowledge gap, we used alchemical and temperature replica exchange molecular dynamics (TREMD) simulations on an RNA duplex and an analogous hairpin to probe the structural effects of modified and/or mutant nucleotides. The simulations successfully predict the modification/mutation-induced relative free energy change for complementary duplex formation, and structural analyses highlight mechanisms driving stability changes. Furthermore, TREMD simulations for a hairpin-forming RNA with and without modification provide reliable estimations of the energy landscape. Illuminating the impact of methylated and/or mutated nucleotides on the structure-function relationship and the folding energy landscape, the simulations provide insights into modification-induced alterations to the folding mechanics of the hairpin. The results here may be biologically significant as hairpins are widespread structure motifs that play critical roles in gene expression and regulation. Specifically, the tetraloop of the probed hairpin is phylogenetically abundant, and the stem mirrors a miRNA seed region whose modification has been implicated in epilepsy pathogenesis.

Transcription termination: nucleotide sequence at 3' end of tryptophan operon in Escherichia coli

We have determined the RNA and DNA sequences in the region specifying termination of transcription at the end of the tryptophan (trp) operon of Escherichia coli. A 3'-terminal mRNA fragment of about 150 nucleotides yielded oligonucleotide products that could be assigned to the end of trpA (the last structural gene in the operon) by correlation with the amino acid sequence of the protein product. Analysis of the DNA corresponding to this region served to align the few noncoding RNA oligonucleotide sequences and demonstrated that termination of trp transcription occurs in vivo at a site 36 nucleotides after trpA, with greater than 95% efficiency. In two different strains partially defective in the transcription termination factor rho, the purified transcript is much longer and more complex, suggesting that a significant amount of read-through occurs in these strains. This is consistent with evidence [Guarente, L. P., Mitchell, D. H. & Beckwith, J. (1977) J. Mol. Biol. 112, 423-436] that efficient termination in vivo at the end of the trp operon is a rho-dependent event. The trp terminator (trp t) shares several features with other known sites of transcription termination, including (i) a 3'-terminal RNA sequence of several uridine residues, C-A-U-U-U-U(OH), (ii) a G.C-rich region in the DNA immediately preceding the site of termination, followed by an A.T-rich region, and (iii) a region of dyad symmetry in the DNA which, in the transcript, is capable of forming a stable hairpin containing seven G.C base pairs and one A.U base pair in its stem.