Effect of the SARS-CoV-2 Delta-associated G15U mutation on the s2m element dimerization and its interactions with miR-1307-3p.

The s2m, a highly conserved 41-nt hairpin structure in the SARS-CoV-2 genome, serves as an attractive therapeutic target that may have important roles in the virus life cycle or interactions with the host. However, the conserved s2m in Delta SARS-CoV-2, a previously dominant variant characterized by high infectivity and disease severity, has received relatively less attention than that of the original SARS-CoV-2 virus. The focus of this work is to identify and define the s2m changes between Delta and SARS-CoV-2 and the subsequent impact of those changes upon the s2m dimerization and interactions with the host microRNA miR-1307-3p. Bioinformatics analysis of the GISAID database targeting the s2m element reveals a >99% correlation of a single nucleotide mutation at the 15th position (G15U) in Delta SARS-CoV-2. Based on 1H NMR spectroscopy assignments comparing the imino proton resonance region of s2m and the s2m G15U at 19°C, we show that the U15-A29 base pair closes, resulting in a stabilization of the upper stem without overall secondary structure deviation. Increased stability of the upper stem did not affect the chaperone activity of the viral N protein, as it was still able to convert the kissing dimers formed by s2m G15U into a stable duplex conformation, consistent with the s2m reference. However, we show that the s2m G15U mutation drastically impacts the binding of host miR-1307-3p. These findings demonstrate that the observed G15U mutation alters the secondary structure of s2m with subsequent impact on viral binding of host miR-1307-3p, with potential consequences on immune responses.

Highly conserved s2m element of SARS-CoV-2 dimerizes via a kissing complex and interacts with host miRNA-1307-3p.

The ongoing COVID-19 pandemic highlights the necessity for a more fundamental understanding of the coronavirus life cycle. The causative agent of the disease, SARS-CoV-2, is being studied extensively from a structural standpoint in order to gain insight into key molecular mechanisms required for its survival. Contained within the untranslated regions of the SARS-CoV-2 genome are various conserved stem-loop elements that are believed to function in RNA replication, viral protein translation, and discontinuous transcription. While the majority of these regions are variable in sequence, a 41-nucleotide s2m element within the genome 3' untranslated region is highly conserved among coronaviruses and three other viral families. In this study, we demonstrate that the SARS-CoV-2 s2m element dimerizes by forming an intermediate homodimeric kissing complex structure that is subsequently converted to a thermodynamically stable duplex conformation. This process is aided by the viral nucleocapsid protein, potentially indicating a role in mediating genome dimerization. Furthermore, we demonstrate that the s2m element interacts with multiple copies of host cellular microRNA (miRNA) 1307-3p. Taken together, our results highlight the potential significance of the dimer structures formed by the s2m element in key biological processes and implicate the motif as a possible therapeutic drug target for COVID-19 and other coronavirus-related diseases.

DOTAP: Structure, hydration, and the counterion effect.

The cationic lipid 1,2-dioleoyl-3-trimethylammonium propane (DOTAP) is one of the original synthetic cationic lipids used for the liposomal transfection of oligonucleotides in gene therapy. The key structural element of DOTAP is its quaternary ammonium headgroup that is responsible for interactions with both nucleic acids and target cell membranes. Because these interactions are fundamental to the design of a major class of transfection lipids, it is important to understand the structure of DOTAP and how it interacts with halide counterions. Here, we use x-ray and neutron diffraction techniques to examine the structure of DOTAP and how chloride (Cl-) and iodide (I-) counterions alter the hydration properties of the DOTAP headgroup. A problem of particular interest is the poor solubility of DOTAP/I- in water solutions. Our results show that the poor solubility results from very tight binding of the I- counterion to the headgroup and the consequent expulsion of water. The structural principles we report here are important for assessing the suitability of DOTAP and its quaternary ammonium derivatives for transfection.

Uridylation of the histone mRNA stem-loop weakens binding interactions with SLBP while maintaining interactions with 3'hExo.

Histone mRNA degradation is controlled by the unique 3' stem-loop of histone mRNA and the stem-loop binding protein (SLBP). As part of this process, the 3' stem-loop is trimmed by the histone-specific 3' exonuclease (3'hExo) and uridylated by the terminal uridylyl transferase 7 (TUT7), creating partially degraded intermediates with short uridylations. The role of these uridylations in degradation is not fully understood. Our work examines changes in the stability of the ternary complex created by trimming and uridylation of the stem-loop to better understand the role of this process in the histone mRNA life cycle. In this study, we used fluorescence polarization and electrophoretic mobility shift assays to demonstrate that both SLBP and 3'hExo can bind to uridylated and partially degraded stem-loop intermediates, although with lower affinity. We further characterized this complex by performing 1-µs molecular dynamics simulations using the AMBER force field and Nanoscale Molecular Dynamics (NAMD). These simulations show that while uridylation helps maintain the overall shape of the stem-loop, the combination of uridylation and dephosphorylation of the TPNK motif in SLBP disrupts key RNA-protein interactions. They also demonstrate that uridylation allows 3'hExo to maintain contact with the stem-loop after partial degradation and plays a role in disrupting key base pairs in partially degraded histone mRNA intermediates. Together, these experiments and simulations suggest that trimming by 3'hExo, uridylation, and SLBP dephosphorylation weakens both RNA-protein interactions and the stem-loop itself. Our results further elucidate the role of uridylation and SLBP dephosphorylation in the early stages of histone mRNA degradation.

The host-defense peptide piscidin P1 reorganizes lipid domains in membranes and decreases activation energies in mechanosensitive ion channels.

The host-defense peptide (HDP) piscidin 1 (P1), isolated from the mast cells of striped bass, has potent activities against bacteria, viruses, fungi, and cancer cells and can also modulate the activity of membrane receptors. Given its broad pharmacological potential, here we used several approaches to better understand its interactions with multicomponent bilayers representing models of bacterial (phosphatidylethanolamine (PE)/phosphatidylglycerol) and mammalian (phosphatidylcholine/cholesterol (PC/Chol)) membranes. Using solid-state NMR, we solved the structure of P1 bound to PC/Chol and compared it with that of P3, a less potent homolog. The comparison disclosed that although both peptides are interfacially bound and α-helical, they differ in bilayer orientations and depths of insertion, and these differences depend on bilayer composition. Although Chol is thought to make mammalian membranes less susceptible to HDP-mediated destabilization, we found that Chol does not affect the permeabilization effects of P1. X-ray diffraction experiments revealed that both piscidins produce a demixing effect in PC/Chol membranes by increasing the fraction of the Chol-depleted phase. Furthermore, P1 increased the temperature required for the lamellar-to-hexagonal phase transition in PE bilayers, suggesting that it imposes positive membrane curvature. Patch-clamp measurements on the inner Escherichia coli membrane showed that P1 and P3, at concentrations sufficient for antimicrobial activity, substantially decrease the activating tension for bacterial mechanosensitive channels. This indicated that piscidins can cause lipid redistribution and restructuring in the microenvironment near proteins. We conclude that the mechanism of piscidin's antimicrobial activity extends beyond simple membrane destabilization, helping to rationalize its broader spectrum of pharmacological effects.

Differential Interactions of Piscidins with Phospholipids and Lipopolysaccharides at Membrane Interfaces.

Piscidins 1 and 3 (P1 and P3) are potent antimicrobial peptides isolated from striped bass. Their mechanism of action involves formation of amphipathic α-helices on contact with phospholipids and destabilization of the microbial cytoplasmic membrane. The peptides are active against both Gram-positive and Gram-negative bacteria, suggesting easy passage across the outer membrane. Here, we performed a comparative study of these two piscidins at the air-water interface on lipopolysaccharide (LPS) monolayers modeling the outer bacterial surface of Gram-negative organisms and on phospholipid monolayers, which mimic the inner membrane. The results show that P1 and P3 are highly surface active (log KAW ∼ 6.8) and have similar affinities to phospholipid monolayers (log Klip ≈ 7.7). P1, which is more potent against Gram negatives, exhibits a much stronger partitioning into LPS monolayers (log KLPS = 8.3). Pressure-area isotherms indicate that under increasing lateral pressures, inserted P1 repartitions from phospholipid monolayers back to the subphase or to a more shallow position with in-plane areas of ∼170 Å2 per peptide, corresponding to fully folded amphipathic α-helices. In contrast, peptide expulsion from LPS occurs with areas of ∼35 Å2, suggesting that the peptides may not form the similarly oriented, rigid secondary structures when they avidly intercalate between LPS molecules. Patch-clamp experiments on Escherichia coli spheroplasts show that when P1 and P3 reach the outer surface of the bacterial cytoplasmic membrane, they produce fluctuating conductive structures at voltages above 80 mV. The data suggests that the strong activity of these piscidins against Gram-negative bacteria begins with the preferential accumulation of peptides in the outer LPS layer followed by penetration into the periplasm, where they form stable amphipathic α-helices upon contact with phospholipids and attack the energized inner membrane.

Characterization of a G-Quadruplex Structure in Pre-miRNA-1229 and in Its Alzheimer's Disease-Associated Variant rs2291418: Implications for miRNA-1229 Maturation.

Alzheimer's disease (AD), the most common age-related neurodegenerative disease, is associated with various forms of cognitive and functional impairment that worsen with disease progression. AD is typically characterized as a protein misfolding disease, in which abnormal plaques form due to accumulation of tau and ß-amyloid (Aß) proteins. An assortment of proteins is responsible for the processing and trafficking of Aß, including sortilin-related receptor 1 (SORL1). Recently, a genome-wide association study of microRNA-related variants found that a single nucleotide polymorphism (SNP) rs2291418 within premature microRNA-1229 (pre-miRNA-1229) is significantly associated with AD. Moreover, the levels of the mature miRNA-1229-3p, which has been shown to regulate the SORL1 translation, are increased in the rs2291418 pre-miRNA-1229 variant. In this study we used various biophysical techniques to show that pre-miRNA-1229 forms a G-quadruplex secondary structure that coexists in equilibrium with the canonical hairpin structure, potentially controlling the production of the mature miR-1229-3p, and furthermore, that the AD-associated SNP rs2291418 pre-miR-1229 changes the equilibrium between these structures. Thus, the G-quadruplex structure we identified within pre-miRNA-1229 could potentially act as a novel therapeutic target in AD.

Structure and Function in Antimicrobial Piscidins: Histidine Position, Directionality of Membrane Insertion, and pH-Dependent Permeabilization.

Piscidins are histidine-enriched antimicrobial peptides that interact with lipid bilayers as amphipathic α-helices. Their activity at acidic and basic pH in vivo makes them promising templates for biomedical applications. This study focuses on p1 and p3, both 22-residue-long piscidins with 68% sequence identity. They share three histidines (H3, H4, and H11), but p1, which is significantly more permeabilizing, has a fourth histidine (H17). This study investigates how variations in amphipathic character associated with histidines affect the permeabilization properties of p1 and p3. First, we show that the permeabilization ability of p3, but not p1, is strongly inhibited at pH 6.0 when the conserved histidines are partially charged and H17 is predominantly neutral. Second, our neutron diffraction measurements performed at low water content and neutral pH indicate that the average conformation of p1 is highly tilted, with its C-terminus extending into the opposite leaflet. In contrast, p3 is surface bound with its N-terminal end tilted toward the bilayer interior. The deeper membrane insertion of p1 correlates with its behavior at full hydration: an enhanced ability to tilt, bury its histidines and C-terminus, induce membrane thinning and defects, and alter membrane conductance and viscoelastic properties. Furthermore, its pH-resiliency relates to the neutral state favored by H17. Overall, these results provide mechanistic insights into how differences in the histidine content and amphipathicity of peptides can elicit different directionality of membrane insertion and pH-dependent permeabilization. This work features complementary methods, including dye leakage assays, NMR-monitored titrations, X-ray and neutron diffraction, oriented CD, molecular dynamics, electrochemical impedance spectroscopy, surface plasmon resonance, and quartz crystal microbalance with dissipation.

G quadruplex RNA structures in PSD-95 mRNA: potential regulators of miR-125a seed binding site accessibility.

Fragile X syndrome (FXS) is the most common inherited form of intellectual disability caused by the CGG trinucleotide expansion in the 3'-untranslated region of the FMR1 gene on the X chromosome, that silences the expression of the Fragile X mental retardation protein (FMRP). FMRP has been shown to bind to a G-rich region within the PSD-95 mRNA which encodes for the postsynaptic density protein 95 (PSD-95), and together with the microRNA miR-125a, to play an important role in the reversible inhibition of the PSD-95 mRNA translation in neurons. The loss of FMRP in Fmr1 KO mice disables this translation control in the production of the PSD-95 protein. Interestingly, the miR-125a binding site on PSD-95 mRNA is embedded in the G-rich region bound by FMRP and postulated to adopt one or more G quadruplex structures. In this study, we have used different biophysical techniques to validate and characterize the formation of parallel G quadruplex structures and binding of miR-125a to its complementary sequence located within the 3' UTR of PSD-95 mRNA. Our results indicate that the PSD-95 mRNA G-rich region folds into alternate G quadruplex conformations that coexist in equilibrium. miR-125a forms a stable complex with PSD-95 mRNA, as evident by characteristic Watson-Crick base-pairing that coexists with one of the G quadruplex forms, suggesting a novel mechanism for G quadruplex structures to regulate the access of miR-125a to its binding site.

Structural interactions of a voltage sensor toxin with lipid membranes.

Protein toxins from tarantula venom alter the activity of diverse ion channel proteins, including voltage, stretch, and ligand-activated cation channels. Although tarantula toxins have been shown to partition into membranes, and the membrane is thought to play an important role in their activity, the structural interactions between these toxins and lipid membranes are poorly understood. Here, we use solid-state NMR and neutron diffraction to investigate the interactions between a voltage sensor toxin (VSTx1) and lipid membranes, with the goal of localizing the toxin in the membrane and determining its influence on membrane structure. Our results demonstrate that VSTx1 localizes to the headgroup region of lipid membranes and produces a thinning of the bilayer. The toxin orients such that many basic residues are in the aqueous phase, all three Trp residues adopt interfacial positions, and several hydrophobic residues are within the membrane interior. One remarkable feature of this preferred orientation is that the surface of the toxin that mediates binding to voltage sensors is ideally positioned within the lipid bilayer to favor complex formation between the toxin and the voltage sensor.

Pore hydration states of KcsA potassium channels in membranes.

Water-filled hydrophobic cavities in channel proteins serve as gateways for transfer of ions across membranes, but their properties are largely unknown. We determined water distributions along the conduction pores in two tetrameric channels embedded in lipid bilayers using neutron diffraction: potassium channel KcsA and the transmembrane domain of M2 protein of influenza A virus. For the KcsA channel in the closed state, the distribution of water is peaked in the middle of the membrane, showing water in the central cavity adjacent to the selectivity filter. This water is displaced by the channel blocker tetrabutyl-ammonium. The amount of water associated with the channel was quantified, using neutron diffraction and solid state NMR. In contrast, the M2 proton channel shows a V-shaped water profile across the membrane, with a narrow constriction at the center, like the hourglass shape of its internal surface. These two types of water distribution are therefore very different in their connectivity to the bulk water. The water and protein profiles determined here provide important evidence concerning conformation and hydration of channels in membranes and the potential role of pore hydration in channel gating.

Biophysical characterization of G-quadruplex forming FMR1 mRNA and of its interactions with different fragile X mental retardation protein isoforms.

Fragile X syndrome, the most common form of inherited mental impairment in humans, is caused by the absence of the fragile X mental retardation protein (FMRP) due to a CGG trinucleotide repeat expansion in the 5'-untranslated region (UTR) and subsequent translational silencing of the fragile x mental retardation-1 (FMR1) gene. FMRP, which is proposed to be involved in the translational regulation of specific neuronal messenger RNA (mRNA) targets, contains an arginine-glycine-glycine (RGG) box RNA binding domain that has been shown to bind with high affinity to G-quadruplex forming mRNA structures. FMRP undergoes alternative splicing, and the binding of FMRP to a proposed G-quadruplex structure in the coding region of its mRNA (named FBS) has been proposed to affect the mRNA splicing events at exon 15. In this study, we used biophysical methods to directly demonstrate the folding of FMR1 FBS into a secondary structure that contains two specific G-quadruplexes and analyze its interactions with several FMRP isoforms. Our results show that minor splice isoforms, ISO2 and ISO3, created by the usage of the second and third acceptor sites at exon 15, bind with higher affinity to FBS than FMRP ISO1, which is created by the usage of the first acceptor site. FMRP ISO2 and ISO3 cannot undergo phosphorylation, an FMRP post-translational modification shown to modulate the protein translation regulation. Thus, their expression has to be tightly regulated, and this might be accomplished by a feedback mechanism involving the FMRP interactions with the G-quadruplex structures formed within FMR1 mRNA.

Structure and hydration of membranes embedded with voltage-sensing domains.

Despite the growing number of atomic-resolution membrane protein structures, direct structural information about proteins in their native membrane environment is scarce. This problem is particularly relevant in the case of the highly charged S1-S4 voltage-sensing domains responsible for nerve impulses, where interactions with the lipid bilayer are critical for the function of voltage-activated ion channels. Here we use neutron diffraction, solid-state nuclear magnetic resonance (NMR) spectroscopy and molecular dynamics simulations to investigate the structure and hydration of bilayer membranes containing S1-S4 voltage-sensing domains. Our results show that voltage sensors adopt transmembrane orientations and cause a modest reshaping of the surrounding lipid bilayer, and that water molecules intimately interact with the protein within the membrane. These structural findings indicate that voltage sensors have evolved to interact with the lipid membrane while keeping energetic and structural perturbations to a minimum, and that water penetrates the membrane, to hydrate charged residues and shape the transmembrane electric field.

Hepatitis C virus RNA: molecular switches mediated by long-range RNA-RNA interactions?

Multiple conserved structural cis-acting regulatory elements have been recognized both in the coding and untranslated regions (UTRs) of the hepatitis C virus (HCV) genome. For example, the cis-element 5BSL3.2 in the HCV-coding region has been predicted to use both its apical and internal loops to interact with the X RNA in the 3'-UTR, with the IIId domain in the 5'-UTR and with the Alt sequence in the coding region. Additionally, the X RNA region uses a palindromic sequence that overlaps the sequence required for the interaction with 5BSL3.2, to dimerize with another HCV genome. The ability of the 5BSL3.2 and X RNA regions to engage in multi-interactions suggests the existence of one or more molecular RNA switches which may regulate different steps of the HCV life cycle. In this study, we used biophysical methods to characterize the essential interactions of these HCV cis-elements at the molecular level. Our results indicate that X RNA interacts with 5BSL3.2 and another X RNA molecule by adopting two different conformations and that 5BSL3.2 engages simultaneously in kissing interactions using its apical and internal loops. Based on these results, we propose a mode of action for possible molecular switches involving the HCV RNA.

FMRP interacts with G-quadruplex structures in the 3'-UTR of its dendritic target Shank1 mRNA.

Fragile X syndrome (FXS), the most common cause of inherited intellectual disability, is caused by the loss of expression of the fragile X mental retardation protein (FMRP). FMRP, which regulates the transport and translation of specific mRNAs, uses its RGG box domain to bind mRNA targets that form G-quadruplex structures. One of the FMRP in vivo targets, Shank1 mRNA, encodes the master scaffold proteins of the postsynaptic density (PSD) which regulate the size and shape of dendritic spines because of their capacity to interact with many different PSD components. Due to their effect on spine morphology, altered translational regulation of Shank1 transcripts may contribute to the FXS pathology. We hypothesized that the FMRP interactions with Shank1 mRNA are mediated by the recognition of the G quadruplex structure, which has not been previously demonstrated. In this study we used biophysical techniques to analyze the Shank1 mRNA 3'-UTR and its interactions with FMRP and its phosphorylated mimic FMRP S500D. We found that the Shank1 mRNA 3 ' -UTR adopts two very stable intramolecular G-quadruplexes which are bound specifically and with high affinity by FMRP both in vitro and in vivo. These results suggest a role of G-quadruplex RNA motif as a structural element in the common mechanism of FMRP regulation of its dendritic mRNA targets.

A membrane-translocating peptide penetrates into bilayers without significant bilayer perturbations.

Using a high throughput screen, we have identified a family of 12-residue long peptides that spontaneously translocate across membranes. These peptides function by a poorly understood mechanism that is very different from that of the well-known, highly cationic cell penetrating peptides such as the tat peptide from HIV. The newly discovered translocating peptides can carry polar cargoes across synthetic bilayers and across cellular membranes quickly and spontaneously without disrupting the membrane. Here we report on the biophysical characterization of a representative translocating peptide from the selected family, TP2, as well as a negative control peptide, ONEG, from the same library. We measured the binding of the two peptides to lipid bilayers, their secondary structure propensities, their dispositions in bilayers by neutron diffraction, and the response of the bilayer to the peptides. Compared to the negative control, TP2 has a greater propensity for membrane partitioning, although it still binds only weakly, and a higher propensity for secondary structure. Perhaps most revealing, TP2 has the ability to penetrate deep into the bilayer without causing significant bilayer perturbations, a property that may help explain its ability to translocate without bilayer permeabilization.

Host microRNA interactions with the SARS-CoV-2 viral genome 3'-untranslated region.

The 2019 pandemic, caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has marked the spread of a novel human coronavirus. While the viral life cycle is well understood, most of the interactions at the virus-host interface remain elusive. Furthermore, the molecular mechanisms behind disease severity and immune evasion are still largely unknown. Conserved elements of the viral genome such as secondary structures within the 5'- and 3'-untranslated regions (UTRs) serve as attractive targets of interest and could prove crucial in furthering our understanding of virus-host interactions. It has been proposed that microRNA (miR) interactions with viral components could be used by both the virus and host for their own benefit. Analysis of the SARS-CoV-2 viral genome 3'-UTR has revealed the potential for host cellular miR binding sites, providing sites for specific interactions with the virus. In this study, we demonstrate that the SARS-CoV-2 genome 3'-UTR binds the host cellular miRNAs miR-760-3p, miR-34a-5p, and miR-34b-5p, which have been shown to influence translation of interleukin-6 (IL-6), the IL-6 receptor (IL-6R), as well as progranulin (PGRN), respectively, proteins that have roles in the host immune response and inflammatory pathways. Furthermore, recent work suggests the potential of miR-34a-5p and miR-34b-5p to target and inhibit translation of viral proteins. Native gel electrophoresis and steady-state fluorescence spectroscopy were utilized to characterize the binding of these miRs to their predicted sites within the SARS-CoV-2 genome 3'-UTR. Additionally, we investigated 2'-fluoro-D-arabinonucleic acid (FANA) analogs of these miRNAs as competitive binding inhibitors for these miR binding interactions. The mechanisms detailed in this study have the potential to drive the development of antiviral treatments for SARS-CoV-2 infection, and provide a potential molecular basis for cytokine release syndrome and immune evasion which could implicate the host-virus interface.

Parameterization of the miniPEG-Modified γPNA Backbone: Toward Induced γPNA Duplex Dissociation.

γ-Modified peptide nucleic acids (γPNAs) serve as potential therapeutic agents against genetic diseases. Miniature poly(ethylene glycol) (miniPEG) has been reported to increase solubility and binding affinity toward genetic targets, yet details of γPNA structure and dynamics are not understood. Within our work, we parameterized missing torsional and electrostatic terms for the miniPEG substituent on the γ-carbon atom of the γPNA backbone in the CHARMM force field. Microsecond timescale molecular dynamics simulations were carried out on six miniPEG-modified γPNA duplexes from NMR structures (PDB ID: 2KVJ). Three NMR models for the γPNA duplex (PDB ID: 2KVJ) were simulated as a reference for structural and dynamic changes captured for the miniPEG-modified γPNA duplex. Principal component analysis performed on the γPNA backbone atoms identified a single isotropic conformational substate (CS) for the NMR simulations, whereas four anisotropic CSs were identified for the ensemble of miniPEG-modified γPNA simulations. The NMR structures were found to have a 23° helical bend toward the major groove, consistent with our simulated CS structure of 19.0°. However, a significant difference between simulated methyl- and miniPEG-modified γPNAs involved the opportunistic invasion of miniPEG through the minor and major groves. Specifically, hydrogen bond fractional analysis showed that the invasion was particularly prone to affect the second G-C base pair, reducing the Watson-Crick base pair hydrogen bond by 60% over the six simulations, whereas the A-T base pairs decreased by only 20%. Ultimately, the invasion led to base stack reshuffling, where the well-ordered base stacking was reduced to segmented nucleobase stacking interactions. Our 6 µs timescale simulations indicate that duplex dissociation suggests the onset toward γPNA single strands, consistent with the experimental observation of decreased aggregation. To complement the insight of miniPEG-modified γPNA structure and dynamics, the new miniPEG force field parameters allow for further exploration of such modified γPNA single strands as potential therapeutic agents against genetic diseases.

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.