Delta SARS-CoV-2 s2m Structure, Dynamics, and Entropy: Consequences of the G15U Mutation

Bioinformatic analysis of the Delta SARS-CoV-2 genome reveals a single nucleotide mutation (G15U) in the stem-loop II motif (s2m) relative to ancestral SARS-CoV-2. Despite sequence similarity, unexpected differences between SARS-CoV-2 and Delta SARS-CoV-2 s2m homodimerization experiments require the discovery of unknown structural and thermodynamic changes necessary to rationalize the data. Using our reported SARS-CoV-2 s2m model, we induced the G15U substitution and performed 3.5 microseconds of unbiased molecular dynamics simulation at 283 and 310 K. The resultant Delta s2m adopted a secondary structure consistent with our reported NMR data, resulting in significant deviations in the tertiary structure and dynamics from our SARS-CoV-2 s2m model. First, we find differences in the overall three-dimensional structure, where the characteristic 90° L-shaped kink of the SARS-CoV-2 s2m did not form in the Delta s2m resulting in a “linear” hairpin with limited bending dynamics. Delta s2m helical parameters are calculated to align closely with A-form RNA, effectively eliminating a hinge point to form the L-shape kink by correcting an upper stem defect in SARS-CoV-2 induced by a noncanonical and dynamic G:A base pair. Ultimately, the shape difference rationalizes the migration differences in reported electrophoresis experiments. Second, increased fluctuation of the Delta s2m palindromic sequence, within the terminal loop, compared to SARS-CoV-2 s2m results in an estimated increase of entropy of 6.8 kcal/mol at 310 K relative to the SARS-CoV-2 s2m. The entropic difference offers a unique perspective on why the Delta s2m homodimerizes less spontaneously, forming fewer kissing dimers and extended duplexes compared to SARS-CoV-2. In this work, both the L-shape reduction and palindromic entropic penalty provides an explanation of our reported in vitro electrophoresis homodimerization results. Ultimately, the structural, dynamical, and entropic differences between the SARS-CoV-2 s2m and Delta s2m serve to establish a foundation for future studies of the s2m function in the viral lifecycle.

[1]  J. Evanseck,et al.  Effect of the SARS-CoV-2 Delta-associated G15U mutation on the s2m element dimerization and its interactions with miR-1307-3p , 2023, bioRxiv.

[2]  J. Evanseck,et al.  Structural, Dynamical, and Entropic Differences between SARS-CoV and SARS-CoV-2 s2m Elements Using Molecular Dynamics Simulations , 2022, ACS physical chemistry Au.

[3]  J. Evanseck,et al.  Bioinformatics analysis of the s2m mutations within the SARS‐CoV‐2 Omicron lineages , 2022, Journal of medical virology.

[4]  R. Granek,et al.  Managing an evolving pandemic: Cryptic circulation of the Delta variant during the Omicron rise , 2022, Science of The Total Environment.

[5]  A. Torres,et al.  SARS-CoV-2: Evolution and Emergence of New Viral Variants , 2022, Viruses.

[6]  P. Fournier,et al.  Culture and identification of a “Deltamicron” SARS‐CoV‐2 in a three cases cluster in southern France , 2022, medRxiv.

[7]  Kang Zhang,et al.  SARS‐CoV‐2 Omicron variant: Immune escape and vaccine development , 2022, MedComm.

[8]  N. Daneman,et al.  Estimates of SARS-CoV-2 Omicron Variant Severity in Ontario, Canada. , 2022, JAMA.

[9]  H. Fan,et al.  SARS-CoV-2 variant Omicron: currently the most complete “escapee” from neutralization by antibodies and vaccines , 2022, Signal Transduction and Targeted Therapy.

[10]  S. Weaver,et al.  Neutralization against Omicron SARS-CoV-2 from previous non-Omicron infection , 2021, bioRxiv.

[11]  J. Bhiman,et al.  Early assessment of the clinical severity of the SARS-CoV-2 Omicron variant in South Africa , 2021, medRxiv.

[12]  Q. Ye,et al.  The global epidemic of SARS‐CoV‐2 variants and their mutational immune escape , 2021, Journal of medical virology.

[13]  S. Madhi,et al.  Reduced neutralization of SARS-CoV-2 B.1.617 by vaccine and convalescent serum , 2021, Cell.

[14]  W. Scott,et al.  Targeting the Conserved Stem Loop 2 Motif in the SARS-CoV-2 Genome , 2021, Journal of virology.

[15]  Shaghayegh Haghjooy Javanmard,et al.  Novel and emerging mutations of SARS-CoV-2: Biomedical implications , 2021, Biomedicine & Pharmacotherapy.

[16]  F. Massi,et al.  Analysis of Emerging Variants in Structured Regions of the SARS-CoV-2 Genome , 2021, Evolutionary bioinformatics online.

[17]  Robert G. Brinson,et al.  SARS-CoV-2 highly conserved s2m element dimerizes via a kissing complex and interacts with host miRNA-1307-3p , 2020, bioRxiv.

[18]  K. Saxena,et al.  Secondary structure determination of conserved SARS-CoV-2 RNA elements by NMR spectroscopy , 2020, Nucleic acids research.

[19]  K. Crandall,et al.  Genome-wide analysis of SARS-CoV-2 virus strains circulating worldwide implicates heterogeneity , 2020, Scientific Reports.

[20]  T. Yeh,et al.  Emerging viral mutants in Australia suggest RNA recombination event in the SARS‐CoV‐2 genome , 2020, The Medical journal of Australia.

[21]  Xiang-Jun Lu DSSR-enabled innovative schematics of 3D nucleic acid structures with PyMOL , 2020, Nucleic acids research.

[22]  A. McLellan,et al.  Implications of SARS-CoV-2 Mutations for Genomic RNA Structure and Host microRNA Targeting , 2020, bioRxiv.

[23]  Katja Petzold,et al.  RNA Dynamics by NMR Spectroscopy , 2019, Chembiochem : a European journal of chemical biology.

[24]  Xiang-Jun Lu,et al.  Web 3DNA 2.0 for the analysis, visualization, and modeling of 3D nucleic acid structures , 2019, Nucleic Acids Res..

[25]  Sandro Bottaro,et al.  Barnaba: software for analysis of nucleic acid structures and trajectories , 2018, bioRxiv.

[26]  Richard A. Cunha,et al.  RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview , 2018, Chemical reviews.

[27]  Stefan Elbe,et al.  Data, disease and diplomacy: GISAID's innovative contribution to global health , 2017, Global challenges.

[28]  Thomas J Lane,et al.  MDTraj: a modern, open library for the analysis of molecular dynamics trajectories , 2014, bioRxiv.

[29]  Ivo L. Hofacker,et al.  Forna (force-directed RNA): Simple and effective online RNA secondary structure diagrams , 2015, Bioinform..

[30]  Robert T. McGibbon,et al.  MDTraj: a modern, open library for the analysis of molecular dynamics trajectories , 2014, bioRxiv.

[31]  Michal Otyepka,et al.  Simulations of A-RNA duplexes. The effect of sequence, solute force field, water model, and salt concentration. , 2012, The journal of physical chemistry. B.

[32]  J. Šponer,et al.  Refinement of the Cornell et al. Nucleic Acids Force Field Based on Reference Quantum Chemical Calculations of Glycosidic Torsion Profiles , 2011, Journal of chemical theory and computation.

[33]  Gaël Varoquaux,et al.  Scikit-learn: Machine Learning in Python , 2011, J. Mach. Learn. Res..

[34]  J. Šponer,et al.  Molecular Dynamics and Quantum Mechanics of RNA: Conformational and Chemical Change We Can Believe In , 2009, Accounts of chemical research.

[35]  John D. Hunter,et al.  Matplotlib: A 2D Graphics Environment , 2007, Computing in Science & Engineering.

[36]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[37]  D. Haussler,et al.  The Structure of a Rigorously Conserved RNA Element within the SARS Virus Genome , 2004, PLoS biology.

[38]  K. Réblová,et al.  Long-residency hydration, cation binding, and dynamics of loop E/helix IV rRNA-L25 protein complex. , 2004, Biophysical journal.

[39]  W. V. van Gunsteren,et al.  Estimating entropies from molecular dynamics simulations. , 2004, The Journal of chemical physics.

[40]  H. Schwalbe,et al.  NMR Spectroscopy of RNA , 2003, Chembiochem : a European journal of chemical biology.

[41]  Ioan Andricioaei,et al.  On the calculation of entropy from covariance matrices of the atomic fluctuations , 2001 .

[42]  E. Westhof,et al.  Geometric nomenclature and classification of RNA base pairs. , 2001, RNA.

[43]  W. Olson,et al.  A-form conformational motifs in ligand-bound DNA structures. , 2000, Journal of molecular biology.

[44]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[45]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[46]  J. Schlitter Estimation of absolute and relative entropies of macromolecules using the covariance matrix , 1993 .

[47]  H. Berendsen,et al.  Essential dynamics of proteins , 1993, Proteins.

[48]  R D Young,et al.  Protein states and proteinquakes. , 1985, Proceedings of the National Academy of Sciences of the United States of America.

[49]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[50]  OUP accepted manuscript , 2021, Nucleic Acids Research.

[51]  Charles C. David,et al.  Principal component analysis: a method for determining the essential dynamics of proteins. , 2014, Methods in molecular biology.

[52]  Jeffrey D. Evanseck,et al.  Chapter 13 Principal Components Analysis: A Review of its Application on Molecular Dynamics Data , 2006 .

[53]  C. Brooks Computer simulation of liquids , 1989 .

[54]  J. MacQueen Some methods for classification and analysis of multivariate observations , 1967 .