Contact model for the prediction of NMR N-H order parameters in globular proteins

An analytical relationship is presented for the estimation of NMR S2 order parameters of N-HN vectors of the protein backbone from high-resolution protein structures. The relationship solely depends on close contacts of the peptide plane to the rest of the protein. Application of the relationship to a number of proteins with high-resolution X-ray and NMR structures yields S2 values that are in good agreement with the ones determined from experimental relaxation data.     

Conformational flexibility of a microcrystalline globular protein: order parameters by solid-state NMR spectroscopy

The majority of protein structures are determined in the crystalline state, yet few methods exist for the characterization of dynamics for crystalline biomolecules. Solid-state NMR can be used to probe detailed dynamic information in crystalline biomolecules. Recent advances in high-resolution solid-state NMR have enabled the site-specific assignment of (13)C and (15)N nuclei in proteins. With the use of multidimensional separated-local-field experiments, we report the backbone and side chain conformational dynamics of ubiquitin, a globular microcrystalline protein. The measurements of molecular conformational order parameters are based on heteronuclear dipolar couplings, and they are correlated to assigned chemical shifts, to obtain a global perspective on the sub-microsecond dynamics in microcrystalline ubiquitin. A total of 38 Calpha, 35 Cbeta and multiple side chain unique order parameters are collected, and they reveal the high mobility of ubiquitin in the microcrystalline state. In general the side chains show elevated motion in comparison with the backbone sites. The data are compared to solution NMR order parameter measurements on ubiquitin. The SSNMR measurements are sensitive to motions on a broader time scale (low microsecond and faster) than solution NMR measurements (low nanosecond and faster), and the SSNMR order parameters are generally lower than the corresponding solution values. Unlike solution NMR relaxation-based order parameters, order parameters for (13)C(1)H(2) spin systems are readily measured from the powder line shape data. These results illustrate the potential for detailed, extensive, and site-specific dynamic studies of biopolymers by solid-state NMR.     

Prediction of NMR order parameters in proteins using weighted protein contact-number model

In the NMR experiment, the protein backbone motion can be described by the N–H order parameters. Though protein dynamics is determined by a complex network of atomic interactions, we show that the order parameter of residues can be determined using a very simple method, the weighted protein contact number model. We computed for each Cα atom the number of neighboring Cα atoms weighted by the inverse distance squared between them. We show that the weighted contact number of each residue is directly related to its order parameter. Despite the simplicity of this model, it performs better than the other method. Since we can compute the order parameters directly from the topological properties (such as protein contact number) of protein structures, our study underscores a very direct link between protein topological structure and its dynamics.     

Quantifying millisecond exchange dynamics in proteins by CPMG relaxation dispersion NMR using side-chain 1H probes

A Carr-Purcell-Meiboom-Gill relaxation dispersion experiment is presented for quantifying millisecond time-scale chemical exchange at side-chain (1)H positions in proteins. Such experiments are not possible in a fully protonated molecule because of magnetization evolution from homonuclear scalar couplings that interferes with the extraction of accurate transverse relaxation rates. It is shown, however, that by using a labeling strategy whereby proteins are produced using {(13)C,(1)H}-glucose and D(2)O a significant number of 'isolated' side-chain (1)H spins are generated, eliminating such effects. It thus becomes possible to record (1)H dispersion profiles at the β positions of Asx, Cys, Ser, His, Phe, Tyr, and Trp as well as the γ positions of Glx, in addition to the methyl side-chain moieties. This brings the total of amino acid side-chain positions that can be simultaneously probed using a single (1)H dispersion experiment to 16. The utility of the approach is demonstrated with an application to the four-helix bundle colicin E7 immunity protein, Im7, which folds via a partially structured low populated intermediate that interconverts with the folded, ground state on the millisecond time-scale. The extracted (1)H chemical shift differences at side-chain positions provide valuable restraints in structural studies of invisible, excited states, complementing backbone chemical shifts that are available from existing relaxation dispersion experiments.     

A 2H NMR relaxation experiment for the measurement of the time scale of methyl side-chain dynamics in large proteins

An NMR experiment is presented for the measurement of the time scale of methyl side-chain dynamics in proteins that are labeled with methyl groups of the (13)CHD(2) variety. The measurement is accomplished by selecting a magnetization mode that to excellent approximation relaxes in a single-exponential manner with a T(1)-like rate. The combination of R(1)((13)CHD(2)) and R(2)((13)CHD(2)) (2)H relaxation rates facilitates the extraction of motional parameters from (13)CHD(2)-labeled proteins exclusively. The utility of the methodology is demonstrated with applications to proteins with tumbling times ranging from 2 ns (protein L, 7.5 kDa, 45 degrees C) to 54 ns (malate synthase G, 82 kDa, 37 degrees C); dynamics parameters are shown to be in excellent agreement with those obtained in (2)H NMR studies of other methyl isotopomers. A consistency relationship is found to exist between R(1)((13)CHD(2)) and the relaxation rates of pure longitudinal and quadrupolar order modes in (13)CH(2)D-labeled methyl groups, and experimental rates measured for a number of proteins are shown to be in excellent agreement with expectations based on theory. The present methodology extends the applicability of (2)H relaxation methods for the quantification of side-chain dynamics in high molecular weight proteins.     

Temperature dependence of NMR order parameters and protein dynamics

The helical subdomain, HP36, of the F-actin-binding headpiece domain of chicken villin, is the smallest naturally occurring polypeptide that folds to a thermostable compact structure. Unconstrained molecular dynamics simulations and constrained molecular dynamics simulations using umbrella sampling are used to study the temperature dependence of internal motions of the backbone amide moieties of HP36. The potential of mean force (PMF) for the N-H bond vector, determined from the constrained simulations, is found to be temperature dependent. A simple analytical expression is derived that describes the temperature dependence of the PMF. The parameters of this model are obtained from the PMF, from the unconstrained molecular dynamics simulations, or from experimental values of the generalized order parameter. The results provide a linkage between experimental and theoretical measures of the temperature dependence of protein motions.     

On relationships between surfactant type and globular proteins interactions in solution

The binding of sodium perfluorooctanoate (C8FONa), sodium octanoate (C8HONa), lithium perfluorooctanoate (C8FOLi), and sodium dodecanoate (C12HONa) onto myoglobin, ovalbumin, and catalase in water has been characterized using electrophoretic mobility. The tendency of the protein-surfactant complexes to change their charge in the order catalase < ovalbumin < myoglobin was observed which was related to the contents of alpha-helices in the proteins. alpha-Helices are more hydrophobic than beta-sheets. The effect of surfactant on the zeta potentials follows C8HONa < C8FONa < C8FOLi < C12HONa for catalase and ovalbumin; and C8HONa < C8FOLi < C8FONa < C12HONa for myoglobin. The numbers of binding sites on the proteins were determined from the observed increases of the zeta-potential as a function of surfactant concentration in the regions where the binding was a consequence of the hydrophobic effect. The Gibbs energies of binding of the surfactants onto the proteins were evaluated. For all systems, Gibbs energies are negative and large at low concentrations (where binding to the high energy sites takes place) and become less negative at higher ones. This fact suggests a saturation process. Changes in Gibbs energies with the different proteins and surfactants under study have been found to follow same sequence than that found for the charge. The role of hydrophobic interactions in these systems has been demonstrated to be the predominant.     

Protein side-chain dynamics observed by solution- and solid-state NMR: comparative analysis of methyl 2H relaxation data

Rapid advances in solid-state MAS NMR made it possible to probe protein dynamics on a per-residue basis, similar to solution experiments. In this work we compare methyl 2H relaxation rates measured in the solid and liquid samples of alpha-spectrin SH3 domain. The solution data are treated using a model-free approach to separate the contributions from the overall molecular tumbling and fast internal motion. The latter part forms the basis for comparison with the solid-state data. Although the accuracy of solid-state measurements is limited by deuterium spin diffusion, the results suggest a significant similarity between methyl dynamics in the two samples. This is a potentially important observation, preparing the ground for combined analysis of the dynamics data by solid- and solution-state NMR.     

Verification of the correlation between geometric parameters of molecules and the parameters of their 1H NMR spectra

The dihedral and bond direction angles between all pairs of vicinal protons of the arabinofuranose residue were calculated from the coordinatees of the hydrogen atoms found by an X-ray study of 3-O-acetyl-β-L -arabinofuranose 1,2,5-orthobenzoate. The values found were compared with those calculated with the help of correlation equations previously proposed by Karplus and recently by the authors, linking the values of those angles with the spin-spin coupling constants of vicinal protons [³J(H,H′)]. It has been found that the best agreement between the angles found crystallographically and calculated from the ¹H NMR data can be achieved using the equation which includes bond direction angles and the sum of the chemical shifts of the protons involved.     

The use of selective deuteration for the sequence specific 1H NMR assignment of larger proteins

The possibility of extending NMR methods for structure determination to larger proteins (MW > 10 kD) depends on the development of isotopic labeling protocols for the simplification of their NMR spectra (isotopic spectral editing). We describe here the successful use of selective deuteration to obtain sequence specific assignments for (thus far) more than 50% of the residues of the trp repressor protein (25 kD). This is the largest protein for which detailed sequence specific assignments have been attempted to-date.     

Conformational interactions between a phenyl ring and a side-chain halide substituent: A 1H NMR and molecular mechanics study of some 2-aryl-1-halopropanes

¹H NMR data are reported for a series of 2-aryl-1-halopropanes. Vicinal coupling constants in the CH₂CH—fragment show that the rotamer populations about the CC bond are sensitive to para substituents. The ratio of anti:gauche aryl/halide conformers is greatest when the para substituent is the electron-donating ethyl group and least when it is the strongly electron-withdrawing nitro group. This points to a non-steric conformational interaction involving the ring and the sidechain heteroatom. Comparison of the empirical results with conformational preferences predicted from molecular mechanics calculations using the COSMIC force field suggests that the interaction serves to enhance the population of the anti arrangement.     

Strong correlation between SHAPE chemistry and the generalized NMR order parameter (S2) in RNA

The functions of most RNA molecules are critically dependent on the distinct local dynamics that characterize secondary structure and tertiary interactions and on structural changes that occur upon binding by proteins and small molecule ligands. Measurements of RNA dynamics at nucleotide resolution set the foundation for understanding the roles of individual residues in folding, catalysis, and ligand recognition. In favorable cases, local order in small RNAs can be quantitatively analyzed by NMR in terms of a generalized order parameter, S2. Alternatively, SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) chemistry measures local nucleotide flexibility in RNAs of any size using structure-sensitive reagents that acylate the 2'-hydroxyl position. In this work, we compare per-residue RNA dynamics, analyzed by both S2 and SHAPE, for three RNAs: the HIV-1 TAR element, the U1A protein binding site, and the Tetrahymena telomerase stem loop 4. We find a very strong correlation between the two measurements: nucleotides with high SHAPE reactivities consistently have low S2 values. We conclude that SHAPE chemistry quantitatively reports local nucleotide dynamics and can be used with confidence to analyze dynamics in large RNAs, RNA-protein complexes, and RNAs in vivo.     

Deuterium spin probes of side-chain dynamics in proteins. 2. Spectral density mapping and identification of nanosecond time-scale side-chain motions

In the previous paper in this issue we have demonstrated that it is possible to measure the five different relaxation rates of a deuteron in (13)CH(2)D methyl groups of (13)C-labeled, fractionally deuterated proteins. The extensive set of data acquired in these experiments provides an opportunity to investigate side-chain dynamics in proteins at a level of detail that heretofore was not possible. The data, acquired on the B1 domain of peptostreptococcal protein L, include 16 (9) relaxation measurements at 4 (2) different magnetic field strengths, 25 degrees C (5 degrees C). These data are shown to be self-consistent and are analyzed using a spectral density mapping procedure which allows extraction of values of the spectral density function at a number of frequencies with no assumptions about the underlying dynamics. Dynamics data from 31 of 35 methyls in the protein for which data could be obtained were well-fitted using the two-parameter Lipari-Szabo model (Lipari, G.; Szabo, A. J. Am. Chem. Soc. 1982, 104, 4546). The data from the remaining 4 methyls can be fitted using a three-parameter version of the Lipari-Szabo model that takes into account, in a simple manner, additional nanosecond time-scale local dynamics. This interpretation is supported by analysis of a molecular dynamics trajectory where spectral density profiles calculated for side-chain methyl sites reflect the influence of slower (nanosecond) time-scale motions involving jumps between rotameric wells. A discussion of the minimum number of relaxation measurements that are necessary to extract the full complement of dynamics information is presented along with an interpretation of the extracted dynamics parameters.     

Specific ion effects in solutions of globular proteins: comparison between analytical models and simulation

Monte Carlo simulations have been performed for ion distributions outside a single globular macroion and for a pair of macroions, in different salt solutions. The model that we use includes both electrostatic and van der Waals interactions between ions and between ions and macroions. Simulation results are compared with the predictions of the Ornstein-Zernike equation with the hypernetted chain closure approximation and the nonlinear Poisson-Boltzmann equation, both augmented by pertinent van der Waals terms. Ion distributions from analytical approximations are generally very close to the simulation results. This demonstrates that properties that are related to ion distributions in the double layer outside a single interface can to a good approximation be obtained from the Poisson-Boltzmann equation. We also present simulation and integral equation results for the mean force between two globular macroions (with properties corresponding to those of hen-egg-white lysozyme protein at pH 4.3) in different salt solutions. The mean force and potential of mean force between the macroions become more attractive upon increasing the polarizability of the counterions (anions), in qualitative agreement with experiments. We finally show that the deduced second virial coefficients agree quite well with experimental results.     

Decharging of globular proteins and protein complexes in electrospray

Electrospray ionization mass spectrometry (ESI-MS) is a valuable tool in structural biology for investigating globular proteins and their biomolecular interactions. During the electrospray ionization process, proteins become desolvated and multiply charged, which may influence their structure. Reducing the net charge obtained during the electrospray process may be relevant for studying globular proteins. In this report we demonstrate the effect of a series of inorganic and organic gas-phase bases on the number of charges that proteins and protein complexes attain. Solution additives with very strong gas-phase basicities (GB) were identified among the so-called "proton sponges". The gas-phase proton affinities (PA) of the compounds that were added to the aqueous protein solutions ranged from 700 to 1050 kJ mol(-1). Circular dichroism studies showed that in these solutions the proteins retain their globular structures. The size of the proteins investigated ranged from the 14.3 kDa lysozyme up to the 800 kDa tetradecameric chaperone complex GroEL. Decharging of the proteins in the electrospray process by up to 60 % could be achieved by adding the most basic compounds rather than the more commonly used ammonium acetate additive. This decharging process probably results from proton competition events between the multiply protonated protein ions and the basic additives just prior to the final desolvation. We hypothesize that such globular protein species, which attain relatively few charges during the ionization event, obtain a gas-phase structure that more closely resembles their solution-phase structure. Thus, these basic additives can be useful in the study of the biologically relevant properties of globular proteins by using mass spectrometry.     

Reorientational eigenmode dynamics: a combined MD/NMR relaxation analysis method for flexible parts in globular proteins

An approach is presented for the interpretation of heteronuclear NMR spin relaxation data in mobile protein parts in terms of reorientational eigenmode dynamics. The method is based on the covariance matrix of the spatial functions of the nuclear spin interactions that cause relaxation expressed as spherical harmonics of rank 2. The approach was applied to characterize the dynamics of a loop region of ubiquitin. The covariance matrix was determined from a conformational ensemble generated by a 5 ns molecular dynamics simulation. It was found that the time correlation functions of the dominant eigenmodes decay in good approximation with a single correlation time. From the reorientational eigenmodes, their eigenvalues, and correlation times, NMR relaxation data were calculated in accordance with Bloch-Wangsness-Redfield relaxation theory and directly compared with experimental (15)N relaxation parameters. Using a fitting procedure, agreement between calculated and experimental data was improved significantly by adjusting eigenvalues and correlation times of the dominant modes. The presented procedure provides detailed information on correlated reorientational dynamics of flexible parts in globular proteins. The covariance matrix was linked to the covariance matrix of backbone dihedral angle fluctuations, allowing one to study the motional behavior of these degrees of freedom on nano- and subnanosecond time scales.     

Carbon-13 NMR distinction between categories of molecular order and disorder in cellulose

Differences between values of proton rotating-frame spin relaxation time constants can be exploited to separate a solid-state¹³C NMR spectrum of cellulose into subspectra of crystalline and noncrystalline regions. Variations in chemical shifts and¹³C spin-lattice relaxation time constants can then be used to study variations in molecular order and disorder within each of the two broader categories. Mechanical damage during Wiley milling increases the content of noncrystalline cellulose and changes the nature of molecular disorder within that category. Resolution enhancement of the subspectrum assigned to crystalline cellulose reveals pairs of signals at 83.9 and 84.9 ppm (cellulose I) or 86.8 and 88.3 ppm (cellulose II) assigned to C-4 on well-ordered crystal surfaces. A broader peak in the subspectrum of crystalline cellulose I is assigned to poorly-ordered surfaces. Relative proportions in Avicel microcrystalline cellulose were estimated as: 54% in crystal interiors, 22% on well-ordered surfaces, 8% on poorly-ordered surfaces, 16% in domains of disorder extending more than a few nanometres.     

Non-Exponential 1H and 2H NMR Relaxation and Self-Diffusion in Asphaltene-Maltene Solutions

The distribution of NMR relaxation times and diffusion coefficients in crude oils results from the vast number of different chemical species. In addition, the presence of asphaltenes provides different relaxation environments for the maltenes, generated by steric hindrance in the asphaltene aggregates and possibly by the spatial distribution of radicals. Since the dynamics of the maltenes is further modified by the interactions between maltenes and asphaltenes, these interactions—either through steric hindrances or promoted by aromatic-aromatic interactions—are of particular interest. Here, we aim at investigating the interaction between individual protonic and deuterated maltene species of different molecular size and aromaticity and the asphaltene macroaggregates by comparing the maltenes’ NMR relaxation (T1 and T2) and translational diffusion (D) properties in the absence and presence of the asphaltene in model solutions. The ratio of the average transverse and longitudinal relaxation rates, describing the non-exponential relaxation of the maltenes in the presence of the asphaltene, and its variation with respect to the asphaltene-free solutions are discussed. The relaxation experiments reveal an apparent slowing down of the maltenes’ dynamics in the presence of asphaltenes, which differs between the individual maltenes. While for single-chained alkylbenzenes, a plateau of the relaxation rate ratio was found for long aliphatic chains, no impact of the maltenes’ aromaticity on the maltene–asphaltene interaction was unambiguously found. In contrast, the reduced diffusion coefficients of the maltenes in presence of the asphaltenes differ little and are attributed to the overall increased viscosity.