RNA dynamics from experimental and computational approaches.

Conformational dynamics is crucial for the biological function of RNA molecules and for their potential as therapeutic targets. This meeting report outlines key "take-home" messages that emerged from the presentations and discussions during the CECAM workshop "RNA dynamics from experimental and computational approaches" in Paris, June 26-28, 2023.

Molecular insights into the interaction between a disordered protein and a folded RNA.

Intrinsically disordered protein regions (IDRs) are well-established as contributors to intermolecular interactions and the formation of biomolecular condensates. In particular, RNA-binding proteins (RBPs) often harbor IDRs in addition to folded RNA-binding domains that contribute to RBP function. To understand the dynamic interactions of an IDR-RNA complex, we characterized the RNA-binding features of a small (68 residues), positively charged IDR-containing protein, SERF. At high concentrations, SERF and RNA undergo charge-driven associative phase separation to form a protein- and RNA-rich dense phase. A key advantage of this model system is that this threshold for demixing is sufficiently high that we could use solution-state biophysical methods to interrogate the stoichiometric complexes of SERF with RNA in the one-phase regime. Herein, we describe our comprehensive characterization of SERF alone and in complex with a small fragment of the HIV-1 TAR RNA (TAR) with complementary biophysical methods and molecular simulations. We find that this binding event is not accompanied by the acquisition of structure by either molecule; however, we see evidence for a modest global compaction of the SERF ensemble when bound to RNA. This behavior likely reflects attenuated charge repulsion within SERF via binding to the polyanionic RNA and provides a rationale for the higher-order assembly of SERF in the context of RNA. We envision that the SERF-RNA system will lower the barrier to accessing the details that support IDR-RNA interactions and likewise deepen our understanding of the role of IDR-RNA contacts in complex formation and liquid-liquid phase separation.

Initial Primer Synthesis of a DNA Primase Monitored by Real-Time NMR Spectroscopy.

Primases are crucial enzymes for DNA replication, as they synthesize a short primer required for initiating DNA replication. We herein present time-resolved nuclear magnetic resonance (NMR) spectroscopy in solution and in the solid state to study the initial dinucleotide formation reaction of archaeal pRN1 primase. Our findings show that the helix-bundle domain (HBD) of pRN1 primase prepares the two substrates and then hands them over to the catalytic domain to initiate the reaction. By using nucleotide triphosphate analogues, the reaction is substantially slowed down, allowing us to study the initial dinucleotide formation in real time. We show that the sedimented protein-DNA complex remains active in the solid-state NMR rotor and that time-resolved 31P-detected cross-polarization experiments allow monitoring the kinetics of dinucleotide formation. The kinetics in the sedimented protein sample are comparable to those determined by solution-state NMR. Protein conformational changes during primer synthesis are observed in time-resolved 1H-detected experiments at fast magic-angle spinning frequencies (100 kHz). A significant number of spectral changes cluster in the HBD pointing to the importance of the HBD for positioning the nucleotides and the dinucleotide.

Kinetic Resolution of the Atomic 3D Structures Formed by Ground and Excited Conformational States in an RNA Dynamic Ensemble.

Knowing the 3D structures formed by the various conformations populating the RNA free-energy landscape, their relative abundance, and kinetic interconversion rates is required to obtain a quantitative and predictive understanding of how RNAs fold and function at the atomic level. While methods integrating ensemble-averaged experimental data with computational modeling are helping define the most abundant conformations in RNA ensembles, elucidating their kinetic rates of interconversion and determining the 3D structures of sparsely populated short-lived RNA excited conformational states (ESs) remains challenging. Here, we developed an approach integrating Rosetta-FARFAR RNA structure prediction with NMR residual dipolar couplings and relaxation dispersion that simultaneously determines the 3D structures formed by the ground-state (GS) and ES subensembles, their relative abundance, and kinetic rates of interconversion. The approach is demonstrated on HIV-1 TAR, whose six-nucleotide apical loop was previously shown to form a sparsely populated (∼13%) short-lived (lifetime ∼ 45 µs) ES. In the GS, the apical loop forms a broad distribution of open conformations interconverting on the pico-to-nanosecond time scale. Most residues are unpaired and preorganized to bind the Tat-superelongation protein complex. The apical loop zips up in the ES, forming a narrow distribution of closed conformations, which sequester critical residues required for protein recognition. Our work introduces an approach for determining the 3D ensemble models formed by sparsely populated RNA conformational states, provides a rare atomic view of an RNA ES, and kinetically resolves the atomic 3D structures of RNA conformational substates, interchanging on time scales spanning 6 orders of magnitude, from picoseconds to microseconds.

Conformational Plasticity of Parallel G-Quadruplex─Implications on Duplex-Quadruplex Motifs.

DNA G-quadruplexes are essential motifs in molecular biology performing a wide range of functions enabled by their unique and diverse structures. In this study, we focus on the conformational plasticity of the most abundant and biologically relevant parallel G-quadruplex topology. A multipronged approach of structure survey, solution-state NMR spectroscopy, and molecular dynamics simulations unravels subtle yet essential features of the parallel G-quadruplex topology. Stark differences in flexibility are observed for the nucleotides depending upon their positioning in the tetrad planes that are intricately correlated with the conformational sampling of the propeller loop. Importantly, the terminal nucleotides in the 5'-end versus the 3'-end of the parallel quadruplex display differential dynamics that manifests their ability to accommodate a duplex on either end of the G-quadruplex. The conformational plasticity characterized in this study provides essential cues toward biomolecular processes such as small molecular binding, intermolecular quadruplex stacking, and implications on how a duplex influences the structure of a neighboring quadruplex.

Elongated Bacterial Pili as a Versatile Alignment Medium for NMR Spectroscopy.

In NMR spectroscopy, residual dipolar couplings (RDCs) have emerged as one of the most exquisite probes of biological structure and dynamics. The measurement of RDCs relies on the partial alignment of the molecule of interest, for example by using a liquid crystal as a solvent. Here, we establish bacterial type 1 pili as an alternative liquid-crystalline alignment medium for the measurement of RDCs. To achieve alignment at pilus concentrations that allow for efficient NMR sample preparation, we elongated wild-type pili by recombinant overproduction of the main structural pilus subunit. Building on the extraordinary stability of type 1 pili against spontaneous dissociation and unfolding, we show that the medium is compatible with challenging experimental conditions such as high temperature, the presence of detergents, organic solvents or very acidic pH, setting it apart from most established alignment media. Using human ubiquitin, HIV-1 TAR RNA and camphor as spectroscopic probes, we demonstrate the applicability of the medium for the determination of RDCs of proteins, nucleic acids and small molecules. Our results show that type 1 pili represent a very useful alternative to existing alignment media and may readily assist the characterization of molecular structure and dynamics by NMR.

Adenosine-Dependent Activation Mechanism of Prodrugs Targeting an Aminoacyl-tRNA Synthetase.

Prodrugs have little or no pharmacological activity and are converted to active drugs in the body by enzymes, metabolic reactions, or through human-controlled actions. However, prodrugs promoting their chemical bioconversion without any of these processes have not been reported before. Here, we present an enzyme-independent prodrug activation mechanism by boron-based compounds (benzoxaboroles) targeting leucyl-tRNA synthetase (LeuRS), including an antibiotic that recently has completed phase II clinical trials to cure tuberculosis. We combine nuclear magnetic resonance spectroscopy and X-ray crystallography with isothermal titration calorimetry to show that these benzoxaboroles do not bind directly to their drug target LeuRS, instead they are prodrugs that activate their bioconversion by forming a highly specific and reversible LeuRS inhibition adduct with ATP, AMP, or the terminal adenosine of the tRNALeu. We demonstrate how the oxaborole group of the prodrugs cyclizes with the adenosine ribose at physiological concentrations to form the active molecule. This bioconversion mechanism explains the remarkably good druglike properties of benzoxaboroles showing efficacy against radically different human pathogens and fully explains the mechanism of action of these compounds. Thus, this adenosine-dependent activation mechanism represents a novel concept in prodrug chemistry that can be applied to improve the solubility, permeability and metabolic stability of challenging drugs.

Multimodal Response to Copper Binding in Superoxide Dismutase Dynamics.

Copper/zinc superoxide dismutase (SOD) is a homodimeric metalloenzyme that has been extensively studied as a benchmark for structure-function relationships in proteins, in particular because of its implication in the familial form of the neurodegenerative disease amyotrophic lateral sclerosis. Here, we investigate microcrystalline preparations of two differently metalated forms of SOD, namely, the fully mature functional Cu,Zn state and the E,Zn-SOD state in which the Cu site is empty. By using solid-state NMR with fast magic-angle spinning (MAS) at high magnetic fields (1H Larmor frequency of 800-1000 MHz), we quantify motions spanning a dynamic range from ns to ms. We determine that metal ion uptake does not act as a rigidification element but as a switch redistributing motional processes on different time scales, with coupling of the dynamics of histidine side chains and those of remote key backbone elements of the protein.

Protein Side-Chain-DNA Contacts Probed by Fast Magic-Angle Spinning NMR.

Protein-nucleic acid interactions are essential in a variety of biological events ranging from the replication of genomic DNA to the synthesis of proteins. Noncovalent interactions guide such molecular recognition events, and protons are often at the center of them, particularly due to their capability of forming hydrogen bonds to the nucleic acid phosphate groups. Fast magic-angle spinning experiments (100 kHz) reduce the proton NMR line width in solid-state NMR of fully protonated protein-DNA complexes to such an extent that resolved proton signals from side-chains coordinating the DNA can be detected. We describe a set of NMR experiments focusing on the detection of protein side-chains from lysine, arginine, and aromatic amino acids and discuss the conclusions that can be obtained on their role in DNA coordination. We studied the 39 kDa enzyme of the archaeal pRN1 primase complexed with DNA and characterize protein-DNA contacts in the presence and absence of bound ATP molecules.

Identifying dynamic, partially occupied residues using anomalous scattering.

Although often presented as taking single `snapshots' of the conformation of a protein, X-ray crystallography provides an averaged structure over time and space within the crystal. The important but difficult task of characterizing structural ensembles in crystals is typically limited to small conformational changes, such as multiple side-chain conformations. A crystallographic method was recently introduced that utilizes residual electron and anomalous density (READ) to characterize structural ensembles encompassing large-scale structural changes. Key to this method is an ability to accurately measure anomalous signals and distinguish them from noise or other anomalous scatterers. This report presents an optimized data-collection and analysis strategy for partially occupied iodine anomalous signals. Using the long-wavelength-optimized beamline I23 at Diamond Light Source, the ability to accurately distinguish the positions of anomalous scatterers with occupancies as low as ∼12% is demonstrated. The number and positions of these anomalous scatterers are consistent with previous biophysical, kinetic and structural data that suggest that the protein Im7 binds to the chaperone Spy in multiple partially occupied conformations. Finally, READ selections demonstrate that re-measured data using the new protocols are consistent with the previously characterized structural ensemble of the chaperone Spy with its client Im7. This study shows that a long-wavelength beamline results in easily validated anomalous signals that are strong enough to be used to detect and characterize highly disordered sections of crystal structures.

A Small Helical Bundle Prepares Primer Synthesis by Binding Two Nucleotides that Enhance Sequence-Specific Recognition of the DNA Template.

Primases have a fundamental role in DNA replication. They synthesize a primer that is then extended by DNA polymerases. Archaeoeukaryotic primases require for synthesis a catalytic and an accessory domain, the exact contribution of the latter being unresolved. For the pRN1 archaeal primase, this domain is a 115-amino acid helix bundle domain (HBD). Our structural investigations of this small HBD by liquid- and solid-state nuclear magnetic resonance (NMR) revealed that only the HBD binds the DNA template. DNA binding becomes sequence-specific after a major allosteric change in the HBD, triggered by the binding of two nucleotide triphosphates. The spatial proximity of the two nucleotides and the DNA template in the quaternary structure of the HBD strongly suggests that this small domain brings together the substrates to prepare the first catalytic step of primer synthesis. This efficient mechanism is likely general for all archaeoeukaryotic primases.

Selecting Conformational Ensembles Using Residual Electron and Anomalous Density (READ).

Heterogeneous and dynamic biomolecular complexes play a central role in many cellular processes but are poorly understood due to experimental challenges in characterizing their structural ensembles. To address these difficulties, we developed a hybrid methodology that combines X-ray crystallography with ensemble selections typically used in NMR studies to determine structural ensembles of heterogeneous biomolecular complexes. The method, termed READ, for residual electron and anomalous density, enables the visualization of heterogeneous conformational ensembles of complexes within crystals. Here we present a detailed protocol for performing the ensemble selections to construct READ ensembles. From a diverse pool of binding poses, a selection scheme is used to determine a subset of conformations that maximizes agreement with the X-ray data. Overall, READ is a general approach for obtaining a high-resolution view of dynamic protein-protein complexes.

Activity of Tumor Necrosis Factor α Is Modulated by Dynamic Conformational Rearrangements.

The homotrimeric ligand tumor necrosis factor α (TNFα) is a key cytokine and immune regulator; however, when deregulated, it leads to several major chronic inflammatory diseases. Perturbation of the protein-protein interface has proven to be an efficient strategy to inactivate TNFα, but the atomic-resolution mechanism of its inactivation remains poorly understood. Here, we probe the solution structure and dynamics of active and inactive TNFα using NMR spectroscopy. The data reveal that TNFα undergoes motions on different time scales. Furthermore, by site-directed mutagenesis of residues at the trimerization interface and by targeting the interface with a low molecular weight inhibitor, we show that TNFα retains its overall structure and trimeric state. However, upon perturbation, TNFα exhibits increased conformational dynamics spanning from the trimerization interface to the regions mediating receptor binding. These findings provide novel insights into the inactivation mechanism of TNFα and the basis for strategies to target TNFα activity.

The Mechanism of HdeA Unfolding and Chaperone Activation.

HdeA is a periplasmic chaperone that is rapidly activated upon shifting the pH to acidic conditions. This activation is thought to involve monomerization of HdeA. There is evidence that monomerization and partial unfolding allow the chaperone to bind to proteins denatured by low pH, thereby protecting them from aggregation. We analyzed the acid-induced unfolding of HdeA using NMR spectroscopy and fluorescence measurements, and obtained experimental evidence suggesting a complex mechanism in HdeA's acid-induced unfolding pathway, as previously postulated from molecular dynamics simulations. Counterintuitively, dissociation constant measurements show a stabilization of the HdeA dimer upon exposure to mildly acidic conditions. We provide experimental evidence that protonation of Glu37, a glutamate residue embedded in a hydrophobic pocket of HdeA, is important in controlling HdeA stabilization and thus the acid activation of this chaperone. Our data also reveal a sharp transition from folded dimer to unfolded monomer between pH3 and pH 2, and suggest the existence of a low-populated, partially folded intermediate that could assist in chaperone activation or function. Overall, this study provides a detailed experimental investigation into the mechanism by which HdeA unfolds and activates.

Dynamic Descriptions of Highly Flexible Molecules from NMR Dipolar Couplings: Physical Basis and Limitations.

Biomolecules that control physiological function by changing their conformation play key roles in biology and remain poorly characterized. NMR dipolar couplings (DCs) depend intrinsically on both molecular shape and structural fluctuations, thereby providing the enticing prospect of tracking these conformational changes at atomic detail. Although this dual dependence has until now severely complicated analysis of DCs from highly dynamic systems, general approaches have recently been proposed that simplify interpretation of experimental DCs, by entirely eliminating molecular alignment from the analysis. Using simple and intuitive simulation of target ensembles, we investigate the impact of such approaches on the resulting descriptions of the conformational energy landscape. We find that ensemble descriptions of highly flexible systems derived from DCs without explicit consideration of the alignment properties of the constituent conformations can be compromised and inaccurate, despite exhibiting high correlation with experimental measurement.

Capturing a Dynamic Chaperone-Substrate Interaction Using NMR-Informed Molecular Modeling.

Chaperones maintain a healthy proteome by preventing aggregation and by aiding in protein folding. Precisely how chaperones influence the conformational properties of their substrates, however, remains unclear. To achieve a detailed description of dynamic chaperone-substrate interactions, we fused site-specific NMR information with coarse-grained simulations. Our model system is the binding and folding of a chaperone substrate, immunity protein 7 (Im7), with the chaperone Spy. We first used an automated procedure in which NMR chemical shifts inform the construction of system-specific force fields that describe each partner individually. The models of the two binding partners are then combined to perform simulations on the chaperone-substrate complex. The binding simulations show excellent agreement with experimental data from multiple biophysical measurements. Upon binding, Im7 interacts with a mixture of hydrophobic and hydrophilic residues on Spy's surface, causing conformational exchange within Im7 to slow down as Im7 folds. Meanwhile, the motion of Spy's flexible loop region increases, allowing for better interaction with different substrate conformations, and helping offset losses in Im7 conformational dynamics that occur upon binding and folding. Spy then preferentially releases Im7 into a well-folded state. Our strategy has enabled a residue-level description of a dynamic chaperone-substrate interaction, improving our understanding of how chaperones facilitate substrate folding. More broadly, we validate our approach using two other binding partners, showing that this approach provides a general platform from which to investigate other flexible biomolecular complexes through the integration of NMR data with efficient computational models.

Visualizing chaperone-assisted protein folding.

Challenges in determining the structures of heterogeneous and dynamic protein complexes have greatly hampered past efforts to obtain a mechanistic understanding of many important biological processes. One such process is chaperone-assisted protein folding. Obtaining structural ensembles of chaperone-substrate complexes would ultimately reveal how chaperones help proteins fold into their native state. To address this problem, we devised a new structural biology approach based on X-ray crystallography, termed residual electron and anomalous density (READ). READ enabled us to visualize even sparsely populated conformations of the substrate protein immunity protein 7 (Im7) in complex with the Escherichia coli chaperone Spy, and to capture a series of snapshots depicting the various folding states of Im7 bound to Spy. The ensemble shows that Spy-associated Im7 samples conformations ranging from unfolded to partially folded to native-like states and reveals how a substrate can explore its folding landscape while being bound to a chaperone.

Inter-helical conformational preferences of HIV-1 TAR-RNA from maximum occurrence analysis of NMR data and molecular dynamics simulations.

Detecting conformational heterogeneity in biological macromolecules is a key for the understanding of their biological function. We here provide a comparison between two independent approaches to assess conformational heterogeneity: molecular dynamics simulations, performed without inclusion of any experimental data, and maximum occurrence (MaxOcc) distribution over the topologically available conformational space. The latter only reflects the extent of the averaging and identifies regions which are most compliant with the experimentally measured NMR Residual Dipolar Couplings (RDCs). The analysis was performed for the HIV-1 TAR RNA, consisting of two helical domains connected by a flexible bulge junction, for which four sets of RDCs were available as well as an 8.2 µs all-atom molecular dynamics simulation. A sample and select approach was previously applied to extract from the molecular dynamics trajectory conformational ensembles in agreement with the four sets of RDCs. The MaxOcc analysis performed here identifies the most likely sampled region in the conformational space of the system which, strikingly, overlaps well with the structures independently sampled in the molecular dynamics calculations and even better with the RDC selected ensemble.

Investigating protein conformational energy landscapes and atomic resolution dynamics from NMR dipolar couplings: a review.

Nuclear magnetic resonance spectroscopy is exquisitely sensitive to protein dynamics. In particular inter-nuclear dipolar couplings, that become measurable in solution when the protein is dissolved in a dilute liquid crystalline solution, report on all conformations sampled up to millisecond timescales. As such they provide the opportunity to describe the Boltzmann distribution present in solution at atomic resolution, and thereby to map the conformational energy landscape in unprecedented detail. The development of analytical methods and approaches based on numerical simulation and their application to numerous biologically important systems is presented.