The disordered N-terminal tail of SARS-CoV-2 Nucleocapsid protein forms a dynamic complex with RNA.

The SARS-CoV-2 Nucleocapsid (N) protein is responsible for condensation of the viral genome. Characterizing the mechanisms controlling nucleic acid binding is a key step in understanding how condensation is realized. Here, we focus on the role of the RNA binding domain (RBD) and its flanking disordered N-terminal domain (NTD) tail, using single-molecule Förster Resonance Energy Transfer and coarse-grained simulations. We quantified contact site size and binding affinity for nucleic acids and concomitant conformational changes occurring in the disordered region. We found that the disordered NTD increases the affinity of the RBD for RNA by about 50-fold. Binding of both nonspecific and specific RNA results in a modulation of the tail configurations, which respond in an RNA length-dependent manner. Not only does the disordered NTD increase affinity for RNA, but mutations that occur in the Omicron variant modulate the interactions, indicating a functional role of the disordered tail. Finally, we found that the NTD-RBD preferentially interacts with single-stranded RNA and that the resulting protein:RNA complexes are flexible and dynamic. We speculate that this mechanism of interaction enables the Nucleocapsid protein to search the viral genome for and bind to high-affinity motifs.

Reweighting of molecular simulations with explicit-solvent SAXS restraints elucidates ion-dependent RNA ensembles.

Small-angle X-ray scattering (SAXS) experiments are increasingly used to probe RNA structure. A number of forward models that relate measured SAXS intensities and structural features, and that are suitable to model either explicit-solvent effects or solute dynamics, have been proposed in the past years. Here, we introduce an approach that integrates atomistic molecular dynamics simulations and SAXS experiments to reconstruct RNA structural ensembles while simultaneously accounting for both RNA conformational dynamics and explicit-solvent effects. Our protocol exploits SAXS pure-solute forward models and enhanced sampling methods to sample an heterogenous ensemble of structures, with no information towards the experiments provided on-the-fly. The generated structural ensemble is then reweighted through the maximum entropy principle so as to match reference SAXS experimental data at multiple ionic conditions. Importantly, accurate explicit-solvent forward models are used at this reweighting stage. We apply this framework to the GTPase-associated center, a relevant RNA molecule involved in protein translation, in order to elucidate its ion-dependent conformational ensembles. We show that (a) both solvent and dynamics are crucial to reproduce experimental SAXS data and (b) the resulting dynamical ensembles contain an ion-dependent fraction of extended structures.

Divalent ions tune the kinetics of a bacterial GTPase center rRNA folding transition from secondary to tertiary structure.

Folding of an RNA from secondary to tertiary structure often depends on divalent ions for efficient electrostatic charge screening (nonspecific association) or binding (specific association). To measure how different divalent cations modify folding kinetics of the 60 nucleotide Ecoli rRNA GTPase center, we combined stopped-flow fluorescence in the presence of Mg2+, Ca2+, or Sr2+ together with time-resolved small angle X-ray scattering (SAXS) in the presence of Mg2+ to observe the folding process. Immediately upon addition of each divalent ion, the RNA undergoes a transition from an extended state with secondary structure to a more compact structure. Subsequently, specific divalent ions modulate populations of intermediates in conformational ensembles along the folding pathway with transition times longer than 10 msec. Rate constants for the five folding transitions act on timescales from submillisecond to tens of seconds. The sensitivity of RNA tertiary structure to divalent cation identity affects all but the fastest events in RNA folding, and allowed us to identify those states that prefer Mg2+ The GTPase center RNA appears to have optimized its folding trajectory to specifically utilize this most abundant intracellular divalent ion.

Climbing the vertebrate branch of U1A/U2B″ protein evolution.

In the vertebrate lineage of the U1A/U2B″/SNF protein family, the U1A and U2B″ proteins bind to RNA stem-loops in the U1 or U2 snRNPs, respectively. However, their specialization is fairly recent, as they evolved from a single ancestral protein. The progress of their specialization (subfunctionalization) can be monitored by the amino acid sequence changes that give rise to their modern RNA-binding specificity. Using ancestral sequence reconstruction to predict the intermediates on the evolutionary branch, a probable path of sequential changes is defined for U1A and U2B″. The RNA-binding affinity for U1A/U2B″ protein ancestors was measured using modern U1 and U2 snRNA stem-loops and RNA stem-loop variants to understand how the proteins' RNA specificities evolved.

Effect of loop composition on the stability and folding kinetics of RNA hairpins with large loops.

RNA hairpins are ubiquitous structural elements in biological RNAs, where they have the potential to regulate RNA folding and interactions with other molecules. There are established methods for predicting the thermodynamic stability of an RNA hairpin, but there are still relatively few detailed examinations of the kinetics of folding. Nonetheless, several recent studies indicate that hairpin folding does not proceed via a simple two-state model. Here, we monitor fluorescence from hairpins constructed as molecular beacons in ensemble, fluorescence correlation spectroscopy, and stopped-flow experiments to describe the folding of RNA hairpins with long (15 nucleotide) loops. Our results show that folding of these hairpins occurs through more than two states and that the mechanism of folding includes a fast intermediate phase observed on the tens of microseconds time scale and a slow phase, attributed to formation of the native folded hairpin loop and stem, observed on the milliseconds time scale. The composition of the RNA loop determines the time scale of intermediate and native folded states. Hairpins with a polyuracil loop sequence exhibit slower relaxation of the intermediate state and faster relaxation of the native folded state when compared to that of hairpins with cytosine or adenine in the loop. We hypothesize this composition dependence could be attributed to nucleobase stacking in cytosine and adenine containing regions of the loop, which would be absent in hairpins containing polyuracil loops. Such base stacking could destabilize the intermediate folds, thereby speeding the relaxation of the intermediate relative to similar sized hairpins with no base stacking in the loop. Likewise, the lower intermediate stability could prolong the relaxation of the native folded state.

A compare-and-contrast NMR dynamics study of two related RRMs: U1A and SNF.

The U1A/U2B″/SNF family of small nuclear ribonucleoproteins uses a phylogenetically conserved RNA recognition motif (RRM1) to bind RNA stemloops in U1 and/or U2 small nuclear RNA (snRNA). RRMs are characterized by their α/ß sandwich topology, and these RRMs use their ß-sheet as the RNA binding surface. Unique to this RRM family is the tyrosine-glutamine-phenylalanine (YQF) triad of solvent-exposed residues that are displayed on the ß-sheet surface; the aromatic residues form a platform for RNA nucleobases to stack. U1A, U2B″, and SNF have very different patterns of RNA binding affinity and specificity, however, so here we ask how YQF in Drosophila SNF RRM1 contributes to RNA binding, as well as to domain stability and dynamics. Thermodynamic double-mutant cycles using tyrosine and phenylalanine substitutions probe the communication between those two residues in the free and bound states of the RRM. NMR experiments follow corresponding changes in the glutamine side-chain amide in both U1A and SNF, providing a physical picture of the RRM1 ß-sheet surface. NMR relaxation and dispersion experiments compare fast (picosecond to nanosecond) and intermediate (microsecond-to-millisecond) dynamics of U1A and SNF RRM1. We conclude that there is a network of amino acid interactions involving Tyr-Gln-Phe in both SNF and U1A RRM1, but whereas mutations of the Tyr-Gln-Phe triad result in small local responses in U1A, they produce extensive microsecond-to-millisecond global motions throughout SNF that alter the conformational states of the RRM.

Linkage and allostery in snRNP protein/RNA complexes.

Drosophila SNF is a member of the U1A/U2B″/SNF protein family that is found in U1 and U2 snRNPs, where it binds to Stemloop II and Stemloop IV of U1 and U2 snRNA, respectively. SNF also binds to the U2A' protein, but only in the U2 snRNP. Although previous reports have implicated U2A' as a necessary auxiliary protein for the binding of SNF to Stemloop IV, there are no mechanisms that explain the partitioning of U2A' to the U2 snRNP and its absence from the U1 snRNP. Using in vitro RNA binding isotherms and isothermal titration calorimetry, the thermodynamics of SNF/RNA/U2A' ternary complex formation have now been characterized. There is a very large binding cooperativity unique to Stemloop IV that favors formation of the SLIV/SNF/U2A' complex. The binding cooperativity, or heterotropic linkage, is interpreted with respect to linked conformational equilibria of both SNF and its RNA ligand and so represents an example of protein-RNA allostery.

Binding affinity and cooperativity control U2B″/snRNA/U2A' RNP formation.

The U1A and U2B″ proteins are components of the U1 and U2 snRNPs, respectively, where they bind to snRNA stemloops. While localization of U1A and U2B″ to their respective snRNP is a well-known phenomenon, binding of U2B″ to U2 snRNA is typically thought to be accompanied by the U2A' protein. The molecular mechanisms that lead to formation of the RNA/U2B″/U2A' complex and its localization to the U2 snRNP are investigated here, using a combination of in vitro RNA-protein and protein-protein fluorescence and isothermal titration calorimetry binding experiments. We find that U2A' protein binds to U2B″ with nanomolar affinity but binds to U1A with only micromolar affinity. In addition, there is RNA-dependent cooperativity (linkage) between protein-protein and protein-RNA binding. The unique combination of tight binding and cooperativity ensures that the U2A'/U2B″ complex is partitioned only to the U2 snRNP.

Intrinsic flexibility of snRNA hairpin loops facilitates protein binding.

Stem-loop II of U1 snRNA and Stem-loop IV of U2 snRNA typically have 10 or 11 nucleotides in their loops. The fluorescent nucleobase 2-aminopurine was used as a substitute for the adenines in each loop to probe the local and global structures and dynamics of these unusually long loops. Using steady-state and time-resolved fluorescence, we find that, while the bases in the loops are stacked, they are able to undergo significant local motion on the picosecond/nanosecond timescale. In addition, the loops have a global conformational change at low temperatures that occurs on the microsecond timescale, as determined using laser T-jump experiments. Nucleobase and loop motions are present at temperatures far below the melting temperature of the hairpin stem, which may facilitate the conformational change required for specific protein binding to these RNA loops.

The dimerization domain of SARS CoV 2 Nucleocapsid protein is partially disordered as a monomer and forms a high affinity dynamic complex.

The SARS-CoV-2 Nucleocapsid (N) is a 419 amino acids protein that drives the compaction and packaging of the viral genome. This compaction is aided not only by protein-RNA interactions, but also by protein-protein interactions that contribute to increasing the valence of the nucleocapsid protein. Here, we focused on quantifying the mechanisms that control dimer formation. Single-molecule Förster Resonance Energy Transfer enabled us to investigate the conformations of the dimerization domain in the context of the full-length protein as well as the energetics associated with dimerization. Under monomeric conditions, we observed significantly expanded configurations of the dimerization domain (compared to the folded dimer structure), which are consistent with a dynamic conformational ensemble. The addition of unlabeled protein stabilizes a folded dimer configuration with a high mean transfer efficiency, in agreement with predictions based on known structures. Dimerization is characterized by a dissociation constant of ∼ 12 nM at 23 O C and is driven by strong enthalpic interactions between the two protein subunits, which originate from the coupled folding and binding. Interestingly, the dimer structure retains some of the conformational heterogeneity of the monomeric units, and the addition of denaturant reveals that the dimer domain can significantly expand before being completely destabilized. Our findings suggest that the inherent flexibility of the monomer form is required to adopt the specific fold of the dimer domain, where the two subunits interlock with one another. We proposed that the retained flexibility of the dimer form may favor the capture and interactions with RNA, and that the temperature dependence of dimerization may explain some of the previous observations regarding the phase separation propensity of the N protein.

Excluded Volume and Weak Interactions in Crowded Solutions Modulate Conformations and RNA Binding of an Intrinsically Disordered Tail.

The cellular milieu is a solution crowded with a significant concentration of different components (proteins, nucleic acids, metabolites, etc.). Such a crowded environment affects protein conformations, dynamics, and interactions. Intrinsically disordered proteins and regions are particularly sensitive to these effects. Here, we investigate the impact on an intrinsically disordered tail that flanks a folded domain, the N-terminal domain, and the RNA-binding domain of the SARS-CoV-2 nucleocapsid protein. We mimic the crowded environment of the cell using polyethylene glycol (PEG) and study its impact on protein conformations using single-molecule Förster resonance energy transfer. We found that high-molecular-weight PEG induces a collapse of the disordered N-terminal tail, whereas low-molecular-weight PEG induces a chain expansion. Our data can be explained by accounting for two opposing contributions: favorable interactions between the protein and crowder molecules and screening of excluded volume interactions. We further characterized the interaction between protein and RNA in the presence of crowding agents. While for all PEG molecules tested, we observed an increase in the binding affinity, the trend is not monotonic as a function of the degree of PEG polymerization. This points to the role of nonspecific protein-PEG interactions on binding in addition to the entropic effects due to crowding. To separate the enthalpic and entropic components of the effects, we investigated the temperature dependence of the association constants in the absence and presence of crowders. Finally, we compared the effects of crowding across mutations in the disordered region and found that the threefold difference in association constants for two naturally occurring variants of the SARS-CoV-2 nucleocapsid protein is reduced to almost identical affinities in the presence of crowders. Overall, our data provide new insights into understanding and modeling the contribution of crowding effects on disordered regions, including the impact of interactions between proteins and crowders and their interplay when binding a ligand.

Millisecond time-scale folding and unfolding of DNA hairpins using rapid-mixing stopped-flow kinetics.

We report stopped-flow kinetics experiments to study the folding and unfolding of 5 base-pair stem and 21 nucleotide polythymidine loop DNA hairpins over various concentrations of NaCl. The reactions occurred on a time scale of milliseconds, considerably longer than the microsecond time scale suggested by previous kinetics studies of similar-sized hairpins. In comparison to a recent fluorescence correlation spectroscopy study (J. Am. Chem. Soc. 2006, 128, 1240-1249), we suggest the microsecond time-scale reactions are due to intermediate states and the millisecond time-scale reactions reported here are due to the formation of the fully folded DNA hairpin. These results support our view that DNA hairpin folding occurs via a minimum three-state mechanism.

The SARS-CoV-2 nucleocapsid protein is dynamic, disordered, and phase separates with RNA.

The SARS-CoV-2 nucleocapsid (N) protein is an abundant RNA-binding protein critical for viral genome packaging, yet the molecular details that underlie this process are poorly understood. Here we combine single-molecule spectroscopy with all-atom simulations to uncover the molecular details that contribute to N protein function. N protein contains three dynamic disordered regions that house putative transiently-helical binding motifs. The two folded domains interact minimally such that full-length N protein is a flexible and multivalent RNA-binding protein. N protein also undergoes liquid-liquid phase separation when mixed with RNA, and polymer theory predicts that the same multivalent interactions that drive phase separation also engender RNA compaction. We offer a simple symmetry-breaking model that provides a plausible route through which single-genome condensation preferentially occurs over phase separation, suggesting that phase separation offers a convenient macroscopic readout of a key nanoscopic interaction.

Coevolution of Drosophila snf protein and its snRNA targets.

SNF is a protein that is found in the U1 and U2 snRNPs (small nuclear ribonucleoproteins) of Drosophila. Its mammalian counterparts are two homologous proteins, U1A and U2B''. In vivo, these proteins segregate to the U1 and U2 snRNPs, respectively, where they bind distinct RNA hairpins. The RNA binding properties and mechanism of U1A have been studied extensively, but much less is known about SNF and U2B'' binding to their RNA targets. By comparing thermodynamic aspects of SNF-RNA interactions with those of U1A-RNA interactions, we find that SNF binds its RNA targets in a manner that is distinct from that of U1A. In vitro, SNF is able to bind both Drosophila U1 stem-loop II and U2 stem-loop IV with high affinity, although it binds stem-loop II more tightly than it binds stem-loop IV. Intriguingly, SNF is unable to bind human U2 stem-loop IV, which suggests that both the protein and RNAs have coevolved to interact with each other such that a single protein can bind RNAs that are more commonly bound by two distinct proteins.

A natural missing link between activated and downhill protein folding scenarios.

We propose protein PTB1 : 4W as a good candidate for engineering into a downhill folder. PTB1 : 4W has a probe-dependent thermal unfolding curve and sub-millisecond T-jump relaxation kinetics on more than one time scale. Its refolding rate in denaturant is a non-linear function of denaturant concentration (curved chevron plot). Yet at high denaturant concentration its unfolding is probe-independent, and the folding kinetics can be fitted to a single exponential decay. The domain appears to fold via a mechanism between downhill folding and activated folding over several small barriers, and when denaturant is added, one of these barriers greatly increases and simplifies the observed folding to apparent two-state kinetics. We predict the simplest free energy function consistent with the thermal denaturation and kinetics experiments by using the singular value Smoluchowski dynamics (SVSD) model. PTB1 : 4W is a natural 'missing link' between downhill and activated folding. We suggest mutations that could move the protein into the downhill folding limit.

Ribosomal Protein L11 Selectively Stabilizes a Tertiary Structure of the GTPase Center rRNA Domain.

The GTPase Center (GAC) RNA domain in bacterial 23S rRNA is directly bound by ribosomal protein L11, and this complex is essential to ribosome function. Previous cocrystal structures of the 58-nucleotide GAC RNA bound to L11 revealed the intricate tertiary fold of the RNA domain, with one monovalent and several divalent ions located in specific sites within the structure. Here, we report a new crystal structure of the free GAC that is essentially identical to the L11-bound structure, which retains many common sites of divalent ion occupation. This new structure demonstrates that RNA alone folds into its tertiary structure with bound divalent ions. In solution, we find that this tertiary structure is not static, but rather is best described as an ensemble of states. While L11 protein cannot bind to the GAC until the RNA has adopted its tertiary structure, new experimental data show that L11 binds to Mg2+-dependent folded states, which we suggest lie along the folding pathway of the RNA. We propose that L11 stabilizes a specific GAC RNA tertiary state, corresponding to the crystal structure, and that this structure reflects the functionally critical conformation of the rRNA domain in the fully assembled ribosome.