The structure of human ubiquitin in 2-methyl-2,4-pentanediol: a new conformational switch

Ubiquitin's Conformational Heterogeneity as Discerned by Nuclear Magnetic Resonance Spectroscopy

Visualizing a protein's molecular motions has been a long standing topic of research in the biophysics community. Largely this has been done by exploiting nuclear magnetic resonance spectroscopy (NMR), and arguably no protein's molecular motions have been better characterized by NMR than that of ubiquitin (Ub), a 76 amino acid polypeptide essential in ubiquitination-a key regulatory system within cells. Herein, we discuss ubiquitin's conformational plasticity as visualized, at atomic resolution, by more than 35 years of NMR work. In our discussions we point out the differences between data acquired in vitro, ex vivo, as well as in vivo and stress the need to investigate Ub's conformational plasticity in more biologically representative backgrounds.

Near atmospheric carbon dioxide activates plant ubiquitin cross-linking

Background Identifying CO₂-binding proteins is vital for our knowledge of CO₂-regulated molecular processes. The carbamate post-translational modification is a reversible CO₂-mediated adduct that can form on neutral N-terminal α-amino or lysine ε-amino groups. Methods We have developed triethyloxonium ion (TEO) as a chemical proteomics tool to trap the carbamate post-translational modification on protein covalently. We use ¹³C-NMR and TEO and identify ubiquitin as a plant CO₂-binding protein. Results We observe the carbamate post-translational modification on the Arabidopsis thaliana ubiquitin ε-amino groups of lysines 6, 33, and 48. We show that biologically relevant near atmospheric PCO₂ levels increase ubiquitin conjugation dependent on lysine 6. We further demonstrate that CO₂ increases the ubiquitin E2 ligase (AtUBC5) charging step via the transthioesterification reaction in which Ub is transferred from the E1 ligase active site to the E2 active site. Conclusions and general significance Therefore, plant ubiquitin is a CO₂-binding protein, and the carbamate post-translational modification represents a potential mechanism through which plant cells can respond to fluctuating PCO₂.

Aromatic ring flips in differently packed ubiquitin protein crystals from MAS NMR and MD

Probing the dynamics of aromatic side chains provides important insights into the behavior of a protein because flips of aromatic rings in a protein's hydrophobic core report on breathing motion involving a large part of the protein. Inherently invisible to crystallography, aromatic motions have been primarily studied by solution NMR. The question how packing of proteins in crystals affects ring flips has, thus, remained largely unexplored. Here we apply magic-angle spinning NMR, advanced phenylalanine 1H-13C/2H isotope labeling and MD simulation to a protein in three different crystal packing environments to shed light onto possible impact of packing on ring flips. The flips of the two Phe residues in ubiquitin, both surface exposed, appear remarkably conserved in the different crystal forms, even though the intermolecular packing is quite different: Phe4 flips on a ca. 10-20 ns time scale, and Phe45 are broadened in all crystals, presumably due to µs motion. Our findings suggest that intramolecular influences are more important for ring flips than intermolecular (packing) effects.

A litmus test for classifying recognition mechanisms of transiently binding proteins

Partner recognition in protein binding is critical for all biological functions, and yet, delineating its mechanism is challenging, especially when recognition happens within microseconds. We present a theoretical and experimental framework based on straight-forward nuclear magnetic resonance relaxation dispersion measurements to investigate protein binding mechanisms on sub-millisecond timescales, which are beyond the reach of standard rapid-mixing experiments. This framework predicts that conformational selection prevails on ubiquitin’s paradigmatic interaction with an SH3 (Src-homology 3) domain. By contrast, the SH3 domain recognizes ubiquitin in a two-state binding process. Subsequent molecular dynamics simulations and Markov state modeling reveal that the ubiquitin conformation selected for binding exhibits a characteristically extended C-terminus. Our framework is robust and expandable for implementation in other binding scenarios with the potential to show that conformational selection might be the design principle of the hubs in protein interaction networks.

The authors provide a litmus test for the recognition mechanism of transiently binding proteins based on nuclear magnetic resonance and find a conformational selection binding mechanism through concentration-dependent kinetics of ubiquitin and SH3.

Selective 1H-1H recoupling via symmetry sequences in fully protonated samples at fast magic angle spinning

Proton-detected solid-state NMR at fast Magic Angle Spinning (MAS) is becoming the norm to characterize molecules. Routinely ¹H-¹H and ¹H-X dipolar couplings are used to characterize the structure and dynamics of molecules. Selective proton recoupling techniques are emerging as a method for structural characterization via estimation of qualitative and quantitative distances. In the present study, we demonstrate through numerical simulations and experiments that the well-characterized CN^v_n sequences can also be tailored for selective recoupling of proton spins by employing C elements of the type (β)_Φ(4β)_Φ+π(3β)_Φ. Herein, several CN^v_n sequences were examined through numerical simulations and experiments. C6¹₄ recoupling sequence with a modified POST-element ((β)_Φ(4β)_Φ+π(3β)_Φ) shows selective polarization transfer efficiencies on the order of 40-50% between various proton spin pairs in fully protonated samples at rf amplitudes ranging from 0.3 to 0.8 times the MAS frequency. These selective recoupling sequences have been labeled as frequency-selective-CN^v_n sequences. The extent of selectivity, polarization transfer efficiency and the feasibility of experimentally measuring proton-proton distances in fully protonated samples are explored here. The development of efficient and robust selective ¹H-¹H recoupling experiments is required to structurally characterize molecules without artificial isotope enrichment or the need for diffracting crystals.

Temperature-Dependent Solid-State NMR Proton Chemical-Shift Values and Hydrogen Bonding

Temperature-dependent NMR experiments are often complicated by rather long magnetic-field equilibration times, for example, occurring upon a change of sample temperature. We demonstrate that the fast temporal stabilization of a magnetic field can be achieved by actively stabilizing the temperature of the magnet bore, which allows quantification of the weak temperature dependence of a proton chemical shift, which can be diagnostic for the presence of hydrogen bonds. Hydrogen bonding plays a central role in molecular recognition events from both fields, chemistry and biology. Their direct detection by standard structure-determination techniques, such as X-ray crystallography or cryo-electron microscopy, remains challenging due to the difficulties of approaching the required resolution, on the order of 1 Å. We, herein, explore a spectroscopic approach using solid-state NMR to identify protons engaged in hydrogen bonds and explore the measurement of proton chemical-shift temperature coefficients. Using the examples of a phosphorylated amino acid and the protein ubiquitin, we show that fast magic-angle spinning (MAS) experiments at 100 kHz yield sufficient resolution in proton-detected spectra to quantify the rather small chemical-shift changes upon temperature variations.

Protein Loop Conformational Free Energy Changes via an Alchemical Path without Reaction Coordinates

We introduce a method called restrain-free energy perturbation-release 2.0 (R-FEP-R 2.0) to estimate conformational free energy changes of protein loops via an alchemical path. R-FEP-R 2.0 is a generalization of the method called restrain-free energy perturbation-release (R-FEP-R) that can only estimate conformational free energy changes of protein side chains but not loops. The reorganization of protein loops is a central feature of many biological processes. Unlike other advanced sampling algorithms such as umbrella sampling and metadynamics, R-FEP-R and R-FEP-R 2.0 do not require predetermined collective coordinates and transition pathways that connect the two endpoint conformational states. The R-FEP-R 2.0 method was applied to estimate the conformational free energy change of a β-turn flip in the protein ubiquitin. The result obtained by R-FEP-R 2.0 agrees with the benchmarks very well. We also comment on problems commonly encountered when applying umbrella sampling to calculate protein conformational free energy changes.

Scaled recoupling of chemical shift anisotropies at high magnetic fields under MAS with interspersed C-elements

The power of chemical shift anisotropy (CSA) measurements for probing structure and dynamics of molecules has been long recognized. NMR pulse sequences that allow measurement of CSA values in an indirect dimension of a protein correlation spectrum have been employed for aliphatic groups, but for practical reasons, carbonyl functional groups have been little studied, despite the fact that carbonyls are expected to give particularly varied and informative CSA values. Specifically, the wide spectral widths of carbonyl tensors make their measurements difficult with typically attainable spectrometer settings. We present here an extended family of experiments that enable the recovery of static CSA lineshapes in an indirect dimension of magic angle spinning (MAS) solid-state NMR experiments, except for various real valued scaling factors. The experiment is suitable for uniformly labeled material, at moderate MAS rates (10 kHz-30 kHz) and at higher magnetic fields (ν0H > 600 MHz). Specifically, the experiments are based on pulse sequence elements from a previous commonly used pulse sequence for CSA measurement, recoupling of chemical shift anisotropy (ROCSA), while modification of scaling factors is achieved by interspersing different blocks of C-elements of the same Cnn 1 cycle. Using experimental conditions similar to the parent ROCSA sequence, a CSA scaling factor between 0 and 0.272 can be obtained, thus allowing a useful practical range of possibilities in experimental conditions for measurement of larger CSA values. Using these blocks, it is also possible to make a constant-time CSA recoupling sequence. The effectiveness of this approach, fROCSA, is shown on model compounds 1-13C-Gly, U-13C,15N-l-His, and microcrystalline U-13C,15N-Ubiquitin.

Accuracy of 1H-1H distances measured using frequency selective recoupling and fast magic-angle spinning

Selective recoupling of protons (SERP) is a method to selectively and quantitatively measure magnetic dipole-dipole interaction between protons and, in turn, the proton-proton distance in solid-state samples at fast magic-angle spinning. We present a bimodal operator-based Floquet approach to describe the numerically optimized SERP recoupling sequence. The description calculates the allowed terms in the first-order effective Hamiltonian, explains the origin of selectivity during recoupling, and shows how different terms are modulated as a function of the radio frequency amplitude and the phase of the sequence. Analytical and numerical simulations have been used to evaluate the effect of higher-order terms and offsets on the polarization transfer efficiency and quantitative distance measurement. The experimentally measured ¹H-¹H distances on a fully protonated thymol sample are ∼10%-15% shorter than those reported from diffraction studies. A semi-quantitative model combined with extensive numerical simulations is used to rationalize the effect of the third-spin and the role of different parameters in the experimentally observed shorter distances. Measurements at high magnetic fields improve the match between experimental and diffraction distances. The measurement of ¹H-¹H couplings at offsets different from the SERP-offset has also been explored. Experiments were also performed on a perdeuterated ubiquitin sample to demonstrate the feasibility of simultaneously measuring multiple quantitative distances and to evaluate the accuracy of the measured distance in the absence of multispin effects. The estimation of proton-proton distances provides a boost to structural characterization of small pharmaceuticals and biomolecules, given that the positions of protons are generally not well defined in x-ray structures.

Quantifying proton NMR coherent linewidth in proteins under fast MAS conditions: a second moment approach

Proton detected solid-state NMR under fast magic-angle-spinning (MAS) conditions is currently redefining the applications of solid-state NMR, in particular in structural biology. Understanding the contributions to the spectral linewidth is thereby of paramount importance. When disregarding the sample-dependent inhomogeneous contributions, the NMR proton linewidth is defined by homogeneous broadening, which has incoherent and coherent contributions. Understanding and disentangling these different contributions in multi-spin systems like proteins is still an open issue. The coherent contribution is mainly caused by the dipolar interaction under MAS and is determined by the molecular structure and the proton chemical shifts. Numerical simulation approaches based on numerically exact direct integration of the Liouville-von Neumann equation can give valuable information about the lineshape, but are limited to small spin systems (<12 spins). We present an alternative simulation method for the coherent contributions based on the rapid and partially analytic calculation of the second moments of large spin systems. We first validate the method on a simple system by predicting the ¹⁹F linewidth in CaF₂ under MAS. We compare simulation results to experimental data for microcrystalline ubiquitin (deuterated 100% back-exchanged at 110 kHz and fully-protonated at 125 kHz). Our results quantitatively explain the observed linewidth per-residue basis for the vast majority of residues.

Protein conformational dynamics studied by 15N and 1H R1ρ relaxation dispersion: Application to wild-type and G53A ubiquitin crystals

Solid-state NMR spectroscopy can provide site-resolved information about protein dynamics over many time scales. Here we combine protein deuteration, fast magic-angle spinning (~45-60kHz) and proton detection to study dynamics of ubiquitin in microcrystals, and in particular a mutant in a region that undergoes microsecond motions in a β-turn region in the wild-type protein. We use ¹⁵N R_1ρ relaxation measurements as a function of the radio-frequency (RF) field strength, i.e. relaxation dispersion, to probe how the G53A mutation alters these dynamics. We report a population-inversion of conformational states: the conformation that in the wild-type protein is populated only sparsely becomes the predominant state. We furthermore explore the potential to use amide-¹H R_1ρ relaxation to obtain insight into dynamics. We show that while quantitative interpretation of ¹H relaxation remains beyond reach under the experimental conditions, due to coherent contributions to decay, one may extract qualitative information about flexibility.

Slow conformational exchange and overall rocking motion in ubiquitin protein crystals

Proteins perform their functions in solution but their structures are most frequently studied inside crystals. Here we probe how the crystal packing alters microsecond dynamics, using solid-state NMR measurements and multi-microsecond MD simulations of different crystal forms of ubiquitin. In particular, near-rotary-resonance relaxation dispersion (NERRD) experiments probe angular backbone motion, while Bloch–McConnell relaxation dispersion data report on fluctuations of the local electronic environment. These experiments and simulations reveal that the packing of the protein can significantly alter the thermodynamics and kinetics of local conformational exchange. Moreover, we report small-amplitude reorientational motion of protein molecules in the crystal lattice with an ~3–5° amplitude on a tens-of-microseconds time scale in one of the crystals, but not in others. An intriguing possibility arises that overall motion is to some extent coupled to local dynamics. Our study highlights the importance of considering the packing when analyzing dynamics of crystalline proteins.

X-ray crystallography is the main method for protein structure determination. Here the authors combine solid-state NMR measurements and molecular dynamics simulations and show that crystal packing alters the thermodynamics and kinetics of local conformational exchange as well as overall rocking motion of protein molecules in the crystal lattice.

Flexibility and Design: Conformational Heterogeneity along the Evolutionary Trajectory of a Redesigned Ubiquitin

Although protein design has been used to introduce new functions, designed variants generally only function as well as natural proteins after rounds of laboratory evolution. One possibility for this pattern is that designed mutants frequently sample nonfunctional conformations. To test this idea, we exploited advances in multiconformer modeling of room-temperature X-ray data collection on redesigned ubiquitin variants selected for increasing binding affinity to the deubiquitinase USP7. Initial core mutations disrupt natural packing and lead to increased flexibility. Additional, experimentally selected mutations quenched conformational heterogeneity through new stabilizing interactions. Stabilizing interactions, such as cation-pi stacking and ordered waters, which are not included in standard protein design energy functions, can create specific interactions that have long-range effects on flexibility across the protein. Our results suggest that increasing flexibility may be a useful strategy to escape local minima during initial directed evolution and protein design steps when creating new functions.

Accelerating proton spin diffusion in perdeuterated proteins at 100 kHz MAS

Fast magic-angle spinning (>60 kHz) has many advantages but makes spin-diffusion-type proton-proton long-range polarization transfer inefficient and highly dependent on chemical-shift offset. Using 100%-HN-[²H,¹³C,¹⁵N]-ubiquitin as a model substance, we quantify the influence of the chemical-shift difference on the spin diffusion between proton spins and compare two experiments which lead to an improved chemical-shift compensation of the transfer: rotating-frame spin diffusion and a new experiment, reverse amplitude-modulated MIRROR. Both approaches enable broadband spin diffusion, but the application of the first variant is limited due to fast spin relaxation in the rotating frame. The reverse MIRROR experiment, in contrast, is a promising candidate for the determination of structurally relevant distance restraints. The applied tailored rf-irradiation schemes allow full control over the range of recoupled chemical shifts and efficiently drive spin diffusion. Here, the relevant relaxation time is the larger longitudinal relaxation time, which leads to a higher signal-to-noise ratio in the spectra.

Allosteric switch regulates protein-protein binding through collective motion

Many biological processes depend on allosteric communication between different parts of a protein, but the role of internal protein motion in propagating signals through the structure remains largely unknown. Through an experimental and computational analysis of the ground state dynamics in ubiquitin, we identify a collective global motion that is specifically linked to a conformational switch distant from the binding interface. This allosteric coupling is also present in crystal structures and is found to facilitate multispecificity, particularly binding to the ubiquitin-specific protease (USP) family of deubiquitinases. The collective motion that enables this allosteric communication does not affect binding through localized changes but, instead, depends on expansion and contraction of the entire protein domain. The characterization of these collective motions represents a promising avenue for finding and manipulating allosteric networks.

Protein stabilizer, NDSB-195, enhances the dynamics of the β4 -α2 loop of ubiquitin

Non-detergent sulfobetaines (NDSBs) are a new group of small, synthetic protein stabilizers, which have advantages over classical compatible osmolytes, such as polyol, amines, and amino acids: they do not increase solution viscosity, unlike polyols, and they are zwitterionic at all pH ranges, unlike amines and amino acids. NDSBs also facilitate the crystallization and refolding of proteins. The mechanism whereby NDSBs exhibit such activities, however, remains elusive. To gain insight into this mechanism, we studied, using nuclear magnetic resonance (NMR), the effects of dimethylethylammonium propane sulfonate (NDSB-195) on the dynamics of ubiquitin, on which a wealth of information has been accumulated. By analyzing the line width of amide proton resonances and the transverse relaxation rates of nitrogen atoms, we found that NDSB-195 enhances the microsecond-millisecond dynamics of a β4 -α2 loop of ubiquitin. Although those compounds that enhance protein dynamics are generally considered to destabilize protein molecules, NDSB-195 enhanced the stability of ubiquitin against guanidium chloride denaturation. Thus, the simultaneous enhancement of stability and flexibility by a single compound can be attained.

Ensembles of a small number of conformations with relative populations

In our previous work, we proposed a new way to represent protein native states, using ensembles of a small number of conformations with relative Populations, or ESP in short. Using Ubiquitin as an example, we showed that using a small number of conformations could greatly reduce the potential of overfitting and assigning relative populations to protein ensembles could significantly improve their quality. To demonstrate that ESP indeed is an excellent alternative to represent protein native states, in this work we compare the quality of two ESP ensembles of Ubiquitin with several well-known regular ensembles or average structure representations. Extensive amount of significant experimental data are employed to achieve a thorough assessment. Our results demonstrate that ESP ensembles, though much smaller in size comparing to regular ensembles, perform equally or even better sometimes in all four different types of experimental data used in the assessment, namely, the residual dipolar couplings, residual chemical shift anisotropy, hydrogen exchange rates, and solution scattering profiles. This work further underlines the significance of having relative populations in describing the native states.

Observing the overall rocking motion of a protein in a crystal

The large majority of three-dimensional structures of biological macromolecules have been determined by X-ray diffraction of crystalline samples. High-resolution structure determination crucially depends on the homogeneity of the protein crystal. Overall ‘rocking' motion of molecules in the crystal is expected to influence diffraction quality, and such motion may therefore affect the process of solving crystal structures. Yet, so far overall molecular motion has not directly been observed in protein crystals, and the timescale of such dynamics remains unclear. Here we use solid-state NMR, X-ray diffraction methods and μs-long molecular dynamics simulations to directly characterize the rigid-body motion of a protein in different crystal forms. For ubiquitin crystals investigated in this study we determine the range of possible correlation times of rocking motion, 0.1–100 μs. The amplitude of rocking varies from one crystal form to another and is correlated with the resolution obtainable in X-ray diffraction experiments.

Small-amplitude overall motion of molecules in crystals limits the achievable resolution in X-ray diffraction, yet little is known about its exact nature. Here, the authors obtain NMR, XRD and MD data from three different crystal forms of a protein (ubiquitin) to gain insight into amplitude and timescale of such motions.

Structural heterogeneity in microcrystalline ubiquitin studied by solid-state NMR

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.

Ensemble MD simulations restrained via crystallographic data: accurate structure leads to accurate dynamics

Currently, the best existing molecular dynamics (MD) force fields cannot accurately reproduce the global free-energy minimum which realizes the experimental protein structure. As a result, long MD trajectories tend to drift away from the starting coordinates (e.g., crystallographic structures). To address this problem, we have devised a new simulation strategy aimed at protein crystals. An MD simulation of protein crystal is essentially an ensemble simulation involving multiple protein molecules in a crystal unit cell (or a block of unit cells). To ensure that average protein coordinates remain correct during the simulation, we introduced crystallography-based restraints into the MD protocol. Because these restraints are aimed at the ensemble-average structure, they have only minimal impact on conformational dynamics of the individual protein molecules. So long as the average structure remains reasonable, the proteins move in a native-like fashion as dictated by the original force field. To validate this approach, we have used the data from solid-state NMR spectroscopy, which is the orthogonal experimental technique uniquely sensitive to protein local dynamics. The new method has been tested on the well-established model protein, ubiquitin. The ensemble-restrained MD simulations produced lower crystallographic R factors than conventional simulations; they also led to more accurate predictions for crystallographic temperature factors, solid-state chemical shifts, and backbone order parameters. The predictions for (15) N R1 relaxation rates are at least as accurate as those obtained from conventional simulations. Taken together, these results suggest that the presented trajectories may be among the most realistic protein MD simulations ever reported. In this context, the ensemble restraints based on high-resolution crystallographic data can be viewed as protein-specific empirical corrections to the standard force fields.

Enhancing the quality of protein conformation ensembles with relative populations

The function and dynamics of many proteins are best understood not from a single structure but from an ensemble. A high quality ensemble is necessary for accurately delineating protein dynamics. However, conformations in an ensemble are generally given equal weights. Few attempts were made to assign relative populations to the conformations, mainly due to the lack of right experimental data. Here we propose a method for assigning relative populations to ensembles using experimental residue dipolar couplings (RDC) as constraints, and show that relative populations can significantly enhance an ensemble's ability in representing the native states and dynamics. The method works by identifying conformation states within an ensemble and assigning appropriate relative populations to them. Each of these conformation states is represented by a sub-ensemble consisting of a subset of the conformations. Application to the ubiquitin X-ray ensemble clearly identifies two key conformation states, with relative populations in excellent agreement with previous work. We then apply the method to a reprotonated ERNST ensemble that is enhanced with a switched conformation, and show that as a result of population reweighting, not only the reproduction of RDCs is significantly improved, but common conformational features (particularly the dihedral angle distributions of ϕ 53 and ψ 52) also emerge for both the X-ray ensemble and the reprotonated ERNST ensemble.

Solubility and supersaturation-dependent protein misfolding revealed by ultrasonication

Although alcohols are useful cosolvents for producing amyloid fibrils, the underlying mechanism of alcohol-dependent fibrillation is unclear. We studied the alcohol-induced fibrillation of hen egg-white lysozyme at various concentrations of ethanol, 2,2,2-trifluoroethanol (TFE), and 1,1,1,3,3,3-hexafluoro-2-propanol (HFIP). Under the conditions where the alcohol-denatured lysozyme retained metastability, ultrasonication effectively triggered fibrillation. The optimal alcohol concentration depended on the alcohol species. HFIP showed a sharp maximum at 12-16%. For TFE, a broad maximum at 40-80% was observed. Ethanol exhibited only an increase in fibrillation above 60%. These profiles were opposite to the equilibrium solubility of lysozyme in water/alcohol mixtures. The results indicate that although fibrillation is determined by solubility, supersaturation prevents conformational transitions and ultrasonication is highly effective in minimizing an effect of supersaturation. We propose an alcohol-dependent protein misfolding funnel useful for examining amyloidogenicity. This misfolding funnel will apply to fibrillation under physiological conditions where biological environments play important roles in decreasing the solubility.

Practical considerations over spectral quality in solid state NMR spectroscopy of soluble proteins

Great theoretical and methodological advances are pushing the limits of resolution and sensitivity in solid state NMR (SSNMR). However, sample preparation remains a critical issue for the success of an experiment. The factors affecting spectral quality in SSNMR samples are discussed, examining cases encountered in the literature and presenting new experimental data. A discussion on resolution and sensitivity in sedimented solutes is framed in this context.

Motional timescale predictions by molecular dynamics simulations: case study using proline and hydroxyproline sidechain dynamics

We propose a new approach for force field optimizations which aims at reproducing dynamics characteristics using biomolecular MD simulations, in addition to improved prediction of motionally averaged structural properties available from experiment. As the source of experimental data for dynamics fittings, we use ¹³C NMR spin-lattice relaxation times T₁ of backbone and sidechain carbons, which allow to determine correlation times of both overall molecular and intramolecular motions. For structural fittings, we use motionally averaged experimental values of NMR J couplings. The proline residue and its derivative 4-hydroxyproline with relatively simple cyclic structure and sidechain dynamics were chosen for the assessment of the new approach in this work. Initially, grid search and simplexed MD simulations identified large number of parameter sets which fit equally well experimental J couplings. Using the Arrhenius-type relationship between the force constant and the correlation time, the available MD data for a series of parameter sets were analyzed to predict the value of the force constant that best reproduces experimental timescale of the sidechain dynamics. Verification of the new force-field (termed as AMBER99SB-ILDNP) against NMR J couplings and correlation times showed consistent and significant improvements compared to the original force field in reproducing both structural and dynamics properties. The results suggest that matching experimental timescales of motions together with motionally averaged characteristics is the valid approach for force field parameter optimization. Such a comprehensive approach is not restricted to cyclic residues and can be extended to other amino acid residues, as well as to the backbone. Proteins 2014; 82:195–215. © 2013 Wiley Periodicals, Inc.

Protein backbone and sidechain torsion angles predicted from NMR chemical shifts using artificial neural networks

A new program, TALOS-N, is introduced for predicting protein backbone torsion angles from NMR chemical shifts. The program relies far more extensively on the use of trained artificial neural networks than its predecessor, TALOS+. Validation on an independent set of proteins indicates that backbone torsion angles can be predicted for a larger, ≥90 % fraction of the residues, with an error rate smaller than ca 3.5 %, using an acceptance criterion that is nearly two-fold tighter than that used previously, and a root mean square difference between predicted and crystallographically observed (ϕ, ψ) torsion angles of ca 12º. TALOS-N also reports sidechain χ(1) rotameric states for about 50 % of the residues, and a consistency with reference structures of 89 %. The program includes a neural network trained to identify secondary structure from residue sequence and chemical shifts.

Atomic-level description of ubiquitin folding

Equilibrium molecular dynamics simulations, in which proteins spontaneously and repeatedly fold and unfold, have recently been used to help elucidate the mechanistic principles that underlie the folding of fast-folding proteins. The extent to which the conclusions drawn from the analysis of such proteins, which fold on the microsecond timescale, apply to the millisecond or slower folding of naturally occurring proteins is, however, unclear. As a first attempt to address this outstanding issue, we examine here the folding of ubiquitin, a 76-residue-long protein found in all eukaryotes that is known experimentally to fold on a millisecond timescale. Ubiquitin folding has been the subject of many experimental studies, but its slow folding rate has made it difficult to observe and characterize the folding process through all-atom molecular dynamics simulations. Here we determine the mechanism, thermodynamics, and kinetics of ubiquitin folding through equilibrium atomistic simulations. The picture emerging from the simulations is in agreement with a view of ubiquitin folding suggested from previous experiments. Our findings related to the folding of ubiquitin are also consistent, for the most part, with the folding principles derived from the simulation of fast-folding proteins, suggesting that these principles may be applicable to a wider range of proteins.

Conformational and dynamic changes at the interface contribute to ligand binding by ubiquitin

Ubiquitin interacts with numerous domains and motifs in its lifetime that vary in structure but bind the same hydrophobic patch. To identify the structural features of ubiquitin that make it an exceptional protein-protein interaction partner, we have studied the interaction of ubiquitin with the signal transducing adaptor molecule-1 ubiquitin interacting motif (UIM) using nuclear magnetic resonance. Our studies bring to light the role of the inherent backbone flexibility of ubiquitin in its interactions with a large array of binding partners, revealed from the changes in C(α) chemical shifts, backbone dynamics, and hydrogen bond lengths upon UIM binding. The crystal structures of ubiquitin complexes lend further support to our findings, underscoring the importance of the unique and flexible hydrogen bond network within ubiquitin and simultaneously providing insights into the nature of the slow motions. Taken together, our studies provide an in-depth view of the molecular changes associated with ligand recognition by ubiquitin.

Site-resolved measurement of microsecond-to-millisecond conformational-exchange processes in proteins by solid-state NMR spectroscopy

We demonstrate that conformational exchange processes in proteins on microsecond-to-millisecond time scales can be detected and quantified by solid-state NMR spectroscopy. We show two independent approaches that measure the effect of conformational exchange on transverse relaxation parameters, namely Carr-Purcell-Meiboom-Gill relaxation-dispersion experiments and measurement of differential multiple-quantum coherence decay. Long coherence lifetimes, as required for these experiments, are achieved by the use of highly deuterated samples and fast magic-angle spinning. The usefulness of the approaches is demonstrated by application to microcrystalline ubiquitin. We detect a conformational exchange process in a region of the protein for which dynamics have also been observed in solution. Interestingly, quantitative analysis of the data reveals that the exchange process is more than 1 order of magnitude slower than in solution, and this points to the impact of the crystalline environment on free energy barriers.

A hydrogen bond regulates slow motions in ubiquitin by modulating a β-turn flip

Proteins exist as conformational ensembles composed of multiple interchanging substates separated by kinetic barriers. Interconverting conformations are often difficult to probe, owing to their sparse population and transient nature. Here, we report the identification and characterization of a subset of conformations in ubiquitin that participate in microsecond-to-millisecond motions in the amides of Ile23, Asn25, and Thr55. A novel side chain to the backbone hydrogen bond that regulates these motions has also been identified. Combining our NMR studies with the available X-ray data, we have unearthed the physical process underlying slow motions-the interconversion of a type I into a type II β-turn flip at residues Glu51 through Arg54. Interestingly, the dominant conformer of wild-type ubiquitin observed in solution near neutral pH is only represented by about 22% of the crystal structures. The conformers generated as a result of the dynamics of the hydrogen bond appear to be correlated to ligand recognition by ubiquitin.