NMR residual dipolar couplings as probes of biomolecular dynamics

Variable helix elongation as a tool to modulate RNA alignment and motional couplings

The application of residual dipolar couplings (RDCs) in studies of RNA structure and dynamics can be complicated by the presence of couplings between collective helix motions and overall alignment and by the inability to modulate overall alignment of the molecule by changing the ordering medium. Here, we show for a 27-nt TAR RNA construct that variable levels of helix elongation can be used to alter both overall alignment and couplings to collective helix motions in a semi-predictable manner. In the absence of elongation, a four base-pair helix II capped by a UUCG apical loop exhibits a higher degree of order compared to a six base-pair helix I (theta(I)/theta(II)=0.56+/-0.1). The principal S(zz) direction is nearly parallel to the axis of helix II but deviates by approximately 40 degrees relative to the axis of helix I. Elongating helix I by three base-pairs equalizes the alignment of the two helices and pushes the RNA into the motional coupling limit such that the two helices have comparable degrees of order (theta(I)/theta(II)=0.92+/-0.04) and orientations relative to S(zz) ( approximately 17 degrees ). Increasing the length of elongation further to 22 base-pairs pushes the RNA into the motional decoupling limit in which helix I dominates alignment (theta(II)/theta(I)=0.45+/-0.05), with S(zz) orientated nearly parallel to its helix axis. Many of these trends can be rationalized using PALES simulations that employ a previously proposed three-state dynamic ensemble of TAR. Our results provide new insights into motional couplings, offer guidelines for assessing their extent, and suggest that variable degrees of helix elongation can allow access to independent sets of RDCs for characterizing RNA structural dynamics.

Analysis of the amide (15)N chemical shift tensor of the C(alpha) tetrasubstituted constituent of membrane-active peptaibols, the alpha-aminoisobutyric acid residue, compared to those of di- and tri-substituted proteinogenic amino acid residues

In protein NMR spectroscopy the chemical shift provides important information for the assignment of residues and a first structural evaluation of dihedral angles. Furthermore, angular restraints are obtained from oriented samples by solution and solid-state NMR spectroscopic approaches. Whereas the anisotropy of chemical shifts, quadrupolar couplings and dipolar interactions have been used to determine the structure, dynamics and topology of oriented membrane polypeptides using solid-state NMR spectroscopy similar concepts have been introduced to solution NMR through the measurements of residual dipolar couplings. The analysis of (15)N chemical shift spectra depends on the accuracy of the chemical shift tensors. When investigating alamethicin and other peptaibols, i.e. polypeptides rich in alpha-aminoisobutyric acid (Aib), the (15)N chemical shift tensor of this C(alpha)-tetrasubstituted amino acid exhibits pronounced differences when compared to glycine, alanine and other proteinogenic residues. Here we present an experimental investigation on the (15)N amide Aib tensor of N-acetyl-Aib-OH and for the Aib residues within peptaibols. Furthermore, a statistical analysis of the tensors published for di- (glycine) and tri-substituted residues has been performed, where for the first time the published data sets are compiled using a common reference. The size of the isotropic chemical shift and main tensor elements follows the order di- < tri- < tetra-substituted amino acids. A (15)N chemical shift-(1)H-(15)N dipolar coupling correlation NMR spectrum of alamethicin is used to evaluate the consequences of variations in the main tensor elements for the structural analysis of this membrane peptide.

Observing biological dynamics at atomic resolution using NMR

Biological macromolecules are highly flexible and continually undergo conformational fluctuations on a broad spectrum of timescales. It has long been recognized that dynamics have an important role in the action of these molecules. However, the relationship between molecular function and motion is extremely challenging to delineate, because the conformational space available to macromolecules is vast and the relevant excursions can be infrequent and short-lived. Recent advances in solution nuclear magnetic resonance (NMR) spectroscopy permit biomolecular dynamics to be observed with unprecedented detail. Applications of these new NMR techniques to the study of fundamental processes such as binding and catalysis have provided new insights into how living systems operate at an atomic level.

Functional dynamics of the folded ankyrin repeats of I kappa B alpha revealed by nuclear magnetic resonance

Inhibition of nuclear factor kappaB (NF-kappaB) is mainly accomplished by IkappaB alpha, which consists of a signal response sequence at the N-terminus, a six-ankyrin repeat domain (ARD) that binds NF-kappaB, and a C-terminal PEST sequence. Previous studies with the ARD revealed that the fifth and sixth repeats are only partially folded in the absence of NF-kappaB. Here we report NMR studies of a truncated version of IkappaB alpha, containing only the first four ankyrin repeats, IkappaB alpha(67-206). This four-repeat segment is well-structured in the free state, enabling full resonance assignments to be made. H-D exchange, backbone dynamics, and residual dipolar coupling (RDC) experiments reveal regions of flexibility. In addition, regions consistent with the presence of micro- to millisecond motions occur periodically throughout the repeat structure. Comparison of the RDCs with the crystal structure gave only moderate agreement, but an ensemble of structures generated by accelerated molecular dynamics gave much better agreement with the measured RDCs. The regions showing flexibility correspond to those implicated in entropic compensation for the loss of flexibility in ankyrin repeats 5 and 6 upon binding to NF-kappaB. The regions showing micro- to millisecond motions in the free protein are the ends of the beta-hairpins that directly interact with NF-kappaB in the complex.

Tilt and azimuthal angles of a transmembrane peptide: a comparison between molecular dynamics calculations and solid-state NMR data of sarcolipin in lipid membranes

We report molecular dynamics simulations in the explicit membrane environment of a small membrane-embedded protein, sarcolipin, which regulates the sarcoplasmic reticulum Ca-ATPase activity in both cardiac and skeletal muscle. In its monomeric form, we found that sarcolipin adopts a helical conformation, with a computed average tilt angle of 28 +/- 6 degrees and azymuthal angles of 66 +/- 22 degrees, in reasonable accord with angles determined experimentally (23 +/- 2 degrees and 50 +/- 4 degrees, respectively) using solid-state NMR with separated-local-field experiments. The effects of time and spatial averaging on both (15)N chemical shift anisotropy and (1)H/(15)N dipolar couplings have been analyzed using short-time averages of fast amide out-of-plane motions and following principal component dynamic trajectories. We found that it is possible to reproduce the regular oscillatory patterns observed for the anisotropic NMR parameters (i.e., PISA wheels) employing average amide vectors. This work highlights the role of molecular dynamics simulations as a tool for the analysis and interpretation of solid-state NMR data.

Effect of pseudorepeat rearrangement on alpha-synuclein misfolding, vesicle binding, and micelle binding

The pathological and physiological hallmarks of the protein alpha-synuclein (aS) are its misfolding into cytotoxic aggregates and its binding to synaptic vesicles, respectively. Both events are mediated by seven 11-residue amphiphilic pseudorepeats and, most generally, involve a transition from intrinsically unstructured conformations to structured conformations. Based on aS interactions with aggregation-inhibiting small molecules, an aS variant termed shuffled alpha-synuclein (SaS), wherein the first six pseudorepeats had been rearranged, was introduced. Here, the effects of this rearrangement on misfolding, vesicle binding, and micelle binding are examined in reference to aS and beta-synuclein to study the sequence characteristics underlying these processes. Fibrillization correlates with the distinct clustering of residues with high beta-sheet propensities, while vesicle affinities depend on the mode of pseudorepeat interchange and loss. In the presence of micelles, the pseudorepeat region of SaS adopts an essentially continuous helix, whereas aS and beta-synuclein encounter a distinct helix break, indicating that a more homogeneous distribution of surfactant affinities in SaS prevented the formation of an extensive helix break in the micelle-bound state. By demonstrating the importance of the distribution of beta-sheet propensities and by revealing inhomogeneous aS surfactant affinities, the present study provides novel insights into two central themes of synuclein biology.

Ultrahigh resolution characterization of domain motions and correlations by multialignment and multireference residual dipolar coupling NMR

Nuclear magnetic resonance (NMR) residual dipolar couplings (RDCs) provide a unique opportunity for spatially characterizing complex motions in biomolecules with time scale sensitivity extending up to milliseconds. Up to five motionally averaged Wigner rotation elements, (D(0k)2(alphaalpha)), can be determined experimentally using RDCs measured in five linearly independent alignment conditions and applied to define motions of axially symmetric bond vectors. Here, we show that up to 25 motionally averaged Wigner rotation elements, (D(mk)2(alphabetagamma)), can be determined experimentally from multialignment RDCs and used to characterize rigid-body motions of chiral domains. The 25 (D(mk)2(alphabetagamma)) elements form a basis set that allows one to measure motions of a domain relative to an isotropic distribution of reference frames anchored on a second domain (and vice versa), thus expanding the 3D spatial resolution with which motions can be characterized. The 25 (D(mk)2(alphabetagamma)) elements can also be used to fit an ensemble consisting of up to eight equally or six unequally populated states. For more than two domains, changing the identity of the domain governing alignment allows access to new information regarding the correlated nature of the domain fluctuations. Example simulations are provided that validate the theoretical derivation and illustrate the high spatial resolution with which rigid-body domain motions can be characterized using multialignment and multireference RDCs. Our results further motivate the development of experimental approaches for both modulating alignment and anchoring it on specifically targeted domains.

Intermolecular dynamics studied by paramagnetic tagging

Yeast cytochrome c and bovine adrenodoxin form a dynamic electron transfer complex, which is a pure encounter complex. It is demonstrated that the dynamic nature of the interaction can readily be probed by using a rigid lanthanide tag attached to cytochrome c. The tag, Caged Lanthanide NMR Probe 5, induces pseudocontact shifts and residual dipolar couplings and does not perturb the binding interface. Due to the dynamics in the complex, residual dipolar couplings in adrenodoxin are very small. Simulation shows that cytochrome c needs to sample a large part of the surface of adrenodoxin to explain the small degree of alignment observed for adrenodoxin. The applied method provides a simple and straightforward way to observe dynamics in protein complexes or domain-domain mobility without the need for external alignment media.

Improved accuracy of 15N-1H scalar and residual dipolar couplings from gradient-enhanced IPAP-HSQC experiments on protonated proteins

The presence of dipole-dipole cross-correlated relaxation as well as unresolved E.COSY effects adversely impacts the accuracy of (1)J(NH) splittings measured from gradient-enhanced IPAP-HSQC spectra. For isotropic samples, the size of the systematic errors caused by these effects depends on the values of (2)J(NHalpha), (3)J(NHbeta) and (3)J(HNHalpha). Insertion of band-selective (1)H decoupling pulses in the IPAP-HSQC experiment eliminates these systematic errors and for the protein GB3 yields (1)J(NH) splittings that agree to within a root-mean-square difference of 0.04 Hz with values measured for perdeuterated GB3. Accuracy of the method is also highlighted by a good fit to the GB3 structure of the (1)H-(15)N RDCs extracted from the minute differences in (1)J(NH) splitting measured at 500 and 750 MHz (1)H frequencies, resulting from magnetic susceptibility anisotropy. A nearly complete set of (2)J(NHalpha) couplings was measured in GB3 in order to evaluate whether the impact of cross-correlated relaxation is dominated by the (15)N-(1)H(alpha) or (15)N-(1)H(beta) dipolar interaction. As expected, we find that (2)J(NHalpha) < or = 2 Hz, with values in the alpha-helix (0.86 +/- 0.52 Hz) slightly larger than in beta-sheet (0.66 +/- 0.26 Hz). Results indicate that under isotropic conditions, N-H(N)/N-H(beta) cross-correlated relaxation often dominates. Unresolved E.COSY effects under isotropic conditions involve (3)J(HNHalpha) and J(NHalpha), but when weakly aligned any aliphatic proton proximate to both N and H(N) can contribute.

A tool for the prediction of structures of complex sugars

In two recent back to back articles(Xia et al., J Chem Theory Comput 3:1620-1628 and 1629-1643, 2007a, b) we have started to address the problem of complex oligosaccharide conformation and folding. The scheme previously presented was based on exhaustive searches in configuration space in conjunction with Nuclear Overhauser Effect (NOE) calculations and the use of a complex rotameric library that takes branching into account. NOEs are extremely useful for structural determination but only provide information about short range interactions and ordering. Instead, the measurement of residual dipolar couplings (RDC), yields information about molecular ordering or folding that is long range in nature. In this article we show the results obtained by incorporation RDC calculations into our prediction scheme. Using this new approach we are able to accurately predict the structure of six human milk sugars: LNF-1, LND-1, LNF-2, LNF-3, LNnT and LNT. Our exhaustive search in dihedral configuration space combined with RDC and NOE calculations allows for highly accurate structural predictions that, because of the non-ergodic nature of these molecules on a time scale compatible with molecular dynamics simulations, are extremely hard to obtain otherwise (Almond et al., Biochemistry 43:5853-5863, 2004). Molecular dynamics simulations in explicit solvent using as initial configurations the structures predicted by our algorithm show that the histo-blood group epitopes in these sugars are relatively rigid and that the whole family of oligosaccharides derives its conformational variability almost exclusively from their common linkage (beta-D: -GlcNAc-(1-->3)-beta-D: -Gal) which can exist in two distinct conformational states. A population analysis based on the conformational variability of this flexible glycosidic link indicates that the relative population of the two distinct states varies for different human milk oligosaccharides.

16-fold degeneracy of peptide plane orientations from residual dipolar couplings: analytical treatment and implications for protein structure determination

Residual dipolar couplings (RDCs) measured for internally rigid molecular fragments provide important information about the relative orientations of these fragments. Dependent on the symmetry of the alignment tensor and the symmetry of the molecular fragment, however, there generally exist more than one solution for the fragment orientation consistent with the measured RDCs. Analytical solutions are presented that describe the complete set of orientations of internally rigid fragments that are consistent with multiple dipolar couplings measured in a single alignment medium that is rhombic. For the first time, it is shown that, for a planar fragment such as the peptide plane, there generally exist 16 different solutions with their analytical expressions presented explicitly. The presence of these solutions is shown to be highly relevant for standard structure determination protocols using RDCs to refine molecular structures. In particular, when using standard protein structure refinement with RDCs that were measured in a single alignment medium as constraints, it is found that often more than one of the peptide plane solutions is physically viable; i.e., despite being consistent with measured RDCs, the local backbone structure can be incorrect. On the basis of experimental and simulated examples, it is rationalized why protein structures that are refined against RDCs measured in a single medium can have lower resolution (precision) than one would expect on the basis of the experimental accuracy of the RDCs. Conditions are discussed under which the correct solution can be identified.

Slow motions in chicken villin headpiece subdomain probed by cross-correlated NMR relaxation of amide NH bonds in successive residues

The villin headpiece subdomain (HP36) is a widely used system for protein-folding studies. Nuclear magnetic resonance cross-correlated relaxation rates arising from correlated fluctuations of two N-H(N) dipole-dipole interactions involving successive residues were measured at two temperatures at which HP36 is at least 99% folded. The experiment revealed the presence of motions slower than overall tumbling of the molecule. Based on the theoretical analysis of the spectral densities we show that the structural and dynamic contributions to the experimental cross-correlated relaxation rate can be separated under certain conditions. As a result, dynamic cross-correlated order parameters describing slow microsecond-to-millisecond motions of N-H bonds in neighboring residues can be introduced for any extent of correlations in the fluctuations of the two bond vectors. These dynamic cross-correlated order parameters have been extracted for HP36. The comparison of their values at two different temperatures indicates that when the temperature is raised, slow motions increase in amplitude. The increased amplitude of these fluctuations may reflect the presence of processes directly preceding the unfolding of the protein.

Self-consistent residual dipolar coupling based model-free analysis for the robust determination of nanosecond to microsecond protein dynamics

Residual dipolar couplings (RDCs) provide information about the dynamic average orientation of inter-nuclear vectors and amplitudes of motion up to milliseconds. They complement relaxation methods, especially on a time-scale window that we have called supra-tau(c) (tau(c) < supra-tau(c) < 50 micros). Here we present a robust approach called Self-Consistent RDC-based Model-free analysis (SCRM) that delivers RDC-based order parameters-independent of the details of the structure used for alignment tensor calculation-as well as the dynamic average orientation of the inter-nuclear vectors in the protein structure in a self-consistent manner. For ubiquitin, the SCRM analysis yields an average RDC-derived order parameter of the NH vectors <S2(rdc)>0.72 +/- 0.02 compared to <S2(LS)> = 0.778 +/- 0.003 for the Lipari-Szabo order parameters, indicating that the inclusion of the supra-tau(c) window increases the averaged amplitude of mobility observed in the sub-supra-tau(c) window by about 34%. For the beta-strand spanned by residues Lys48 to Leu50, an alternating pattern of backbone NH RDC order parameter S2(rdc)(NH) = (0.59, 0.72, 0.59) was extracted. The backbone of Lys48, whose side chain is known to be involved in the poly-ubiquitylation process that leads to protein degradation, is very mobile on the supra-tau(c) time scale (S2(rdc)(NH) = 0.59 +/- 0.03), while it is inconspicuous (S2(LS)(NH)= 0.82) on the sub-tau(c) as well as on micros-ms relaxation dispersion time scales. The results of this work differ from previous RDC dynamics studies of ubiquitin in the sense that the results are essentially independent of structural noise providing a much more robust assessment of dynamic effects that underlie the RDC data.

De novo determination of internuclear vector orientations from residual dipolar couplings measured in three independent alignment media

The straightforward interpretation of solution state residual dipolar couplings (RDCs) in terms of internuclear vector orientations generally requires prior knowledge of the alignment tensor, which in turn is normally estimated using a structural model. We have developed a protocol which allows the requirement for prior structural knowledge to be dispensed with as long as RDC measurements can be made in three independent alignment media. This approach, called Rigid Structure from Dipolar Couplings (RSDC), allows vector orientations and alignment tensors to be determined de novo from just three independent sets of RDCs. It is shown that complications arising from the existence of multiple solutions can be overcome by careful consideration of alignment tensor magnitudes in addition to the agreement between measured and calculated RDCs. Extensive simulations as well applications to the proteins ubiquitin and Staphylococcal protein GB1 demonstrate that this method can provide robust determinations of alignment tensors and amide N-H bond orientations often with better than 10 degrees accuracy, even in the presence of modest levels of internal dynamics.

Solid-state nuclear magnetic resonance spectroscopy studies of furanose ring dynamics in the DNA HhaI binding site

The dynamics of the furanose rings in the GCGC moiety of the DNA oligomer [d(G 1A 2T 3A 4 G 5 C 6 G 7 C 8T 9A 10T 11C 12)] 2 are studied by using deuterium solid-state NMR (SSNMR). SSNMR spectra obtained from DNAs selectively deuterated on the furanose rings of nucleotides within the 5'-GCGC-3' moiety indicated that all of these positions are structurally flexible. The furanose ring within the deoxycytidine that is the methylation target displays the largest-amplitude structural changes according to the observed deuterium NMR line shapes, whereas the furanose rings of nucleotides more remote from the methylation site have less-mobile furanose rings (i.e., with puckering amplitudes < 0.3 A). Previous work has shown that methylation reduces the amplitude of motion in the phosphodiester backbone of the same DNA, and our observations indicate that methylation perturbs backbone dynamics through the furanose ring. These NMR data indicate that the 5'-GCGC-3' is dynamic, with the largest-amplitude motions occurring nearest the methylation site. The inherent flexibility of this moiety in DNA makes the molecule more amenable to the large-amplitude structural rearrangements that must occur when the DNA binds to the HhaI methyltransferase.

Protein backbone motions viewed by intraresidue and sequential HN-Halpha residual dipolar couplings

Triple resonance E.COSY-based techniques were used to measure intra-residue and sequential H(N)-H(alpha) residual dipolar couplings (RDCs) for the third IgG-binding domain of protein G (GB3), aligned in Pf1 medium. Measurements closely correlate with values predicted on the basis of an NMR structure, previously determined on the basis of a large number of one-bond backbone RDCs measured in five alignment media. However, in particular the sequential H(N)-H(alpha) RDCs are smaller than predicted for a static structure, suggesting a degree of motion for these internuclear vectors that exceeds that of the backbone amide N-H vectors. Of all experimentally determined GB3 structures available, the best correlation between experimental (1)H-(1)H couplings is observed for a GB3 ensemble, previously derived to generate a realistic picture of the conformational space sampled by GB3 (Clore and Schwieters, J Mol Biol 355:879-886, 2006). However, for both NMR and X-ray-derived structures the (1)H-(1)H couplings are found to be systematically smaller than expected on the basis of alignment tensors derived from (15)N-(1)H amide RDCs, assuming librationally corrected N-H bond lengths of 1.041 A.

Simultaneous NMR study of protein structure and dynamics using conservative mutagenesis

A novel iterative procedure is described that allows both the orientation and dynamics of internuclear bond vectors to be determined from direct interpretation of NMR dipolar couplings, measured under at least three orthogonal alignment conditions. If five orthogonal alignments are available, the approach also yields information on the degree of motional anisotropy and the direction in which the largest amplitude internal motion of each bond vector takes place. The method is demonstrated for the backbone (15)N-(1)H, (13)C(alpha)-(1)H(alpha), and (13)C(alpha)-13C' interactions in the previously well-studied protein domain GB3, dissolved in a liquid crystalline suspension of filamentous phage Pf1. Alignment variation is achieved by using conservative mutations of charged surface residues. Results indicate remarkably uniform backbone dynamics, with amplitudes that agree well with those of previous (15)N relaxation studies for most residues involved in elements of secondary structure, but larger amplitude dynamics than those found by (15)N relaxation for residues in loop and turn regions. In agreement with a previous analysis of dipolar couplings, the N-H bonds in the second beta-strand, which is involved in antibody recognition, show elevated dynamics with largest amplitudes orthogonal to the chain direction.

Structure and dynamics of Ca2+-binding domain 1 of the Na+/Ca2+ exchanger in the presence and in the absence of Ca2+

The Na(+)/Ca(2+) exchanger is the major exporter of Ca(2+) across the cell membrane of cardiomyocytes. The activity of the exchanger is regulated by a large intracellular loop that contains two Ca(2+)-binding domains, calcium-binding domain (CBD) 1 and CBD2. CBD1 binds Ca(2+) with much higher affinity than CBD2 and is considered to be the primary Ca(2+) sensor. The effect of Ca(2+) on the structure and dynamics of CBD1 has been characterized by NMR spectroscopy using chemical shifts, residual dipolar couplings, and spin relaxation. Residual dipolar couplings are used in a new way for residue selection in the determination of the anisotropic rotational diffusion tensor from spin relaxation data. The results provide a highly consistent description across these complementary data sets and show that Ca(2+) binding is accompanied by a selective conformational change among the binding site residues. Residues that exhibit a significant conformational change are also sites of altered dynamics. In particular, Ca(2+) binding restricts the mobility of the major acidic segment and affects the dynamics of several nearby binding loops. These observations indicate that Ca(2+) elicits a local transition to a well-ordered coordination geometry in the CBD1-binding site.

Lanthanide-binding peptides for NMR measurements of residual dipolar couplings and paramagnetic effects from multiple angles

Lanthanide-binding peptide tags (LBTs) containing a single cysteine residue can be attached to proteins via a disulfide bond, presenting a flexible means of tagging proteins site-specifically with a lanthanide ion. Here we show that cysteine residues placed in different positions of the LBT can be used to expose the protein to different orientations of the magnetic susceptibility anisotropy (delta chi) tensor and to generate different molecular alignments in a magnetic field. Delta chi tensors determined by nuclear magnetic resonance (NMR) spectroscopy for LBT complexes with Yb3+, Tm3+, and Er3+ suggest a rational way of producing alignment tensors with different orientations. In addition, knowledge of the delta chi tensor of LBT allows modeling of the protein-LBT structures. Despite evidence for residual mobility of the LBTs with respect to the protein, the pseudocontact shifts and residual dipolar couplings displayed by proteins disulfide-bonded to LBTs are greater than those achievable with most other lanthanide binding tags.

Influence of the fluctuations of the alignment tensor on the analysis of the structure and dynamics of proteins using residual dipolar couplings

It has been suggested that the fluctuations of the alignment tensor can affect the results of procedures for characterizing the structure and the dynamics of proteins using residual dipolar couplings. We show here that the very significant fluctuations of the steric alignment tensor caused by the dynamics of proteins can be safely ignored when they do not correlate with those of the bond vectors. A detailed analysis of these correlations in the protein ubiquitin reveals that their effects are negligible for the analysis of backbone motions within secondary structure elements, but also that they may be significant in turns, loops and side chains, especially for bond vectors that have small residual dipolar couplings. Our results suggest that methods that explicitly consider the motions of the alignment tensor will be needed to study the large-scale structural fluctuations that take place on the millisecond timescale, which are often important for the biological function of proteins, from residual dipolar coupling measurements.

Application of structural dynamic approaches provide novel insights into the enzymatic mechanism of the tumor necrosis factor-alpha-converting enzyme

Zinc dependent metalloproteinases comprise a large family of structurally homologous enzymes with a wide variety of biological roles. Originally described as proteinases involved in extracellular matrix (ECM) catabolism, these enzymes were later found to serve major roles as initiators of signaling pathways in many aspects of biology, ranging from cell proliferation, differentiation and communication, to pathological states associated with tumor metastasis, inflammation, tissue degeneration and cell death. From these enzymes, the tumor necrosis factor-α converting enzyme (TACE) stands out as a central shedding activity mediating the regulated release of a host of cytokines, receptors and other cell surface molecules. Selective drugs targeted at blocking TACE for treatment of rheumatoid arthritis and other disease indications are highly sought. Yet, the structural and chemical knowledge underlying its enzymatic activity is very limited. This is in part due to the fact that the catalytic zinc atom of metalloproteinases is usually spectroscopically silent and hence difficult to study using conventional spectroscopic and analytical tools. Most structural and biochemical studies, as well as medicinal chemistry efforts carried out so far were limited to non-dynamic structure/function characterization. Thus, to date, our mechanistic knowledge comes from theoretical calculations derived from static crystal structures from family members that are highly similar in their amino acid sequence and three-dimensional structure.

This review introduces the importance of real-time quantification of biophysical properties and structural kinetic behavior applied to the study of TACE and other zinc metalloproteinases to dissect their molecular mechanisms. The molecular details that link the catalytic chemistry to key kinetic, electronic and structural events have remained elusive because of the difficulties associated with probing time-dependent structure-function aspects of enzymatic reactions. Here we discuss the use of conventional and real-time structural-spectroscopic tools to study the reactive metal site during catalysis, and initial lessons on the enzymatic mechanism that we are learning. Approaches such as the ones presented here may be useful in the design of specific inhibitors as drug candidates.

A large data set comparison of protein structures determined by crystallography and NMR: statistical test for structural differences and the effect of crystal packing

The existence of a large number of proteins for which both nuclear magnetic resonance (NMR) and X-ray crystallographic coordinates have been deposited into the Protein Data Bank (PDB) makes the statistical comparison of the corresponding crystal and NMR structural models over a large data set possible, and facilitates the study of the effect of the crystal environment and other factors on structure. We present an approach for detecting statistically significant structural differences between crystal and NMR structural models which is based on structural superposition and the analysis of the distributions of atomic positions relative to a mean structure. We apply this to a set of 148 protein structure pairs (crystal vs NMR), and analyze the results in terms of methodological and physical sources of structural difference. For every one of the 148 structure pairs, the backbone root-mean-square distance (RMSD) over core atoms of the crystal structure to the mean NMR structure is larger than the average RMSD of the members of the NMR ensemble to the mean, with 76% of the structure pairs having an RMSD of the crystal structure to the mean more than a factor of two larger than the average RMSD of the NMR ensemble. On average, the backbone RMSD over core atoms of crystal structure to the mean NMR is approximately 1 A. If non-core atoms are included, this increases to 1.4 A due to the presence of variability in loops and similar regions of the protein. The observed structural differences are only weakly correlated with the age and quality of the structural model and differences in conditions under which the models were determined. We examine steric clashes when a putative crystalline lattice is constructed using a representative NMR structure, and find that repulsive crystal packing plays a minor role in the observed differences between crystal and NMR structures. The observed structural differences likely have a combination of physical and methodological causes. Stabilizing attractive interactions arising from intermolecular crystal contacts which shift the equilibrium of the crystal structure relative to the NMR structure is a likely physical source which can account for some of the observed differences. Methodological sources of apparent structural difference include insufficient sampling or other issues which could give rise to errors in the estimates of the precision and/or accuracy.

NMR identification of transient complexes critical to adenylate kinase catalysis

A fundamental question in protein chemistry is how the native energy landscape of enzymes enables efficient catalysis of chemical reactions. Adenylate kinase is a small monomeric enzyme that catalyzes the reversible conversion of AMP and ATP into two ADP molecules. Previous structural studies have revealed that substrate binding is accompanied by large rate-limiting spatial displacements of both the ATP and AMP binding motifs. In this report a solution-state NMR approach was used to probe the native energy landscape of adenylate kinase in its free form, in complex with its natural substrates, and in the presence of a tight binding inhibitor. Binding of ATP induces a dynamic equilibrium in which the ATP binding motif populates both the open and the closed conformations with almost equal populations. A similar scenario is observed for AMP binding, which induces an equilibrium between open and closed conformations of the AMP binding motif. These ATP- and AMP-bound structural ensembles represent complexes that exist transiently during catalysis. Simultaneous binding of AMP and ATP is required to force both substrate binding motifs to close cooperatively. In addition, a previously unknown unidirectional energetic coupling between the ATP and AMP binding sites was discovered. On the basis of these and previous results, we propose that adenylate kinase belongs to a group of enzymes whose substrates act to shift pre-existing equilibria toward catalytically active states.

Magnetic field induced residual dipolar couplings of imino groups in nucleic acids from measurements at a single magnetic field

For base-paired nucleic acids, variations in (1) J (NH) and the imino (1)H chemical shift are both dominated by hydrogen bond length. In the absence of molecular alignment, the (1) J (NH) coupling for the imino proton then can be approximated by (1) J (NH) = (1.21Hz/ppm)delta(H) - 103.5 +/- 0.6 Hz, where delta(H) represents the chemical shift of the imino proton in ppm. This relation permits imino residual dipolar couplings (RDCs) resulting from magnetic susceptibility anisotropy (MSA) to be extracted from measurement of ((1) J (NH) + RDC) splittings at a single magnetic field strength. Magnetic field-induced RDCs were measured for tRNA(Val) and the alignment tensor determined from magnetic-field alignment of tRNA(Val) agrees well with the tensor calculated by summation of the MSA tensors of the individual nucleobases.

NMR studies of RNA dynamics and structural plasticity using NMR residual dipolar couplings

An increasing number of RNAs are being discovered that perform their functions by undergoing large changes in conformation in response to a variety of cellular signals, including recognition of proteins and small molecular targets, changes in temperature, and RNA synthesis itself. The measurement of NMR residual dipolar couplings (RDCs) in partially aligned systems is providing new insights into the structural plasticity of RNA through combined characterization of large-amplitude collective helix motions and local flexibility in noncanonical regions over a wide window of biologically relevant timescales (<milliseconds). Here, we review RDC methodology for studying RNA structural dynamics and survey what has been learnt thus far from application of these methods. Future methodological challenges are also identified.

NMR studies of dynamics in RNA and DNA by 13C relaxation

RNA and DNA molecules experience motions on a wide range of time scales, ranging from rapid localized motions to much slower collective motions of entire helical domains. The many functions of RNA in biology very often require this molecule to change its conformation in response to biological signals in the form of small molecules, proteins or other nucleic acids, whereas local motions in DNA may facilitate protein recognition and allow enzymes acting on DNA to access functional groups on the bases that would otherwise be buried in Watson-Crick base pairs. Although these statements make a compelling case to study the sequence dependent dynamics in nucleic acids, there are few residue-specific studies of nucleic acid dynamics. Fortunately, NMR studies of dynamics of nucleic acids and nucleic acids-protein complexes are gaining increased attention. The aim of this review is to provide an update of the recent progress in studies of nucleic acid dynamics by NMR based on the application of solution relaxation techniques.

Molecular crowding inhibits intramolecular breathing motions in proteins

In aqueous solution some proteins undergo large-scale movements of secondary structures, subunits or domains, referred to as protein "breathing", that define a native-state ensemble of structures. These fluctuations are sensitive to the nature and concentration of solutes and other proteins and are thereby expected to be different in the crowded interior of a cell than in dilute solution. Here we use a combination of wide angle X-ray scattering (WAXS) and computational modeling to derive a quantitative measure of the spatial scale of conformational fluctuations in a protein solution. Concentration-dependent changes in the observed scattering intensities are consistent with a model of structural fluctuations in which secondary structures undergo rigid-body motions relative to one another. This motion increases with decreasing protein concentration or increasing temperature. Analysis of a set of five structurally and functionally diverse proteins reveals a diversity of kinetic behaviors. Proteins with multiple disulfide bonds exhibit little or no increase in breathing in dilute solutions. The spatial extent of structural fluctuations appears highly dependent on both protein structure and concentration and is universally suppressed at very high protein concentrations.

Characterizing the relative orientation and dynamics of RNA A-form helices using NMR residual dipolar couplings

We present a protocol for determining the relative orientation and dynamics of A-form helices in 13C/15N isotopically enriched RNA samples using NMR residual dipolar couplings (RDCs). Non-terminal Watson-Crick base pairs in helical stems are experimentally identified using NOE and trans-hydrogen bond connectivity and modeled using the idealized A-form helix geometry. RDCs measured in the partially aligned RNA are used to compute order tensors describing average alignment of each helix relative to the applied magnetic field. The order tensors are translated into Euler angles defining the average relative orientation of helices and order parameters describing the amplitude and asymmetry of interhelix motions. The protocol does not require complete resonance assignments and therefore can be implemented rapidly to RNAs much larger than those for which complete high-resolution NMR structure determination is feasible. The protocol is particularly valuable for exploring adaptive changes in RNA conformation that occur in response to biologically relevant signals. Following resonance assignments, the procedure is expected to take no more than 2 weeks of acquisition and data analysis time.

Top-down approach in protein RDC data analysis: de novo estimation of the alignment tensor

In solution NMR spectroscopy the residual dipolar coupling (RDC) is invaluable in improving both the precision and accuracy of NMR structures during their structural refinement. The RDC also provides a potential to determine protein structure de novo. These procedures are only effective when an accurate estimate of the alignment tensor has already been made. Here we present a top-down approach, starting from the secondary structure elements and finishing at the residue level, for RDC data analysis in order to obtain a better estimate of the alignment tensor. Using only the RDCs from N-H bonds of residues in alpha-helices and CA-CO bonds in beta-strands, we are able to determine the offset and the approximate amplitude of the RDC modulation-curve for each secondary structure element, which are subsequently used as targets for global minimization. The alignment order parameters and the orientation of the major principal axis of individual helix or strand, with respect to the alignment frame, can be determined in each of the eight quadrants of a sphere. The following minimization against RDC of all residues within the helix or strand segment can be carried out with fixed alignment order parameters to improve the accuracy of the orientation. For a helical protein Bax, the three components A(xx), A(yy) and A(zz), of the alignment order can be determined with this method in average to within 2.3% deviation from the values calculated with the available atomic coordinates. Similarly for beta-sheet protein Ubiquitin they agree in average to within 8.5%. The larger discrepancy in beta-strand parameters comes from both the diversity of the beta-sheet structure and the lower precision of CA-CO RDCs. This top-down approach is a robust method for alignment tensor estimation and also holds a promise for providing a protein topological fold using limited sets of RDCs.

Explicit treatment of force contribution from alignment tensor using overdetermined linear equations and its application in NMR structure determination

Residual dipolar coupling (RDC) provides valuable information about the orientation of each internuclear vector in a macromolecule with respect to the static magnetic field. However, structure determination utilizing RDC still remains challenging without additional restraints such as NOE. In this context, a novel approach has been developed to efficiently extract structural information from RDC by successive application of singular value decomposition (SVD) method in the course of NMR structure determination. Force contribution from the alignment tensor is rigorously formulated in the context of SVD, and assessments have been made to verify its numerical accuracy. The efficacy of this approach is illustrated by showing that RDC restraints alone can restore a distorted beta-hairpin to native-like structure using the replica-exchange molecular dynamics simulations.

SESAME-HSQC for simultaneous measurement of NH and CH scalar and residual dipolar couplings

We present a novel pulse sequence, SESAME-HSQC, for the simultaneous measurement of several NH and CH scalar and residual dipolar couplings in double labeled proteins. The proposed Spin-statE Selective All Multiplicity Edited (SESAME)-HSQC combines gradient selected and sensitivity enhanced (15)N- and constant-time (13)C-HSQC experiments with the recently introduced spin-state selective method (Nolis et al., J. Magn. Reson. 180 (2006) 39-50) for measuring couplings simultaneously at amide and aliphatic regions. Excellent resolution and high sensitivity is warranted by removing all coupling interactions during the indirectly detected t(1) period, and by employing pulsed field gradients for coherence selection and utilizing coherence order selective spin-state selection. The scalar and residual dipolar couplings can be readily measured from a two-dimensional (15)N/(13)C-HSQC spectrum without additional spectral crowding. SESAME-HSQC can be used for epitope mapping by observing chemical shift changes in both amide and aliphatic regions. Simultaneously, potential conversion in protein conformation can be probed by analyzing changes in residual dipolar couplings induced by ligand binding. The pulse sequence is experimentally verified with a sample of (15)N/(13)C enriched human ubiquitin. The internuclear vector directions determined from the residual dipolar couplings are found to be in excellent correlation with those predicted from ubiquitin's refined solution structure.

Geometry, energetics, and dynamics of hydrogen bonds in proteins: structural information derived from NMR scalar couplings

An accurate description of hydrogen bonds is essential to identify the determinants of protein stability and function as well as folding and misfolding behavior. We describe a method of using J couplings through hydrogen bonds as ensemble-averaged restraints in molecular dynamics simulations. Applications to the cases of ubiquitin and protein G show that these scalar couplings provide powerful structural information that, when used through the methodology that we present here, enables the description of the geometry and energetics of hydrogen bonds with an accuracy approaching that of high-resolution X-ray structures.

Computational Approach to Nuclear Magnetic Resonance in 1-Alkyl-3-methylimidazolium Ionic Liquids

A quantum-chemical computational approach to accurately predict the nuclear magnetic resonance (NMR) properties of 1-alkyl-3-methylimidazolium ionic liquids has been performed by the gauge-including atomic orbitals method at the B3LYP/6-31++G** level using different simulated ionic liquid environments. The first molecular model chosen to describe the ionic liquid system includes the gas-phase optimized structures of ion pairs and separated ions of a series of imidazolium salts containing methyl, butyl, and octyl substituents and PF6-, BF4-, and Br- anions. In addition, a continuum polarizable model of solvation has been applied to predict the effects of the medium polarity on the molecular properties of 1,3-dimethylimidazolium hexafluorophosphate (MmimPF6). Furthermore, the specific acidic and basic solute-solvent interactions have been simulated by a discrete solvation model based on molecular clusters formed by MmimPF6 species and a discrete number of water molecules. The computational prediction of the NMR spectra allows a consistent interpretation of the dispersed experimental evidence in the literature. The following are main contributions of this work: (a) Theoretical results state the presence of a chemical equilibrium between ion-pair aggregates and solvent-separated counterions of 1-alkyl-3-methylimidazolium salts which is tuned by the solvent environment; thus, strong specific (acidic and basic) and nonspecific (polarity and polarizability) solvent interactions are predicted favoring the dissociated ionic species. (b) The calculated 1H and 13C NMR properties of these ionic liquids are revealed as highly dependent on the nature of solute-solvent interactions. Thus, the chemical shift of the hydrogen atom in position two of the imidazolium ring is deviated to high values by the specific interactions with water molecules, whereas nonspecific interaction with water (as a solvent) affects, in the opposite direction, this 1H NMR parameter. (c) Last, current calculations support the presence of hydrogen bonding between counterions, suggesting the importance of this interaction in the properties of the solvent in the 1-alkyl-3-methylimidazolium ionic liquids.