Protein dynamics from solution NMR

Solution nuclear magnetic resonance (NMR) spectroscopy is unique in its ability to elucidate the details of atomic-level structural and dynamical properties of biological macromolecules under native-like conditions. Recent advances in NMR techniques and protein sample preparation now allow comprehensive investigation of protein dynamics over timescales ranging 14 orders of magnitude at nearly every atomic site. Thus, solution NMR is poised to reveal aspects of the physico-chemical properties that govern the ensemble distribution of protein conformers and the dynamics of their interconversion. We review these advances as well as their recent application to the study of proteins.     

Protein dynamics from NMR

This review surveys recent investigations of conformational fluctuations of proteins in solution using NMR techniques. Advances in experimental methods have provided more accurate means of characterizing fast and slow internal motions as well as overall diffusion. The information obtained from NMR dynamics experiments provides insights into specific structural changes or configurational energetics associated with function. A variety of applications illustrate that studies of protein dynamics provide insights into protein-protein interactions, target recognition, ligand binding, and enzyme function.     

Protein dynamics from chemical shift and dipolar rotational spin-echo 15N NMR

The partial collapse of dipolar and chemical shift tensors for peptide NH and for the amide NH at cross-link sites in cell wall peptidoglycan, of intact lyophilized cells of Aerococcus viridans, indicates NH vector root-mean-square fluctuations of 23 degrees. This result is consistent with the local mobility calculated in typical picosecond regime computer simulations of protein dynamics in the solid state. The experimental root-mean-square angular fluctuations for both types of NH vectors increase to 37 degrees for viable wet cells at 10 degrees C. The similarity in mobilities for both general protein and cell wall peptidoglycan suggests that one additional motion in wet cells involves cooperative fluctuations of segments of cell walls, attached proteins, and associated cytoplasmic proteins.     

NMR solution structure and dynamics of motilin in isotropic phospholipid bicellar solution

The structure and dynamics of the gastrointestinal peptide hormone motilin, consisting of 22 amino acid residues, have been studied in the presence of isotropic q=0.5 phospholipid bicelles. The NMR solution structure of the peptide in acidic bicelle solution was determined from 203 NOE-derived distance constraints and six backbone torsion angle constraints. Dynamic properties for the ¹³C-¹H vector in Leu10 were determined for motilin specifically labeled with ¹³C at this position by analysis of multiple-field relaxation data. The structure reveals an ordered -helical conformation between Glu9 and Lys20. The N-terminus is also well structured with a turn resembling that of a classical -turn. The ¹³C dynamics clearly show that motilin tumbles slowly in solution, with a correlation time characteristic of a large object. It was also found that motilin has a large degree of local flexibility as compared with what has previously been reported in SDS micelles. The results show that motilin interacts with the bicelle, displaying motional properties of a peptide bound to a membrane. In comparison, motilin in neutral bicelles seems less structured and more flexible. This study shows that the small isotropic bicelles are well suited for use as membrane-mimetic for structural as well as dynamical investigations of membrane-bound peptides by high-resolution NMR.     

Protein dynamics in solution and in a crystalline environment: a molecular dynamics study

The effect of a solvent and a crystalline environment on the dynamics of proteins is investigated by the method of computer simulation. Three 25-ps molecular dynamics simulations at 300 K of the bovine pancreatic trypsin inhibitor (BPTI), consisting of 454 heavy atoms, are compared: one of BPTI in vacuo, one of BPTI in a box with 2647 spherical nonpolar solvent atoms, and one of BPTI surrounded by fixed crystal image atoms. Both average and time-dependent molecular properties are examined to determine the effect of the environment on the behavior of the protein. The dynamics of BPTI in solution or in the crystal environment are found to be very similar to that found in the vacuum calculation. The primary difference in the average properties is that the equilibrium structure in the presence of solvent or the crystal field is significantly closer to the X-ray structure than is the vacuum result; concomitantly, the more realistic environment leads to a number density closer to experiment. The presence of solvent has a negligible effect on the overall magnitude of the positional or dihedral angle fluctuations in the interior of the protein; however, there are changes in the decay times of the fluctuations of interior atoms. For surface residues, both the magnitude and the time course of the motions are significantly altered by the solvent. There tends to be an increase in the displacements of long side chains and the flexible parts of the main chain that protrude into the solvent. Further, these motions tend to have a more diffusive character with longer relaxation times than in vacuo. The crystal environment has a specific effect on a number of side chains which are held in relatively fixed positions through hydrogen-bond and electric interactions with the neighboring protein atoms. Most of the effects of the solution environment seem to be sufficiently nonspecific that it may be possible to model them by applying a mean field and stochastic dynamic methods.     

Protein structure determination in solution by NMR spectroscopy

The introduction of nuclear magnetic resonance (NMR) spectroscopy as a second method for protein structure determination at atomic resolution, in addition to x-ray diffraction in single crystals, has already led to a significant increase in the number of known protein structures. The NMR method provides data that are in many ways complementary to those obtained from x-ray crystallography and thus promises to widen our view of protein molecules, giving a clearer insight into the relation between structure and function.     

Structure and dynamics of a membrane protein in micelles from three solution NMR experiments

Three solution NMR experiments on a uniformly ¹⁵N labeled membrane protein in micelles provide sufficient information to describe the structure, topology, and dynamics of its helices, as well as additional information that characterizes the principal features of residues in terminal and inter-helical loop regions. The backbone amide resonances are assigned with an HMQC-NOESY experiment and the backbone dynamics are characterized by a ¹H-¹⁵N heteronuclear NOE experiment, which clearly distinguishes between the structured helical residues and the more mobile residues in the terminal and interhelical loop regions of the protein. The structure and topology of the helices are described by Dipolar waves and PISA wheels derived from experimental measurements of residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs). The results show that the membrane-bound form of Pf1 coat protein has a 20-residue trans-membrane hydrophobic helix with an orientation that differs by about 90° from that of an 8-residue amphipathic helix. This combination of three-experiments that yields Dipolar waves and PISA wheels has the potential to contribute to high-throughput structural characterizations of membrane proteins.     

A theory of protein dynamics to predict NMR relaxation

We present a theoretical, site-specific, approach to predict protein subunit correlation times, as measured by NMR experiments of ¹H-¹⁵N nuclear Overhauser effect, spin-lattice relaxation, and spin-spin relaxation. Molecular dynamics simulations are input to our equation of motion for protein dynamics, which is solved analytically to produce the eigenvalues and the eigenvectors that specify the NMR parameters. We directly compare our theoretical predictions to experiments and to simulation data for the signal transduction chemotaxis protein Y (CheY), which regulates the swimming response of motile bacteria. Our theoretical results are in good agreement with both simulations and experiments, without recourse to adjustable parameters. The theory is general, since it allows calculations of any dynamical property of interest. As an example, we present theoretical calculations of NMR order parameters and x-ray Debye-Waller temperature factors; both quantities show good agreement with experimental data.     

Protein dynamics and enzyme catalysis: insights from simulations

The role of protein dynamics in enzyme catalysis is one of the most active and controversial areas in enzymology today. Some researchers claim that protein dynamics are at the heart of enzyme catalytic efficiency, while others state that dynamics make no significant contribution to catalysis. What is the biochemist - or student - to make of the ferocious arguments in this area? Protein dynamics are complex and fascinating, as molecular dynamics simulations and experiments have shown. The essential question is: do these complex motions have functional significance? In particular, how do they affect or relate to chemical reactions within enzymes, and how are chemical and conformational changes coupled together? Biomolecular simulations can analyse enzyme reactions and dynamics in atomic detail, beyond that achievable in experiments: accurate atomistic modelling has an essential part to play in clarifying these issues. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.     

ProDy: protein dynamics inferred from theory and experiments

SUMMARY: We developed a Python package, ProDy, for structure-based analysis of protein dynamics. ProDy allows for quantitative characterization of structural variations in heterogeneous datasets of structures experimentally resolved for a given biomolecular system, and for comparison of these variations with the theoretically predicted equilibrium dynamics. Datasets include structural ensembles for a given family or subfamily of proteins, their mutants and sequence homologues, in the presence/absence of their substrates, ligands or inhibitors. Numerous helper functions enable comparative analysis of experimental and theoretical data, and visualization of the principal changes in conformations that are accessible in different functional states. ProDy application programming interface (API) has been designed so that users can easily extend the software and implement new methods. AVAILABILITY: ProDy is open source and freely available under GNU General Public License from http://www.csb.pitt.edu/ProDy/.     

Transmembrane helix orientation and dynamics: insights from ensemble dynamics with solid-state NMR observables

As the major component of membrane proteins, transmembrane helices embedded in anisotropic bilayer environments adopt preferential orientations that are characteristic or related to their functional states. Recent developments in solid-state nuclear magnetic resonance (SSNMR) spectroscopy have made it possible to measure NMR observables that can be used to determine such orientations in a native bilayer environment. A quasistatic single conformer model is frequently used to interpret the SSNMR observables, but important motional information can be missing or misinterpreted in the model. In this work, we have investigated the orientation of the single-pass transmembrane domain of viral protein ”u“ (VpuTM) from HIV-1 by determining an ensemble of structures using multiple conformer models based on the SSNMR ensemble dynamics technique. The resulting structure ensemble shows significantly larger orientational fluctuations while the ensemble-averaged orientation is compatible with the orientation based on the quasistatic model. This observation is further corroborated by comparison with the VpuTM orientation from comparative molecular dynamics simulations in explicit bilayer membranes. SSNMR ensemble dynamics not only reveals the importance of transmembrane helix dynamics in interpretation of SSNMR observables, but also provides a means to simultaneously extract both transmembrane helix orientation and dynamics information from the SSNMR measurements.     

Protein solution structure calculations in solution: Solvated molecular dynamics refinement of calbindin D9k

The three-dimensional solution structures of proteins determinedwith NMR-derived constraints are almost always calculated in vacuo. Thesolution structure of (Ca²⁺)_{_2}-calbindinD_9k has been redetermined by new restrained molecular dynamics(MD) calculations that include Ca²⁺ ions and explicit solventmolecules. Four parallel sets of MD refinements were run to provide accuratecomparisons of structures produced in vacuo, in vacuo withCa²⁺ ions, and with two different protocols in a solvent bathwith Ca²⁺ ions. The structural ensembles were analyzed interms of structural definition, molecular energies, packing density,solvent-accessible surface, hydrogen bonds, and the coordination of calciumions in the two binding loops. Refinement including Ca²⁺ ionsand explicit solvent results in significant improvements in the precisionand accuracy of the structure, particularly in the binding loops. Theseresults are consistent with results previously obtained in free MDsimulations of proteins in solution and show that the rMD refinedNMR-derived solution structures of proteins, especially metalloproteins, canbe significantly improved by these strategies.     

A solution NMR view of protein dynamics in the biological membrane

Structure determination of membrane-associated proteins (MPs) represents a frontier of structural biology that is characterized by unique challenges in sample preparation and data acquisition. No less important is our ability to study the dynamics of MPs, since MP flexibility and characteristic motions often make sizeable contributions to their function. This review focuses on solution state NMR methods to characterize dynamics of MPs in the membrane environment. NMR approaches to study molecular motions on a wide range of time-scales and their application to membrane proteins are described. Studies of polytopic and bitopic MPs demonstrating the power of such methods to characterize the dynamic behavior of MPs and their interaction with the membrane-mimicking surroundings are presented. Attempts are made to place the dynamic conclusions into a biological context. The importance and limitations of such investigations guarantee that further developments in this field will be actively pursued.     

Protein structure and dynamics in nonaqueous solvents: insights from molecular dynamics simulation studies

Protein structure and dynamics in nonaqueous solvents are here investigated using molecular dynamics simulation studies, by considering two model proteins (ubiquitin and cutinase) in hexane, under varying hydration conditions. Ionization of the protein groups is treated assuming "pH memory," i.e., using the ionization states characteristic of aqueous solution. Neutralization of charged groups by counterions is done by considering a counterion for each charged group that cannot be made neutral by establishing a salt bridge with another charged group; this treatment is more physically reasonable for the nonaqueous situation, contrasting with the usual procedures. Our studies show that hydration has a profound effect on protein stability and flexibility in nonaqueous solvents. The structure becomes more nativelike with increasing values of hydration, up to a certain point, when further increases render it unstable and unfolding starts to occur. There is an optimal amount of water, approximately 10% (w/w), where the protein structure and flexibility are closer to the ones found in aqueous solution. This behavior can explain the experimentally known bell-shaped dependence of enzyme catalysis on hydration, and the molecular reasons for it are examined here. Water and counterions play a fundamental and dynamic role on protein stabilization, but they also seem to be important for protein unfolding at high percentages of bound water.