Conformational Dynamics of the Hepatitis C Virus 3'X RNA.

The 3' end of the hepatitis C virus genome is terminated by a highly conserved, 98-nucleotide sequence called 3'X. This untranslated structural element is thought to regulate several essential RNA-dependent processes associated with infection. 3'X has two proposed conformations comprised of either three- or two stem-loop structures that result from different base pairing interactions within the first 55 nucleotides. Here, we used single-molecule FRET spectroscopy to monitor the conformational status of fluorescently labeled constructs that isolate this region of the RNA (3'X55). We observed that 3'X55 can adopt both proposed conformations and the relative abundance of them can be modulated by either solution conditions or nucleotide deletions. Furthermore, interconversion between the two conformations is slow and takes place over the course of several hours. The simultaneous existence of two slowly interconverting conformations may help prime individual copies of the viral genome for either viral protein or RNA synthesis, thereby minimizing conflicts between these two competing processes.

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.

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).

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.

The SARS-CoV-2 Conserved Macrodomain Is a Mono-ADP-Ribosylhydrolase.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and other SARS-related CoVs encode 3 tandem macrodomains within nonstructural 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 antiviral 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 Middle East respiratory syndrome coronavirus (MERS-CoV) Mac1 domains exhibit similar structural folds, and all 3 proteins bound to ADP-ribose with affinities in the low micromolar range. Importantly, using ADP-ribose-detecting binding reagents in both a gel-based assay and novel enzyme-linked immunosorbent assays (ELISAs), 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 than 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 more than 1.2 million deaths worldwide. With no currently approved treatments, novel therapeutic strategies are desperately needed. All coronaviruses encode a highly conserved macrodomain (Mac1) that binds to and removes ADP-ribose adducts from proteins in a dynamic posttranslational process that is increasingly being 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 those of 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.

Real-time compaction of nanoconfined DNA by an intrinsically disordered macromolecular counterion.

We demonstrate how a recently developed nanofluidic device can be used to study protein-induced compaction of genome-length DNA freely suspended in solution. The protein we use in this study is the hepatitis C virus core protein (HCVcp), which is a positively charged, intrinsically disordered protein. Using nanofluidic devices in combination with fluorescence microscopy, we observe that protein-induced compaction preferentially begins at the ends of linear DNA. This observation would be difficult to make with many other single-molecule techniques, which generally require the DNA ends to be anchored to a substrate. We also demonstrate that this protein-induced compaction is reversible and can be dynamically modulated by exposing the confined DNA molecules to solutions containing either HCVcp (to promote compaction) or Proteinase K (to disassemble the compact nucleo-protein complex). Although the natural binding partner for HCVcp is genomic viral RNA, the general biophysical principles governing protein-induced compaction of DNA are likely relevant for a broad range of nucleic acid-binding proteins and their targets.

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.

Disordered RNA chaperones can enhance nucleic acid folding via local charge screening.

RNA chaperones are proteins that aid in the folding of nucleic acids, but remarkably, many of these proteins are intrinsically disordered. How can these proteins function without a well-defined three-dimensional structure? Here, we address this question by studying the hepatitis C virus core protein, a chaperone that promotes viral genome dimerization. Using single-molecule fluorescence spectroscopy, we find that this positively charged disordered protein facilitates the formation of compact nucleic acid conformations by acting as a flexible macromolecular counterion that locally screens repulsive electrostatic interactions with an efficiency equivalent to molar salt concentrations. The resulting compaction can bias unfolded nucleic acids towards folding, resulting in faster folding kinetics. This potentially widespread mechanism is supported by molecular simulations that rationalize the experimental findings by describing the chaperone as an unstructured polyelectrolyte.

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 dT19, dA19, and rA19 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.

Tests of Kramers' Theory at the Single-Molecule Level: Evidence for Folding of an Isolated RNA Tertiary Interaction at the Viscous Speed Limit.

Dissipation and friction influence the conformational dynamics of biological polymers as they traverse barriers on rugged free energy surfaces. It is well established that the "speed limit" for macromolecular folding is dictated by a combination of (i) solvent friction, which depends on solvent viscosity, η, and (ii) internal friction, which is independent of solvent and depends solely on the molecular folding pathway. In this work, single-molecule Förster resonance energy transfer (FRET) confocal spectroscopy is used to study viscosity-dependent folding kinetics of an isolated RNA tertiary motif, that of the GAAA tetraloop receptor, allowing both solvent and internal frictional contributions to be investigated and extracted independently for both flexible PEG- and RNA-based (rU7, rA7) linkers in the unimolecular construct. Specifically, our single-molecule data reveal that (i) folding rate constants scale linearly with the inverse solvent viscosity (η), which supports Kramers'/Grote-Hynes' rate theory for η-dependent RNA folding and that (ii) they provide quantitative upper limits for the intrinsic viscosity, [ηint ≈ 0.1(2) cP], arising from internal friction associated with folding/unfolding of an isolated RNA tertiary interaction. Furthermore, in contrast to strong viscosity-induced shifts in the folding/unfolding rate constants, temperature-dependent studies demonstrate that the enthalpic, entropic, and free energy contributions to the transition state barrier are largely insensitive to the solvent viscosity. This supports a very simple picture for the conformational kinetics of isolated RNA tertiary interactions wherein rate constants for folding/unfolding are both inversely dependent on viscosity and limited by diffusional access to the transition state region on a multidimensional free energy surface. Particularly under cellular conditions, where ηsolv > 1 cp, this suggests that RNAs fold/unfold at a "speed limit" dictated by solvent viscosity and transition-state barrier thermodynamics rather than internal molecular friction.

Phase-transition assisted mechanical behavior of TiZrHfTax high-entropy alloys.

Recent developments of high-entropy alloys with high strength and high ductility draw attention to the metastability-engineering strategy. Using first-principle theory, here we demonstrate that reducing the Ta level in the refractory TiZrHfTax system destabilizes the body-centered cubic (bcc) phase and leads to the appearance of the hexagonal close-packed (hcp) phase embedded in the bcc matrix. The alloying-induced features of the elastic parameters for the cubic and hexagonal structures are mapped out in details, and strong sensitivity to the crystal lattice and chemistry is revealed. Results show softening of the bcc matrix with decreasing Ta concentration which ensures ductile behavior. However, the elastically nearly isotropic hcp precipitates possess enhanced resistance against shear which promotes strengthening of the TiZrHfTax dual-phase system. The present atomic-level insight provides strong evidence to the experimental observation, and emphasizes the significance of quantum-design for advanced multi-phase high-entropy alloys with excellent strength-ductility combinations.

Twinning in metastable high-entropy alloys.

Twinning is a fundamental mechanism behind the simultaneous increase of strength and ductility in medium- and high-entropy alloys, but its operation is not yet well understood, which limits their exploitation. Since many high-entropy alloys showing outstanding mechanical properties are actually thermodynamically unstable at ambient and cryogenic conditions, the observed twinning challenges the existing phenomenological and theoretical plasticity models. Here, we adopt a transparent approach based on effective energy barriers in combination with first-principle calculations to shed light on the origin of twinning in high-entropy alloys. We demonstrate that twinning can be the primary deformation mode in metastable face-centered cubic alloys with a fraction that surpasses the previously established upper limit. The present advance in plasticity of metals opens opportunities for tailoring the mechanical response in engineering materials by optimizing metastable twinning in high-entropy alloys.

Incorporation of isotopic, fluorescent, and heavy-atom-modified nucleotides into RNAs by position-selective labeling of RNA.

Site-specific incorporation of labeled nucleotides is an extremely useful synthetic tool for many structural studies (e.g., NMR, electron paramagnetic resonance (EPR), fluorescence resonance energy transfer (FRET), and X-ray crystallography) of RNA. However, specific-position-labeled RNAs >60 nt are not commercially available on a milligram scale. Position-selective labeling of RNA (PLOR) has been applied to prepare large RNAs labeled at desired positions, and all the required reagents are commercially available. Here, we present a step-by-step protocol for the solid-liquid hybrid phase method PLOR to synthesize 71-nt RNA samples with three different modification applications, containing (i) a 13C15N-labeled segment; (ii) discrete residues modified with Cy3, Cy5, or biotin; or (iii) two iodo-U residues. The flexible procedure enables a wide range of downstream biophysical analyses using precisely localized functionalized nucleotides. All three RNAs were obtained in <2 d, excluding time for preparing reagents and optimizing experimental conditions. With optimization, the protocol can be applied to other RNAs with various labeling schemes, such as ligation of segmentally labeled fragments.

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.

First Principles Theory of the hcp-fcc Phase Transition in Cobalt.

Identifying the forces that drive a phase transition is always challenging. The hcp-fcc phase transition that occurs in cobalt at ~700 K has not yet been fully understood, although early theoretical studies have suggested that magnetism plays a main role in the stabilization of the fcc phase at high temperatures. Here, we perform a first principles study of the free energies of these two phases, which we break into contributions arising from the vibration of the lattice, electronic and magnetic systems and volume expansion. Our analysis of the energy of the phases shows that magnetic effects alone cannot drive the fcc-hcp transition in Co and that the largest contribution to the stabilization of the fcc phase comes from the vibration of the ionic lattice. By including all the contributions to the free energy considered here we obtain a theoretical transition temperature of 825 K.