Mechanism of Ligand Binding to Theophylline RNA Aptamer

Studying RNA-ligand interactions and quantifying their binding thermodynamics and kinetics are of particular relevance in the field of drug discovery. Here, we combined biochemical binding assays and accelerated molecular simulations to investigate ligand binding and dissociation in RNA using the theophylline-binding RNA as a model system. All-atom simulations using a Ligand Gaussian accelerated Molecular Dynamics method (LiGaMD) have captured repetitive binding and dissociation of theophylline and caffeine to RNA. Theophylline's binding free energy and kinetic rate constants align with our experimental data, while caffeine's binding affinity is over 10,000 times weaker, and its kinetics could not be determined. LiGaMD simulations allowed us to identify distinct low-energy conformations and multiple ligand binding pathways to RNA. Simulations revealed a "conformational selection" mechanism for ligand binding to the flexible RNA aptamer, which provides important mechanistic insights into ligand binding to the theophylline-binding model. Our findings suggest that compound docking using a structural ensemble of representative RNA conformations would be necessary for structure-based drug design of flexible RNA.

Single-molecule-binding studies of antivirals targeting the hepatitis C virus core protein

IMPORTANCE

The hepatitis C virus is associated with nearly 300,000 deaths annually. At the core of the virus is an RNA-protein complex called the nucleocapsid, which consists of the viral genome and many copies of the core protein. Because the assembly of the nucleocapsid is a critical step in viral replication, a considerable amount of effort has been devoted to identifying antiviral therapeutics that can bind to the core protein and disrupt assembly. Although several candidates have been identified, little is known about how they interact with the core protein or how those interactions alter the structure and thus the function of this viral protein. Our work biochemically characterizes several of these binding interactions, highlighting both similarities and differences as well as strengths and weaknesses. These insights bolster the notion that this viral protein is a viable target for novel therapeutics and will help to guide future developments of these candidate antivirals.

Apolar chemical environments compact unfolded RNAs and can promote folding

It is well documented that the structure, and thus function, of nucleic acids depends on the chemical environment surrounding them, which often includes potential proteinaceous binding partners. The nonpolar amino acid side chains of these proteins will invariably alter the polarity of the local chemical environment around the nucleic acid. However, we are only beginning to understand how environmental polarity generally influences the structural and energetic properties of RNA folding. Here, we use a series of aqueous-organic cosolvent mixtures to systematically modulate the solvent polarity around two different RNA folding constructs that can form either secondary or tertiary structural elements. Using single-molecule Förster resonance energy transfer spectroscopy to simultaneously monitor the structural and energetic properties of these RNAs, we show that the unfolded conformations of both model RNAs become more compact in apolar environments characterized by dielectric constants less than that of pure water. In the case of tertiary structure formation, this compaction also gives rise to more energetically favorable folding. We propose that these physical changes arise from an enhanced accumulation of counterions in the low dielectric environment surrounding the unfolded RNA.

Recognition of single-stranded nucleic acids by small-molecule splicing modulators

Risdiplam is the first approved small-molecule splicing modulator for the treatment of spinal muscular atrophy (SMA). Previous studies demonstrated that risdiplam analogues have two separate binding sites in exon 7 of the SMN2 pre-mRNA: (i) the 5'-splice site and (ii) an upstream purine (GA)-rich binding site. Importantly, the sequence of this GA-rich binding site significantly enhanced the potency of risdiplam analogues. In this report, we unambiguously determined that a known risdiplam analogue, SMN-C2, binds to single-stranded GA-rich RNA in a sequence-specific manner. The minimum required binding sequence for SMN-C2 was identified as GAAGGAAGG. We performed all-atom simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method, which captured spontaneous binding of a risdiplam analogue to the target nucleic acids. We uncovered, for the first time, a ligand-binding pocket formed by two sequential GAAG loop-like structures. The simulation findings were highly consistent with experimental data obtained from saturation transfer difference (STD) NMR and structure-affinity-relationship studies of the risdiplam analogues. Together, these studies illuminate us to understand the molecular basis of single-stranded purine-rich RNA recognition by small-molecule splicing modulators with an unprecedented binding mode.

Integrating single-molecule FRET and biomolecular simulations to study diverse interactions between nucleic acids and proteins

The conformations of biological macromolecules are intimately related to their cellular functions. Conveniently, the well-characterized dipole-dipole distance-dependence of Förster resonance energy transfer (FRET) makes it possible to measure and monitor the nanoscale spatial dimensions of these conformations using fluorescence spectroscopy. For this reason, FRET is often used in conjunction with single-molecule detection to study a wide range of conformationally dynamic biochemical processes. Written for those not yet familiar with the subject, this review aims to introduce biochemists to the methodology associated with single-molecule FRET, with a particular emphasis on how it can be combined with biomolecular simulations to study diverse interactions between nucleic acids and proteins. In the first section, we highlight several conceptual and practical considerations related to this integrative approach. In the second section, we review a few recent research efforts wherein various combinations of single-molecule FRET and biomolecular simulations were used to study the structural and dynamic properties of biochemical systems involving different types of nucleic acids (e.g., DNA and RNA) and proteins (e.g., folded and disordered).

The SARS-CoV-2 conserved macrodomain is a mono-ADP-ribosylhydrolase

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and other SARS-like-CoVs encode 3 tandem macrodomains within non-structural protein 3 (nsp3). The first macrodomain, Mac1, is conserved throughout CoVs, and binds to and hydrolyzes mono-ADP-ribose (MAR) from target proteins. Mac1 likely counters host-mediated anti-viral ADP-ribosylation, a posttranslational modification that is part of the host response to viral infections. Mac1 is essential for pathogenesis in multiple animal models of CoV infection, implicating it as a virulence factor and potential therapeutic target. Here we report the crystal structure of SARS-CoV-2 Mac1 in complex with ADP-ribose. SARS-CoV-2, SARS-CoV and MERS-CoV Mac1 exhibit similar structural folds and all 3 proteins bound to ADP-ribose with low μM affinities. Importantly, using ADP-ribose detecting binding reagents in both a gel-based assay and novel ELISA assays, we demonstrated de-MARylating activity for all 3 CoV Mac1 proteins, with the SARS-CoV-2 Mac1 protein leading to a more rapid loss of substrate compared to the others. In addition, none of these enzymes could hydrolyze poly-ADP-ribose. We conclude that the SARS-CoV-2 and other CoV Mac1 proteins are MAR-hydrolases with similar functions, indicating that compounds targeting CoV Mac1 proteins may have broad anti-CoV activity.

IMPORTANCE

SARS-CoV-2 has recently emerged into the human population and has led to a worldwide pandemic of COVID-19 that has caused greater than 900 thousand deaths worldwide. With, no currently approved treatments, novel therapeutic strategies are desperately needed. All coronaviruses encode for a highly conserved macrodomain (Mac1) that binds to and removes ADP-ribose adducts from proteins in a dynamic post-translational process increasingly recognized as an important factor that regulates viral infection. The macrodomain is essential for CoV pathogenesis and may be a novel therapeutic target. Thus, understanding its biochemistry and enzyme activity are critical first steps for these efforts. Here we report the crystal structure of SARS-CoV-2 Mac1 in complex with ADP-ribose, and describe its ADP-ribose binding and hydrolysis activities in direct comparison to SARS-CoV and MERS-CoV Mac1 proteins. These results are an important first step for the design and testing of potential therapies targeting this unique protein domain.

Binding without folding - the biomolecular function of disordered polyelectrolyte complexes

Recent evidence shows that oppositely charged intrinsically disordered proteins (IDPs) can form high-affinity complexes that involve neither the formation of secondary or tertiary structure nor site-specific interactions between individual residues. Similar electrostatically dominated interactions have also been identified for positively charged IDPs binding to nucleic acids. These highly disordered polyelectrolyte complexes constitute an extreme case within the spectrum of biomolecular interactions involving disorder. Such interactions are likely to be widespread, since sequence analysis predicts proteins with highly charged disordered regions to be surprisingly numerous. Here, we summarize the insights that have emerged from the highly disordered polyelectrolyte complexes identified so far, and we highlight recent developments and future challenges in (i) establishing models for the underlying highly dynamic structural ensembles, (ii) understanding the novel binding mechanisms associated with them, and (iii) identifying the functional consequences.

Accurate Transfer Efficiencies, Distance Distributions, and Ensembles of Unfolded and Intrinsically Disordered Proteins From Single-Molecule FRET

Intrinsically disordered proteins (IDPs) sample structurally diverse ensembles. Characterizing the underlying distributions of conformations is a key step toward understanding the structural and functional properties of IDPs. One increasingly popular method for obtaining quantitative information on intramolecular distances and distributions is single-molecule Förster resonance energy transfer (FRET). Here we describe two essential elements of the quantitative analysis of single-molecule FRET data of IDPs: the sample-specific calibration of the single-molecule instrument that is required for determining accurate transfer efficiencies, and the use of state-of-the-art methods for inferring accurate distance distributions from these transfer efficiencies. First, we illustrate how to quantify the correction factors for instrument calibration with alternating donor and acceptor excitation measurements of labeled samples spanning a wide range of transfer efficiencies. Second, we show how to infer distance distributions based on suitably parameterized simple polymer models, and how to obtain conformational ensembles from Bayesian reweighting of molecular simulations or from parameter optimization in simplified coarse-grained models.

Transition path times of coupled folding and binding reveal the formation of an encounter complex

The association of biomolecules is the elementary event of communication in biology. Most mechanistic information of how the interactions between binding partners form or break is, however, hidden in the transition paths, the very short parts of the molecular trajectories from the encounter of the two molecules to the formation of a stable complex. Here we use single-molecule spectroscopy to measure the transition path times for the association of two intrinsically disordered proteins that form a folded dimer upon binding. The results reveal the formation of a metastable encounter complex that is electrostatically favored and transits to the final bound state within tens of microseconds. Such measurements thus open a new window into the microscopic events governing biomolecular interactions.

Dispersion Correction Alleviates Dye Stacking of Single-Stranded DNA and RNA in Simulations of Single-Molecule Fluorescence Experiments

We combine single-molecule Förster resonance energy transfer (single-molecule FRET) experiments with extensive all-atom molecular dynamics (MD) simulations (>100 μs) to characterize the conformational ensembles of single-stranded (ss) DNA and RNA in solution. From MD simulations with explicit dyes attached to single-stranded nucleic acids via flexible linkers, we calculate FRET efficiencies and fluorescence anisotropy decays. We find that dispersion-corrected water models alleviate the problem of overly abundant interactions between fluorescent dyes and the aromatic ring systems of nucleobases. To model dye motions in a computationally efficient and conformationally exhaustive manner, we introduce a dye-conformer library, built from simulations of dinucleotides with covalently attached dye molecules. We use this library to calculate FRET efficiencies for dT₁₉, dA₁₉, and rA₁₉ simulated without explicit labels over a wide range of salt concentrations. For end-labeled homopolymeric pyrimidine ssDNA, MD simulations with the parmBSC1 force field capture the overall trend in salt-dependence of single-molecule FRET based distance measurements. For homopolymeric purine ssRNA and ssDNA, the DESRES and parmBSC1 force fields, respectively, provide useful starting points, even though our comparison also identifies clear deviations from experiment.

Biophysical Insights from Temperature-Dependent Single-Molecule Förster Resonance Energy Transfer

Single-molecule fluorescence microscopy techniques can be used in combination with micrometer length-scale temperature control and Förster resonance energy transfer (FRET) in order to gain detailed information about fundamental biophysical phenomena. In particular, this combination of techniques has helped foster the development of remarkable quantitative tools for studying both time- and temperature-dependent structural kinetics of biopolymers. Over the past decade, multiple research efforts have successfully incorporated precise spatial and temporal control of temperature into single-molecule FRET (smFRET)-based experiments, which have uncovered critical thermodynamic information on a wide range of biological systems such as conformational dynamics of nucleic acids. This review provides an overview of various temperature-dependent smFRET approaches from our laboratory and others, highlighting efforts in which such methods have been successfully applied to studies of single-molecule nucleic acid folding.

Conformational Plasticity of Hepatitis C Virus Core Protein Enables RNA-Induced Formation of Nucleocapsid-like Particles

Many of the unanswered questions associated with hepatitis C virus assembly are related to the core protein (HCVcp), which forms an oligomeric nucleocapsid encompassing the viral genome. The structural properties of HCVcp have been difficult to quantify, at least in part because it is an intrinsically disordered protein. We have used single-molecule Förster Resonance Energy Transfer techniques to study the conformational dimensions and dynamics of the HCVcp nucleocapsid domain (HCVncd) at various stages during the RNA-induced formation of nucleocapsid-like particles. Our results indicate that HCVncd is a typical intrinsically disordered protein. When it forms small ribonucleoprotein complexes with various RNA hairpins from the 3' end of the HCV genome, it compacts but remains intrinsically disordered and conformationally dynamic. Above a critical RNA concentration, these ribonucleoprotein complexes rapidly and cooperatively assemble into large nucleocapsid-like particles, wherein the individual HCVncd subunits become substantially more extended.

Mechanistic Insights into Cofactor-Dependent Coupling of RNA Folding and mRNA Transcription/Translation by a Cobalamin Riboswitch

Riboswitches are mRNA elements regulating gene expression in response to direct binding of a metabolite. While these RNAs are increasingly well understood with respect to interactions between receptor domains and their cognate effector molecules, little is known about the specific mechanistic relationship between metabolite binding and gene regulation by the downstream regulatory domain. Using a combination of cell-based, biochemical, and biophysical techniques, we reveal the specific RNA architectural features enabling a cobalamin-dependent hairpin loop docking interaction between receptor and regulatory domains. Furthermore, these data demonstrate that docking kinetics dictate a regulatory response involving the coupling of translation initiation to general mechanisms that control mRNA abundance. These results yield a comprehensive picture of how RNA structure in the riboswitch regulatory domain enables kinetically constrained ligand-dependent regulation of gene expression.

Single-Molecule FRET Reveals Three Conformations for the TLS Domain of Brome Mosaic Virus Genome

Metabolite-dependent conformational switching in RNA riboswitches is now widely accepted as a critical regulatory mechanism for gene expression in bacterial systems. More recently, similar gene regulation mechanisms have been found to be important for viral systems as well. One of the most abundant and best-studied systems is the tRNA-like structure (TLS) domain, which has been found to occur in many plant viruses spread across numerous genera. In this work, folding dynamics for the TLS domain of Brome Mosaic Virus have been investigated using single-molecule fluorescence resonance energy transfer techniques. In particular, burst fluorescence methods are exploited to observe metal-ion ([M(n+)])-induced folding in freely diffusing RNA constructs resembling the minimal TLS element of brome mosaic virus RNA3. The results of these experiments reveal a complex equilibrium of at least three distinct populations. A stepwise, or consecutive, thermodynamic model for TLS folding is developed, which is in good agreement with the [M(n+)]-dependent evolution of conformational populations and existing structural information in the literature. Specifically, this folding pathway explains the metal-ion dependent formation of a functional TLS domain from unfolded RNAs via two consecutive steps: 1) hybridization of a long-range stem interaction, followed by 2) formation of a 3'-terminal pseudoknot. These two conformational transitions are well described by stepwise dissociation constants for [Mg(2+)] (K1 = 328 ± 30 μM and K2 = 1092 ± 183 μM) and [Na(+)] (K1 = 74 ± 6 mM and K2 = 243 ± 52 mM)-induced folding. The proposed thermodynamic model is further supported by inhibition studies of the long-range stem interaction using a complementary DNA oligomer, which effectively shifts the dynamic equilibrium toward the unfolded conformation. Implications of this multistep conformational folding mechanism are discussed with regard to regulation of virus replication.

Kinetic and thermodynamic origins of osmolyte-influenced nucleic acid folding

The influential role of monovalent and divalent metal cations in facilitating conformational transitions in both RNA and DNA has been a target of intense biophysical research efforts. However, organic neutrally charged cosolutes can also significantly alter nucleic acid conformational transitions. For example, highly soluble small molecules such as trimethylamine N-oxide (TMAO) and urea are occasionally utilized by organisms to regulate cellular osmotic pressure. Ensemble studies have revealed that these so-called osmolytes can substantially influence the thermodynamics of nucleic acid conformational transitions. In the present work, we exploit single-molecule FRET (smFRET) techniques to measure, for first time, the kinetic origins of these osmolyte-induced changes to the folding free energy. In particular, we focus on smFRET RNA and DNA constructs designed as model systems for secondary and tertiary structure formation. These findings reveal that TMAO preferentially stabilizes both secondary and tertiary interactions by increasing kfold and decreasing kunfold, whereas urea destabilizes both conformational transitions, resulting in the exact opposite shift in kinetic rate constants (i.e., decreasing kfold and increasing kunfold). Complementary temperature-dependent smFRET experiments highlight a thermodynamic distinction between the two different mechanisms responsible for TMAO-facilitated conformational transitions, while only a single mechanism is seen for the destabilizing osmolyte urea. Finally, these results are interpreted in the context of preferential interactions between osmolytes, and the solvent accessible surface area (SASA) associated with the (i) nucleobase, (ii) sugar, and (iii) phosphate groups of nucleic acids in order to map out structural changes that occur during the conformational transitions.

Single-molecule conformational dynamics of a biologically functional hydroxocobalamin riboswitch

Riboswitches represent a family of highly structured regulatory elements found primarily in the leader sequences of bacterial mRNAs. They function as molecular switches capable of altering gene expression; commonly, this occurs via a conformational change in a regulatory element of a riboswitch that results from ligand binding in the aptamer domain. Numerous studies have investigated the ligand binding process, but little is known about the structural changes in the regulatory element. A mechanistic description of both processes is essential for deeply understanding how riboswitches modulate gene expression. This task is greatly facilitated by studying all aspects of riboswitch structure/dynamics/function in the same model system. To this end, single-molecule fluorescence resonance energy transfer (smFRET) techniques have been used to directly observe the conformational dynamics of a hydroxocobalamin (HyCbl) binding riboswitch (env8HyCbl) with a known crystallographic structure.1 The single-molecule RNA construct studied in this work is unique in that it contains all of the structural elements both necessary and sufficient for regulation of gene expression in a biological context. The results of this investigation reveal that the undocking rate constant associated with the disruption of a long-range kissing-loop (KL) interaction is substantially decreased when the ligand is bound to the RNA, resulting in a preferential stabilization of the docked conformation. Notably, the formation of this tertiary KL interaction directly sequesters the Shine-Dalgarno sequence (i.e., the ribosome binding site) via base-pairing, thus preventing translation initiation. These results reveal that the conformational dynamics of this regulatory switch are quantitatively described by a four-state kinetic model, whereby ligand binding promotes formation of the KL interaction. The results of complementary cell-based gene expression experiments conducted in Escherichia coli are highly correlated with the smFRET results, suggesting that KL formation is directly responsible for regulating gene expression.

Single-molecule fluorescence resonance energy transfer studies of the human telomerase RNA pseudoknot: temperature-/urea-dependent folding kinetics and thermodynamics

The ribonucleoprotein telomerase is an RNA-dependent DNA polymerase that catalyzes the repetitive addition of a short, species-specific, DNA sequence to the ends of linear eukaryotic chromosomes. The single RNA component of telomerase contains both the template sequence for DNA synthesis and a functionally critical pseudoknot motif, which can also exist as a less stable hairpin. Here we use a minimal version of the human telomerase RNA pseudoknot to study this hairpin-pseudoknot structural equilibrium using temperature-controlled single-molecule fluorescence resonance energy transfer (smFRET) experiments. The urea dependence of these experiments aids in determination of the folding kinetics and thermodynamics. The wild-type pseudoknot behavior is compared and contrasted to a mutant pseudoknot sequence implicated in a genetic disorder-dyskeratosis congenita. These findings clearly identify that this 2nt noncomplementary mutation destabilizes the folding of the wild-type pseudoknot by substantially reducing the folding rate constant (≈ 400-fold) while only nominally increasing the unfolding rate constant (≈ 5-fold). Furthermore, the urea dependence of the equilibrium and rate constants is used to develop a free energy landscape for this unimolecular equilibrium and propose details about the structure of the transition state. Finally, the urea-dependent folding experiments provide valuable physical insights into the mechanism for destabilization of RNA pseudoknots by such chemical denaturants.

Single-molecule kinetics reveal cation-promoted DNA duplex formation through ordering of single-stranded helices

In this work, the kinetics of short, fully complementary oligonucleotides are investigated at the single-molecule level. Constructs 6-9 bp in length exhibit single exponential kinetics over 2 orders of magnitude time for both forward (kon, association) and reverse (koff, dissociation) processes. Bimolecular rate constants for association are weakly sensitive to the number of basepairs in the duplex, with a 2.5-fold increase between 9 bp (k'on = 2.1(1) × 10(6) M(-1) s(-1)) and 6 bp (k'on = 5.0(1) × 10(6) M(-1) s(-1)) sequences. In sharp contrast, however, dissociation rate constants prove to be exponentially sensitive to sequence length, varying by nearly 600-fold over the same 9 bp (koff = 0.024 s(-1)) to 6 bp (koff = 14 s(-1)) range. The 8 bp sequence is explored in more detail, and the NaCl dependence of kon and koff is measured. Interestingly, kon increases by >40-fold (kon = 0.10(1) s(-1) to 4.0(4) s(-1) between [NaCl] = 25 mM and 1 M), whereas in contrast, koff decreases by fourfold (0.72(3) s(-1) to 0.17(7) s(-1)) over the same range of conditions. Thus, the equilibrium constant (Keq) increases by ≈160, largely due to changes in the association rate, kon. Finally, temperature-dependent measurements reveal that increased [NaCl] reduces the overall exothermicity (ΔΔH° > 0) of duplex formation, albeit by an amount smaller than the reduction in entropic penalty (-TΔΔS° < 0). This reduced entropic cost is attributed to a cation-facilitated preordering of the two single-stranded species, which lowers the association free-energy barrier and in turn accelerates the rate of duplex formation.

Thermodynamic origins of monovalent facilitated RNA folding

Cations have long been associated with formation of native RNA structure and are commonly thought to stabilize the formation of tertiary contacts by favorably interacting with the electrostatic potential of the RNA, giving rise to an "ion atmosphere". A significant amount of information regarding the thermodynamics of structural transitions in the presence of an ion atmosphere has accumulated and suggests stabilization is dominated by entropic terms. This work provides an analysis of how RNA-cation interactions affect the entropy and enthalpy associated with an RNA tertiary transition. Specifically, temperature-dependent single-molecule fluorescence resonance energy transfer studies have been exploited to determine the free energy (ΔG°), enthalpy (ΔH°), and entropy (ΔS°) of folding for an isolated tetraloop-receptor tertiary interaction as a function of Na(+) concentration. Somewhat unexpectedly, increasing the Na(+) concentration changes the folding enthalpy from a strongly exothermic process [e.g., ΔH° = -26(2) kcal/mol at 180 mM] to a weakly exothermic process [e.g., ΔH° = -4(1) kcal/mol at 630 mM]. As a direct corollary, it is the strong increase in folding entropy [Δ(ΔS°) > 0] that compensates for this loss of exothermicity for the achievement of more favorable folding [Δ(ΔG°) < 0] at higher Na(+) concentrations. In conjunction with corresponding measurements of the thermodynamics of the transition state barrier, these data provide a detailed description of the folding pathway associated with the GAAA tetraloop-receptor interaction as a function of Na(+) concentration. The results support a potentially universal mechanism for monovalent facilitated RNA folding, whereby an increasing monovalent concentration stabilizes tertiary structure by reducing the entropic penalty for folding.

Real-Time Infrared Overtone Laser Control of Temperature in Picoliter H(2)O Samples: "Nanobathtubs" for Single Molecule Microscopy

An approach for high spatiotemporal control of aqueous sample temperatures in confocal microscopy is reported. This technique exploits near-IR diode-laser illumination to locally heat picoliter volumes of water via first-overtone excitation in the OH-stretch manifold. A thin water cell after the objective resonantly removes any residual IR light from the detection system, allowing for continuous observation of single-molecule fluorescence throughout the heating event. This technique is tested quantitatively by reproducing single-molecule RNA folding results obtained from "bulk" stage heating measurements. Calibration of sample temperatures is obtained from time-correlated single-photon counting studies of Rhodamine B fluorescence decay. We obtain an upper limit to the heating response time (τ(heat) < 20 ms) consistent with even faster estimates (τ(heat) ≈ 0.25 ms) based on laser spot size, H(2)O heat capacit,y and absorption cross section. This combination of fast, noncontact heating of picoliter volumes provides new opportunities for real-time thermodynamic/kinetic studies at the single-molecule level.