Topological Constraints and Their Conformational Entropic Penalties on RNA Folds

SELEX based aptamers with diagnostic and entry inhibitor therapeutic potential for SARS-CoV-2

Frequent mutation and variable immunological protection against vaccination is a common feature for COVID-19 pandemic. Early detection and confinement remain key to controlling further spread of infection. In response, we have developed an aptamer-based system that possesses both diagnostic and therapeutic potential towards the virus. A random aptamer library (~ 1017 molecules) was screened using systematic evolution of ligands by exponential enrichment (SELEX) and aptamer R was identified as a potent binder for the SARS-CoV-2 spike receptor binding domain (RBD) using in vitro binding assay. Using a pseudotyped viral entry assay we have shown that aptamer R specifically inhibited the entry of a SARS-CoV-2 pseudotyped virus in HEK293T-ACE2 cells but did not inhibit the entry of a Vesicular Stomatitis Virus (VSV) glycoprotein (G) pseudotyped virus, hence establishing its specificity towards SARS-CoV-2 spike protein. The antiviral potential of aptamers R and J (same central sequence as R but lacking flanked primer regions) was tested and showed 95.4% and 82.5% inhibition, respectively, against the SARS-CoV-2 virus. Finally, intermolecular interactions between the aptamers and the RBD domain were analyzed using in silico docking and molecular dynamics simulations that provided additional insight into the binding and inhibitory action of aptamers R and J.

Thermo-Programmed Synthetic DNA-Based Receptors

Herein, we present a generalizable and versatile strategy to engineer synthetic DNA ligand-binding devices that can be programmed to load and release a specific ligand at a defined temperature. We do so by re-engineering two model DNA-based receptors: a triplex-forming bivalent DNA-based receptor that recognizes a specific DNA sequence and an ATP-binding aptamer. The temperature at which these receptors load/release their ligands can be finely modulated by controlling the entropy associated with the linker connecting the two ligand-binding domains. The availability of a set of receptors with tunable and reversible temperature dependence allows achieving complex load/release behavior such as sustained ligand release over a wide temperature range. Similar programmable thermo-responsive synthetic ligand-binding devices can be of utility in applications such as drug delivery and production of smart materials.

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.

Thermodynamics of unfolding mechanisms of mouse mammary tumor virus pseudoknot from a coarse-grained loop-entropy model

Pseudoknotted RNA molecules play important biological roles that depend on their folded structure. To understand the underlying principles that determine their thermodynamics and folding/unfolding mechanisms, we carried out a study on a variant of the mouse mammary tumor virus pseudoknotted RNA (VPK), a widely studied model system for RNA pseudoknots. Our method is based on a coarse-grained discrete-state model and the algorithm of PK3D (pseudoknot structure predictor in three-dimensional space), with RNA loops explicitly constructed and their conformational entropic effects incorporated. Our loop entropy calculations are validated by accurately capturing previously measured melting temperatures of RNA hairpins with varying loop lengths. For each of the hairpins that constitutes the VPK, we identified alternative conformations that are more stable than the hairpin structures at low temperatures and predicted their populations at different temperatures. Our predictions were validated by thermodynamic experiments on these hairpins. We further computed the heat capacity profiles of VPK, which are in excellent agreement with available experimental data. Notably, our model provides detailed information on the unfolding mechanisms of pseudoknotted RNA. Analysis of the distribution of base-pairing probability of VPK reveals a cooperative unfolding mechanism instead of a simple sequential unfolding of first one stem and then the other. Specifically, we find a simultaneous "loosening" of both stems as the temperature is raised, whereby both stems become partially melted and co-exist during the unfolding process.

Vfold2D-MC: A Physics-Based Hybrid Model for Predicting RNA Secondary Structure Folding

Accurate prediction of RNA structure and folding stability has a far-reaching impact on our understanding of RNA functions. Here we develop Vfold2D-MC, a new physics-based model, to predict RNA structure and folding thermodynamics from the sequence. The model employs virtual bond-based coarse-graining of RNA backbone conformation and generates RNA conformations through Monte Carlo sampling of the bond angles and torsional angles of the virtual bonds. Using a coarse-grained statistical potential derived from the known structures, we assign each conformation with a statistical weight. The weighted average over the conformational ensemble gives the entropy and free energy parameters for the hairpin, bulge, and internal loops, and multiway junctions. From the thermodynamic parameters, we predict RNA structures, melting curves, and structural changes from the sequence. Theory-experiment comparisons indicate that Vfold2D-MC not only gives improved structure predictions but also enables the interpretation of thermodynamic results for different RNA structures, including multibranched junctions. This new model sets a promising framework to treat more complicated RNA structures, such as pseudoknotted and intramolecular kissing loops, for which experimental thermodynamic parameters are often unavailable.

Diagrammatic approaches to RNA structures with trinucleotide repeats

Trinucleotide repeat expansion disorders are associated with the overexpansion of (CNG) repeats on the genome. Messenger RNA transcripts of sequences with greater than 60–100 (CNG) tandem units have been implicated in trinucleotide repeat expansion disorder pathogenesis. In this work, we develop a diagrammatic theory to study the structural diversity of these (CNG)_n RNA sequences. Representing structural elements on the chain’s conformation by a set of graphs and employing elementary diagrammatic methods, we have formulated a renormalization procedure to re-sum these graphs and arrive at a closed-form expression for the ensemble partition function. With a simple approximation for the renormalization and applied to extended (CNG)_n sequences, this theory can comprehensively capture an infinite set of conformations with any number and any combination of duplexes, hairpins, multiway junctions, and quadruplexes. To quantify the diversity of different (CNG)_n ensembles, the analytical equations derived from the diagrammatic theory were solved numerically to derive equilibrium estimates for the secondary structural contents of the chains. The results suggest that the structural ensembles of (CNG)_n repeat sequence with n ∼60 are surprisingly diverse, and the distribution is sensitive to the ability of the N nucleotide to make noncanonical pairs and whether the (CNG)_n sequence can sustain stable quadruplexes. The results show how perturbations in the form of biases on the stabilities of the various structural motifs, duplexes, junctions, helices, and quadruplexes could affect the secondary structures of the chains and how these structures may switch when they are perturbed.

Theoretical Model for Solvent-Induced Base Stacking Interactions in Solvent-Free DNA Simulations

Ultrahigh-throughput conformational sampling of biopolymers like nucleic acids are most effectively carried out without explicit solvents, but the physical origins of almost all inter- and intramolecular interactions controlling nucleic acid structures are rooted in water. Single-stranded (ss) DNAs or RNAs in water are characterized by ensembles of diverse conformations. To properly capture solvent-induced nucleobase stacking interactions in an otherwise solvent-free Monte Carlo algorithm, theoretical models are developed here to describe the solvent entropy and dispersion terms in base stacking free energies. To validate these models, equilibrium ensembles of ss (dA) _n and (dT) _n sequences ( n = 30, 40, and 50) were simulated, and they quantitatively reproduced experimental small-angle X-ray scattering (SAXS) data. Simulated dA ensembles show substantial stacking. While less prevalent, stacking in dT chains is not negligible. Analysis of SAXS profiles suggests that excess features between wavevector 0.03 and 0.18 Å^-1 correlate with stacking, and stacking in dA versus dT chains is chain length-dependent, where (dT)₃₀ and (dA)₃₀ chains have more similar structures, but longer dA chains show more stacking over dT. The average stack length in ss-dA chains is 5-10 nucleotides, yielding an estimate for the overall A|A stacking free energy at ∼1 kcal/mol.

DNA-Based Nanodevices Controlled by Purely Entropic Linker Domains

We demonstrate here the rational design of purely entropic domains as a versatile approach to achieve control of the input/output response of synthetic molecular receptors. To do so and to highlight the versatility and generality of this approach, we have rationally re-engineered two model DNA-based receptors: a clamp-like DNA-based switch that recognizes a specific DNA sequence and an ATP-binding aptamer. We show that, by varying the length of the linker domain that connects the two recognition portions of these receptors, it is possible to finely control their affinity for their specific ligand. Through mathematical modeling and thermodynamic characterization, we also demonstrate for both systems that entropy changes associated with changes in linker length are responsible for affinity modulation and that the linker we have designed behaves as a disordered random-coil polymer. The approach also allows us to regulate the ligand concentration range at which the receptors respond and show optimal specificity. Given these attributes, the use of purely entropic domains appears as a versatile and general approach to finely control the activity of synthetic receptors in a highly predictable and controlled fashion.