Search

Nuclear magnetic resonance (NMR) spectroscopy allows studying proteins in solution and under physiological temperatures. Frequently, either the amide groups of the protein backbone or the methyl groups in side chains are used as reporters of structural dynamics in proteins. A structural dynamics study of the protein backbone of globular proteins on 15 N labeled and fully protonated samples usually works well for proteins with a molecular weight of up to 50 kDa. When side chain deuteration in combination with transverse relaxation optimized spectroscopy (TROSY) is applied, this limit can be extended up to 200 kDa for globular proteins and up to 1 MDa when the focus is on the side chains. When intrinsically disordered proteins (IDPs) or proteins with intrinsically disordered regions (IDRs) are investigated, these weight limitations do not apply but can go well beyond. The reason is that IDPs or IDRs, characterized by high internal flexibility, are frequently dynamically decoupled. Various NMR methods offer atomic-resolution insights into structural protein dynamics across a wide range of time scales, from picoseconds up to hours. Standard 15 N relaxation measurements overview a protein's internal flexibility and characterize the protein backbone dynamics experienced on the fast pico- to nanosecond timescale. This article presents a hands-on protocol for setting up and recording NMR 15 N R1, R2, and heteronuclear Overhauser effect (hetNOE) experiments. We show exemplary data and explain how to interpret them simply qualitatively before any more sophisticated analysis.

NMR relaxation experiments provide residue-specific insights into the structural dynamics of proteins. Here, we present an optimized set of sensitivity-enhanced 15 N R 1 and R 1ρ relaxation experiments applicable to fully protonated proteins. The NMR pulse sequences are conceptually similar to the set of TROSY-based sequences and their HSQC counterpart (Lakomek et al., J. Biomol. NMR 2012). Instead of the TROSY read-out scheme, a sensitivity-enhanced HSQC read-out scheme is used, with improved and easier optimized water suppression. The presented pulse sequences are applied on the cytoplasmic domain of the SNARE protein Synpatobrevin-2 (Syb-2), which is intrinsically disordered in its monomeric pre-fusion state. A two-fold increase in the obtained signal-to-noise ratio is observed for this intrinsically disordered protein, therefore offering a four-fold reduction of measurement time compared to the TROSY-detected version. The inter-scan recovery delay can be shortened to two seconds. Pulse sequences were tested at 600 MHz and 1200 MHz 1 H Larmor frequency, thus applicable over a wide magnetic field range. A comparison between protonated and deuterated protein samples reveals high agreement, indicating that reliable 15 N R 1 and R 1ρ rate constants can be extracted for fully protonated and deuterated samples. The presented pulse sequences will benefit not only for IDPs but also for an entire range of low and medium-sized proteins., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Introduction: Misfolding of amyloidogenic proteins is a molecular hallmark of neurodegenerative diseases in humans. A detailed understanding of the underlying molecular mechanisms is mandatory for developing innovative therapeutic approaches. The bovine PI3K-SH3 domain has been a model system for aggregation and fibril formation. Methods: We monitored the fibril formation kinetics of low pH-denatured recombinantly expressed [U- 13 C, 15 N] labeled bovine PI3K-SH3 by a combination of solution NMR, high-resolution magic angle spinning (HR-MAS) NMR and solid-state NMR spectra. Solution NMR offers the highest sensitivity and, therefore, allows for the recording of two-dimensional NMR spectra with residue-specific resolution for individual time points of the time series. However, it can only follow the decay of the aggregating monomeric species. In solution NMR, aggregation occurs under quiescent experimental conditions. Solid-state NMR has lower sensitivity and allows only for the recording of one-dimensional spectra during the time series. Conversely, solid-state NMR is the only technique to detect disappearing monomers and aggregated species in the same sample by alternatingly recoding scalar coupling and dipolar coupling (CP)-based spectra. HR-MAS NMR is used here as a hybrid method bridging solution and solid-state NMR. In solid-state NMR and HR-MAS NMR the sample is agitated due to magic angle spinning. Results: Good agreement of the decay rate constants of monomeric SH3, measured by the three different NMR methods, is observed. Moderate MAS up to 8 kHz seems to influence the aggregation kinetics of seeded fibril formation only slightly. Therefore, under sufficient seeding (1% seeds used here), quiescent conditions (solution NMR), and agitated conditions deliver similar results, arguing against primary nucleation induced by MAS as a major contributor. Using solid-state NMR, we find that the amount of disappeared monomer corresponds approximately to the amount of aggregated species under the applied experimental conditions (250 µM PI3K-SH3, pH 2.5, 298 K, 1% seeds) and within the experimental error range. Data can be fitted by simple mono-exponential conversion kinetics, with lifetimes τ in the 14-38 h range. Atomic force microscopy confirms that fibrils substantially grew in length during the aggregation experiment. This argues for fibril elongation as the dominant growth mechanism in fibril mass (followed by the CP-based solid-state NMR signal). Conclusion: We suggest a combined approach employing both solution NMR and solid-state NMR, back-to-back, on two aliquots of the same sample under seeding conditions as an additional approach to follow monomer depletion and growth of fibril mass simultaneously. Atomic force microscopy images confirm fibril elongation as a major contributor to the increase in fibril mass., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Gardon, Becker, Rähse, Hölbling, Apostolidis, Schulz, Bochinsky, Gremer, Heise and Lakomek.)

The neuronal SNARE protein SNAP25a (isoform 2) forms part of the SNARE complex eliciting synaptic vesicle fusion during neuronal exocytosis. While the post-fusion cis-SNARE complex has been studied extensively, little is known about the pre-fusion conformation of SNAP25a. Here we analyze monomeric SNAP25a by NMR spectroscopy, further supported by small-angle X-ray scattering (SAXS) experiments. SAXS data indicate that monomeric SNAP25 is more compact than a Gaussian chain but still a random coil. NMR shows that for monomeric SNAP25a, before SNAP25a interacts with its SNARE partners to drive membrane fusion, only the N-terminal part (region A5 to V36) of the first SNARE motif, SN1 (L11 - L81), is helical, comprising two α-helices (ranging from A5 to Q20 and S25 toV36). From E37 onwards, SNAP25a is mostly disordered and displays high internal flexibility, including the C-terminal part of SN1, almost the entire second SNARE motif (SN2, N144-A199), and the connecting loop region. Apart from the N-terminal helices, only the C-termini of both SN1 (E73 - K79) and SN2 (region T190 - A199), as well as two short regions in the connecting loop (D99 - K102 and E123 - M127) show a weak α-helical propensity (α-helical population < 25%). We speculate that the N-terminal helices (A5 to Q20 and S25 to V36) which constitute the N-terminus of SN1 act as a nucleation site for initiating SNARE zippering., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

During the replication process of SARS-CoV-2, the main protease of the virus [3-chymotrypsin-like protease (3CL pro )] plays a pivotal role and is essential for the life cycle of the pathogen. Numerous studies have been conducted so far, which have confirmed 3CL pro as an attractive drug target to combat COVID-19. We describe a novel and efficient next-generation sequencing (NGS) supported phage display selection strategy for the identification of a set of SARS-CoV-2 3CL pro targeting peptide ligands that inhibit the 3CL protease, in a competitive or noncompetitive mode, in the low μM range. From the most efficient l-peptides obtained from the phage display, we designed all-d-peptides based on the retro-inverso (ri) principle. They had IC 50 values also in the low μM range and in combination, even in the sub-micromolar range. Additionally, the combination with Rutinprivir decreases 10-fold the IC 50 value of the competitive inhibitor. The inhibition modes of these d-ri peptides were the same as their respective l-peptide versions. Our results demonstrate that retro-inverso obtained all-d-peptides interact with high affinity and inhibit the SARS-CoV-2 3CL protease, thus reinforcing their potential for further development toward therapeutic agents. The here described d-ri peptides address limitations associated with current l-peptide inhibitors and are promising lead compounds. Further optimization regarding pharmacokinetic properties will allow the development of even more potent d-peptides to be used for the prevention and treatment of COVID-19.

In the last two decades, solid-state nuclear magnetic resonance (ssNMR) spectroscopy has transformed from a spectroscopic technique investigating small molecules and industrial polymers to a potent tool decrypting structure and underlying dynamics of complex biological systems, such as membrane proteins, fibrils, and assemblies, in near-physiological environments and temperatures. This transformation can be ascribed to improvements in hardware design, sample preparation, pulsed methods, isotope labeling strategies, resolution, and sensitivity. The fundamental engagement between nuclear spins and radio-frequency pulses in the presence of a strong static magnetic field is identical between solution and ssNMR, but the experimental procedures vastly differ because of the absence of molecular tumbling in solids. This review discusses routinely employed state-of-the-art static and MAS pulsed NMR methods relevant for biological samples with rotational correlation times exceeding 100's of nanoseconds. Recent developments in signal filtering approaches, proton methodologies, and multiple acquisition techniques to boost sensitivity and speed up data acquisition at fast MAS are also discussed. Several examples of protein structures (globular, membrane, fibrils, and assemblies) solved with ssNMR spectroscopy have been considered. We also discuss integrated approaches to structurally characterize challenging biological systems and some newly emanating subdisciplines in ssNMR spectroscopy.

Severe respiratory syndrome coronavirus-2 (SARS-CoV-2) is a highly contagious beta-class coronavirus. Although vaccinations have shown high efficacy, the emergence of novel variants of concern (VOCs) has already exhibited traits of immune evasion. Thus, the development of tailored antiviral medications for patients with incomplete, inefficient, or non-existent immunization, is essential. The attachment of viral surface proteins to the cell surface is the first crucial step in the viral replication cycle, which for SARS-CoV-2 is mediated by the high affinity interaction of the viral trimeric spike with the host cell surface-located human angiotensin converting enzyme-2 (hACE2). Here, we used a novel and efficient next generation sequencing (NGS) supported phage display strategy for the selection of a set of SARS-CoV-2 receptor binding domain (RBD)-targeting peptide ligands that bind to the target protein with low µM to nM dissociation constants. Compound CVRBDL-3 inhibits the SARS-CoV-2 spike protein association to hACE2 in a concentration-dependent manner for pre- as well as post-complex formation conditions. Further rational optimization yielded a CVRBDL-3 based divalent compound, which demonstrated inhibitory efficacy with an IC 50 value of 47 nM. The obtained compounds were not only efficient for the different spike constructs from the originally isolated "wt" SARS-CoV-2, but also for B.1.1.7 mutant trimeric spike protein. Our work demonstrates that phage display-derived peptide ligands are potential fusion inhibitors of viral cell entry. Moreover, we show that rational optimization of a combination of peptide sequences is a potential strategy in the further development of therapeutics for the treatment of acute COVID-19.

Chronic mental illnesses (CMIs) pose a significant challenge to global health due to their complex and poorly understood etiologies and hence, absence of causal therapies. Research of the past two decades has revealed dysfunction of the disrupted in schizophrenia 1 (DISC1) protein as a predisposing factor involved in several psychiatric disorders. DISC1 is a multifaceted protein that serves myriads of functions in mammalian cells, for instance, influencing neuronal development and synapse maintenance. It serves as a scaffold hub forming complexes with a variety (~300) of partners that constitute its interactome. Herein, using combinations of structural and biophysical tools, we demonstrate that the C-region of the DISC1 protein is highly polymorphic, with important consequences for its physiological role. Results from solid-state NMR spectroscopy and electron microscopy indicate that the protein not only forms symmetric oligomers but also gives rise to fibrils closely resembling those found in certain established amyloid proteinopathies. Furthermore, its aggregation as studied by isothermal titration calorimetry (ITC) is an exergonic process, involving a negative enthalpy change that drives the formation of oligomeric (presumably tetrameric) species as well as β-fibrils. We have been able to narrow down the β-core region participating in fibrillization to residues 716-761 of full-length human DISC1. This region is absent in the DISC1 Δ22aa splice variant, resulting in reduced association with proteins from the dynein motor complex, viz., NDE-like 1 (NDEL1) and lissencephaly 1 (LIS1), which are crucial during mitosis. By employing surface plasmon resonance, we show that the oligomeric DISC1 C-region has an increased affinity and shows cooperativity in binding to LIS1 and NDEL1, in contrast to the noncooperative binding mode exhibited by the monomeric version. Based on the derived structural models, we propose that the association between the binding partners involves two neighboring subunits of DISC1 C-region oligomers. Altogether, our findings highlight the significance of the DISC1 C-region as a crucial factor governing the balance between its physiological role as a multifunctional scaffold protein and aggregation-related aberrations with potential significance for disease., (© 2021. The Author(s).)

The Hepatitis C virus nonstructural protein 5A (NS5A) is a membrane-associated protein involved in multiple steps of the viral life cycle. Direct-acting antivirals (DAAs) targeting NS5A are a cornerstone of antiviral therapy, but the mode-of-action of these drugs is poorly understood. This is due to the lack of information on the membrane-bound NS5A structure. Herein, we present the structural model of an NS5A AH-linker-D1 protein reconstituted as proteoliposomes. We use highly sensitive proton-detected solid-state NMR methods suitable to study samples generated through synthetic biology approaches. Spectra analyses disclose that both the AH membrane anchor and the linker are highly flexible. Paramagnetic relaxation enhancements (PRE) reveal that the dimer organization in lipids requires a new type of NS5A self-interaction not reflected in previous crystal structures. In conclusion, we provide the first characterization of NS5A AH-linker-D1 in a lipidic environment shedding light onto the mode-of-action of clinically used NS5A inhibitors., (© 2020 The Authors. Angewandte Chemie International Edition published by Wiley-VCH GmbH.)

Proton-detected 100 kHz magic-angle-spinning (MAS) solid-state NMR is an emerging analysis method for proteins with only hundreds of microgram quantities, and thus allows structural investigation of eukaryotic membrane proteins. This is the case for the cell-free synthesized hepatitis C virus (HCV) nonstructural membrane protein 4B (NS4B). We demonstrate NS4B sample optimization using fast reconstitution schemes that enable lipid-environment screening directly by NMR. 2D spectra and relaxation properties guide the choice of the best sample preparation to record 2D 1 H-detected 1 H, 15 N and 3D 1 H, 13 C, 15 N correlation experiments with linewidths and sensitivity suitable to initiate sequential assignments. Amino-acid-selectively labeled NS4B can be readily obtained using cell-free synthesis, opening the door to combinatorial labeling approaches which should enable structural studies., (© 2019 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA.)

Intrinsically disordered proteins (IDPs) and their conformational transitions play an important role in neurotransmitter release at the neuronal synapse. Here, the SNARE proteins are essential by forming the SNARE complex that drives vesicular membrane fusion. While it is widely accepted that the SNARE proteins are intrinsically disordered in their monomeric prefusion form, important mechanistic aspects of this prefusion conformation and its lipid interactions, before forming the SNARE complex, are not fully understood at the molecular level and remain controversial. Here, by a combination of NMR and fluorescence spectroscopy methods, we find that vesicular synaptobrevin-2 (syb-2) in its monomeric prefusion conformation shows high flexibility, characteristic for an IDP, but also a high dynamic range and increasing rigidity from the N to C terminus. The gradual increase in rigidity correlates with an increase in lipid binding affinity from the N to C terminus. It could also explain the increased rate for C-terminal SNARE zippering, known to be faster than N-terminal SNARE zippering. Also, the syb-2 SNARE motif and, in particular, the linker domain show transient and weak membrane binding, characterized by a high off-rate and low (millimolar) affinity. The transient membrane binding of syb-2 may compensate for the repulsive forces between the two membranes and/or the SNARE motifs and the membranes, helping to destabilize the hydrophilic-hydrophobic boundary in the bilayer. Therefore, we propose that optimum flexibility and membrane binding of syb-2 regulate SNARE assembly and minimize repulsive forces during membrane fusion., Competing Interests: The authors declare no conflict of interest., (Copyright © 2019 the Author(s). Published by PNAS.)

15 N R 2 relaxation measurements are key for the elucidation of the dynamics of both folded and intrinsically disordered proteins (IDPs). Here we show, on the example of the intrinsically disordered protein α-synuclein and the folded domain PDZ2, that at physiological pH and near physiological temperatures amide-water exchange can severely skew Hahn-echo based 15 N R 2 relaxation measurements as well as low frequency data points in CPMG relaxation dispersion experiments. The nature thereof is the solvent exchange with deuterium in the sample buffer, which modulates the 15 N chemical shift tensor via the deuterium isotope effect, adding to the apparent relaxation decay which leads to systematic errors in the relaxation data. This results in an artificial increase of the measured apparent 15 N R 2 rate constants-which should not be mistaken with protein inherent chemical exchange contributions, R ex , to 15 N R 2 . For measurements of 15 N R 2 rate constants of IDPs and folded proteins at physiological temperatures and pH, we recommend therefore the use of a very low D 2 O molar fraction in the sample buffer, as low as 1%, or the use of an external D 2 O reference along with a modified 15 N R 2 Hahn-echo based experiment. This combination allows for the measurement of R ex contributions to 15 N R 2 originating from conformational exchange in a time window from µs to ms.

Through-bond J-coupling based experiments in solid-state NMR spectroscopy are challenging because the J couplings are typically much smaller than the dipolar couplings. This often leads to a lower transfer efficiency compared to dipolar-coupling based sequences. One of the reasons for the low transfer efficiency are the second-order cross terms involving the strong heteronuclear dipolar couplings leading to fast magnetization decay. Here, we show that by employing a symmetry-based C9 sequence, which was carefully selected to suppress second-order terms, efficient polarization transfers of up to 80% can be achieved without decoupling on fully protonated two-spin model systems at a MAS frequency of 55.5 kHz with rf-field amplitudes of about 25 kHz. In addition, we analyse the effects of rf inhomogeneity and crystallites selection due to the polarization preparation method on the TOBSY transfer efficiency. We demonstrate on small model substances as well as on deuterated and 100% back-exchanged ubiquitin that C9 39 1 and C9 48 1 are efficient and practical TOBSY sequences at experimental conditions ranging from proton Larmor frequencies of 400-850 MHz, and MAS frequencies ranging from 55.5 to 111.1 kHz., (Copyright © 2017 Elsevier Inc. All rights reserved.)

The structural and dynamical characterization of membrane proteins in a lipid bilayer at physiological pH and temperature and free of crystal constraints is crucial for the elucidation of a structure/dynamics-activity relationship. Toward this aim, we explore here the properties of the outer-membrane protein OmpX embedded in lipid bilayer nanodiscs using proton-detected magic angle spinning (MAS) solid-state NMR at 60 and 110 kHz. [ 1 H, 15 N]-correlation spectra overlay well with the corresponding solution-state NMR spectra. Line widths as well as line intensities in solid and solution both depend critically on the sample temperature and, in particular, on the crossing of the lipid phase transition temperature. MAS (110 kHz) experiments yield well-resolved NMR spectra also for fully protonated OmpX and both below and above the lipid phase transition temperature.

15 N R 1ρ relaxation experiments in solid-state NMR spectroscopy are sensitive to timescales and amplitudes of internal protein motions in the hundreds of nano- to microsecond time window, which is difficult to probe by solution-state NMR spectroscopy. By using 15 N R 1ρ relaxation experiments, a simplified approach to detect low microsecond protein dynamics is described and residue-specific correlation times are determined from the ratio of 15 N R 1ρ rate constants at different magic angle spinning frequencies. Microcrystalline ubiquitin exhibits small-amplitude dynamics on a timescale of about 1 μs across the entire protein, and larger amplitude motions, also on the 1 μs timescale, for several sites, including the β 1 -β 2 turn and the N terminus of the α helix. According to the analysis, the microsecond protein backbone dynamics are of lower amplitude than that concluded in previous solid-state NMR spectroscopy studies, but persist across the entire protein with a rather uniform timescale of 1 μs., (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Detergents are often used to investigate the structure and dynamics of membrane proteins. Whereas the structural integrity seems to be preserved in detergents for many membrane proteins, their functional activity is frequently compromised, but can be restored in a lipid environment. Herein we show with per-residue resolution that while OmpX forms a stable β-barrel in DPC detergent micelles, DHPC/DMPC bicelles, and DMPC nanodiscs, the pico- to nanosecond and micro- to millisecond motions differ substantially between the detergent and lipid environment. In particular for the β-strands, there is pronounced dynamic variability in the lipid environment, which appears to be suppressed in micelles. This unexpected complex and membrane-mimetic-dependent dynamic behavior indicates that the frequent loss of membrane protein activity in detergents might be related to reduced internal dynamics and that membrane protein activity correlates with lipid flexibility., (© 2016 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA.)

Proteins are inserted into the bacterial plasma membrane cotranslationally after translating ribosomes are targeted to the translocon in the membrane via the signal recognition particle (SRP) pathway. The targeting pathway involves an interaction between SRP and the SRP receptor, FtsY. Here we focus on the role of FtsY and its interaction with the translocon in controlling targeting. We show that in unbound FtsY the NG and A domains interact with one another. The interaction involves the membrane-targeting region at the junction between A and N domain. The closed form of FtsY is impaired in binding to SRP. Upon binding to the phospholipid-embedded translocon the domains of FtsY move apart. This enhances the docking of the FtsY NG domain to the homologous NG domain of the SRP protein Ffh. Thus, FtsY binding to the translocon has a central role in orchestrating the formation of a quaternary transfer complex in which the nascent peptide is transferred to the translocon. We propose that FtsY activation at the translocon ensures that ribosome-SRP complexes are directed to available translocons. This way sequestering SRP in futile complexes with unbound FtsY can be avoided and efficient targeting to the translocon achieved., (© 2016 The Authors. Molecular Microbiology Published by John Wiley & Sons Ltd.)

Integral membrane proteins in bacteria are co-translationally targeted to the SecYEG translocon for membrane insertion via the signal recognition particle (SRP) pathway. The SRP receptor FtsY and its N-terminal A domain, which is lacking in any structural model of FtsY, were studied using NMR and fluorescence spectroscopy. The A domain is mainly disordered and highly flexible; it binds to lipids via its N terminus and the C-terminal membrane targeting sequence. The central A domain binds to the translocon non-specifically and maintains disorder. Translocon targeting and binding of the A domain is driven by electrostatic interactions. The intrinsically disordered A domain tethers FtsY to the translocon, and because of its flexibility, allows the FtsY NG domain to scan a large area for binding to the NG domain of ribosome-bound SRP, thereby promoting the formation of the quaternary transfer complex at the membrane., (© 2016 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA.)

By applying [1-(13) C]- and [2-(13) C]-glucose labeling schemes to the folded globular protein ubiquitin, a strong reduction of spectral crowding and increase in resolution in solid-state NMR (ssNMR) spectra could be achieved. This allowed spectral resonance assignment in a straightforward manner and the collection of a wealth of long-range distance information. A high precision solid-state NMR structure of microcrystalline ubiquitin was calculated with a backbone rmsd of 1.57 to the X-ray structure and 1.32 Å to the solution NMR structure. Interestingly, we can resolve structural heterogeneity as the presence of three slightly different conformations. Structural heterogeneity is most significant for the loop region β1-β2 but also for β-strands β1, β2, β3, and β5 as well as for the loop connecting α1 and β3. This structural polymorphism observed in the solid-state NMR spectra coincides with regions that showed dynamics in solution NMR experiments on different timescales., (© 2015 The Protein Society.)

Human immunodeficiency viral (HIV-1) fusion is mediated by the viral envelope gp120/gp41 complex (ENVelope glycoprotein). After gp120 shedding, gp41 is exposed and elicits membrane fusion via a cascade of conformational changes. In contrast to prefusion and postfusion conformation, little is known about any intermediate conformation. We report on a solution NMR investigation of homotrimeric HIV-1 gp41(27-194), comprising the transmembrane region and reconstituted in dodecyl phosphocholine (DPC) micelles. The protein is mainly α-helical, but experiences internal dynamics on the nanosecond and micro to millisecond time scale and transient α-helical behavior for certain residues in the N-terminal heptad repeat (NHR). Strong lipid interactions are observed, in particular for C-terminal residues of the NHR and imunodominant loop region connecting NHR and C-terminal heptad repeat (CHR). Our data indicate an extended conformation with features anticipated for a prefusion intermediate, presumably in exchange with a lowly populated postfusion six-helical bundle conformation., (Copyright © 2014 Elsevier Ltd. All rights reserved.)

While extracting dynamics parameters from backbone (15)N relaxation measurements in proteins has become routine over the past two decades, it is increasingly recognized that accurate quantitative analysis can remain limited by the potential presence of systematic errors associated with the measurement of (15)N R(1) and R(2) or R(1ρ) relaxation rates as well as heteronuclear (15)N-{(1)H} NOE values. We show that systematic errors in such measurements can be far larger than the statistical error derived from either the observed signal-to-noise ratio, or from the reproducibility of the measurement. Unless special precautions are taken, the problem of systematic errors is shown to be particularly acute in perdeuterated systems, and even more so when TROSY instead of HSQC elements are used to read out the (15)N magnetization through the NMR-sensitive (1)H nucleus. A discussion of the most common sources of systematic errors is presented, as well as TROSY-based pulse schemes that appear free of systematic errors to the level of <1 %. Application to the small perdeuterated protein GB3, which yields exceptionally high S/N and therefore is an ideal test molecule for detection of systematic errors, yields relaxation rates that show considerably less residue by residue variation than previous measurements. Measured R(2)'/R(1)' ratios fit an axially symmetric diffusion tensor with a Pearson's correlation coefficient of 0.97, comparable to fits obtained for backbone amide RDCs to the Saupe matrix.

Long-range correlated motions in proteins are candidate mechanisms for processes that require information transfer across protein structures, such as allostery and signal transduction. However, the observation of backbone correlations between distant residues has remained elusive, and only local correlations have been revealed using residual dipolar couplings measured by NMR spectroscopy. In this work, we experimentally identified and characterized collective motions spanning four β-strands separated by up to 15 Å in ubiquitin. The observed correlations link molecular recognition sites and result from concerted conformational changes that are in part mediated by the hydrogen-bonding network.

This study presents the first application of the model-free analysis (MFA) (Meiler in J Am Chem Soc 123:6098-6107, 2001; Lakomek in J Biomol NMR 34:101-115, 2006) to methyl group RDCs measured in 13 different alignment media in order to describe their supra-tau (c) dynamics in ubiquitin. Our results indicate that methyl groups vary from rigid to very mobile with good correlation to residue type, distance to backbone and solvent exposure, and that considerable additional dynamics are effective at rates slower than the correlation time tau (c). In fact, the average amplitude of motion expressed in terms of order parameters S (2) associated with the supra-tau (c) window brings evidence to the existence of fluctuations contributing as much additional mobility as those already present in the faster ps-ns time scale measured from relaxation data. Comparison to previous results on ubiquitin demonstrates that the RDC-derived order parameters are dominated both by rotameric interconversions and faster libration-type motions around equilibrium positions. They match best with those derived from a combined J-coupling and residual dipolar coupling approach (Chou in J Am Chem Soc 125:8959-8966, 2003) taking backbone motion into account. In order to appreciate the dynamic scale of side chains over the entire protein, the methyl group order parameters are compared to existing dynamic ensembles of ubiquitin. Of those recently published, the broadest one, namely the EROS ensemble (Lange in Science 320:1471-1475, 2008), fits the collection of methyl group order parameters presented here best. Last, we used the MFA-derived averaged spherical harmonics to perform highly-parameterized rotameric searches of the side chains conformation and find expanded rotamer distributions with excellent fit to our data. These rotamer distributions suggest the presence of concerted motions along the side chains.

Conformational ensembles are increasingly recognized as a useful representation to describe fundamental relationships between protein structure, dynamics and function. Here we present an ensemble of ubiquitin in solution that is created by sampling conformational space without experimental information using "Backrub" motions inspired by alternative conformations observed in sub-Angstrom resolution crystal structures. Backrub-generated structures are then selected to produce an ensemble that optimizes agreement with nuclear magnetic resonance (NMR) Residual Dipolar Couplings (RDCs). Using this ensemble, we probe two proposed relationships between properties of protein ensembles: (i) a link between native-state dynamics and the conformational heterogeneity observed in crystal structures, and (ii) a relation between dynamics of an individual protein and the conformational variability explored by its natural family. We show that the Backrub motional mechanism can simultaneously explore protein native-state dynamics measured by RDCs, encompass the conformational variability present in ubiquitin complex structures and facilitate sampling of conformational and sequence variability matching those occurring in the ubiquitin protein family. Our results thus support an overall relation between protein dynamics and conformational changes enabling sequence changes in evolution. More practically, the presented method can be applied to improve protein design predictions by accounting for intrinsic native-state dynamics.

RDCs (residual dipolar couplings) in NMR spectroscopy provide information about protein dynamics complementary to NMR relaxation methods, especially in the previously inaccessible time window between the protein correlation time tau(c) and 50 micros. For ubiquitin, new modes of motion of the protein backbone could be detected using RDC-based techniques. An ensemble of ubiquitin based on these RDC values is found to comprise all different conformations that ubiquitin adopts upon binding to different recognition proteins. These conformations in protein-protein complexes had been derived from 46 X-ray structures. Thus, for ubiquitin recognition by other proteins, conformational selection rather than induced fit seems to be the dominant mechanism.

Residual dipolar couplings (RDCs) provide information about the dynamic average orientation of inter-nuclear vectors and amplitudes of motion up to milliseconds. They complement relaxation methods, especially on a time-scale window that we have called supra-tau(c) (tau(c) < supra-tau(c) < 50 micros). Here we present a robust approach called Self-Consistent RDC-based Model-free analysis (SCRM) that delivers RDC-based order parameters-independent of the details of the structure used for alignment tensor calculation-as well as the dynamic average orientation of the inter-nuclear vectors in the protein structure in a self-consistent manner. For ubiquitin, the SCRM analysis yields an average RDC-derived order parameter of the NH vectors 0.72 +/- 0.02 compared to = 0.778 +/- 0.003 for the Lipari-Szabo order parameters, indicating that the inclusion of the supra-tau(c) window increases the averaged amplitude of mobility observed in the sub-supra-tau(c) window by about 34%. For the beta-strand spanned by residues Lys48 to Leu50, an alternating pattern of backbone NH RDC order parameter S2(rdc)(NH) = (0.59, 0.72, 0.59) was extracted. The backbone of Lys48, whose side chain is known to be involved in the poly-ubiquitylation process that leads to protein degradation, is very mobile on the supra-tau(c) time scale (S2(rdc)(NH) = 0.59 +/- 0.03), while it is inconspicuous (S2(LS)(NH)= 0.82) on the sub-tau(c) as well as on micros-ms relaxation dispersion time scales. The results of this work differ from previous RDC dynamics studies of ubiquitin in the sense that the results are essentially independent of structural noise providing a much more robust assessment of dynamic effects that underlie the RDC data.

Protein dynamics are essential for protein function, and yet it has been challenging to access the underlying atomic motions in solution on nanosecond-to-microsecond time scales. We present a structural ensemble of ubiquitin, refined against residual dipolar couplings (RDCs), comprising solution dynamics up to microseconds. The ensemble covers the complete structural heterogeneity observed in 46 ubiquitin crystal structures, most of which are complexes with other proteins. Conformational selection, rather than induced-fit motion, thus suffices to explain the molecular recognition dynamics of ubiquitin. Marked correlations are seen between the flexibility of the ensemble and contacts formed in ubiquitin complexes. A large part of the solution dynamics is concentrated in one concerted mode, which accounts for most of ubiquitin's molecular recognition heterogeneity and ensures a low entropic complex formation cost.

The presence of slow motions with large amplitudes, as detected by measurements based on residual dipolar couplings [Peti, W., Meiler, J., Brueschweiler, R. and Griesinger, C. (2002) J. Am. Chem. Soc., 124, 5822-5833], has stirred up much discussion in recent years. Based on ubiquitin NH residual dipolar couplings (rdcs) measured in 31 different alignment conditions, a model-free analysis of structure and dynamics [Meiler, J., Peti, W., Prompers, J., Griesinger, C. and Brueschweiler, R. (2001) J. Am. Chem. Soc., 123, 6098-6107] is presented. Starting from this broad experimental basis, rdc-based order parameters with so far unattained accuracy were determined. These rdc-based order parameters underpin the presence of new modes of motion slower than the inverse overall tumbling correlation time. Amplitudes and anisotropies of the motion were derived. The effect of structural noise on the results was proven to be negligible.