Search

A transverse relaxation optimized spectroscopy (TROSY) approach is described for the optimal detection of NH 2 groups in asparagine and glutamine side chains of proteins. Specifically, we have developed NMR experiments for isolating the slow-relaxing 15 N and 1 H components of NH 2 multiplets. Although even modest sensitivity gains in 2D NH 2 -TROSY correlation maps compared to their decoupled NH 2 -HSQC counterparts can be achieved only occasionally, substantial improvements in resolution of the NMR spectra are demonstrated for asparagine and glutamine NH 2 sites of a buried cavity mutant, L99A, of T4 lysozyme at 5 ºC. The NH 2 -TROSY approach is applied to CPMG relaxation dispersion measurements at the side chain NH 2 positions of the L99A T4 lysozyme mutant - a model system for studies of the role of protein dynamics in ligand binding., Competing Interests: Declarations. Competing interests: The authors have no competing interests to declare that are relevant to the content of this article., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

We describe the utility of small nutation angle (acute; <90°) 1 H radiofrequency pulses for efficient manipulation of magnetization in selectively [ 13 CH 3 ]-labeled methyl groups of otherwise deuterated proteins. Focusing primarily on NMR applications that target either fast (pico-to-nanosecond) motions of the methyl group three-fold rotation axis, or slow (micro-to-millisecond) processes associated with chemical exchange, we show that significant simplification of the 13 CH 3 spin-system and, as a consequence, of NMR pulse schemes, may be achieved in certain cases by the proper choice of the flip-angle of the 1 H acute-angle pulse. In other instances, appropriate adjustment of acute-angle 1 H pulses permits optimization of the sensitivity of NMR experiments. The results of acute-angle pulse based NMR experiments are validated by comparison with well-established NMR techniques for the characterization of fast dynamics of methyl-containing side-chains and chemical exchange processes., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Published by Elsevier B.V.)

The present work describes an update to the protein covalent geometry and atomic radii parameters in the Xplor-NIH biomolecular structure determination package. In combination with an improved treatment of selected non-bonded interactions between atoms three bonds apart, such as those involving methyl hydrogens, and a previously developed term that affects the system's gyration volume, the new parameters are tested using structure calculations on 30 proteins with restraints derived from nuclear magnetic resonance data. Using modern structure validation criteria, including several formally adopted by the Protein Data Bank, and a clear measure of structural accuracy, the results show superior performance relative to previous Xplor-NIH implementations. Additionally, the Xplor-NIH structures compare favorably against originally determined NMR models., (Published 2024. This article is a U.S. Government work and is in the public domain in the USA.)

An NMR experiment is designed for accurate and robust measurement of transverse relaxation rates of degenerate 1 H transitions in selectively 13 CH 3 -labeled, deuterated small proteins. The measurement is based on the use of acute (<90°) angle 1 H radio-frequency pulses and relies on selection of the slow- and fast-relaxing components of methyl magnetization following the relaxation period in separate experiments. The R 2 decay series recorded with selection of the fast-relaxing components serves as a useful complement to the R 2 series acquired with selection of the slow-relaxing part, and permits the extension of the range of relative contributions of the fast- and slow-relaxing parts to apparent signal decay. The approach is experimentally verified on 13 CH 3 methyl groups of the ILV-{ 13 CH 3 }-labeled protein ubiquitin at 10 °C and 25 °C. The obtained methyl 1 H relaxation rates are in remarkably good agreement with the values obtained from well-established NMR techniques., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023. Published by Elsevier Inc.)

A methyl Transverse Relaxation Optimized Spectroscopy (methyl-TROSY) based, multiple quantum (MQ) 13 C Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion NMR experiment is described. The experiment is derived from the previously developed MQ 13 C- 1 H CPMG scheme (Korzhnev in J Am Chem Soc 126: 3964-73, 2004) supplemented with a CPMG train of refocusing 1 H pulses applied with constant frequency and synchronized with the 13 C CPMG pulse train. The optimal 1 H 'decoupling' scheme that minimizes the amount of fast-relaxing methyl MQ magnetization present during CPMG intervals, makes use of an XY-4 phase cycling of the refocusing composite 1 H pulses. For small-to-medium sized proteins, the MQ 13 C CPMG experiment has the advantage over its single quantum (SQ) 13 C counterpart of significantly reducing intrinsic, exchange-free relaxation rates of methyl coherences. For high molecular weight proteins, the MQ 13 C CPMG experiment eliminates complications in the interpretation of MQ 13 C- 1 H CPMG relaxation dispersion profiles arising from contributions to exchange from differences in methyl 1 H chemical shifts between ground and excited states. The MQ 13 C CPMG experiment is tested on two protein systems: (1) a triple mutant of the Fyn SH3 domain that interconverts slowly on the chemical shift time scale between the major folded state and an excited state folding intermediate; and (2) the 82-kDa enzyme Malate Synthase G (MSG), where chemical exchange at individual Ile δ1 methyl positions occurs on a much faster time-scale., (© 2023. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

A survey, primarily based on work in the authors' laboratory during the last 10 years, is provided of recent developments in NMR studies of exchange processes involving protein-ligand and protein-protein interactions. We start with a brief overview of the theoretical background of Dark state Exchange Saturation Transfer (DEST) and lifetime line-broadening (ΔR 2 ) NMR methodology. Some limitations of the DEST/ΔR 2 methodology in applications to molecular systems with intermediate molecular weights are discussed, along with the means of overcoming these limitations with the help of closely related exchange NMR techniques, such as the measurements of Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion, exchange-induced chemical shifts or rapidly-relaxing components of relaxation decays. Some theoretical underpinnings of the quantitative description of global dynamics of proteins on the surface of very high molecular weight particles (nanoparticles) are discussed. Subsequently, several applications of DEST/ΔR 2 methodology are described from a methodological perspective with an emphasis on providing examples of how kinetic and relaxation parameters for exchanging systems can be reliably extracted from NMR data for each particular model of exchange. Among exchanging systems that are not associated with high molecular weight species, we describe several exchange NMR-based studies that focus on kinetic modelling of transient pre-nucleation oligomerization of huntingtin peptides that precedes aggregation and fibril formation., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Published by Elsevier B.V.)

A simple NMR experiment for the measurement of transverse relaxation rates of degenerate 1 H spins in 13 CH 3 methyl groups of deuterated proteins, is described. The experiment relies on the use of acute-angle 1 H radio-frequency pulses, whereby a series of methyl 1 H R 2 decays is acquired with different values of the 1 H pulse flip-angle, which modulates the relative contributions of the slow- and fast-relaxing components of 1 H magnetization during the relaxation delay. These are subsequently analyzed simultaneously to extract the transverse relaxation rates of both components (R S and R F ), along with the rate of cross-relaxation between these components, δ. The dipolar 1 H- 1 H cross-correlated relaxation rate η = (R F - R S )/2 and the cross-relaxation rate δ, extracted from such acute-angle-pulse R 2 decay series are compared with those derived from well-established methodology that uses relaxation-violated methyl 1 H coherence transfer (so-called 'forbidden' experiments). Good agreement is achieved for the rates η derived from the two techniques, while slightly more accurate values of δ are obtained from analysis of acute-angle-pulse R 2 series. Recording of acute-angle-pulse R 2 series represents an attractive alternative to existing methodologies for quantifying methyl 1 H relaxation and dynamics of the methyl three-fold symmetry axis in selectively [ 13 CH 3 ]-methyl-labeled, deuterated proteins - particularly so for very high-molecular-weight proteins, where the measurements of R F rates or the build-up of methyl 1 H magnetization in relaxation-violated coherence transfer experiments, are problematic due to low sensitivity., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Published by Elsevier Inc.)

The dynamics of methyl-bearing side chains in proteins were probed by 13 C relaxation measurements of a number of 13 C magnetization modes in selectively 13 CH 3 -labeled methyl groups of proteins. We first show how 13 C magnetization modes in a 13 CH 3 spin-system can be isolated using acute-angle 1 H radio-frequency pulses. The parameters of methyl-axis dynamics, a measure of methyl-axis ordering ( S axis 2 ) and the correlation time of fast local methyl-axis motions (τ f ), derived from 13 C relaxation in 13 CH 3 groups are compared with their counterparts obtained from 13 methyl isotopomers. We show that in high-molecular-weight proteins, excellent correlations are obtained between the [ 13 CHD 2 methyl isotopomers. We show that in high-molecular-weight proteins, excellent correlations are obtained between the [ 13 CHD 2 ]-derived S axis 2 methyls. In smaller proteins, a certain degree of anticorrelation is observed between the 13 C magnetization of the I = 1/2 manifold in 13 CH 3 methyls. In smaller proteins, a certain degree of anticorrelation is observed between the S axis 2 and τ f values obtained from 13 C relaxation of the I = 1/2 manifold magnetization in 13 CH 3 methyls. These parameters can be partially decorrelated by inclusion in the analysis of relaxation data of the I = 3/2 manifold 13 C magnetization.

Optimized selection of the slow-relaxing components of single-quantum 13 C magnetization in 13 CH 3 methyl groups of proteins using acute (< 90°) angle 1 H radio-frequency pulses, is described. The optimal selection scheme is more relaxation-tolerant and provides sensitivity gains in comparison to the experiment where the undesired (fast-relaxing) components of 13 C magnetization are simply 'filtered-out' and only 90° 1 H pulses are employed for magnetization transfer to and from 13 C nuclei. When applied to methyl 13 C single-quantum Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiments for studies of chemical exchange, the selection of the slow-relaxing 13 C transitions results in a significant decrease in intrinsic (exchange-free) transverse spin relaxation rates of all exchanging species. For exchanging systems involving high-molecular-weight species, the lower transverse relaxation rates translate into an increase in the information content of the resulting relaxation dispersion profiles.

Dynamics of protein side chains is one of the principal determinants of conformational entropy in protein structures and molecular recognition events. We describe NMR experiments that rely on the use of magic-angle pulses for efficient isolation of degenerate 1 H transitions of the I=3/2 manifold of 13 CH 3 methyl groups, and serve as 'building blocks' for the measurement of transverse spin relaxation rates of the fast- and slow-relaxing 1 H transitions - the primary quantitative reporters of methyl axis dynamics in selectively { 13 CH 3 }-methyl-labelled, highly deuterated proteins. The magic-angle-pulse driven experiments are technically simpler and, in the absence of relaxation, predicted to be 2.3-fold more sensitive than previously developed analogous schemes. Validation of the methodology on a sample of { 13 CH 3 }-labeled ubiquitin demonstrates quantitative agreement between order parameters of methyl three-fold symmetry axis obtained with magic-angle-pulse driven experiments and other established NMR techniques, paving the way for studies of methyl axis dynamics in human DNAJB6b chaperone, a protein that undergoes exchange with high-molecular-weight oligomeric species., (© 2020 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

We introduce a new hydrogen bonding potential of mean force generated from high-quality crystal structures for use in Xplor-NIH structure calculations. This term applies to hydrogen bonds involving both backbone and sidechain atoms. When used in structure refinement calculations of 10 example protein systems with experimental distance, dihedral and residual dipolar coupling restraints, we demonstrate that the new term has superior performance to the previously developed hydrogen bonding potential of mean force used in Xplor-NIH., (Published 2019. This article is a U.S. Government work and is in the public domain in the USA.)

Over the past approximately three decades, isotope-directed NMR spectroscopy has become a powerful method for determining 3D structures of biological macromolecules and their complexes in solution. From a structural perspective NMR provides an invaluable tool for studying systems that are not amenable to crystallization, including intrinsically disordered proteins and weak complexes. In contrast to both X-ray crystallography and cryo-electron microscopy which afford a largely static view of the systems under consideration, the great power of NMR lies in its ability to quantitatively probe exchange dynamics between interconverting states, and to reveal and characterize at atomic resolution the existence of transient states that may be populated at levels as low as 1%. Such "excited" states play a key role in macromolecular recognition, allostery, signal transduction and macromolecular assembly, including the initial events involved in aggregation and amyloid formation. Optimal application of NMR to such systems of fundamental biological interest requires a sound footing of the physical underpinnings of today's and tomorrow's sophisticated NMR experiments., (Published by Elsevier Inc.)

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.

The global motions and exchange kinetics of a model protein, ubiquitin, bound to the surface of negatively charged lipid-based nanoparticles (liposomes) are derived from combined analysis of exchange lifetime broadening arising from binding to nanoparticles of differing size. The relative contributions of residence time and rotational tumbling to the total effective correlation time of the bound protein are modulated by nanoparticle size, thereby permitting the various motional and exchange parameters to be determined. The residence time of ubiquitin bound to the surface of both large and small unilamellar liposomes is ∼20 μs. Bound ubiquitin undergoes internal rotation about an axis approximately perpendicular to the lipid surface on a low microsecond time scale (∼2 μs), while simultaneously wobbling in a cone of semiangle 30-55° centered about the internal rotation axis on the nanosecond time scale. The binding interface of ubiquitin with liposomes is mapped by intermolecular paramagnetic relaxation enhancement using Gd(3+)-tagged vesicles, to a predominantly positively charged surface orthogonal to the internal rotation axis.

The energetic and volumetric properties of a three-state protein folding system, comprising a metastable triple mutant of the Fyn SH3 domain, have been investigated using pressure-dependent (15) N-relaxation dispersion NMR from 1 to 2500 bar. Changes in partial molar volumes (ΔV) and isothermal compressibilities (ΔκT ) between all the states along the folding pathway have been determined to reasonable accuracy. The partial volume and isothermal compressibility of the folded state are 100 mL mol(-1) and 40 μL mol(-1)  bar(-1) , respectively, higher than those of the unfolded ensemble. Of particular interest are the findings related to the energetic and volumetric properties of the on-pathway folding intermediate. While the latter is energetically close to the unfolded state, its volumetric properties are similar to those of the folded protein. The compressibility of the intermediate is larger than that of the folded state reflecting the less rigid nature of the former relative to the latter., (© 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Pulsed double electron-electron resonance (DEER) provides pairwise P(r) distance distributions in doubly spin labeled proteins. We report that in protonated proteins, P(r) is dependent on the length of the second echo period T owing to local environmental effects on the spin-label phase memory relaxation time Tm . For the protein ABD, this effect results in a 1.4 Å increase in the P(r) maximum from T=6 to 20 μs. Protein A has a bimodal P(r) distribution, and the relative height of the shorter distance peak at T=10 μs, the shortest value required to obtain a reliable P(r), is reduced by 40 % relative to that found by extrapolation to T=0. Our results indicate that data at a series of T values are essential for quantitative interpretation of DEER to determine the extent of the T dependence and to extrapolate the results to T=0. Complete deuteration (99 %) of the protein was accompanied by a significant increase in Tm and effectively abolished the P(r) dependence on T., (© 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

In this brief review, we summarize various aspects of NMR paramagnetic relaxation enhancement (PRE). We discuss the types of spin labels used in NMR studies, describe the relevant theory used to accurately calculate PREs from coordinates, including how to take into account the fact that paramagnetic labels tend to be highly mobile and sample a wide range of conformational space, and outline methods to refine structures or ensembles of structures directly against PRE data using simulated annealing. Finally, we show how the PRE can be used to detect, characterize, and visualize sparsely populated states of proteins and their complexes that are invisible to all other biophysical techniques., (© 2015 Elsevier Inc. All rights reserved.)

An outstanding challenge has been to understand the mechanism whereby proteins associate. We report here the results of exhaustively sampling the conformational space in protein-protein association using a physics-based energy function. The agreement between experimental intermolecular paramagnetic relaxation enhancement (PRE) data and the PRE profiles calculated from the docked structures shows that the method captures both specific and non-specific encounter complexes. To explore the energy landscape in the vicinity of the native structure, the nonlinear manifold describing the relative orientation of two solid bodies is projected onto a Euclidean space in which the shape of low energy regions is studied by principal component analysis. Results show that the energy surface is canyon-like, with a smooth funnel within a two dimensional subspace capturing over 75% of the total motion. Thus, proteins tend to associate along preferred pathways, similar to sliding of a protein along DNA in the process of protein-DNA recognition. DOI: http://dx.doi.org/10.7554/eLife.01370.001.

A methyl-detected 'out-and-back' NMR experiment for obtaining simultaneous correlations of methyl resonances of valine and isoleucine/leucine residues with backbone carbonyl chemical shifts, SIM-HMCM(CGCBCA)CO, is described. The developed pulse-scheme serves the purpose of convenience in recording a single data set for all Ile(δ1), Leu(δ) and Val(γ) (ILV) methyl positions instead of acquiring two separate spectra selective for valine or leucine/isoleucine residues. The SIM-HMCM(CGCBCA)CO experiment can be used for ILV methyl assignments in moderately sized protein systems (up to ~100 kDa) where the backbone chemical shifts of (13)C(α), (13)Cβ and (13)CO are known from prior NMR studies and where some losses in sensitivity can be tolerated for the sake of an overall reduction in NMR acquisition time.

Quantitative studies in molecular and structural biology generally require accurate and precise determination of protein concentrations, preferably via a method that is both quick and straightforward to perform. The measurement of ultraviolet absorbance at 280 nm has proven especially useful, since the molar absorptivity (extinction coefficient) at 280 nm can be predicted directly from a protein sequence. This method, however, is only applicable to proteins that contain tryptophan or tyrosine residues. Absorbance at 205 nm, among other wavelengths, has been used as an alternative, although generally using absorptivity values that have to be uniquely calibrated for each protein, or otherwise only roughly estimated. Here, we propose and validate a method for predicting the molar absorptivity of a protein or peptide at 205 nm directly from its amino acid sequence, allowing one to accurately determine the concentrations of proteins that do not contain tyrosine or tryptophan residues. This method is simple to implement, requires no calibration, and should be suitable for a wide range of proteins and peptides., (© 2013 The Protein Society.)

Statistical potentials that embody torsion angle probability densities in databases of high-quality X-ray protein structures supplement the incomplete structural information of experimental nuclear magnetic resonance (NMR) datasets. By biasing the conformational search during the course of structure calculation toward highly populated regions in the database, the resulting protein structures display better validation criteria and accuracy. Here, a new statistical torsion angle potential is developed using adaptive kernel density estimation to extract probability densities from a large database of more than 10⁶ quality-filtered amino acid residues. Incorporated into the Xplor-NIH software package, the new implementation clearly outperforms an older potential, widely used in NMR structure elucidation, in that it exhibits simultaneously smoother and sharper energy surfaces, and results in protein structures with improved conformation, nonbonded atomic interactions, and accuracy., (Copyright © 2012 The Protein Society.)

While an extensive body of knowledge has accumulated on the structures of transcription factors, DNA and their complexes from both NMR and crystallography, much less is known at a molecular level regarding the mechanisms whereby transcription factors locate their specific DNA target site within an overwhelming sea of non-specific DNA sites. Indirect kinetic data suggested that three processes are involved in the search procedure: jumping by dissociation of the protein from the DNA followed by re-association at another site, direct transfer from one DNA molecule or segment to another, and one-dimensional sliding. In this brief perspective I summarize recent NMR developments from our laboratory that have permitted direct characterization of the species and molecular mechanisms involved in the target search process, including the detection of highly transient sparsely-populated states. The main tool in these studies involves the application of paramagnetic relaxation enhancement, supplemented by z-exchange spectroscopy, lineshape analysis and residual dipolar couplings. These studies led to the first direct demonstration of rotation-coupled sliding of a protein along the DNA and the direct transfer of a protein from one DNA molecule to another without dissociating into free solution., (© Springer Science+Business Media B.V. (outside the USA) 2011)

Methyl-transverse relaxation optimized spectroscopy is rapidly becoming the preferred NMR technique for probing structure and dynamics of very large proteins up to ~1 MDa in molecular size. Data interpretation, however, necessitates assignment of methyl groups which still presents a very challenging and time-consuming process. Here we demonstrate that, in combination with a known 3D structure, paramagnetic relaxation enhancement (PRE), induced by nitroxide spin-labels incorporated at only a few surface-exposed engineered cysteines, provides fast, straightforward and robust access to methyl group resonance assignments, including stereoassignments for the methyl groups of leucine and valine. Neither prior assignments, including backbone assignments, for the protein, nor experiments that transfer magnetization between methyl groups and the protein backbone, are required. PRE-derived assignments are refined by 4D methyl-methyl nuclear Overhauser enhancement data, eliminating ambiguities and errors that may arise due to the high sensitivity of PREs to the potential presence of sparsely-populated transient states., (© Springer Science+Business Media B.V. (outside the USA) 2011)

(15)N R(2)/R(1) relaxation data contain information on molecular shape and size as well as on bond vector orientations relative to the diffusion tensor. Since the diffusion tensor can be directly calculated from the molecular coordinates, direct inclusion of (15)N R(2)/R(1) restraints in NMR structure calculations without any a priori assumptions is possible. Here we show that (15)N R(2)/R(1) restraints are particularly valuable when only sparse distance restraints are available. Using three examples of proteins of varying size, namely, GB3 (56 residues), ubiquitin (76 residues), and the N-terminal domain of enzyme I (EIN, 249 residues), we show that incorporation of (15)N R(2)/R(1) restraints results in large and significant increases in coordinate accuracy that can make the difference between being able or unable to determine an approximate global fold. For GB3 and ubiquitin, good coordinate accuracy was obtained using only backbone hydrogen-bond restraints supplemented by (15)N R(2)/R(1) relaxation restraints. For EIN, the global fold could be determined using sparse nuclear Overhauser enhancement (NOE) distance restraints involving only NH and methyl groups in conjunction with (15)N R(2)/R(1) restraints. These results are of practical significance in the study of larger and more complex systems, where the increasing spectral complexity and number of chemical shift degeneracies reduce the number of unambiguous NOE assignments that can be readily obtained, resulting in progressively reduced NOE coverage as the size of the protein increases.

Sparsely populated states of macromolecules, characterized by short lifetimes and high free-energies relative to the predominant ground state, often play a key role in many biological, chemical, and biophysical processes. In this review, we briefly summarize various new developments in NMR spectroscopy that permit these heretofore invisible, sparsely populated states to be detected, characterized, and in some instances visualized. Relaxation dispersion spectroscopy yields detailed kinetic information on processes involving species characterized by distinct chemical shifts with lifetimes in the ∼50 μs-10 ms range and populations as low as 0.5%. In the fast exchange regime (time scale less than ∼250-500 μs), the footprint of sparsely populated states can be observed on paramagnetic relaxation enhancement profiles measured on the resonances of the major species, thereby yielding structural information that is directly related to paramagnetic center-nuclei distances from which it is possible, under suitable circumstances, to compute a structure or ensemble of structures for the minor species. Finally, differential transverse relaxation measurements can be used to detect lifetime broadening effects that directly reflect the unidirectional rates for the conversion of NMR-visible into high-molecular weight NMR-invisible species. Examples of these various approaches are presented., (Copyright © 2010 The Protein Society.)

(15)N relaxation rates contain information on overall molecular shape and size, as well as residue specific orientations of N-H bond vectors relative to the axes of the diffusion tensor. Here we describe a pseudopotential E(relax) that permits direct use of (15)N relaxation rates, in the form of R(2)/R(1) ratios, as experimental restraints in structure calculations without requiring prior information to be extracted from a known molecular structure. The elements of the rotational diffusion tensor are calculated from the atomic coordinates at each step of the structure calculation and then used together with the N-H bond vector orientations to compute the (15)N R(2)/R(1) ratios. We show that the E(relax) term can be reliably used for protein-protein docking of complexes and illustrate its applicability to the 40 kDa complex of the N-terminal domain of enzyme I and the histidine phosphocarrier protein HPr and to the symmetric HIV-1 protease dimer.

A pseudopotential has been introduced into NMR protein structure determination that effectively restrains molecular volume based on the observation that rigid, well-packed protein structures are approximately ellipsoidal. Allowing an ellipsoidal shape is more general than the single approximately spherical shape imposed by the radius of gyration pseudopotential introduced previously. We demonstrate that this new gyration volume term improves structures both during the calculation of an initial unknown fold and during final refinement.

In this paper, we present a series of heteronuclear NMR experiments for the direct observation and characterization of lysine NH3 groups in proteins. In the context of the HoxD9 homeodomain bound specifically to DNA we were able to directly observe three cross-peaks, arising from lysine NH3 groups, with 15N chemical shifts around approximately 33 ppm at pH 5.8 and 35 degrees C. Measurement of water-exchange rates and various types of 15N transverse relaxation rates for these NH3 groups, reveals that rapid water exchange dominates the 15N relaxation for antiphase coherence with respect to 1H through scalar relaxation of the second kind. As a consequence of this phenomenon, 15N line shapes of NH3 signals in a conventional 1H-15N heteronuclear single quantum coherence (HSQC) correlation experiment are much broader than those of backbone amide groups. A 2D 1H-15N correlation experiment that exclusively observes in-phase 15N transverse coherence (termed HISQC for heteronuclear in-phase single quantum coherence spectroscopy) is independent of scalar relaxation in the t(1) (15N) time domain and as a result exhibits strikingly sharper 15N line shapes and higher intensities for NH3 cross-peaks than either HSQC or heteronuclear multiple quantum coherence (HMQC) correlation experiments. Coherence transfer through the relatively small J-coupling between 15Nzeta and 13Cepsilon (4.7-5.0 Hz) can be achieved with high efficiency by maintaining in-phase 15N coherence owing to its slow relaxation. With the use of a suite of triple resonance experiments based on the same design principles as the HISQC, all the NH3 cross-peaks observed in the HISQC spectrum could be assigned to lysines that directly interact with DNA phosphate groups. Selective observation of functional NH3 groups is feasible because of hydrogen bonding or salt bridges that protect them from rapid water exchange. Finally, we consider the potential use of lysine NH3 groups as an alternative probe for larger systems as illustrated by data obtained on the 128-kDa enzyme I dimer.

Due to practical limitations in available 15N rf field strength, imperfections in 15N 180 degrees pulses arising from off-resonance effects can result in significant sensitivity loss, even if the chemical shift offset is relatively small. Indeed, in multi-dimensional NMR experiments optimized for protein backbone amide groups, cross-peaks arising from the Arg guanidino 15Nepsilon (approximately 85 ppm) are highly attenuated by the presence of multiple INEPT transfer steps. To improve the sensitivity for correlations involving Arg Nepsilon-Hepsilon groups, we have incorporated 15N broadband 180 degrees pulses into 3D 15N-separated NOE-HSQC and HNCACB experiments. Two 15N-WURST pulses incorporated at the INEPT transfer steps of the 3D 15N-separated NOE-HSQC pulse sequence resulted in a approximately 1.5-fold increase in sensitivity for the Arg Nepsilon-Hepsilon signals at 800 MHz. For the 3D HNCACB experiment, five 15N Abramovich-Vega pulses were incorporated for broadband inversion and refocusing, and the sensitivity of Arg1Hepsilon-15Nepsilon-13Cgamma/13Cdelta correlation peaks was enhanced by a factor of approximately 1.7 at 500 MHz. These experiments eliminate the necessity for additional experiments to assign Arg 1Hepsilon and 15Nepsilon resonances. In addition, the increased sensitivity afforded for the detection of NOE cross-peaks involving correlations with the 15Nepsilon/1Hepsilon of Arg in 3D 15N-separated NOE experiments should prove to be very useful for structural analysis of interactions involving Arg side-chains.

A large fraction of hydrogens in proteins and nucleic acids is of the methylene type. Their detailed study, however, in terms of structure and dynamics by NMR spectroscopy is hampered by their fast relaxation properties, which give rise to low sensitivity and resolution. It is demonstrated that six different relaxation interference processes, involving 1H-13C and 1H-1H dipolar interactions and 1H and 13C chemical shift anisotropy, can be used simultaneously to mitigate these problems effectively. The approach is applicable to the majority of NMR experiments commonly used to study side chain and backbone conformation. For proteins, its efficiency is evaluated quantitatively for two samples: the third IgG-binding domain from Streptococcal Protein G and the protein calmodulin complexed with a 26-residue target peptide. Gains in both resolution and sensitivity by up to factors of 3.2 and 2.0, respectively, are observed for Gly residues at high magnetic field strengths, but even at much lower fields gains remain substantial. The resolution enhancement obtained for methylene groups makes possible a detailed analysis of spectral regions commonly considered inaccessible due to spectral crowding. For DNA, the high resolution now obtainable for C5' sites permits an H5'/H5''-based sequential NOE assignment procedure, complementary to the conventional base-H1'/H2'/H2'' pathway.

Backbone residual dipolar coupling (N-H, Calpha-Halpha, N-C', and Calpha-C') data collected in five different media on the B3 IgG binding domain of streptococcal protein G (GB3) have been analyzed by simultaneous refinement of the coordinates and optimization of the magnitudes and orientations of the alignment tensors using single and multiple structure representations. We show, using appropriate error analysis, that agreement between observed and calculated dipolar couplings at the level of experimental uncertainty is obtained with a two-structure (N(e) = 2) ensemble representation which represents the simplest equilibrium description of anisotropic motions. The data permit one to determine the magnitude of the anisotropic motions along the four different backbone bond vectors in terms of order parameters. The order parameters, , for the N-H bond vectors are in qualitative agreement with the generalized order parameters, S(2)NH(relaxation), derived from (15)N relaxation measurements, with a correlation coefficient of 0.84. S(2)NH(relaxation) can be regarded as the product of an anisotropic order parameter, corresponding to derived from the residual dipolar couplings, and an axially symmetric order parameter, S(2)NH(axial), corresponding to bond librations which are expected to be essentially uniform along the polypeptide chain. The current data indicate that the average value of S(2)NH(axial) is approximately 0.9. The close correspondence of and S(2)NH(relaxation) indicates that any large-scale displacements from the mean coordinate positions on time scales longer than the rotational correlation time are rare and hence do not perturb the observed dipolar couplings. Analysis of a set of 100 N(e) = 2 ensembles reveals the presence of some long-range correlated motions of N-H and Calpha-Halpha vectors involving residues far apart in the sequence but close together in space. In addition, direct evidence is obtained for ubiquitous crankshaft motions along the entire length of the polypeptide backbone manifested by the anticorrelation of the backbone torsion angles phi(i) and psi(i-1).

The major rate-limiting step in high-throughput NMR protein structure determination involves the calculation of a reliable initial fold, the elimination of incorrect nuclear Overhauser enhancement (NOE) assignments, and the resolution of NOE assignment ambiguities. We present a robust approach to automatically calculate structures with a backbone coordinate accuracy of 1.0-1.5 A from datasets in which as much as 80% of the long-range NOE information (i.e., between residues separated by more than five positions in the sequence) is incorrect. The current algorithm differs from previously published methods in that it has been expressly designed to ensure that the results from successive cycles are not biased by the global fold of structures generated in preceding cycles. Consequently, the method is highly error tolerant and is not easily funnelled down an incorrect path in either three-dimensional structure or NOE assignment space. The algorithm incorporates three main features: a linear energy function representation of the NOE restraints to allow maximization of the number of simultaneously satisfied restraints during the course of simulated annealing; a method for handling the presence of multiple possible assignments for each NOE cross-peak which avoids local minima by treating each possible assignment as if it were an independent restraint; and a probabilistic method to permit both inactivation and reactivation of all NOE restraints on the fly during the course of simulated annealing. NOE restraints are never removed permanently, thereby significantly reducing the likelihood of becoming trapped in a false minimum of NOE assignment space. The effectiveness of the algorithm is demonstrated using completely automatically peak-picked experimental NOE data from two proteins: interleukin-4 (136 residues) and cyanovirin-N (101 residues). The limits of the method are explored using simulated data on the 56-residue B1 domain of Streptococcal protein G.

It is often the case that a substantial number of torsion angles (both backbone and sidechain) in structures of proteins and nucleic acids determined by NMR are found in physically unlikely and energetically unfavorable conformations. We have previously proposed a database-derived potential of mean force comprising one-, two-, three-, and four-dimensional potential surfaces which describe the likelihood of various torsion angle combinations to bias conformational sampling during simulated annealing refinement toward those regions that are populated in very high resolution (< or =1.75 A) crystal structures. We now note a shortcoming of our original implementation of this approach: namely, the forces it places on atoms are very rough. When the density of experimental restraints is low, this roughness can both hinder convergence to commonly populated regions of torsion angle space and reduce overall conformational sampling. In this paper we describe a modification that completely eliminates these problems by replacing the original potential surfaces by a sum of multidimensional Gaussian functions. Structures refined with the new Gaussian implementation now simultaneously enjoy excellent global sampling and excellent local choices of torsion angles.

A simple and rapid method is presented for solving the three-dimensional structures of protein-protein complexes in solution on the basis of experimental NMR restraints that provide the requisite translational (i.e., intermolecular nuclear Overhauser enhancement, NOE, data) and orientational (i.e., backbone (1)H-(15)N dipolar couplings and intermolecular NOEs) information. Providing high-resolution structures of the proteins in the unbound state are available and no significant backbone conformational changes occur upon complexation (which can readily be assessed by analysis of dipolar couplings measured on the complex), accurate and rapid docking of the two proteins can be achieved. The method, which is demonstrated for the 40-kDa complex of enzyme I and the histidine phosphocarrier protein, involves the application of rigid body minimization using a target function comprising only three terms, namely experimental NOE-derived intermolecular interproton distance and dipolar coupling restraints, and a simple intermolecular van der Waals repulsion potential. This approach promises to dramatically reduce the amount of time and effort required to solve the structures of protein-protein complexes by NMR, and to extend the capabilities of NMR to larger protein-protein complexes, possibly up to molecular masses of 100 kDa or more.

Recent advances in multidimensional nuclear magnetic resonance methodology to obtain 1H, 15N and 13C resonance assignments, interproton distance and torsion angle restraints, and restraints that characterize long-range order, coupled with new methods of structure refinement and novel methods for reducing linewidths, have permitted three-dimensional solution structures of single chain proteins in excess of 250 residues and multimeric protein in excess of 40 kDa to be solved. These developments may permit the determination by nuclear magnetic resonance of macromolecular structures up to molecular weights in the 50-60 kDa range, thereby bringing into reach numerous systems of considerable biological interest, including a large variety of protein-protein and protein-nucleic acid complexes.

It has recently been shown that the degree of alignment of macromolecules in an aqueous dilute liquid crystalline medium of bicelles is sufficient to permit accurate values of residual 15N-1H, 13C-1H, and 13Calpha-C' dipolar couplings to be obtained on a routine basis, thereby providing potentially unique long-range structural information. To make use of this information in macromolecular structure determination, the magnitude of the axial and rhombic components of the molecular alignment tensor must be determined. This can be achieved by taking advantage of the fact that different, fixed-distance internuclear vector types are differently distributed relative to the alignment tensor. A histogram of the ensemble of normalized residual dipolar couplings for several such vector types approximates a powder pattern from which the magnitude of the axial and rhombic components are readily extracted in the absence of any prior structural information. The applicability of this method is demonstrated using synthetic data derived from four proteins representative of different sizes, topologies, and secondary structures, and experimental data measured on the small protein ubiquitin., (Copyright 1998 Academic Press.)

Residual dipolar couplings arising from small degrees of alignment of molecules in a magnetic field provide unique long-range structural information. The potential of this approach for structure refinement has recently been demonstrated for a protein-DNA complex in which the magnetic susceptibility tensor was axially symmetric. For most macromolecules and macromolecular complexes, however, axial symmetry cannot be assumed. Moreover, the presence of significant rhombicity will clearly affect the accuracy of the resulting coordinates. In this Communication we present a simple calculational strategy that makes use of simulated annealing refinement against the residual dipolar couplings in combination with a grid search, to simultaneously refine the structures and ascertain the magnitude of the axial and rhombic components of the tensor., (Copyright 1998 Academic Press.)

Recent advances in multidimensional NMR to obtain resonance assignments, interproton distance and torsion angle restraints, and restraints that characterize long range order, coupled with new methods of structure refinement, have permitted solution structures of proteins in excess of 250 residues to be solved.

Here we present a comparative study of the stability of several proteins in H2O and D2O as a function of pH/pH*. We show that the substitution of D2O for H2O leads to an increase in the transition temperature and a decrease in the enthalpy of unfolding. The stability of the proteins, however, appears to be largely unchanged as a result of entropic compensation for the decrease in enthalpy. This enthalpy-entropy compensation is attributed to changes in hydration of proteins in D2O compared to H2O. Analysis of thermodynamic data for the transfer of model compounds from H2O to D2O shows that almost all the changes in the enthalpy of unfolding and in the protein-ligand interactions due to water isotopic substitution can be rationalized by changes in hydration of the buried non-polar groups.

With the advent of multidimensional heteronuclear-edited and -filtered NMR experiments, the field of three-dimensional structure determination by NMR has again increased in scope, making it possible to move the technology beyond the approximately 10 kDa limit inherent to conventional two-dimensional NMR to systems up to potentially 35 to 40 kDa. This article outlines the basic strategies for solving three-dimensional structures of larger systems, in particular, protein complexes and multimeric proteins using three- and four-dimensional NMR spectroscopy, summarizes the key experiments, and illustrates the power of these methods using several examples of protein-DNA, protein-peptide complexes, and oligomeric proteins from the authors' laboratories.

We have examined the 13C alpha and 13C beta chemical shifts of a number of proteins and found that their values at the N-terminal end of a helix provide a good predictor for the presence of a capping box. A capping box consists of a hydrogen-bonded cycle of four amino acids in which the side chain of the N-cap residue forms a hydrogen bond with the backbone amide of the N3 residue, whose side chain in turn may accept a hydrogen bond from the amide of the N-cap residue. The N-cap residue exhibits characteristic values for its backbone torsion angles, with phi and psi clustering around 94 +/- 15 degrees and 167 +/- 5 degrees, respectively. This is manifested by a 1-2 ppm upfield shift of the 13 C alpha resonance and a 1-4 ppm downfield shift of the 13C beta resonance, relative to their random coil values, and is mainly associated with the unusually large value of psi. The residues following the N-cap residue exhibit downfield shifts of 1-3 ppm for the 13C alpha resonances and small upfield shifts for the 13C beta ones, typical of an alpha-helix.

The recent development of a whole panoply of multidimensional heteronuclear-edited and -filtered NMR experiments has revolutionized the field of protein structure determination by NMR, making it possible to extend the methodology from the 10-kDa limit of conventional 2-dimensional NMR to systems up to potentially 35-40 kDa. The basic strategy for solving 3-dimensional structures of larger proteins and protein-ligand complexes in solution using 3- and 4-dimensional NMR spectroscopy is summarized, and the power of these methods is illustrated using 3 examples: interleukin-1 beta, the complex of calmodulin with a target peptide, and the specific complex of the transcription factor GATA-1 with its cognate DNA target site.

Structure determination of macromolecules in solution by nuclear magnetic resonance (NMR) spectroscopy involves the fitting of atomic models to the observed nuclear Overhauser effect (NOE) data. Complete cross-validation has been used to define reliable and unbiased criteria for the quality of solution NMR structures. The method is based on the partitioning of NOE data into test sets and the cross-validation of statistical quantities for each of the test sets. A high correlation between cross-validated measures of fit, such as distance bound violations and NMR R values, and the quality of solution NMR structures was observed. Less complete data resulted in poorer satisfaction of the cross-validated measures of fit. Optimization of cross-validated measures of fit will likely produce solution NMR structures with maximal information content.

The effects of the number, precision and accuracy of interproton distance restraints, of direct refinement against nuclear Overhauser enhancement (NOE) intensities and of the description of the non-bonded contacts on the precision and accuracy of a nuclear magnetic resonance (NMR) protein structure determination have been investigated. The model system employed is the 56 residue immunoglobulin G binding domain of streptococcal protein G. This choice was based on the availability of a very high resolution NMR structure (atomic root-mean-square distribution of the ensemble of 60 calculated structures about the mean co-ordinate positions of 0.25 A for the backbone atoms, 0.65 A for all atoms and 0.39 A for all atoms excluding disordered surface side-chains). The experimental NMR data set for this structure determination comprised a total of 1058 experimental restraints of which 854 were approximate interproton distance restraints corresponding to all the structurally useful NOEs observable for this protein. The calculations presented in this paper reveal the following. (1) The number of interproton distance restraints constitutes the single most important determinant of both precision and accuracy. The ensemble precision and accuracy improves significantly as the number of interproton distance restraints is increased to an average of approximately 15 per residue, of which approximately 60% involve unique proton pairs; subsequent additions of interproton distance restraints, however, lead to less dramatic improvements as information redundancy sets in. (2) The ratio of ensemble precision to ensemble accuracy (which ranges from 0.5 to 0.7 for the backbone atoms) is approximately independent both of the number, precision and accuracy of the interproton distance restraints, and of whether the structures are refined against interproton distance restraints or directly against NOE intensities. (3) In an ensemble of structures generated from a large number of loose approximate interproton distance restraints (an average of approximately 15 restraints per residue with approximately 60% involving unique proton pairs), the interproton distance vectors corresponding to the restraints are very well defined with approximately 80% of vectors between unique proton pairs having a standard deviation of < or = 0.1 A. (4) The accuracy of the mean co-ordinates of an ensemble of structures is significantly higher than the average accuracy of the individual structures comprising the ensemble. For an average ensemble precision of > or = 0.6 A, the dependence of the accuracy of the mean co-ordinates on ensemble precision is approximately linear.(ABSTRACT TRUNCATED AT 400 WORDS)

Joint refinement of macromolecules against crystallographic and nuclear magnetic resonance (NMR) observations is presented as a way of combining experimental information from the two methods. The model of interleukin-1 beta derived by the joint x-ray and NMR refinement is shown to be consistent with the experimental observations of both methods and to have crystallographic R value and geometrical parameters that are of the same quality as or better than those of models obtained by conventional crystallographic studies. The few NMR observations that are violated by the model serve as an indicator for genuine differences between the crystal and solution structures. The joint x-ray-NMR refinement can resolve structural ambiguities encountered in studies of multidomain proteins, in which low- to medium-resolution diffraction data can be complemented by higher resolution NMR data obtained for the individual domains.