Deriving quantitative dynamics information for proteins and RNAs using ROTDIF with a graphical user interface

Key arginine residues in R2D2 dsRBD1 and dsRBD2 lead the siRNA recognition in Drosophila melanogaster RNAi pathway

In Drosophila melanogaster, Dcr-2:R2D2 heterodimer binds to the 21 nucleotide siRNA duplex to form the R2D2/Dcr-2 Initiator (RDI) complex, which is critical for the initiation of siRNA-induced silencing complex (RISC) assembly. During RDI complex formation, R2D2, a protein that contains three dsRNA binding domains (dsRBD), senses two aspects of the siRNA: thermodynamically more stable end (asymmetry sensing) and the 5'-phosphate (5'-P) recognition. Despite several detailed studies to date, the molecular determinants arising from R2D2 for performing these two tasks remain elusive. In this study, we have performed structural, biophysical, and biochemical characterization of R2D2 dsRBDs. We found that the solution NMR-derived structure of R2D2 dsRBD1 yielded a canonical α1-β1-β2-β3-α2 fold, wherein two arginine salt bridges provide additional stability to the R2D2 dsRBD1. Furthermore, we show that R2D2 dsRBD1 interacts with thermodynamically asymmetric siRNA duplex independent of its 5'-phosphorylation state, whereas R2D2 dsRBD2 prefers to interact with 5'-P siRNA duplex. The mutation of key arginine residues, R53 and R101, in concatenated dsRBDs of R2D2 results in a significant loss of siRNA duplex recognition. Our study deciphers the active roles of R2D2 dsRBDs by showing that dsRBD1 initiates siRNA recognition, whereas dsRBD2 senses 5'-phosphate as an authentic mark on functional siRNA.

Specifying conformational heterogeneity of multi-domain proteins at atomic resolution

The conformational landscape of multi-domain proteins is inherently linked to their specific functions. This also holds for polyubiquitin chains that are assembled by two or more ubiquitin domains connected by a flexible linker thus showing a large interdomain mobility. However, molecular recognition and signal transduction are associated with particular conformational substates that are populated in solution. Here, we apply high-resolution NMR spectroscopy in combination with dual-scale MD simulations to explore the conformational space of K6-, K29-, and K33-linked diubiquitin molecules. The conformational ensembles are evaluated utilizing a paramagnetic cosolute reporting on solvent exposure plus a set of complementary NMR parameters. This approach unravels a conformational heterogeneity of diubiquitins and explains the diversity of structural models that have been determined for K6-, K29-, and K33-linked diubiquitins in free and ligand-bound states so far. We propose a general application of the approach developed here to demystify multi-domain proteins occurring in nature.

Cytochrome c lysine acetylation regulates cellular respiration and cell death in ischemic skeletal muscle

Skeletal muscle is more resilient to ischemia-reperfusion injury than other organs. Tissue specific post-translational modifications of cytochrome c (Cytc) are involved in ischemia-reperfusion injury by regulating mitochondrial respiration and apoptosis. Here, we describe an acetylation site of Cytc, lysine 39 (K39), which was mapped in ischemic porcine skeletal muscle and removed by sirtuin5 in vitro. Using purified protein and cellular double knockout models, we show that K39 acetylation and acetylmimetic K39Q replacement increases cytochrome c oxidase (COX) activity and ROS scavenging while inhibiting apoptosis via decreased binding to Apaf-1, caspase cleavage and activity, and cardiolipin peroxidase activity. These results are discussed with X-ray crystallography structures of K39 acetylated (1.50 Å) and acetylmimetic K39Q Cytc (1.36 Å) and NMR dynamics. We propose that K39 acetylation is an adaptive response that controls electron transport chain flux, allowing skeletal muscle to meet heightened energy demand while simultaneously providing the tissue with robust resilience to ischemia-reperfusion injury.

Upstream of N-Ras C-terminal cold shock domains mediate poly(A) specificity in a novel RNA recognition mode and bind poly(A) binding protein

RNA binding proteins (RBPs) often engage multiple RNA binding domains (RBDs) to increase target specificity and affinity. However, the complexity of target recognition of multiple RBDs remains largely unexplored. Here we use Upstream of N-Ras (Unr), a multidomain RBP, to demonstrate how multiple RBDs orchestrate target specificity. A crystal structure of the three C-terminal RNA binding cold-shock domains (CSD) of Unr bound to a poly(A) sequence exemplifies how recognition goes beyond the classical ππ-stacking in CSDs. Further structural studies reveal several interaction surfaces between the N-terminal and C-terminal part of Unr with the poly(A)-binding protein (pAbp). All interactions are validated by mutational analyses and the high-resolution structures presented here will guide further studies to understand how both proteins act together in cellular processes.

Conformational Dynamics of Phytoglobin BvPgb1.2 from Beta vulgaris ssp. vulgaris

Plant hemoglobins, often referred to as phytoglobins, play important roles in abiotic stress tolerance. Several essential small physiological metabolites can be bound to these heme proteins. In addition, phytoglobins can catalyze a range of different oxidative reactions in vivo. These proteins are often oligomeric, but the degree and relevance of subunit interactions are largely unknown. In this study, we delineate which residues are involved in dimer formation of a sugar beet phytoglobin type 1.2 (BvPgb1.2) using NMR relaxation experiments. E. coli cells harboring a phytoglobin expression vector were cultivated in isotope-labeled (²H, ¹³C and ¹⁵N) M9 medium. The triple-labeled protein was purified to homogeneity using two chromatographic steps. Two forms of BvPgb1.2 were examined, the oxy-form and the more stable cyanide-form. Using three-dimensional triple-resonance NMR experiments, sequence-specific assignments for CN-bound BvPgb1.2 were achieved for 137 backbone amide cross-peaks in the ¹H-¹⁵N TROSY spectrum, which amounts to 83% of the total number of 165 expected cross-peaks. A large proportion of the non-assigned residues are located in α-helixes G and H, which are proposed to be involved in protein dimerization. Such knowledge around dimer formation will be instrumental for developing a better understanding of phytoglobins' roles in planta.

Dynamic allostery in substrate binding by human thymidylate synthase

Human thymidylate synthase (hTS) is essential for DNA replication and therefore a therapeutic target for cancer. Effective targeting requires knowledge of the mechanism(s) of regulation of this 72 kDa homodimeric enzyme. Here, we investigate the mechanism of binding cooperativity of the nucleotide substrate. We have employed exquisitely sensitive methyl-based CPMG and CEST NMR experiments enabling us to identify residues undergoing bifurcated linear 3-state exchange, including concerted switching between active and inactive conformations in the apo enzyme. The inactive state is populated to only ~1.3%, indicating that conformational selection contributes negligibly to the cooperativity. Instead, methyl rotation axis order parameters, determined by ²H transverse relaxation rates, suggest that rigidification of the enzyme upon substrate binding is responsible for the entropically-driven cooperativity. Lack of the rigidification in product binding and substrate binding to an N-terminally truncated enzyme, both non-cooperative, support this idea. In addition, the lack of this rigidification in the N-terminal truncation indicates that interactions between the flexible N-terminus and the rest of the protein, which are perturbed by substrate binding, play a significant role in the cooperativity—a novel mechanism of dynamic allostery. Together, these findings yield a rare depth of insight into the substrate binding cooperativity of an essential enzyme.

Cross-correlated relaxation rates provide facile exchange signature in selectively labeled RNA

Gerhard Wagner has made numerous contributions to NMR spectroscopy, particularly his developments in the field of spin-relaxation stand out in directly mapping the spectral density functions of proteins. He and his group developed experimental techniques to reveal the importance of dynamics to protein biological function and drug discovery. On his 75th birthday, we take this opportunity to highlight how some of those seminal ideas developed for proteins are being extended to RNAs. The role of dynamics in the structure and function of RNA has been a major interest in drug design and therapeutics. Here we present the use of cross-correlated relaxation rates (η_xy) from anti-TROSY (R_2α) and TROSY (R_2β) to rapidly obtain qualitative information about the chemical exchange taking place within the bacterial and human A-site RNA system while reducing the sets of relaxation experiments required to map dynamics. We show that η_xy correlates with the order parameter which gives information on how flexible or rigid a residue is. We further show R_2β/η_xy can rapidly be used to probe chemical exchange as seen from its agreement with R_ex. In addition, we report the ability of R_2β/η_xy to determine chemical exchange taking place within the bacterial A-site RNA during structural transitions at pH 6.2 and 6.5. Finally, comparison of the R_2β/η_xy ratios indicates bacterial A-site has greater R_2β/η_xy values for G19 (1.34 s^-1), A20 (1.38 s^-1), U23 (1.63 s^-1) and C24 (1.51 s^-1) than human A-site [A19 (0.76 s^-1), A20 (1.01 s^-1), U23 (0.74 s^-1) and C24 (0.71 s^-1)]. Taken together, we have shown that the chemical exchange can quickly be analyzed for RNA systems from cross-correlated relaxation rates.

Conformational Dynamics of the Hepatitis B Virus Pre-genomic RNA on Multiple Time Scales: Implications for Viral Replication

Human hepatitis B virus (HBV) replication is initiated by the binding of the viral polymerase (P) to epsilon (ε), an ≈85-nucleotide (nt) cis-acting regulatory stem-loop RNA located at the 5'-end of the pre-genomic RNA (pgRNA). This interaction triggers P and pgRNA packaging and protein-primed reverse transcription and is therefore an attractive therapeutic target. Our recent nuclear magnetic resonance (NMR) structure of ε provides a useful starting point toward a detailed understanding of HBV replication, and hints at the functional importance of ε dynamics. Here, we present a detailed description of ε motions on the ps to ns and μs to ms time scales by NMR spin relaxation and relaxation dispersion, respectively. We also carried out molecular dynamics simulations to provide additional insight into ε conformational dynamics. These data outline a series of complex motions on multiple time scales within ε. Moreover, these motions occur in mostly conserved nucleotides from structural regions (i.e., priming loop, pseudo-triloop, and U43 bulge) that biochemical and mutational studies have shown to be essential for P binding, P-pgRNA packaging, protein-priming, and DNA synthesis. Taken together, our work implicates RNA dynamics as an integral feature that governs HBV replication.

Structure and dynamics of the quaternary hunchback mRNA translation repression complex

A key regulatory process during Drosophila development is the localized suppression of the hunchback mRNA translation at the posterior, which gives rise to a hunchback gradient governing the formation of the anterior-posterior body axis. This suppression is achieved by a concerted action of Brain Tumour (Brat), Pumilio (Pum) and Nanos. Each protein is necessary for proper Drosophila development. The RNA contacts have been elucidated for the proteins individually in several atomic-resolution structures. However, the interplay of all three proteins during RNA suppression remains a long-standing open question. Here, we characterize the quaternary complex of the RNA-binding domains of Brat, Pum and Nanos with hunchback mRNA by combining NMR spectroscopy, SANS/SAXS, XL/MS with MD simulations and ITC assays. The quaternary hunchback mRNA suppression complex comprising the RNA binding domains is flexible with unoccupied nucleotides functioning as a flexible linker between the Brat and Pum-Nanos moieties of the complex. Moreover, the presence of the Pum-HD/Nanos-ZnF complex has no effect on the equilibrium RNA binding affinity of the Brat RNA binding domain. This is in accordance with previous studies, which showed that Brat can suppress mRNA independently and is distributed uniformly throughout the embryo.

Hinge like domain motion facilitates human RBMS1 protein binding to proto-oncogene c-myc promoter

DNA binding proteins recognize DNA specifically or non-specifically using direct and indirect readout mechanisms like sliding, hopping, and diffusion. However, a common difficulty in explicitly elucidating any particular mechanism of site-specific DNA-protein recognition is the lack of knowledge regarding target sequences and inadequate account of non-specific interactions, in general. Here, we decipher the structural basis of target search performed by the key regulator of expression of c-myc proto-oncogene, the human RBMS1 protein. In this study, we have shown the structural reorganization of this multi-domain protein required for recognizing the specific c-myc promoter sequence. The results suggest that a synergy between structural re-organization and thermodynamics is necessary for the recognition of target sequences. The study presents another perspective of looking at the DNA-protein interactions.

Quantifying the effects of long-range 13C-13C dipolar coupling on measured relaxation rates in RNA

Selective stable isotope labeling has transformed structural and dynamics analysis of RNA by NMR spectroscopy. These methods can remove ¹³C-¹³C dipolar couplings that complicate ¹³C relaxation analyses. While these phenomena are well documented for sites with adjacent ¹³C nuclei (e.g. ribose C1'), less is known about so-called isolated sites (e.g. adenosine C2). To investigate and quantify the effects of long-range (> 2 Å) ¹³C-¹³C dipolar interactions on RNA dynamics, we simulated adenosine C2 relaxation rates in uniformly [U-¹³C/¹⁵N]-ATP or selectively [2-¹³C]-ATP labeled RNAs. Our simulations predict non-negligible ¹³C-¹³C dipolar contributions from adenosine C4, C5, and C6 to C2 longitudinal (R₁) relaxation rates in [U-¹³C/¹⁵N]-ATP labeled RNAs. Moreover, these contributions increase at higher magnetic fields and molecular weights to introduce discrepancies that exceed 50%. This will become increasingly important at GHz fields. Experimental R₁ measurements in the 61 nucleotide human hepatitis B virus encapsidation signal ε RNA labeled with [U-¹³C/¹⁵N]-ATP or [2-¹³C]-ATP corroborate these simulations. Thus, in the absence of selectively labeled samples, long-range ¹³C-¹³C dipolar contributions must be explicitly taken into account when interpreting adenosine C2 R₁ rates in terms of motional models for large RNAs.

Structure analysis suggests Ess1 isomerizes the carboxy-terminal domain of RNA polymerase II via a bivalent anchoring mechanism

Accurate gene transcription in eukaryotes depends on isomerization of serine-proline bonds within the carboxy-terminal domain (CTD) of RNA polymerase II. Isomerization is part of the "CTD code" that regulates recruitment of proteins required for transcription and co-transcriptional RNA processing. Saccharomyces cerevisiae Ess1 and its human ortholog, Pin1, are prolyl isomerases that engage the long heptad repeat (YSPTSPS)₂₆ of the CTD by an unknown mechanism. Here, we used an integrative structural approach to decipher Ess1 interactions with the CTD. Ess1 has a rigid linker between its WW and catalytic domains that enforces a distance constraint for bivalent interaction with the ends of long CTD substrates (≥4-5 heptad repeats). Our binding results suggest that the Ess1 WW domain anchors the proximal end of the CTD substrate during isomerization, and that linker divergence may underlie evolution of substrate specificity.

How wide is the window opened by high-resolution relaxometry on the internal dynamics of proteins in solution?

The dynamics of molecules in solution is usually quantified by the determination of timescale-specific amplitudes of motions. High-resolution nuclear magnetic resonance (NMR) relaxometry experiments-where the sample is transferred to low fields for longitudinal (T₁) relaxation, and back to high field for detection with residue-specific resolution-seeks to increase the ability to distinguish the contributions from motion on timescales slower than a few nanoseconds. However, tumbling of a molecule in solution masks some of these motions. Therefore, we investigate to what extent relaxometry improves timescale resolution, using the "detector" analysis of dynamics. Here, we demonstrate improvements in the characterization of internal dynamics of methyl-bearing side chains by carbon-13 relaxometry in the small protein ubiquitin. We show that relaxometry data leads to better information about nanosecond motions as compared to high-field relaxation data only. Our calculations show that gains from relaxometry are greater with increasing correlation time of rotational diffusion.

Connecting Longitudinal and Transverse Relaxation Rates in Live-Cell NMR

In the cytosolic environment, protein crowding and Brownian motions result in numerous transient encounters. Each such encounter event increases the apparent size of the interacting molecules, leading to slower rotational tumbling. The extent of transient protein complexes formed in live cells can conveniently be quantified by an apparent viscosity, based on NMR-detected spin-relaxation measurements, that is, the longitudinal (T1) and transverse (T2) relaxation. From combined analysis of three different proteins and surface mutations thereof, we find that T2 implies significantly higher apparent viscosity than T1. At first sight, the effect on T1 and T2 seems thus nonunifiable, consistent with previous reports on other proteins. We show here that the T1 and T2 deviation is actually not a inconsistency but an expected feature of a system with fast exchange between free monomers and transient complexes. In this case, the deviation is basically reconciled by a model with fast exchange between the free-tumbling reporter protein and a transient complex with a uniform 143 kDa partner. The analysis is then taken one step further by accounting for the fact that the cytosolic content is by no means uniform but comprises a wide range of molecular sizes. Integrating over the complete size distribution of the cytosolic interaction ensemble enables us to predict both T1 and T2 from a single binding model. The result yields a bound population for each protein variant and provides a quantification of the transient interactions. We finally extend the approach to obtain a correction term for the shape of a database-derived mass distribution of the interactome in the mammalian cytosol, in good accord with the existing data of the cellular composition.

Deleterious effects of carbon-carbon dipolar coupling on RNA NMR dynamics

Many regulatory RNAs undergo dynamic exchanges that are crucial for their biological functions and NMR spectroscopy is a versatile tool for monitoring dynamic motions of biomolecules. Meaningful information on biomolecular dynamics requires an accurate measurement of relaxation parameters such as longitudinal (R₁) rates, transverse (R₂) rates and heteronuclear Overhauser effect (hNOE). However, earlier studies have shown that the large ¹³C-¹³C interactions complicate analysis of the carbon relaxation parameters. To investigate the effect of ¹³C-¹³C interactions on RNA dynamic studies, we performed relaxation measurements on various RNA samples with different labeling patterns and compared these measurements with the computational simulations. For uniformly labeled samples, contributions of the neighboring carbon to R₁ measurements were observed. These contributions increased with increasing magnetic field and overall correlation time ([Formula: see text]) for R₁ rates, necessitating more careful analysis for uniformly labeled large RNAs. In addition, the hNOE measurements were also affected by the adjacent carbon nuclei. Unlike R₁ rates, R_1ρ rates showed relatively good agreement between uniformly- and site-selectively labeled samples, suggesting no dramatic effect from their attached carbon, in agreement with previous observations. Overall, having more accurate rate measurements avoids complex analysis and will be a key for interpreting ¹³C relaxation rates for molecular motion that can provide valuable insights into cellular molecular recognition events.

13C NMR Relaxation Analysis of Protein GB3 for the Assessment of Side Chain Dynamics Predictions by Current AMBER and CHARMM Force Fields

Molecular simulations with seven current AMBER- and CHARMM-based force fields yield markedly differing internal bond vector autocorrelation function predictions for many of the 223 methine and methylene H-C bonds of the 56-residue protein GB3. To enable quantification of accuracy, ¹³C R₁, R₂, and heteronuclear NOE relaxation rates have been determined for the methine and stereochemically assigned methylene C^α and C^β positions. With only three experimental relaxation values for each bond vector, central to this analysis is the accuracy with which MD-derived autocorrelation curves can be represented by a 3-parameter equation which, in turn, maps onto the NMR relaxation values. In contrast to the more widely used extended Lipari-Szabo order parameter representation, 95% of these MD-derived internal autocorrelation curves for GB3 can be fitted to within 1.0% rmsd over the time frame from 30 ps to 4 ns by a biexponential Larmor frequency-selective representation (LF-S²). Applying the LF-S² representation to the experimental relaxation rates and uncertainties serves to determine the boundary range for the autocorrelation function of each bond vector consistent with the experimental data. Not surprisingly, all seven force fields predict the autocorrelation functions for the more motionally restricted ¹H^α-¹³C^α and ¹H^β-¹³C^β bond vectors with reasonable accuracy. However, for the ¹H^β-¹³C^β bond vectors exhibiting aggregate order parameter S² values less than 0.85, only 1% of the MD-derived predictions lie with 1 σ of the experimentally determined autocorrelation functions and only 7% within 2 σ. On the other hand, substantial residue type-specific improvements in predictive performance were observed among the recent AMBER force fields. This analysis indicates considerable potential for the use of ¹³C relaxation measurements in guiding the optimization of the side chain dynamics characteristics of protein molecular simulations.

Ultrafast Dynamics of Water-Protein Coupled Motions around the Surface of Eye Crystallin

Water dynamics on the protein surface mediate both protein structure and function. However, many questions remain about the role of the protein hydration layers in protein fluctuations and how the dynamics of these layers relate to specific protein properties. The fish eye lens protein γM7-crystallin (γM7) is found in vivo at extremely high concentrations nearing the packing limit, corresponding to only a few water layers between adjacent proteins. In this study, we conducted a site-specific probing of hydration water motions and side-chain dynamics at nine selected sites around the surface of γM7 using a tryptophan scan with femtosecond spectroscopy and NMR nuclear spin relaxation (NSR). We observed correlated fluctuations between hydration water and protein side chains on the time scales of a few picoseconds and hundreds of picoseconds, corresponding to local reorientations and network restructuring, respectively. These motions are heterogeneous over the protein surface and relate to the various steric and chemical properties of the local protein environment. Overall, we found that γM7 has relatively slower water dynamics within the hydration shell than a similar β-sheet protein, which may contribute to the high packing limit of this unique protein.

RNA Dynamics by NMR Spectroscopy

An ever-increasing number of functional RNAs require a mechanistic understanding. RNA function relies on changes in its structure, so-called dynamics. To reveal dynamic processes and higher energy structures, new NMR methods have been developed to elucidate these dynamics in RNA with atomic resolution. In this Review, we provide an introduction to dynamics novices and an overview of methods that access most dynamic timescales, from picoseconds to hours. Examples are provided as well as insight into theory, data acquisition and analysis for these different methods. Using this broad spectrum of methodology, unprecedented detail and invisible structures have been obtained and are reviewed here. RNA, though often more complicated and therefore neglected, also provides a great system to study structural changes, as these RNA structural changes are more easily defined-Lego like-than in proteins, hence the numerous revelations of RNA excited states.

Molecular recognition of ubiquitin and Lys63-linked diubiquitin by STAM2 UIM-SH3 dual domain: the effect of its linker length and flexibility

Multidomain proteins represent a broad spectrum of the protein landscape and are involved in various interactions. They could be considered as modular building blocks assembled in distinct fashion and connected by linkers of varying lengths and sequences. Due to their intrinsic flexibility, these linkers provide proteins a subtle way to modulate interactions and explore a wide range of conformational space. In the present study, we are seeking to understand the effect of the flexibility and dynamics of the linker involved in the STAM2 UIM-SH3 dual domain protein with respect to molecular recognition. We have engineered several constructs of UIM-SH3 with different length linkers or domain deletion. By means of SAXS and NMR experiments, we have shown that the modification of the linker modifies the flexibility and the dynamics of UIM-SH3. Indeed, the global tumbling of both the UIM and SH3 domain is different but not independent from each other while the length of the linker has an impact on the ps-ns time scale dynamics of the respective domains. Finally, the modification of the flexibility and dynamics of the linker has a drastic effect on the interaction of UIM-SH3 with Lys63-linked diubiquitin with a roughly eight-time weaker dissociation constant.

Modular type I polyketide synthase acyl carrier protein domains share a common N-terminally extended fold

Acyl carrier protein (ACP) domains act as interaction hubs within modular polyketide synthase (PKS) systems, employing specific protein-protein interactions to present acyl substrates to a series of enzyme active sites. Many domains from the multimodular PKS that generates the toxin mycolactone display an unusually high degree of sequence similarity, implying that the few sites which vary may do so for functional reasons. When domain boundaries based on prior studies were used to prepare two isolated ACP segments from this system for studies of their interaction properties, one fragment adopted the expected tertiary structure, but the other failed to fold, despite sharing a sequence identity of 49%. Secondary structure prediction uncovered a previously undetected helical region (H0) that precedes the canonical helix-bundle ACP topology in both cases. This article reports the NMR solution structures of two N-terminally extended mycolactone mACP constructs, mH0ACPa and mH0ACPb, both of which possess an additional α-helix that behaves like a rigid component of the domain. The interactions of these species with a phosphopantetheinyl transferase and a ketoreductase domain are unaffected by the presence of H0, but a shorter construct that lacks the H0 region is shown to be substantially less thermostable than mH0ACPb. Bioinformatics analysis suggests that the extended H0-ACP motif is present in 98% of type I cis-acyltransferase PKS chain-extension modules. The polypeptide linker that connects an H0-ACP motif to the preceding domain must therefore be ~12 residues shorter than previously thought, imposing strict limits on ACP-mediated substrate delivery within and between PKS modules.

UBC9 Mutant Reveals the Impact of Protein Dynamics on Substrate Selectivity and SUMO Chain Linkages

SUMO, a conserved ubiquitin-like protein, is conjugated to a multitude of cellular proteins to maintain genomic integrity and resist genotoxic stress. Studies of the SUMO E2 conjugating enzyme mutant, UBC9_P123L, suggested that altered substrate specificity enhances cell sensitivity to DNA damaging agents. Using nuclear magnetic resonance chemical shift studies, we confirm that the mutation does not alter the core globular fold of UBC9, while ¹⁵N relaxation measurements demonstrate mutant-induced stabilization of distinct chemical states in residues near the active site cysteine and substrate recognition motifs. We further demonstrate that the P123L substitution induces a switch from the preferential addition of SUMO to lysine residues in unstructured sites to acceptor lysines embedded in secondary structures, thereby also inducing alterations in SUMO chain linkages. Our results provide new insights regarding the impact that structural dynamics of UBC9 have on substrate selection and specifically SUMO chain formation. These findings highlight the potential contribution of nonconsensus SUMO targets and/or alternative SUMO chain linkages on DNA damage response and chemotherapeutic sensitivity.

High-Sensitivity Rheo-NMR Spectroscopy for Protein Studies

Shear stress can induce structural deformation of proteins, which might result in aggregate formation. Rheo-NMR spectroscopy has the potential to monitor structural changes in proteins under shear stress at the atomic level; however, existing Rheo-NMR methodologies have insufficient sensitivity to probe protein structure and dynamics. Here we present a simple and versatile approach to Rheo-NMR, which maximizes sensitivity by using a spectrometer equipped with a cryogenic probe. As a result, the sensitivity of the instrument ranks highest among the Rheo-NMR spectrometers reported so far. We demonstrate that the newly developed Rheo-NMR instrument can acquire high-quality relaxation data for a protein under shear stress and can trace structural changes in a protein during fibril formation in real time. The described approach will facilitate rheological studies on protein structural deformation, thereby aiding a physical understanding of shear-induced amyloid fibril formation.

Fe65-PTB2 Dimerization Mimics Fe65-APP Interaction

Physiological function and pathology of the Alzheimer’s disease causing amyloid precursor protein (APP) are correlated with its cytosolic adaptor Fe65 encompassing a WW and two phosphotyrosine-binding domains (PTBs). The C-terminal Fe65-PTB2 binds a large portion of the APP intracellular domain (AICD) including the GYENPTY internalization sequence fingerprint. AICD binding to Fe65-PTB2 opens an intra-molecular interaction causing a structural change and altering Fe65 activity. Here we show that in the absence of the AICD, Fe65-PTB2 forms a homodimer in solution and determine its crystal structure at 2.6 Å resolution. Dimerization involves the unwinding of a C-terminal α-helix that mimics binding of the AICD internalization sequence, thus shielding the hydrophobic binding pocket. Specific dimer formation is validated by nuclear magnetic resonance (NMR) techniques and cell-based analyses reveal that Fe65-PTB2 together with the WW domain are necessary and sufficient for dimerization. Together, our data demonstrate that Fe65 dimerizes via its APP interaction site, suggesting that besides intra- also intermolecular interactions between Fe65 molecules contribute to homeostatic regulation of APP mediated signaling.

Structure of the Escherichia coli ProQ RNA-binding protein

The protein ProQ has recently been identified as a global small noncoding RNA-binding protein in Salmonella, and a similar role is anticipated for its numerous homologs in divergent bacterial species. We report the solution structure of Escherichia coli ProQ, revealing an N-terminal FinO-like domain, a C-terminal domain that unexpectedly has a Tudor domain fold commonly found in eukaryotes, and an elongated bridging intradomain linker that is flexible but nonetheless incompressible. Structure-based sequence analysis suggests that the Tudor domain was acquired through horizontal gene transfer and gene fusion to the ancestral FinO-like domain. Through a combination of biochemical and biophysical approaches, we have mapped putative RNA-binding surfaces on all three domains of ProQ and modeled the protein's conformation in the apo and RNA-bound forms. Taken together, these data suggest how the FinO, Tudor, and linker domains of ProQ cooperate to recognize complex RNA structures and serve to promote RNA-mediated regulation.

Molecular impact of covalent modifications on nonribosomal peptide synthetase carrier protein communication

Nonribosomal peptide synthesis involves the interplay between covalent protein modifications, conformational fluctuations, catalysis, and transient protein-protein interactions. Delineating the mechanisms involved in orchestrating these various processes will deepen our understanding of domain-domain communication in nonribosomal peptide synthetases (NRPSs) and lay the groundwork for the rational reengineering of NRPSs by swapping domains handling different substrates to generate novel natural products. Although many structural and biochemical studies of NRPSs exist, few studies have focused on the energetics and dynamics governing the interactions in these systems. Here, we present detailed binding studies of an adenylation domain and its partner carrier protein in apo-, holo-, and substrate-loaded forms. Results from fluorescence anisotropy, isothermal titration calorimetry, and NMR titrations indicated that covalent modifications to a carrier protein modulate domain communication, suggesting that chemical modifications to carrier proteins during NRPS synthesis may impart directionality to sequential NRPS domain interactions. Comparison of the structure and dynamics of an apo-aryl carrier protein with those of its modified forms revealed structural fluctuations induced by post-translational modifications and mediated by modulations of protein dynamics. The results provide a comprehensive molecular description of a carrier protein throughout its life cycle and demonstrate how a network of dynamic residues can propagate the molecular impact of chemical modifications throughout a protein and influence its affinity toward partner domains.

Molecular Cross-Talk between Nonribosomal Peptide Synthetase Carrier Proteins and Unstructured Linker Regions

Nonribosomal peptide synthetases (NRPSs) employ multiple domains separated by linker regions to incorporate substrates into natural products. During synthesis, substrates are covalently tethered to carrier proteins that translocate between catalytic partner domains. The molecular parameters that govern translocation and associated linker remodeling remain unknown. Here, we used NMR to characterize the structure, dynamics, and invisible states of a peptidyl carrier protein flanked by its linkers. We showed that the N-terminal linker stabilizes and interacts with the protein core while modulating dynamics at specific sites involved in post-translational modifications and/or domain interactions. The results detail the molecular communication between peptidyl carrier proteins and their linkers and could guide efforts in engineering NRPSs to obtain new pharmaceuticals.

A critical assessment of methods to recover information from averaged data

Conformational heterogeneity is key to the function of many biomacromolecules, but only a few groups have tried to characterize it until recently. Now, thanks to the increased throughput of experimental data and the increased computational power, the problem of the characterization of protein structural variability has become more and more popular. Several groups have devoted their efforts in trying to create quantitative, reliable and accurate protocols for extracting such information from averaged data. We analyze here different approaches, discussing strengths and weaknesses of each. All approaches can roughly be clustered into two groups: those satisfying the maximum entropy principle and those recovering ensembles composed of a restricted number of molecular conformations. In the first case, the solution focuses on the features that are common to all the infinite solutions satisfying the experimental data; in the second case, the reconstructed ensemble shows the conformational regions where a large probability can be placed. The upper limits for conformational probabilities (MaxOcc) can also be calculated. We also give an overview of the mainstream experimental observables, with considerations on the assumptions underlying their usage.

Ubiquitylation Directly Induces Fold Destabilization of Proteins

Ubiquitin is a common post-translational modifier and its conjugation is a key signal for proteolysis by the proteasome. Because the molecular mass of ubiquitin is larger than that of other modifiers such as phosphate, acetyl, or methyl groups, ubiquitylation not only influences biochemical signaling, but also may exert physical effects on its substrate proteins by increasing molecular volume and altering shape anisotropy. Here we show that ubiquitylation destabilizes the fold of two proteins, FKBP12 and FABP4, and that elongation of the conjugated ubiquitin chains further enhances this destabilization effect. Moreover, NMR relaxation analysis shows that ubiquitylation induces characteristic structural fluctuations in the backbone of both proteins. These results suggest that the ubiquitylation-driven structural fluctuations lead to fold destabilization of its substrate proteins. Thus, physical destabilization by ubiquitylation may facilitate protein degradation by the proteasome.

Solution NMR structures of the C-domain of Tetrahymena cytoskeletal protein Tcb2 reveal distinct calcium-induced structural rearrangements

Tcb2 is a calcium-binding protein that localizes to the membrane-associated skeleton of the ciliated protozoan Tetrahymena thermophila with hypothesized roles in ciliary movement, cell cortex signaling, and pronuclear exchange. Tcb2 has also been implicated in a unique calcium-triggered, ATP-independent type of contractility exhibited by filamentous networks isolated from the Tetrahymena cytoskeleton. To gain insight into Tcb2's structure-function relationship and contractile properties, we determined solution NMR structures of its C-terminal domain in the calcium-free and calcium-bound states. The overall architecture is similar to other calcium-binding proteins, with paired EF-hand calcium-binding motifs. Comparison of the two structures reveals that Tcb2-C's calcium-induced conformational transition differs from the prototypical calcium sensor calmodulin, suggesting that the two proteins play distinct functional roles in Tetrahymena and likely have different mechanisms of target recognition. Future studies of the full-length protein and the identification of Tcb2 cellular targets will help establish the molecular basis of Tcb2 function and its unique contractile properties. Proteins 2016; 84:1748-1756. © 2016 Wiley Periodicals, Inc.

Linkage-specific conformational ensembles of non-canonical polyubiquitin chains

Polyubiquitination is a critical protein post-translational modification involved in a variety of processes in eukaryotic cells. The molecular basis for selective recognition of the polyubiquitin signals by cellular receptors is determined by the conformations polyubiquitin chains adopt; this has been demonstrated for K48- and K63-linked chains. Recent studies of the so-called non-canonical chains (linked via K6, K11, K27, K29, or K33) suggest they play important regulatory roles in growth, development, and immune system pathways, but biophysical studies are needed to elucidate the physical/structural basis of their interactions with receptors. A first step towards this goal is characterization of the conformations these chains adopt in solution. We assembled diubiquitins (Ub2) comprised of every lysine linkage. Using solution NMR measurements, small-angle neutron scattering (SANS), and in silico ensemble generation, we determined population-weighted conformational ensembles that shed light on the structure and dynamics of the non-canonical polyubiquitin chains. We found that polyubiquitin is conformationally heterogeneous, and each chain type exhibits unique conformational ensembles. For example, K6-Ub2 and K11-Ub2 (at physiological salt concentration) are in dynamic equilibrium between at least two conformers, where one exhibits a unique Ub/Ub interface, distinct from that observed in K48-Ub2 but similar to crystal structures of these chains. Conformers for K29-Ub2 and K33-Ub2 resemble recent crystal structures in the ligand-bound state. Remarkably, a number of diubiquitins adopt conformers similar to K48-Ub2 or K63-Ub2, suggesting potential overlap of biological function among different lysine linkages. These studies highlight the potential power of determining function from elucidation of conformational states.

Solution Structure of a Nonribosomal Peptide Synthetase Carrier Protein Loaded with Its Substrate Reveals Transient, Well-Defined Contacts

Nonribosomal peptide synthetases (NRPSs) are microbial enzymes that produce a wealth of important natural products by condensing substrates in an assembly line manner. The proper sequence of substrates is obtained by tethering them to phosphopantetheinyl arms of holo carrier proteins (CPs) via a thioester bond. CPs in holo and substrate-loaded forms visit NRPS catalytic domains in a series of transient interactions. A lack of structural information on substrate-loaded carrier proteins has hindered our understanding of NRPS synthesis. Here, we present the first structure of an NRPS aryl carrier protein loaded with its substrate via a native thioester bond, together with the structure of its holo form. We also present the first quantification of NRPS CP backbone dynamics. Our results indicate that prosthetic moieties in both holo and loaded forms are in contact with the protein core, but they also sample states in which they are disordered and extend in solution. We observe that substrate loading induces a large conformational change in the phosphopantetheinyl arm, thereby modulating surfaces accessible for binding to other domains. Our results are discussed in the context of NRPS domain interactions.

The structure of the SOLE element of oskar mRNA

mRNA localization by active transport is a regulated process that requires association of mRNPs with protein motors for transport along either the microtubule or the actin cytoskeleton. oskar mRNA localization at the posterior pole of the Drosophila oocyte requires a specific mRNA sequence, termed the SOLE, which comprises nucleotides of both exon 1 and exon 2 and is assembled upon splicing. The SOLE folds into a stem-loop structure. Both SOLE RNA and the exon junction complex (EJC) are required for oskar mRNA transport along the microtubules by kinesin. The SOLE RNA likely constitutes a recognition element for a yet unknown protein, which either belongs to the EJC or functions as a bridge between the EJC and the mRNA. Here, we determine the solution structure of the SOLE RNA by Nuclear Magnetic Resonance spectroscopy. We show that the SOLE forms a continuous helical structure, including a few noncanonical base pairs, capped by a pentanucleotide loop. The helix displays a widened major groove, which could accommodate a protein partner. In addition, the apical helical segment undergoes complex dynamics, with potential functional significance.

The structure of the SOLE element of oskar mRNA

mRNA localization by active transport is a regulated process that requires association of mRNPs with protein motors for transport along either the microtubule or the actin cytoskeleton. oskar mRNA localization at the posterior pole of the Drosophila oocyte requires a specific mRNA sequence, termed the SOLE, which comprises nucleotides of both exon 1 and exon 2 and is assembled upon splicing. The SOLE folds into a stem-loop structure. Both SOLE RNA and the exon junction complex (EJC) are required for oskar mRNA transport along the microtubules by kinesin. The SOLE RNA likely constitutes a recognition element for a yet unknown protein, which either belongs to the EJC or functions as a bridge between the EJC and the mRNA. Here, we determine the solution structure of the SOLE RNA by Nuclear Magnetic Resonance spectroscopy. We show that the SOLE forms a continuous helical structure, including a few noncanonical base pairs, capped by a pentanucleotide loop. The helix displays a widened major groove, which could accommodate a protein partner. In addition, the apical helical segment undergoes complex dynamics, with potential functional significance.

Structural fidelity and NMR relaxation analysis in a prototype RNA hairpin

RNA hairpins are widespread and very stable motifs that contribute decisively to RNA folding and biological function. The GTP1G2C3A4C5U6U7C8G9G10U11G12C13C14 construct (with a central UUCG tetraloop) has been extensively studied by solution NMR, and offers and excellent opportunity to evaluate the structure and dynamical description afforded by molecular dynamics (MD) simulations. Here, we compare average structural parameters and NMR relaxation rates estimated from a series of multiple independent explicit solvent MD simulations using the two most recent RNA AMBER force fields (ff99 and ff10). Predicted overall tumbling times are ∼20% faster than those inferred from analysis of NMR data and follow the same trend when temperature and ionic strength is varied. The Watson-Crick stem and the "canonical" UUCG loop structure are maintained in most simulations including the characteristic syn conformation along the glycosidic bond of G9, although some key hydrogen bonds in the loop are partially disrupted. Our analysis pinpoints G9-G10 backbone conformations as a locus of discrepancies between experiment and simulation. In general the results for the more recent force-field parameters (ff10) are closer to experiment than those for the older ones (ff99). This work provides a comprehensive and detailed comparison of state of the art MD simulations against a wide variety of solution NMR measurements.

PREdator: a python based GUI for data analysis, evaluation and fitting

The analysis of a series of experimental data is an essential procedure in virtually every field of research. The information contained in the data is extracted by fitting the experimental data to a mathematical model. The type of the mathematical model (linear, exponential, logarithmic, etc.) reflects the physical laws that underlie the experimental data. Here, we aim to provide a readily accessible, user-friendly python script for data analysis, evaluation and fitting. PREdator is presented at the example of NMR paramagnetic relaxation enhancement analysis.

pH dependence of conformational fluctuations of the protein backbone

Proton binding equilibria (pK(a) values) of ionizable groups in proteins are exquisitely sensitive to their microenvironments. Apparent pK(a) values measured for individual ionizable residues with NMR spectroscopy are actually population-weighted averages of the pK(a) in different conformational microstates. NMR spectroscopy experiments with staphylococcal nuclease were used to test the hypothesis that pK(a) values of surface Glu and Asp residues are affected by pH-sensitive fluctuations of the backbone between folded and locally unfolded conformations. (15)N spin relaxation studies showed that as the pH decreases from the neutral into the acidic range the amplitudes of backbone fluctuations in the ps-ns timescale increase near carboxylic residues. Hydrogen exchange experiments suggested that backbone conformational fluctuations promoted by decreasing pH also reflect slower local or sub-global unfolding near carboxylic groups. This study has implications for structure-based pKa calculations: (1) The timescale of the backbone's response to ionization events in proteins can range from ps to ms, and even longer; (2) pH-sensitive fluctuations of the backbone can be localized to both the segment the ionizable residue is attached to or the one that occludes the ionizable group; (3) Structural perturbations are not necessarily propagated through Coulomb interactions; instead, local fluctuations appear to be coupled through the co-operativity inherent to elements of secondary structure and to networks of hydrogen bonds. These results are consistent with the idea that local conformational fluctuations and stabilities are important determinants of apparent pK(a) values of ionizable residues in proteins.

NMR approaches for structural analysis of multidomain proteins and complexes in solution

NMR spectroscopy is a key method for studying the structure and dynamics of (large) multidomain proteins and complexes in solution. It plays a unique role in integrated structural biology approaches as especially information about conformational dynamics can be readily obtained at residue resolution. Here, we review NMR techniques for such studies focusing on state-of-the-art tools and practical aspects. An efficient approach for determining the quaternary structure of multidomain complexes starts from the structures of individual domains or subunits. The arrangement of the domains/subunits within the complex is then defined based on NMR measurements that provide information about the domain interfaces combined with (long-range) distance and orientational restraints. Aspects discussed include sample preparation, specific isotope labeling and spin labeling; determination of binding interfaces and domain/subunit arrangements from chemical shift perturbations (CSP), nuclear Overhauser effects (NOEs), isotope editing/filtering, cross-saturation, and differential line broadening; and based on paramagnetic relaxation enhancements (PRE) using covalent and soluble spin labels. Finally, the utility of complementary methods such as small-angle X-ray or neutron scattering (SAXS, SANS), electron paramagnetic resonance (EPR) or fluorescence spectroscopy techniques is discussed. The applications of NMR techniques are illustrated with studies of challenging (high molecular weight) protein complexes.