NMR experiments for the sign determination of homonuclear scalar and residual dipolar couplings

Structure determination of high-energy states in a dynamic protein ensemble

Macromolecular function frequently requires that proteins change conformation into high-energy states^1-4. However, methods for solving the structures of these functionally essential, lowly populated states are lacking. Here we develop a method for high-resolution structure determination of minorly populated states by coupling NMR spectroscopy-derived pseudocontact shifts⁵ (PCSs) with Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion⁶ (PCS-CPMG). Our approach additionally defines the corresponding kinetics and thermodynamics of high-energy excursions, thereby characterizing the entire free-energy landscape. Using a large set of simulated data for adenylate kinase (Adk), calmodulin and Src kinase, we find that high-energy PCSs accurately determine high-energy structures (with a root mean squared deviation of less than 3.5 angström). Applying our methodology to Adk during catalysis, we find that the high-energy excursion involves surprisingly small openings of the AMP and ATP lids. This previously unresolved high-energy structure solves a longstanding controversy about conformational interconversions that are rate-limiting for catalysis. Primed for either substrate binding or product release, the high-energy structure of Adk suggests a two-step mechanism combining conformational selection to this state, followed by an induced-fit step into a fully closed state for catalysis of the phosphoryl-transfer reaction. Unlike other methods for resolving high-energy states, such as cryo-electron microscopy and X-ray crystallography, our solution PCS-CPMG approach excels in cases involving domain rearrangements of smaller systems (less than 60 kDa) and populations as low as 0.5%, and enables the simultaneous determination of protein structure, kinetics and thermodynamics while proteins perform their function.

Hsp40s play complementary roles in the prevention of tau amyloid formation

The microtubule-associated protein, tau, is the major subunit of neurofibrillary tangles associated with neurodegenerative conditions, such as Alzheimer's disease. In the cell, however, tau aggregation can be prevented by a class of proteins known as molecular chaperones. While numerous chaperones are known to interact with tau, though, little is known regarding the mechanisms by which these prevent tau aggregation. Here, we describe the effects of ATP-independent Hsp40 chaperones, DNAJA2 and DNAJB1, on tau amyloid-fiber formation and compare these to the small heat shock protein HSPB1. We find that the chaperones play complementary roles, with each preventing tau aggregation differently and interacting with distinct sets of tau species. Whereas HSPB1 only binds tau monomers, DNAJB1 and DNAJA2 recognize aggregation-prone conformers and even mature fibers. In addition, we find that both Hsp40s bind tau seeds and fibers via their C-terminal domain II (CTDII), with DNAJA2 being further capable of recognizing tau monomers by a second, distinct site in CTDI. These results lay out the mechanisms by which the diverse members of the Hsp40 family counteract the formation and propagation of toxic tau aggregates and highlight the fact that chaperones from different families/classes play distinct, yet complementary roles in preventing pathological protein aggregation.

Several neurological conditions, such as Alzheimer’s and Parkinson’s disease, are characterized by the build-up of protein clumps known as aggregates. In the case of Alzheimer’s disease, a key protein, called tau, aggregates to form fibers that are harmful to neuronal cells in the brain. One of the ways our cells can prevent this from occurring is through the action of proteins known as molecular chaperones, which can bind to tau proteins and prevent them from sticking together.

Tau can take on many forms. For example, a single tau protein on its own, known as a monomer, is unstructured. In patients with Alzheimer’s, these monomers join together into small clusters, known as seeds, that rapidly aggregate and accumulate into rigid, structured fibers. One chaperone, HSPB1, is known to bind to tau monomers and prevent them from being incorporated into fibers. Recently, another group of chaperones, called J-domain proteins, was also found to interact with tau. However, it was unclear how these chaperones prevent aggregation and whether they bind to tau in a similar manner as HSPB1.

To help answer this question, Irwin, Faust et al. studied the effect of two J-domain proteins, as well as the chaperone HSBP1, on tau aggregation. This revealed that, unlike HSBP1, the two J-domain proteins can bind to multiple forms of tau, including when it has already aggregated in to seeds and fibers. This suggests that these chaperones can stop the accumulation of fibers at several different stages of the aggregation process. Further experiments examining which sections of the J-domain proteins bind to tau, showed that both attach to fibers via the same region. However, the two J-domain proteins are not identical in their interaction with tau. While one of them uses a distinct region to bind to tau monomers, the other does not bind to single tau proteins at all.

These results demonstrate how different cellular chaperones can complement one another in order to inhibit harmful protein aggregation. Further studies will be needed to understand the full role of J-domain proteins in preventing tau from accumulating into fibers, as well as their potential as drug targets for developing new treatments.

Probing Long-Range Anisotropic Interactions: a General and Sign-Sensitive Strategy to Measure 1 H-1 H Residual Dipolar Couplings as a Key Advance for Organic Structure Determination

Residual dipolar couplings (RDCs) are amongst the most powerful NMR parameters for organic structure elucidation. In order to maximize their effectiveness in increasingly complex cases such as flexible compounds, a maximum of RDCs between nuclei sampling a large distribution of orientations is needed, including sign information. For this, the easily accessible one-bond ¹ H-¹³ C RDCs alone often fall short. Long-range ¹ H-¹ H RDCs are both abundant and typically sample highly complementary orientations, but accessing them in a sign-sensitive way has been severely obstructed due to the overflow of ¹ H-¹ H couplings. Here, we present a generally applicable strategy that allows the measurement of a large number of ¹ H-¹ H RDCs, including their signs, which is based on a combination of an improved PSYCHEDELIC method and a new selective constant-time β-COSY experiment. The potential of ¹ H-¹ H RDCs to better determine molecular alignment and to discriminate between enantiomers and diastereomers is demonstrated.

Overview of refinement procedures within REFMAC5: utilizing data from different sources

Refinement is a process that involves bringing into agreement the structural model, available prior knowledge and experimental data. To achieve this, the refinement procedure optimizes a posterior conditional probability distribution of model parameters, including atomic coordinates, atomic displacement parameters (B factors), scale factors, parameters of the solvent model and twin fractions in the case of twinned crystals, given observed data such as observed amplitudes or intensities of structure factors. A library of chemical restraints is typically used to ensure consistency between the model and the prior knowledge of stereochemistry. If the observation-to-parameter ratio is small, for example when diffraction data only extend to low resolution, the Bayesian framework implemented in REFMAC5 uses external restraints to inject additional information extracted from structures of homologous proteins, prior knowledge about secondary-structure formation and even data obtained using different experimental methods, for example NMR. The refinement procedure also generates the `best' weighted electron-density maps, which are useful for further model (re)building. Here, the refinement of macromolecular structures using REFMAC5 and related tools distributed as part of the CCP4 suite is discussed.

NMR resonance splitting of urea in stretched hydrogels: proton exchange and (1)H/(2)H isotopologues

Urea at ∼12 M in concentrated gelatin gel, that was stretched, gave (1)H and (2)H NMR spectral splitting patterns that varied in a predictable way with changes in the relative proportions of (1)H2O and (2)H2O in the medium. This required consideration of the combinatorics of the two amide groups in urea that have a total of four protonation/deuteration sites giving rise to 16 different isotopologues, if all the atoms were separately identifiable. The rate constant that characterized the exchange of the protons with water was estimated by back-transformation analysis of 2D-EXSY spectra. There was no (1)H NMR spectral evidence that the chiral gelatin medium had caused in-equivalence in the protons bonded to each amide nitrogen atom. The spectral splitting patterns in (1)H and (2)H NMR spectra were accounted for by intra-molecular scalar and dipolar interactions, and quadrupolar interactions with the electric field gradients of the gelatin matrix, respectively.

1H, 13C, and 15N backbone and side chain resonance assignments of thermophilic Geobacillus kaustophilus cyclophilin-A

Cyclophilins catalyze the reversible peptidyl-prolyl isomerization of their substrates and are present across all kingdoms of life from humans to bacteria. Although numerous biological roles have now been discovered for cyclophilins, their function was initially ascribed to their chaperone-like activity in protein folding where they catalyze the often rate-limiting step of proline isomerization. This chaperone-like activity may be especially important under extreme conditions where cyclophilins are often over expressed, such as in tumors for human cyclophilins (Lee Archiv Pharm Res 33(2): 181-187, 2010), but also in organisms that thrive under extreme conditions, such as theromophilic bacteria. Moreover, the reversible nature of the peptidyl-prolyl isomerization reaction catalyzed by cyclophilins has allowed these enzymes to serve as model systems for probing the role of conformational changes during catalytic turnover (Eisenmesser et al. Science 295(5559): 1520-1523, 2002; Eisenmesser et al. Nature 438(7064): 117-121, 2005). Thus, we present here the resonance assignments of a thermophilic cyclophilin from Geobacillus kaustophilus derived from deep-sea sediment (Takami et al. Extremophiles 8(5): 351-356, 2004). This thermophilic cyclophilin may now be studied at a variety of temperatures to provide insight into the comparative structure, dynamics, and catalytic mechanism of cyclophilins.

Protein structure validation and identification from unassigned residual dipolar coupling data using 2D-PDPA

More than 90% of protein structures submitted to the PDB each year are homologous to some previously characterized protein structure. The extensive resources that are required for structural characterization of proteins can be justified for the 10% of the novel structures, but not for the remaining 90%. This report presents the 2D-PDPA method, which utilizes unassigned residual dipolar coupling in order to address the economics of structure determination of routine proteins by reducing the data acquisition and processing time. 2D-PDPA has been demonstrated to successfully identify the correct structure of an array of proteins that range from 46 to 445 residues in size from a library of 619 decoy structures by using unassigned simulated RDC data. When using experimental data, 2D-PDPA successfully identified the correct NMR structures from the same library of decoy structures. In addition, the most homologous X-ray structure was also identified as the second best structural candidate. Finally, success of 2D-PDPA in identifying and evaluating the most appropriate structure from a set of computationally predicted structures in the case of a previously uncharacterized protein Pf2048.1 has been demonstrated. This protein exhibits less than 20% sequence identity to any protein with known structure and therefore presents a compelling and practical application of our proposed work.

Molecular conformations in the pentasaccharide LNF-1 derived from NMR spectroscopy and molecular dynamics simulations

The conformational dynamics of the human milk oligosaccharide lacto-N-fucopentaose (LNF-1), α-L-Fucp-(1 → 2)-β-D-Galp-(1 → 3)-β-D-GlcpNAc-(1 → 3)-β-D-Galp-(1 → 4)-D-Glcp, has been analyzed using NMR spectroscopy and molecular dynamics (MD) computer simulations. Employing the Hadamard (13)C-excitation technique and the J-HMBC experiment, (1)H,(13)C trans-glycosidic J coupling constants were obtained, and from one- and two-dimensional (1)H,(1)H T-ROESY experiments, proton-proton cross-relaxation rates were determined in isotropic D(2)O solution. In the lyotropic liquid-crystalline medium consisting of ditetradecylphosphatidylcholine, dihexylphosphatidylcholine, N-cetyl-N,N,N-trimethylammonium bromide, and D(2)O, (1)H, (1)H and one-bond (1)H, (13)C residual dipolar couplings (RDCs), as well as relative sign information on homonuclear RDCs, were determined for the pentasaccharide. Molecular dynamics simulations with explicit water were carried out from which the internal isomerization relaxation time constant, τ(N), was calculated for transitions at the ψ torsion angle of the β-(1 → 3) linkage to the lactosyl group in LNF-1. Compared to the global reorientation time, τ(M), of ∼0.6 ns determined experimentally in D(2)O solution, the time constant for the isomerization relaxation process, τ(N(scaled)), is about one-third as large. The NMR parameters derived from the isotropic solution show very good agreement with those calculated from the MD simulations. The only notable difference occurs at the reducing end, which should be more flexible than observed by the molecular simulation, a conclusion in complete agreement with previous (13)C NMR relaxation data. A hydrogen-bond analysis of the MD simulation revealed that inter-residue hydrogen bonds on the order of ∼30% were present across the glycosidic linkages to sugar ring oxygens. This finding highlights that intramolecular hydrogen bonds might be important in preserving well-defined structures in otherwise flexible molecules. An analysis including generalized order parameters obtained from nuclear spin relaxation experiments was performed and successfully shown to limit the conformational space accessible to the molecule when the number of experimental data are too scarce for a complete conformational analysis.

Determination of the structures of symmetric protein oligomers from NMR chemical shifts and residual dipolar couplings

Symmetric protein dimers, trimers, and higher-order cyclic oligomers play key roles in many biological processes. However, structural studies of oligomeric systems by solution NMR can be difficult due to slow tumbling of the system and the difficulty in identifying NOE interactions across protein interfaces. Here, we present an automated method (RosettaOligomers) for determining the solution structures of oligomeric systems using only chemical shifts, sparse NOEs, and domain orientation restraints from residual dipolar couplings (RDCs) without a need for a previously determined structure of the monomeric subunit. The method integrates previously developed Rosetta protocols for solving the structures of monomeric proteins using sparse NMR data and for predicting the structures of both nonintertwined and intertwined symmetric oligomers. We illustrated the performance of the method using a benchmark set of nine protein dimers, one trimer, and one tetramer with available experimental data and various interface topologies. The final converged structures are found to be in good agreement with both experimental data and previously published high-resolution structures. The new approach is more readily applicable to large oligomeric systems than conventional structure-determination protocols, which often require a large number of NOEs, and will likely become increasingly relevant as more high-molecular weight systems are studied by NMR.

Correlation of (2)J couplings with protein secondary structure

Geminal two-bond couplings ((2)J) in proteins were analyzed in terms of correlation with protein secondary structure. NMR coupling constants measured and evaluated for a total six proteins comprise 3999 values of (2)J(CalphaN'), (2)J(C'HN), (2)J(HNCalpha), (2)J(C'Calpha), (2)J(HalphaC'), (2)J(HalphaCalpha), (2)J(CbetaC'), (2)J(N'Halpha), (2)J(N'Cbeta), and (2)J(N'C'), encompassing an aggregate 969 amino-acid residues. A seamless chain of pattern comparisons across the spectrum datasets recorded allowed the absolute signs of all (2)J coupling constants studied to be retrieved. Grouped by their mediating nucleus, C', N' or C(alpha), (2)J couplings related to C' and N' depend significantly on phi,psi torsion-angle combinations. beta turn types I, I', II and II', especially, can be distinguished on the basis of relative-value patterns of (2)J(CalphaN'), (2)J(HNCalpha), (2)J(C'HN), and (2)J(HalphaC'). These coupling types also depend on planar or tetrahedral bond angles, whereas such dependences seem insignificant for other types. (2)J(HalphaCbeta) appears to depend on amino-acid type only, showing negligible correlation with torsion-angle geometry. Owing to its unusual properties, (2)J(CalphaN') can be considered a "one-bond" rather than two-bond interaction, the allylic analog of (1)J(N'Calpha), as it were. Of all protein J coupling types, (2)J(CalphaN') exhibits the strongest dependence on molecular conformation, and among the (2)J types, (2)J(HNCalpha) comes second in terms of significance, yet was hitherto barely attended to in protein structure work.

E46K Parkinson's-linked mutation enhances C-terminal-to-N-terminal contacts in alpha-synuclein

Parkinson's disease (PD) is associated with the deposition of fibrillar aggregates of the protein alpha-synuclein (alphaS) in neurons. Intramolecular contacts between the acidic C-terminal tail of alphaS and its N-terminal region have been proposed to regulate alphaS aggregation, and two originally described PD mutations, A30P and A53T, reportedly reduce such contacts. We find that the most recently discovered PD-linked alphaS mutation E46K, which also accelerates the aggregation of the protein, does not interfere with C-terminal-to-N-terminal contacts and instead enhances such contacts. Furthermore, we do not observe a substantial reduction in such contacts in the two previously characterized mutants. Our results suggest that C-terminal-to-N-terminal contacts in alphaS are not strongly protective against aggregation, and that the dominant mechanism by which PD-linked mutations facilitate alphaS aggregation may be altering the physicochemical properties of the protein such as net charge (E46K) and secondary structure propensity (A30P and A53T).

Evidence of reciprocal reorientation of the catalytic and hemopexin-like domains of full-length MMP-12

The proteolytic activity of matrix metalloproteinases toward extracellular matrix components (ECM), cytokines, chemokines, and membrane receptors is crucial for several homeostatic and pathological processes. Active MMPs are a family of single-chain enzymes (23 family members in the human genome), most of which constituted by a catalytic domain and by a hemopexin-like domain connected by a linker. The X-ray structures of MMP-1 and MMP-2 suggest a conserved and well-defined spatial relationship between the two domains. Here we present structural data for MMP-12, suitably stabilized against self-hydrolysis, both in solution (NMR and SAXS) and in the solid state (X-ray), showing that the hemopexin-like and the catalytic domains experience conformational freedom with respect to each other on a time scale shorter than 10(-8) s. Hints on the probable conformations are also obtained. This experimental finding opens new perspectives for the often hypothesized active role of the hemopexin-like domain in the enzymatic activity of MMPs.

Structural and functional characterization of the interaction between cyclophilin B and a heparin-derived oligosaccharide

The chemotaxis and integrin-mediated adhesion of T lymphocytes triggered by secreted cyclophilin B (CypB) depend on interactions with both cell surface heparan sulfate proteoglycans (HSPG) and the extracellular domain of the CD147 membrane receptor. Here, we use NMR spectroscopy to characterize the interaction of CypB with heparin-derived oligosaccharides. Chemical shift perturbation experiments allowed the precise definition of the heparan sulfate (HS) binding site of CypB. The N-terminal extremity of CypB, which contains a consensus sequence for heparin-binding proteins was modeled on the basis of our experimental NMR data. Because the HS binding site extends toward the CypB catalytic pocket, we measured its peptidyl-prolyl cis-trans isomerase (PPIase) activity in the absence or presence of a HS oligosaccharide toward a CD147-derived peptide. We report the first direct evidence that CypB is enzymatically active on CD147, as it is able to accelerate the cis/trans isomerization of the Asp(179)-Pro(180) bond in a CD147-derived peptide. However, HS binding has no significant influence on this PPIase activity. We thus conclude that the glycanic moiety of HSPG serves as anchor for CypB at the cell surface, and that the signal could be transduced by CypB via its PPIase activity toward CD147.

Structural studies on the RNA-recognition motif of NELF E, a cellular negative transcription elongation factor involved in the regulation of HIV transcription

The elongation of transcription of HIV RNA at the TAR (transactivation-response element) is highly regulated by positive and negative factors. The cellular negative transcription elongation factor NELF (negative elongation factor) was suggested to be involved in transcriptional regulation of HIV-1 (HIV type 1) by binding to the stem of the viral TAR RNA which is synthesized by cellular RNA polymerase II at the viral long terminal repeat. NELF is a heterotetrameric protein consisting of NELF A, B, C or the splice variant D, and E. In the present study, we determined the solution structure of the RRM (RNA-recognition motif) of the RNA-binding subunit NELF E and studied its interaction with the viral TAR RNA. Our results show that the separately expressed recombinant NELF E RRM has alpha-helical and beta-strand elements adopting a betaalphabetabetaalphabeta fold and is able to bind to TAR RNA. Fluorescence equilibrium titrations with fluorescently labelled double- and single-stranded oligoribonucleotides representing the TAR RNA stem imply that NELF E RRM binds to the single-stranded TAR RNAs with K(d) values in the low-micromolar range.

Molecular conformations of a disaccharide investigated using NMR spectroscopy

The molecular structure of alpha-L-Rhap-(1--> 2)-alpha-L-Rhap-OMe has been investigated using conformation sensitive NMR parameters: cross-relaxation rates, scalar 3J(CH) couplings and residual dipolar couplings obtained in a dilute liquid crystalline phase. The order matrices of the two sugar residues are different, which indicates that the molecule cannot exist in a single conformation. The conformational distribution function, P(phi, psi) related to the two glycosidic linkage torsion angles phi and psi was constructed using the APME method, valid in the low orientational order limit. The APME approach is based on the additive potential (AP) and maximum entropy (ME) models. The analyses of the trajectories generated in molecular dynamics and Langevin dynamics (LD) computer simulations gave support to the distribution functions constructed from the experimental NMR parameters. It is shown that at least two conformational regions are populated on the Ramachandran map and that these regions exhibit very different molecular order.

Cooperation between fixed and low pH-inducible interfaces controls lipoprotein release by the LDL receptor

Low-density lipoprotein (LDL) receptors bind lipoprotein particles at the cell surface and release them in the low pH environment of the endosome. The published structure of the receptor determined at endosomal pH reveals an interdomain interface between its beta propeller and its fourth and fifth ligand binding (LA) repeats, suggesting that the receptor adopts a closed conformation at low pH to release LDL. Here, we combine lipoprotein binding and release assays with NMR spectroscopy to examine structural features of the receptor promoting release of LDL at low pH. These studies lead to a model in which the receptor uses a pH-invariant scaffold as an anchor to restrict conformational search space, combining it with flexible linkers between ligand binding repeats to interconvert between open and closed conformations. This finely tuned balance between interdomain rigidity and flexibility is likely to represent a shared structural feature in proteins of the LDL receptor family.

The absolute sign of J coupling constants determined using the order matrix calculation

A novel methodology using the order matrix calculation to determine the absolute sign of spin-spin couplings based on the structure of organic compounds is presented. The sign of the residual dipolar coupling (RDC) depends on the sign of corresponding scalar spin-spin coupling constant and the sign of the RDC has a dramatic influence on the order matrix calculation. Therefore, the sign of the spin-spin coupling constant can be obtained by an order matrix calculation through the corresponding RDC. Six types of spin-spin coupling constants, including 2J(H,H), 1J(C,F), 2J(C,F), 3J(C,F), 2J(F,H) and 3J(F,H), were obtained simultaneously. Except for 3J(C,F) where the measured RDCs have very small magnitudes, the signs were determined unambiguously.

Measurement of long-range 1H-1H dipolar couplings in weakly aligned proteins

Measurement of 1H-1H dipolar couplings in macromolecules, weakly oriented by a dilute liquid crystalline medium, is generally limited to the largest such interactions. By removing dipolar couplings to nearest neighbors, either by decoupling, deuteration, or both, more remote interactions become accessible. The approach is demonstrated for measurement of amide-amide interactions in the proteins calmodulin and ubiquitin and permits observation of direct dipolar couplings between protons up to 7 A apart. Quantitative evaluation of 1H-1H dipolar couplings measured in ubiquitin shows excellent agreement with its solution structure.

Solution structure of a hydrophobic analogue of the winter flounder antifreeze protein

The solution structure of a synthetic mutant type I antifreeze protein (AFP I) was determined in aqueous solution at pH 7.0 using nuclear magnetic resonance (NMR) spectroscopy. The mutations comprised the replacement of the four Thr residues by Val and the introduction of two additional Lys-Glu salt bridges. The antifreeze activity of this mutant peptide, VVVV2KE, has been previously shown to be similar to that of the wild type protein, HPLC6 (defined here as TTTT). The solution structure reveals an alphahelix bent in the same direction as the more bent conformer of the published crystal structure of TTTT, while the side chain chi1 rotamers of VVVV2KE are similar to those of the straighter conformer in the crystal of TTTT. The Val side chains of VVVV2KE assume the same orientations as the Thr side chains of TTTT, confirming the conservative nature of this mutation. The combined data suggest that AFP I undergoes an equilibrium between straight and bent helices in solution, combined with independent equilibria between different side chain rotamers for some of the amino acid residues. The present study presents the first complete sequence-specific resonance assignments and the first complete solution structure determination by NMR of any AFP I protein.

Measurement of homonuclear proton couplings based on cross-peak nulling in CT-COSY

A method in which 1H-1H scalar and dipolar couplings are obtained from the cross-peak nulling condition in a series of constant-time (CT) COSY spectra, as a function of the duration of the CT period, is described. The method is best suited for measurement of 1H-1H couplings in the range 5-20 Hz. It is shown, however, that results can be sensitive to cross-correlated relaxation effects. Also, artifactual resonances, resulting from strong coupling, can be quite pronounced in CT-COSY spectra, even when /J(AB)/(deltaA-deltaB)/<0.1. The experiments are demonstrated for the DNA dodecamer d(CGCGAATTCGCG)2, both in isotropic solution and in a liquid crystalline phase.

Measurement of homonuclear proton couplings from regular 2D COSY spectra

An interactive computer procedure is described which determines (1)H--(1)H couplings from fitting the cross-peak multiplets in regular phase-sensitive COSY spectra. The robustness and simplicity of the method rely on the fact that a given cross-peak intensity is not an independent variable in the fitting procedure, making it possible to measure couplings accurately even from individual cross peaks with unresolved multiplet structure.