Concordance of Residual Dipolar Couplings, Backbone Order Parameters and Crystallographic B-factors for a Small α/β Protein: A Unified Picture of High Probability, Fast Atomic Motions in Proteins

On the capacity and superposition of minima in neural network loss function landscapes

Minima of the loss function landscape (LFL) of a neural network are locally optimal sets of weights that extract and process information from the input data to make outcome predictions. In underparameterised networks, the capacity of the weights may be insufficient to fit all the relevant information. We demonstrate that different local minima specialise in certain aspects of the learning problem, and process the input information differently. This effect can be exploited using a meta-network in which the predictive power from multiple minima of the LFL is combined to produce a better classifier. With this approach, we can increase the area under the receiver operating characteristic curve by around 20 % for a complex learning problem. We propose a theoretical basis for combining minima and show how a meta-network can be trained to select the representative that is used for classification of a specific data item. Finally, we present an analysis of symmetry-equivalent solutions to machine learning problems, which provides a systematic means to improve the efficiency of this approach.

Flexible Loop in Carbohydrate-Binding Module 48 Allosterically Modulates Substrate Binding of the 1,4-α-Glucan Branching Enzyme

The 1,4-α-glucan branching enzyme (GBE, EC 2.4.1.18) catalyzes the formation of α-1,6 branching points in starch and plays a key role in synthesis. To obtain mechanistic insights into the catalytic action of the enzyme, we first determined the crystal structure of GBE from Rhodothermus obamensis STB05 (RoGBE) to a resolution of 2.39 Å (PDB ID: 6JOY). The structure consists of three domains: domain A, domain C, and the carbohydrate-binding module 48 (CBM48). An engineered truncated mutant lacking the CBM48 domain (ΔCBM48) showed significantly reduced ligand binding affinity and enzyme activity. Comparison of the structures of RoGBE with other GBEs showed that CBM48 of RoGBE had a longer flexible loop. Truncation of the flexible loops resulted in reduced binding affinity and activity, thereby substantiating the importance of the optimum loop structure for catalysis. In essence, our study shows that CBM48, especially the flexible loop, plays an important role in substrate binding and enzymatic activity of RoGBE. Further, based on the structural analysis, kinetics, and activity assays on wild type and mutants, as well as homology modeling, we proposed a mechanistic model (called the "lid model") to illustrate how the flexible loop triggers substrate binding, ultimately leading to catalysis.

BANΔIT: B'-Factor Analysis for Drug Design and Structural Biology

The analysis of B‐factor profiles from X‐ray protein structures can be utilized for structure‐based drug design since protein mobility changes have been associated with the quality of protein‐ligand interactions. With the BANΔIT (B’‐factor analysis and ΔB’ interpretation toolkit), we have developed a JavaScript‐based browser application that provides a graphical user interface for the normalization and analysis of B’‐factor profiles. To emphasize the usability for rational drug design applications, we have analyzed a selection of crystallographic protein‐ligand complexes and have given exemplary conclusions for further drug optimization including the development of a B’‐factor‐supported pharmacophore model for SARS CoV‐2 main protease inhibitors. BANΔIT is available online at https://bandit.uni‐mainz.de. The source code can be downloaded from https://github.com/FBarthels/BANDIT.

Remarkable Rigidity of the Single α-Helical Domain of Myosin-VI As Revealed by NMR Spectroscopy

Although the α-helix has long been recognized as an all-important element of secondary structure, it generally requires stabilization by tertiary interactions with other parts of a protein's structure. Highly charged single α-helical (SAH) domains, consisting of a high percentage (>75%) of Arg, Lys, and Glu residues, are exceptions to this rule but have been difficult to characterize structurally. Our study focuses on the 68-residue medial tail domain of myosin-VI, which is found to contain a highly ordered α-helical structure extending from Glu-6 to Lys-63. High hydrogen exchange protection factors (15-150), small (ca. 4 Hz) 3 JHNHα couplings, and a near-perfect fit to an ideal model α-helix for its residual dipolar couplings (RDCs), measured in a filamentous phage medium, support the high regularity of this helix. Remarkably, the hydrogen exchange rates are far more homogeneous than the protection factors derived from them, suggesting that for these transiently broken helices the intrinsic exchange rates derived from the amino acid sequence are not appropriate reference values. 15N relaxation data indicate a very high degree of rotational diffusion anisotropy ( D∥/ D⊥ ≈ 7.6), consistent with the hydrodynamic behavior predicted for such a long, nearly straight α-helix. Alignment of the helix by a paramagnetic lanthanide ion attached to its N-terminal region shows a decrease in alignment as the distance from the tagging site increases. This decrease yields a precise measure for the persistence length of 224 ± 10 Å at 20 °C, supporting the idea that the role of the SAH helix is to act as an extension of the myosin-VI lever arm.

Mapping Local Conformational Landscapes of Proteins in Solution

The ability of proteins to adopt multiple conformational states is essential to their function, and elucidating the details of such diversity under physiological conditions has been a major challenge. Here we present a generalized method for mapping protein population landscapes by NMR spectroscopy. Experimental NOESY spectra are directly compared with a set of expectation spectra back-calculated across an arbitrary conformational space. Signal decomposition of the experimental spectrum then directly yields the relative populations of local conformational microstates. In this way, averaged descriptions of conformation can be eliminated. As the method quantitatively compares experimental and expectation spectra, it inherently delivers an R factor expressing how well structural models explain the input data. We demonstrate that our method extracts sufficient information from a single 3D NOESY experiment to perform initial model building, refinement, and validation, thus offering a complete de novo structure determination protocol.

Distance-independent Cross-correlated Relaxation and Isotropic Chemical Shift Modulation in Protein Dynamics Studies

Cross-correlated relaxation (CCR) in multiple-quantum coherences differs from other relaxation phenomena in its theoretical ability to be mediated across an infinite distance. The two interfering relaxation mechanisms may be dipolar interactions, chemical shift anisotropies, chemical shift modulations or quadrupolar interactions. These properties make multiple-quantum CCR an attractive probe for structure and dynamics of biomacromolecules not accessible from other measurements. Here, we review the use of multiple-quantum CCR measurements in dynamics studies of proteins. We compile a list of all experiments proposed for CCR rate measurements, provide an overview of the theory with a focus on protein dynamics, and present applications to various protein systems.

The Structure of Melanoregulin Reveals a Role for Cholesterol Recognition in the Protein's Ability to Promote Dynein Function

Melanoregulin (Mreg) is a small, highly charged, multiply palmitoylated protein present on the membrane of melanosomes. Mreg is implicated in the transfer of melanosomes from melanocytes to keratinocytes, and in promoting the microtubule minus end-directed transport of these organelles. The possible molecular function of Mreg was identified by solving its structure using nuclear magnetic resonance (NMR) spectroscopy. Mreg contains six α helices forming a fishhook-like fold in which positive and negative charges occupy opposite sides of the protein's surface and sandwich a putative, cholesterol recognition sequence (CRAC motif). Mreg containing a point mutation within its CRAC motif still targets to late endosomes/lysosomes, but no longer promotes their microtubule minus end-directed transport. Moreover, wild-type Mreg does not promote the microtubule minus end-directed transport of late endosomes/lysosomes in cells transiently depleted of cholesterol. Finally, reversing the charge of three clustered acidic residues partially inhibits Mreg's ability to drive these organelles to microtubule minus ends.

Detection of Correlated Protein Backbone and Side-Chain Angle Fluctuations

NMR methods for the characterization of local protein motions have attained a high level of sophistication. Measurement of the synchronization between those motions, however, poses a serious challenge. Such correlated motions are one of the underlying mechanisms for the propagation of local changes to remote sites and as such for information transfer. Here, we demonstrate the experimental detection of the synchronization of motion over an intermediate range. To that purpose, we designed pulse sequences for the measurement of cross-correlated relaxation between the backbone H^N -N and side-chain H^β -C^β dipoles in Ile, Thr, and Val in the protein GB3. These bonds are related through two and three intervening dihedral angles. We show that the correlated motions inherent in a structural ensemble obtained from a large and diverse array of NMR probes are in excellent agreement with our measurements.

Cross-correlated relaxation rates between protein backbone H-X dipolar interactions

The relaxation interference between dipole-dipole interactions of two separate spin pairs carries structural and dynamics information. In particular, when compared to individual dynamic behavior of those spin pairs, such cross-correlated relaxation (CCR) rates report on the correlation between the spin pairs. We have recently mapped out correlated motion along the backbone of the protein GB3, using CCR rates among and between consecutive H^N-N and H^α-C^α dipole-dipole interactions. Here, we provide a detailed account of the measurement of the four types of CCR rates. All rates were obtained from at least two different pulse sequences, of which the yet unpublished ones are presented. Detailed comparisons between the different methods and corrections for unwanted pathways demonstrate that the averaged CCR rates are highly accurate and precise with errors of 1.5-3% of the entire value ranges.

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.

Protein dynamics and function from solution state NMR spectroscopy

It is well-established that dynamics are central to protein function; their importance is implicitly acknowledged in the principles of the Monod, Wyman and Changeux model of binding cooperativity, which was originally proposed in 1965. Nowadays the concept of protein dynamics is formulated in terms of the energy landscape theory, which can be used to understand protein folding and conformational changes in proteins. Because protein dynamics are so important, a key to understanding protein function at the molecular level is to design experiments that allow their quantitative analysis. Nuclear magnetic resonance (NMR) spectroscopy is uniquely suited for this purpose because major advances in theory, hardware, and experimental methods have made it possible to characterize protein dynamics at an unprecedented level of detail. Unique features of NMR include the ability to quantify dynamics (i) under equilibrium conditions without external perturbations, (ii) using many probes simultaneously, and (iii) over large time intervals. Here we review NMR techniques for quantifying protein dynamics on fast (ps-ns), slow (μs-ms), and very slow (s-min) time scales. These techniques are discussed with reference to some major discoveries in protein science that have been made possible by NMR spectroscopy.

The Exact NOE as an Alternative in Ensemble Structure Determination

The structure-function paradigm is increasingly replaced by the structure-dynamics-function paradigm. All protein activity is steered by the interplay between enthalpy and entropy. Conformational dynamics serves as a proxy of conformational entropy. Therefore, it is essential to study not only the average conformation but also the spatial sampling of a protein on all timescales. To this purpose, we have established a protocol for determining multiple-state ensembles of proteins based on exact nuclear Overhauser effects (eNOEs). We have recently extended our previously reported eNOE data set for the protein GB3 by a very large set of backbone and side-chain residual dipolar couplings and three-bond J couplings. Here, we demonstrate that at least four structural states are required to represent the complete data set by dissecting the contributions to the CYANA target function, which quantifies restraint violations in structure calculation. We present a four-state ensemble of GB3, which largely preserves the characteristics obtained from eNOEs only. Due to the abundance of the input data, the ensemble and χ(1) angles in particular are well suited for cross-validation of the input data and comparison to x-ray structures. Principal component analysis is used to automatically identify and validate relevant states of the ensembles. Overall, our findings suggest that eNOEs are a valuable alternative to traditional NMR probes in spatial elucidation of proteins.

Direct Investigation of Slow Correlated Dynamics in Proteins via Dipolar Interactions

The synchronization of native state motions as they transition between microstates influences catalysis kinetics, mediates allosteric interactions, and reduces the conformational entropy of proteins. However, it has proven difficult to describe native microstates because they are usually minimally frustrated and may interconvert on the micro- to millisecond time scale. Direct observation of concerted equilibrium fluctuations would therefore be an important tool for describing protein native states. Here we propose a strategy that relates NMR cross-correlated relaxation (CCR) rates between dipolar interactions to residual dipolar couplings (RDCs) of individual consecutive H(N)-N and H(α)-C(α) bonds, which act as a proxy for the peptide planes and the side chains, respectively. Using Xplor-NIH ensemble structure calculations restrained with the RDC and CCR data, we observe collective motions on time scales slower than nanoseconds in the backbone for GB3. To directly access the correlations from CCR, we develop a structure-free data analysis. The resulting dynamic correlation map is consistent with the ensemble-restrained simulations and reveals a complex network. In general, we find that the bond motions are on average slightly correlated and that the local environment dominates many observations. Despite this, some patterns are typical over entire secondary structure elements. In the β-sheet, nearly all bonds are weakly correlated, and there is an approximately binary alternation in correlation intensity corresponding to the solvent exposure/shielding alternation of the side chains. For α-helices, there is also a weak correlation in the H(N)-N bonds. The degree of correlation involving H(α)-C(α) bonds is directly affected by side-chain fluctuations, whereas loops show complex and nonuniform behavior.

Hybrid Approaches to Structural Characterization of Conformational Ensembles of Complex Macromolecular Systems Combining NMR Residual Dipolar Couplings and Solution X-ray Scattering

Solving structures or structural ensembles of large macromolecular systems in solution poses a challenging problem. While NMR provides structural information at atomic resolution, increased spectral complexity, chemical shift overlap, and short transverse relaxation times (associated with slow tumbling) render application of the usual techniques that have been so successful for medium sized systems (<50 kDa) difficult. Solution X-ray scattering, on the other hand, is not limited by molecular weight but only provides low resolution structural information related to the overall shape and size of the system under investigation. Here we review how combining atomic resolution structures of smaller domains with sparse experimental data afforded by NMR residual dipolar couplings (which yield both orientational and shape information) and solution X-ray scattering data in rigid-body simulated annealing calculations provides a powerful approach for investigating the structural aspects of conformational dynamics in large multidomain proteins. The application of this hybrid methodology is illustrated for the 128 kDa dimer of bacterial Enzyme I which exists in a variety of open and closed states that are sampled at various points in the catalytic cycles, and for the capsid protein of the human immunodeficiency virus.

NMR Order Parameter Determination from Long Molecular Dynamics Trajectories for Objective Comparison with Experiment

Functional protein motions covering a wide range of time scales can be studied, among other techniques, by NMR and by molecular dynamics (MD) computer simulations. MD simulations of proteins now routinely extend into the hundreds of nanoseconds time scale range exceeding the overall tumbling correlation times of proteins in solution by several orders of magnitude. This provides a unique opportunity to rigorously validate these simulations by quantitative comparison with model-free order parameters derived from NMR relaxation experiments. However, presently there is no consensus on how such a comparison is best done. We address here how this can be accomplished in a way that is both efficient and objective. For this purpose, we analyze (15)N R1 and R2 and heteronuclear {(1)H}-(15)N NOE NMR relaxation parameters computed from 500 ns MD trajectories of 10 different protein systems using the model-free analysis. The resulting model-free S(2) order parameters are then used as targets for S(2) values computed directly from the trajectories by the iRED method by either averaging over blocks of variable lengths or by using exponentially weighted snapshots (wiRED). We find that the iRED results are capable of reproducing the target S(2) values with high accuracy provided that the averaging window is chosen 5 times the length of the overall tumbling correlation time. These results provide useful guidelines for the derivation of NMR order parameters from MD for a meaningful comparison with their experimental counterparts.

Compiled data set of exact NOE distance limits, residual dipolar couplings and scalar couplings for the protein GB3

We compiled an NMR data set consisting of exact nuclear Overhauser enhancement (eNOE) distance limits, residual dipolar couplings (RDCs) and scalar (J) couplings for GB3, which forms one of the largest and most diverse data set for structural characterization of a protein to date. All data have small experimental errors, which are carefully estimated. We use the data in the research article Vogeli et al., 2015, Complementarity and congruence between exact NOEs and traditional NMR probes for spatial decoding of protein dynamics, J. Struct. Biol., 191, 3, 306–317, doi:10.1016/j.jsb.2015.07.008 [1] for cross-validation in multiple-state structural ensemble calculation. We advocate this set to be an ideal test case for molecular dynamics simulations and structure calculations.

Unraveling the complexity of protein backbone dynamics with combined (13)C and (15)N solid-state NMR relaxation measurements

Typically, protein dynamics involve a complex hierarchy of motions occurring on different time scales between conformations separated by a range of different energy barriers. NMR relaxation can in principle provide a site-specific picture of both the time scales and amplitudes of these motions, but independent relaxation rates sensitive to fluctuations in different time scale ranges are required to obtain a faithful representation of the underlying dynamic complexity. This is especially pertinent for relaxation measurements in the solid state, which report on dynamics in a broader window of time scales by more than 3 orders of magnitudes compared to solution NMR relaxation. To aid in unraveling the intricacies of biomolecular dynamics we introduce (13)C spin-lattice relaxation in the rotating frame (R1ρ) as a probe of backbone nanosecond-microsecond motions in proteins in the solid state. We present measurements of (13)C'R1ρ rates in fully protonated crystalline protein GB1 at 600 and 850 MHz (1)H Larmor frequencies and compare them to (13)C'R1, (15)N R1 and R1ρ measured under the same conditions. The addition of carbon relaxation data to the model free analysis of nitrogen relaxation data leads to greatly improved characterization of time scales of protein backbone motions, minimizing the occurrence of fitting artifacts that may be present when (15)N data is used alone. We also discuss how internal motions characterized by different time scales contribute to (15)N and (13)C relaxation rates in the solid state and solution state, leading to fundamental differences between them, as well as phenomena such as underestimation of picosecond-range motions in the solid state and nanosecond-range motions in solution.

Molecular Dynamics of Biomolecules through Direct Analysis of Dipolar Couplings

Residual dipolar couplings (RDCs) are important probes in structural biology, but their analysis is often complicated by the determination of an alignment tensor or its associated assumptions. We here apply the maximum entropy principle to derive a tensor-free formalism which allows for direct, dynamic analysis of RDCs and holds the classic tensor formalism as a special case. Specifically, the framework enables us to robustly analyze data regardless of whether a clear separation of internal and overall dynamics is possible. Such a separation is often difficult in the core subjects of current structural biology, which include multidomain and intrinsically disordered proteins as well as nucleic acids. We demonstrate the method is tractable and self-consistent and generalizes to data sets comprised of observations from multiple different alignment conditions.

Integrative, dynamic structural biology at atomic resolution--it's about time

Biomolecules adopt a dynamic ensemble of conformations, each with the potential to interact with binding partners or perform the chemical reactions required for a multitude of cellular functions. Recent advances in X-ray crystallography, nuclear magnetic resonance (NMR) spectroscopy and other techniques are helping us realize the dream of seeing--in atomic detail--how different parts of biomolecules shift between functional substates using concerted motions. Integrative structural biology has advanced our understanding of the formation of large macromolecular complexes and how their components interact in assemblies by leveraging data from many low-resolution methods. Here, we review the growing opportunities for integrative, dynamic structural biology at the atomic scale, contending there is increasing synergistic potential between X-ray crystallography, NMR and computer simulations to reveal a structural basis for protein conformational dynamics at high resolution.

Visualizing transient dark states by NMR spectroscopy

Myriad biological processes proceed through states that defy characterization by conventional atomic-resolution structural biological methods. The invisibility of these 'dark' states can arise from their transient nature, low equilibrium population, large molecular weight, and/or heterogeneity. Although they are invisible, these dark states underlie a range of processes, acting as encounter complexes between proteins and as intermediates in protein folding and aggregation. New methods have made these states accessible to high-resolution analysis by nuclear magnetic resonance (NMR) spectroscopy, as long as the dark state is in dynamic equilibrium with an NMR-visible species. These methods - paramagnetic NMR, relaxation dispersion, saturation transfer, lifetime line broadening, and hydrogen exchange - allow the exploration of otherwise invisible states in exchange with a visible species over a range of timescales, each taking advantage of some unique property of the dark state to amplify its effect on a particular NMR observable. In this review, we introduce these methods and explore two specific techniques - paramagnetic relaxation enhancement and dark state exchange saturation transfer - in greater detail.

A relationship between the transient structure in the monomeric state and the aggregation propensities of α-synuclein and β-synuclein

α-Synuclein is an intrinsically disordered protein whose aggregation is implicated in Parkinson's disease. A second member of the synuclein family, β-synuclein, shares significant sequence similarity with α-synuclein but is much more resistant to aggregation. β-Synuclein is missing an 11-residue stretch in the central non-β-amyloid component region that forms the core of α-synuclein amyloid fibrils, yet insertion of these residues into β-synuclein to produce the βSHC construct does not markedly increase the aggregation propensity. To investigate the structural basis of these different behaviors, quantitative nuclear magnetic resonance data, in the form of paramagnetic relaxation enhancement-derived interatomic distances, are combined with molecular dynamics simulations to generate ensembles of structures representative of the solution states of α-synuclein, β-synuclein, and βSHC. Comparison of these ensembles reveals that the differing aggregation propensities of α-synuclein and β-synuclein are associated with differences in the degree of residual structure in the C-terminus coupled to the shorter separation between the N- and C-termini in β-synuclein and βSHC, making protective intramolecular contacts more likely.

Protein apparent dielectric constant and its temperature dependence from remote chemical shift effects

A NMR protocol is introduced that permits accurate measurement of minute, remote chemical shift perturbations (CSPs), caused by a mutation-induced change in the electric field. Using protein GB3 as a model system, ¹H^N CSPs in K19A and K19E mutants can be fitted to small changes in the electric field at distal sites in the protein using the Buckingham equation, yielding an apparent dielectric constant ε_a of 8.6 ± 0.8 at 298 K. These CSPs, and their derived ε_a value, scale strongly with temperature. For example, CSPs at 313 K are about ∼30% smaller than those at 278 K, corresponding to an effective ε_a value of about 7.3 at 278 K and 10.5 at 313 K. Molecular dynamics simulations in explicit solvent indicate that solvent water makes a significant contribution to ε_a.

Simultaneous use of solution NMR and X-ray data in REFMAC5 for joint refinement/detection of structural differences

The program REFMAC5 from CCP4 was modified to allow the simultaneous use of X-ray crystallographic data and paramagnetic NMR data (pseudocontact shifts and self-orientation residual dipolar couplings) and/or diamagnetic residual dipolar couplings. Incorporation of these long-range NMR restraints in REFMAC5 can reveal differences between solid-state and solution conformations of molecules or, in their absence, can be used together with X-ray crystallographic data for structural refinement. Since NMR and X-ray data are complementary, when a single structure is consistent with both sets of data and still maintains reasonably `ideal' geometries, the reliability of the derived atomic model is expected to increase. The program was tested on five different proteins: the catalytic domain of matrix metalloproteinase 1, GB3, ubiquitin, free calmodulin and calmodulin complexed with a peptide. In some cases the joint refinement produced a single model consistent with both sets of observations, while in other cases it indicated, outside the experimental uncertainty, the presence of different protein conformations in solution and in the solid state.

Towards a true protein movie: a perspective on the potential impact of the ensemble-based structure determination using exact NOEs

Confined by the Boltzmann distribution of the energies of the states, a multitude of structural states are inherent to biomolecules. For a detailed understanding of a protein's function, its entire structural landscape at atomic resolution and insight into the interconversion between all the structural states (i.e. dynamics) are required. Whereas dedicated trickery with NMR relaxation provides aspects of local dynamics, and 3D structure determination by NMR is well established, only recently have several attempts been made to formulate a more comprehensive description of the dynamics and the structural landscape of a protein. Here, a perspective is given on the use of exact NOEs (eNOEs) for the elucidation of structural ensembles of a protein describing the covered conformational space.

The nuclear Overhauser effect from a quantitative perspective

The nuclear Overhauser enhancement or effect (NOE) is the most important measure in liquid-state NMR with macromolecules. Thus, the NOE is the subject of numerous reviews and books. Here, the NOE is revisited in light of our recently introduced measurements of exact nuclear Overhauser enhancements (eNOEs), which enabled the determination of multiple-state 3D protein structures. This review encompasses all relevant facets from the theoretical considerations to the use of eNOEs in multiple-state structure calculation. Important aspects include a detailed presentation of the relaxation theory relevant for the nuclear Overhauser effect, the estimation of the correction for spin diffusion, the experimental determination of the eNOEs, the conversion of eNOE rates into distances and validation of their quality, the distance-restraint classification and the protocols for calculation of structures and ensembles.

An overview of tools for the validation of protein NMR structures

Biomolecular structures at atomic resolution present a valuable resource for the understanding of biology. NMR spectroscopy accounts for 11% of all structures in the PDB repository. In response to serious problems with the accuracy of some of the NMR-derived structures and in order to facilitate proper analysis of the experimental models, a number of program suites are available. We discuss nine of these tools in this review: PROCHECK-NMR, PSVS, GLM-RMSD, CING, Molprobity, Vivaldi, ResProx, NMR constraints analyzer and QMEAN. We evaluate these programs for their ability to assess the structural quality, restraints and their violations, chemical shifts, peaks and the handling of multi-model NMR ensembles. We document both the input required by the programs and output they generate. To discuss their relative merits we have applied the tools to two representative examples from the PDB: a small, globular monomeric protein (Staphylococcal nuclease from S. aureus, PDB entry 2kq3) and a small, symmetric homodimeric protein (a region of human myosin-X, PDB entry 2lw9).

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.

Integrated description of protein dynamics from room-temperature X-ray crystallography and NMR

Detailed descriptions of atomic coordinates and motions are required for an understanding of protein dynamics and their relation to molecular recognition, catalytic function, and allostery. Historically, NMR relaxation measurements have played a dominant role in the determination of the amplitudes and timescales (picosecond-nanosecond) of bond vector fluctuations, whereas high-resolution X-ray diffraction experiments can reveal the presence of and provide atomic coordinates for multiple, weakly populated substates in the protein conformational ensemble. Here we report a hybrid NMR and X-ray crystallography analysis that provides a more complete dynamic picture and a more quantitative description of the timescale and amplitude of fluctuations in atomic coordinates than is obtainable from the individual methods alone. Order parameters (S(2)) were calculated from single-conformer and multiconformer models fitted to room temperature and cryogenic X-ray diffraction data for dihydrofolate reductase. Backbone and side-chain order parameters derived from NMR relaxation experiments are in excellent agreement with those calculated from the room-temperature single-conformer and multiconformer models, showing that the picosecond timescale motions observed in solution occur also in the crystalline state. These motions are quenched in the crystal at cryogenic temperatures. The combination of NMR and X-ray crystallography in iterative refinement promises to provide an atomic resolution description of the alternate conformational substates that are sampled through picosecond to nanosecond timescale fluctuations of the protein structure. The method also provides insights into the structural heterogeneity of nonmethyl side chains, aromatic residues, and ligands, which are less commonly analyzed by NMR relaxation measurements.

Solution NMR refinement of a metal ion bound protein using metal ion inclusive restrained molecular dynamics methods

Correctly calculating the structure of metal coordination sites in a protein during the process of nuclear magnetic resonance (NMR) structure determination and refinement continues to be a challenging task. In this study, we present an accurate and convenient means by which to include metal ions in the NMR structure determination process using molecular dynamics (MD) simulations constrained by NMR-derived data to obtain a realistic and physically viable description of the metal binding site(s). This method provides the framework to accurately portray the metal ions and its binding residues in a pseudo-bond or dummy-cation like approach, and is validated by quantum mechanical/molecular mechanical (QM/MM) MD calculations constrained by NMR-derived data. To illustrate this approach, we refine the zinc coordination complex structure of the zinc sensing transcriptional repressor protein Staphylococcus aureus CzrA, generating over 130 ns of MD and QM/MM MD NMR-data compliant sampling. In addition to refining the first coordination shell structure of the Zn(II) ion, this protocol benefits from being performed in a periodically replicated solvation environment including long-range electrostatics. We determine that unrestrained (not based on NMR data) MD simulations correlated to the NMR data in a time-averaged ensemble. The accurate solution structure ensemble of the metal-bound protein accurately describes the role of conformational sampling in allosteric regulation of DNA binding by zinc and serves to validate our previous unrestrained MD simulations of CzrA. This methodology has potentially broad applicability in the structure determination of metal ion bound proteins, protein folding and metal template protein-design studies.

Conformational and functional effects induced by D- and L-amino acid epimerization on a single gene encoded peptide from the skin secretion of Hypsiboas punctatus

Skin secretion of Hypsiboas punctatus is the source of a complex mixture of bioactive compounds where peptides and small proteins prevail, similarly to many other amphibians. Among dozens of molecules isolated from H. punctatus in a proteomic based approach, we report here the structural and functional studies of a novel peptide named Phenylseptin (FFFDTLKNLAGKVIGALT-NH2) that was purified as two naturally occurring D- and L-Phes configurations. The amino acid epimerization and C-terminal amidation for both molecules were confirmed by a combination of techniques including reverse-phase UFLC, ion mobility mass spectrometry, high resolution MS/MS experiments, Edman degradation, cDNA sequencing and solid-phase peptide synthesis. RMSD analysis of the twenty lowest-energy (1)H NMR structures of each peptide revealed a major 90° difference between the two backbones at the first four N-terminal residues and substantial orientation changes of their respective side chains. These structural divergences were considered to be the primary cause of the in vitro quantitative differences in antimicrobial activities between the two molecules. Finally, both molecules elicited equally aversive reactions in mice when delivered orally, an effect that depended entirely on peripheral gustatory pathways.

Assessing and refining molecular dynamics simulations of proteins with nuclear magnetic resonance data

The sophistication of the force fields, algorithms and hardware used for molecular dynamics (MD) simulations of proteins is continuously increasing. No matter how advanced the methodology, however, it is essential to evaluate the appropriateness of the structures sampled in a simulation by comparison with quantitative experimental data. Solution nuclear magnetic resonance (NMR) data are particularly useful for checking the quality of protein simulations, as they provide both structural and dynamic information on a variety of temporal and spatial scales. Here, various features and implications of using NMR data to validate and bias MD simulations are outlined, including an overview of the different types of NMR data that report directly on structural properties and of relevant simulation techniques. The focus throughout is on how to properly account for conformational averaging, particularly within the context of the assumptions inherent in the relationships that link NMR data to structural properties.

MaxOcc: a web portal for maximum occurrence analysis

The MaxOcc web portal is presented for the characterization of the conformational heterogeneity of two-domain proteins, through the calculation of the Maximum Occurrence that each protein conformation can have in agreement with experimental data. Whatever the real ensemble of conformations sampled by a protein, the weight of any conformation cannot exceed the calculated corresponding Maximum Occurrence value. The present portal allows users to compute these values using any combination of restraints like pseudocontact shifts, paramagnetism-based residual dipolar couplings, paramagnetic relaxation enhancements and small angle X-ray scattering profiles, given the 3D structure of the two domains as input. MaxOcc is embedded within the NMR grid services of the WeNMR project and is available via the WeNMR gateway at http://py-enmr.cerm.unifi.it/access/index/maxocc . It can be used freely upon registration to the grid with a digital certificate.

The relation of the X-ray B-factor to protein dynamics: insights from recent dynamic solid-state NMR data

In addressing the potential use of B-factors derived from X-ray scattering data of proteins for the understanding the (functional) dynamics of proteins, we present a comparison of B-factors of five different proteins (SH3 domain, Crh, GB1, ubiquitin and thioredoxin) with data from recent solid-state nuclear magnetic resonance experiments reflecting true (rotational) dynamics on well-defined timescales. Apart from trivial correlations involving mobile loop regions and chain termini, we find no significant correlation of B-factors with the dynamic data on any of the investigated timescales, concluding that there is no unique and general correlation of B-factors with the internal reorientational dynamics of proteins.

Understanding biomolecular motion, recognition, and allostery by use of conformational ensembles

We review the role conformational ensembles can play in the analysis of biomolecular dynamics, molecular recognition, and allostery. We introduce currently available methods for generating ensembles of biomolecules and illustrate their application with relevant examples from the literature. We show how, for binding, conformational ensembles provide a way of distinguishing the competing models of induced fit and conformational selection. For allostery we review the classic models and show how conformational ensembles can play a role in unravelling the intricate pathways of communication that enable allostery to occur. Finally, we discuss the limitations of conformational ensembles and highlight some potential applications for the future.

Active site dynamics in NADH oxidase from Thermus thermophilus studied by NMR spin relaxation

We have characterized the backbone dynamics of NADH oxidase from Thermus thermophilus (NOX) using a recently-developed suite of NMR experiments designed to isolate exchange broadening, together with (15)N R (1), R (1ρ ), and {(1)H}-(15)N steady-state NOE relaxation measurements performed at 11.7 and 18.8 T. NOX is a 54 kDa homodimeric enzyme that belongs to a family of structurally homologous flavin reductases and nitroreductases with many potential biotechnology applications. Prior studies have suggested that flexibility is involved in the catalytic mechanism of the enzyme. The active site residue W47 was previously identified as being particularly important, as its level of solvent exposure correlates with enzyme activity, and it was observed to undergo "gating" motions in computer simulations. The NMR data are consistent with these findings. Signals from W47 are dynamically broadened beyond detection and several other residues in the active site have significant R ( ex ) contributions to transverse relaxation rates. In addition, the backbone of S193, whose side chain hydroxyl proton hydrogen bonds directly with the FMN cofactor, exhibits extensive mobility on the ns-ps timescale. We hypothesize that these motions may facilitate structural rearrangements of the active site that allow NOX to accept both FMN and FAD as cofactors.

How uniform is the peptide plane geometry? A high-accuracy NMR study of dipolar Cα-C'/H N-N cross-correlated relaxation

Highly precise and accurate measurements of very small NMR cross-correlated relaxation rates, namely those between protein H (i) (N) -N(i) and C (i-1) (α) -C(i-1)' dipoles, are demonstrated with an error of 0.03 s(-1) for GB3. Because the projection angles between the two dipole vectors are very close to the magic angle the rates range only from -0.2 to +0.2 s(-1). Small changes of the average vector orientations have a dramatic impact on the relative values. The rates suggest deviation from idealized peptide plane geometry caused by twists around the C'-N bonds and/or pyramidalization of the nitrogen atoms. A clear alternating pattern along the sequence is observed in β strands 1, 3 and 4 of GB3, where the side chains of almost all residues with large positive rates are solvent exposed. In the α helix all rates are relatively large and positive. Some of the currently most accurate structures of GB3 determined by both high resolution X-ray crystallography and NMR are in satisfactory agreement with the experimental rates in the helix and β strand 3, but not in the loops and the two central strands of the sheet for which no alternating pattern is predicted.