Solution structure calculations through self-orientation in a magnetic field of a cerium(III) substituted calcium-binding protein

Pseudocontact Shifts in Biomolecular NMR Spectroscopy

Paramagnetic centers in biomolecules, such as specific metal ions that are bound to a protein, affect the nuclei in their surrounding in various ways. One of these effects is the pseudocontact shift (PCS), which leads to strong chemical shift perturbations of nuclear spins, with a remarkably long range of 50 Å and beyond. The PCS in solution NMR is an effect originating from the anisotropic part of the dipole-dipole interaction between the magnetic momentum of unpaired electrons and nuclear spins. The PCS contains spatial information that can be exploited in multiple ways to characterize structure, function, and dynamics of biomacromolecules. It can be used to refine structures, magnify effects of dynamics, help resonance assignments, allows for an intermolecular positioning system, and gives structural information in sensitivity-limited situations where all other methods fail. Here, we review applications of the PCS in biomolecular solution NMR spectroscopy, starting from early works on natural metalloproteins, following the development of non-natural tags to chelate and attach lanthanoid ions to any biomolecular target to advanced applications on large biomolecular complexes and inside living cells. We thus hope to not only highlight past applications but also shed light on the tremendous potential the PCS has in structural biology.

A protocol for the refinement of NMR structures using simultaneously pseudocontact shift restraints from multiple lanthanide ions

The binding of paramagnetic metal ions to proteins produces a number of different effects on the NMR spectra of the system. In particular, when the magnetic susceptibility of the metal ion is anisotropic, pseudocontact shifts (PCSs) arise and can be easily measured. They constitute very useful restraints for the solution structure determination of metal-binding proteins. In this context, there has been great interest in the use of lanthanide(III) ions to induce PCSs in diamagnetic proteins, e.g. through the replacement native calcium(II) ions. By preparing multiple samples in each of which a different ion of the lanthanide series is introduced, it is possible to obtain multiple independent PCS datasets that can be used synergistically to generate protein structure ensembles (typically called bundles). For typical NMR-based determination of protein structure, it is necessary to perform an energetic refinement of such initial bundles to obtain final structures whose geometric quality is suitable for deposition in the PDB. This can be conveniently done by using restrained molecular dynamics simulations (rMD) in explicit solvent. However, there are no available protocols for rMD using multiple PCS datasets as part of the restraints. In this work, we extended the PCS module of the AMBER MD package to handle multiple datasets and tuned a previously developed protocol for NMR structure refinement to achieve consistent convergence with PCS restraints. Test calculations with real experimental data show that this new implementation delivers the expected improvement of protein geometry, resulting in final structures that are of suitable quality for deposition. Furthermore, we observe that also initial structures generated only with traditional restraints can be successfully refined using traditional and PCS restraints simultaneously.

How to tackle protein structural data from solution and solid state: An integrated approach

Long-range NMR restraints, such as diamagnetic residual dipolar couplings and paramagnetic data, can be used to determine 3D structures of macromolecules. They are also used to monitor, and potentially to improve, the accuracy of a macromolecular structure in solution by validating or "correcting" a crystal model. Since crystal structures suffer from crystal packing forces they may not be accurate models for the macromolecular structures in solution. However, the presence of real differences should be tested for by simultaneous refinement of the structure using both crystal and solution NMR data. To achieve this, the program REFMAC5 from CCP4 was modified to allow the simultaneous use of X-ray crystallographic and paramagnetic NMR data and/or diamagnetic residual dipolar couplings. Inconsistencies between crystal structures and solution NMR data, if any, may be due either to structural rearrangements occurring on passing from the solution to solid state, or to a greater degree of conformational heterogeneity in solution with respect to the crystal. In the case of multidomain proteins, paramagnetic restraints can provide the correct mutual orientations and positions of domains in solution, as well as information on the conformational variability experienced by the macromolecule.

Modulating RNA Alignment Using Directional Dynamic Kinks: Application in Determining an Atomic-Resolution Ensemble for a Hairpin using NMR Residual Dipolar Couplings

Approaches that combine experimental data and computational molecular dynamics (MD) to determine atomic resolution ensembles of biomolecules require the measurement of abundant experimental data. NMR residual dipolar couplings (RDCs) carry rich dynamics information, however, difficulties in modulating overall alignment of nucleic acids have limited the ability to fully extract this information. We present a strategy for modulating RNA alignment that is based on introducing variable dynamic kinks in terminal helices. With this strategy, we measured seven sets of RDCs in a cUUCGg apical loop and used this rich data set to test the accuracy of an 0.8 μs MD simulation computed using the Amber ff10 force field as well as to determine an atomic resolution ensemble. The MD-generated ensemble quantitatively reproduces the measured RDCs, but selection of a sub-ensemble was required to satisfy the RDCs within error. The largest discrepancies between the RDC-selected and MD-generated ensembles are observed for the most flexible loop residues and backbone angles connecting the loop to the helix, with the RDC-selected ensemble resulting in more uniform dynamics. Comparison of the RDC-selected ensemble with NMR spin relaxation data suggests that the dynamics occurs on the ps-ns time scales as verified by measurements of R(1ρ) relaxation-dispersion data. The RDC-satisfying ensemble samples many conformations adopted by the hairpin in crystal structures indicating that intrinsic plasticity may play important roles in conformational adaptation. The approach presented here can be applied to test nucleic acid force fields and to characterize dynamics in diverse RNA motifs at atomic resolution.

Information content of long-range NMR data for the characterization of conformational heterogeneity

Long-range NMR data, namely residual dipolar couplings (RDCs) from external alignment and paramagnetic data, are becoming increasingly popular for the characterization of conformational heterogeneity of multidomain biomacromolecules and protein complexes. The question addressed here is how much information is contained in these averaged data. We have analyzed and compared the information content of conformationally averaged RDCs caused by steric alignment and of both RDCs and pseudocontact shifts caused by paramagnetic alignment, and found that, despite the substantial differences, they contain a similar amount of information. Furthermore, using several synthetic tests we find that both sets of data are equally good towards recovering the major state(s) in conformational distributions.

Exploring regions of conformational space occupied by two-domain proteins

The presence of heterogeneity in the interdomain arrangement of several biomolecules is required for their function. Here we present a method to obtain crucial clues to distinguish between different kinds of protein conformational distributions based on experimental NMR data. The method explores subregions of the conformational space and provides both upper and lower bounds of probability for the system to be in each subregion.

Convenient method for resolving degeneracies due to symmetry of the magnetic susceptibility tensor and its application to pseudo contact shift-based protein-protein complex structure determination

Pseudo contact shifts (PCSs) induced by paramagnetic lanthanide ions fixed in a protein frame provide long-range distance and angular information, and are valuable for the structure determination of protein-protein and protein-ligand complexes. We have been developing a lanthanide-binding peptide tag (hereafter LBT) anchored at two points via a peptide bond and a disulfide bond to the target proteins. However, the magnetic susceptibility tensor displays symmetry, which can cause multiple degenerated solutions in a structure calculation based solely on PCSs. Here we show a convenient method for resolving this degeneracy by changing the spacer length between the LBT and target protein. We applied this approach to PCS-based rigid body docking between the FKBP12-rapamycin complex and the mTOR FRB domain, and demonstrated that degeneracy could be resolved using the PCS restraints obtained from two-point anchored LBT with two different spacer lengths. The present strategy will markedly increase the usefulness of two-point anchored LBT for protein complex structure determination.

PCS-based structure determination of protein-protein complexes

A simple and fast nuclear magnetic resonance method for docking proteins using pseudo-contact shift (PCS) and (1)H(N)/(15)N chemical shift perturbation is presented. PCS is induced by a paramagnetic lanthanide ion that is attached to a target protein using a lanthanide binding peptide tag anchored at two points. PCS provides long-range (approximately 40 A) distance and angular restraints between the lanthanide ion and the observed nuclei, while the (1)H(N)/(15)N chemical shift perturbation data provide loose contact-surface information. The usefulness of this method was demonstrated through the structure determination of the p62 PB1-PB1 complex, which forms a front-to-back 20 kDa homo-oligomer. As p62 PB1 does not intrinsically bind metal ions, the lanthanide binding peptide tag was attached to one subunit of the dimer at two anchoring points. Each monomer was treated as a rigid body and was docked based on the backbone PCS and backbone chemical shift perturbation data. Unlike NOE-based structural determination, this method only requires resonance assignments of the backbone (1)H(N)/(15)N signals and the PCS data obtained from several sets of two-dimensional (15)N-heteronuclear single quantum coherence spectra, thus facilitating rapid structure determination of the protein-protein complex.

Two-point anchoring of a lanthanide-binding peptide to a target protein enhances the paramagnetic anisotropic effect

Paramagnetic lanthanide ions fixed in a protein frame induce several paramagnetic effects such as pseudo-contact shifts and residual dipolar couplings. These effects provide long-range distance and angular information for proteins and, therefore, are valuable in protein structural analysis. However, until recently this approach had been restricted to metal-binding proteins, but now it has become applicable to non-metalloproteins through the use of a lanthanide-binding tag. Here we report a lanthanide-binding peptide tag anchored via two points to the target proteins. Compared to conventional single-point attached tags, the two-point linked tag provides two to threefold stronger anisotropic effects. Though there is slight residual mobility of the lanthanide-binding tag, the present tag provides a higher anisotropic paramagnetic effect.

Double-lanthanide-binding tags: design, photophysical properties, and NMR applications

Lanthanide-binding tags (LBTs) are peptide sequences of up to 20 encoded amino acids that tightly and selectively complex lanthanide ions and can sensitize terbium (Tb3+) luminescence. On the basis of these properties, it was predicted that increasing the number of bound lanthanides would improve the capabilities of these tags. Therefore, using a structurally well-characterized single-LBT sequence as a starting point, a "double-LBT" (dLBT), which concatenates two lanthanide-binding motifs, was designed. Herein we report the generation of dLBT peptides and luminescence and NMR studies on a dLBT-tagged ubiquitin fusion protein. These lanthanide-bound constructs are shown to be improved luminescent tags with avid lanthanide binding and up to 3-fold greater luminescence intensity. NMR experiments were conducted on the ubiquitin construct, wherein bound paramagnetic lanthanides were used as alignment-inducing agents to gain residual dipolar couplings, which are valuable restraints for macromolecular structure determination. Together, these results indicate that dLBTs will be valuable chemical tools for biophysical applications leading to new approaches for studying the structure, function, and dynamics of proteins.

Magnetic susceptibility-induced alignment of proteins in reverse micelles

Proteins encapsulated within the aqueous core of reverse micelles are found to partially align in a magnetic field. The degree of alignment is sufficient to result in sizable residual 15N-1H dipolar couplings that can be easily measured. It is found that the magnetic susceptibility of the reverse micelle particle is not dominated by the encapsulated protein. The residual dipolar couplings are found to be structurally meaningful.

NMR Spectroscopy of Paramagnetic Metalloproteins

This article deals with the solution structure determination of paramagnetic metalloproteins by NMR spectroscopy. These proteins were believed not to be suitable for NMR investigations for structure determination until a decade ago, but eventually novel experiments and software protocols were developed, with the aim of making the approach suitable for the goal and as user-friendly and safe as possible. In the article, we also give hints for the optimization of experiments with respect to each particular metal ion, with the aim of also providing a handy tool for nonspecialists. Finally, a section is dedicated to the significant progress made on 13C direct detection, which reduces the negative effects of paramagnetism and may constitute a new chapter in the whole field of NMR spectroscopy.

Backbone-only restraints for fast determination of the protein fold: the role of paramagnetism-based restraints. Cytochrome b562 as an example

CH(alpha) residual dipolar couplings (Deltardc's) were measured for the oxidized cytochrome b562 from Escherichia coli as a result of its partial self-orientation in high magnetic fields due to the anisotropy of the overall magnetic susceptibility tensor. Both the low spin iron (III) heme and the four-helix bundle fold contribute to the magnetic anisotropy tensor. CH(alpha) Deltardc's, which span a larger range than the analogous NH values (already available in the literature) sample large space variations at variance with NH Deltardc's, which are largely isooriented within alpha helices. The whole structure is now significantly refined with the chemical shift index and CH(alpha) Deltardc's. The latter are particularly useful also in defining the molecular magnetic anisotropy parameters. It is shown here that the backbone folding can be conveniently and accurately determined using backbone restraints only, which include NOEs, hydrogen bonds, residual dipolar couplings, pseudocontact shifts, and chemical shift index. All these restraints are easily and quickly determined from the backbone assignment. The calculated backbone structure is comparable to that obtained by using also side chain restraint. Furthermore, the structure obtained with backbone only restraints is, in its whole, very similar to that obtained with the complete set of restraints. The paramagnetism based restraints are shown to be absolutely relevant, especially for Deltardc's.

Energetics and Mechanism of Ca2+ Displacement by Lanthanides in a Calcium Binding Protein

The displacement of Ca(2+) by trivalent lanthanide ions (Tm(3+)) in a protozoan (Entamoeba histolytica) Ca(2+) binding protein has been studied by NMR and isothermal calorimetry (ITC). The study provides a basis for understanding the behavior of lanthanides when used as a substitute for Ca(2+), the pattern of sequential binding, the structural changes involved, the range and magnitude of paramagnetic interaction, and the associated energetics and mechanism. The progressive Ca(2+) displacement from site III first, followed by displacement from site II, I, and IV, as observed during the NMR titration experiments, is interpreted in the light of ITC data to provide a deeper insight into the intradomain and, for the first time, interdomain cooperativity and information about the statistical phenomenon involved in it. A theoretical model governing Ca(2+) displacement is provided. The small structural changes involved in Ca(2+) displacement by a diamagnetic lanthanide (La(3+)) has also been monitored.

Paramagnetism-Based Refinement Strategy for the Solution Structure of Human α-Parvalbumin

In the frame of a research aimed at the detailed structural characterization of human calcium-binding proteins of the EF-hand family, the solution structure of human alpha-parvalbumin has been solved by NMR and refined with the help of substitution of the Ca(2+) ion in the EF site with the paramagnetic Dy(3+) ion. A simple (1)H-(15)N HSQC spectrum allowed the NH assignments based on the properties of Dy(3+). This allowed us to exploit pseudocontact shifts and residual dipolar couplings for solution structure refinement. The backbone and heavy atom RMSD are 0.55 +/- 0.08 and 1.02 +/- 0.08 A, respectively, and decrease to 0.39 +/- 0.05 and 0.90 +/- 0.06 A upon refinement with paramagnetism-based restraints. The RMSD for the metal itself in the EF site in the refined structure is 0.26 +/- 0.12 A. Backbone NH R(1), R(2), and NOE measured at two temperatures show the protein to be relatively rigid. The NH orientations are well determined by the paramagnetism-based restraints. This allows us to detect small but significant local structural differences with the orthologue protein from rat, whose X-ray structure is available at 2.0 A resolution. All differences are related to local changes in the amino acidic composition.

A strategy for the NMR characterization of type II copper(II) proteins: the case of the copper trafficking protein CopC from Pseudomonas Syringae

CopC from Pseudomonas syringae was found to be a protein capable of binding both Cu(I) and Cu(II) at two different sites. The solution structure of the apo protein is available, and structural information has been obtained on the Cu(I) bound form. We attempt here to set the limits for the determination of the solution structure of a Cu(II) protein, such as the Cu(II) bound form of CopC, in which the Cu(II) ion takes a type II coordination. The electron relaxation time is estimated from NMRD measurements to be 3 ns which leads to a correlation time for the nuclear spin-electron spin dipolar interaction of 2 ns. This information allowed us to tailor the NMR experiments and to fully exploit purely heteronuclear spectroscopy to assign as many signals as possible. In this way, 37 (13)C and 11 (15)N signals that completely escape detection with conventional approaches were assigned. Paramagnetic based structural constraints were obtained by measuring paramagnetic longitudinal relaxation enhancements (rho(para)) which allowed us to precisely locate the copper ion within the protein frame. Pseudocontact shifts (pcs's) were also used as constraints for 83 (1)H and 18 (13)C nuclei. With them, together with other standard structural constraints, a structure is obtained (and submitted to PDB) where information is only missing in a sphere with a 6 A radius from the copper ion. If we borrow information from EXAFS data, which show evidence of two copper coordinated histidines, then His 1 and His 91 are unambiguously identified as copper ligands. EXAFS data indicate two more light donor atoms (O/N) which could be from Asp 27 and Glu 89, whereas the NMRD data indicate the presence of a semicoordinated water molecule at 2.8 A (Cu-O distance) roughly orthogonal to the plane identified by the other four ligands. This represents the most extensively characterized structure of a type II Cu(II) protein obtained employing the most advanced NMR methods and with the aid of EXAFS data. The knowledge of the location of the Cu(II) in the protein is important for the copper transfer mechanism.

15N-1H Residual dipolar coupling analysis of native and alkaline-K79A Saccharomyces cerevisiae cytochrome c

Residual dipolar couplings (RDCs) and pseudocontact shifts are experimentally accessible properties in nuclear magnetic resonance that are related to structural parameters and to the magnetic susceptibility anisotropy. We have determined RDCs due to field-induced orientation of oxidized-K79A and reduced cytochrome c at pH 7.0 and oxidized-K79A cytochrome c at pH 11.1 through measurements of amide (15)N-(1)H (1)J couplings at 800 and 500 MHz. The pH 7.0 RDCs for Fe(III)- and Fe(II)-cytochrome c together with available nuclear Overhauser effects were used to recalculate solution structures that were consistent with both sets of constraints. Molecular magnetic susceptibility anisotropy values were calculated for both redox states of the protein. By subtracting the residual dipolar couplings (RDCs) of the reduced form from those of the oxidized form measured at the same magnetic field (800 MHz), we found the RDC contribution of the paramagnetic metal ion in the oxidized protein. The magnetic susceptibility anisotropy, which was calculated from the structure, was found to be the same as that of the paramagnetic metal ion obtained independently from pseudocontact shifts, thereby indicating that the elements of secondary structure either are rigid or display the same mobility in both oxidation states. The residual dipolar coupling values of the alkaline-K79A form are small with respect to those of oxidized native cytochrome, whereas the pseudocontact shifts are essentially of the same magnitude, indicating local mobility. Importantly, this is the first time that mobility has been found through comparison of RDCs with pseudocontact shifts.

Structural basis for sequential displacement of Ca(2+) by Yb(3+) in a protozoan EF-hand calcium binding protein

We have studied the displacement of Ca(2+)by the trivalent lanthanide ions (Yb(3+)) in a protozoan (Entamoeba histolytica) Ca(2+)-binding protein (EhCaBP), by NMR and thermodynamics. We have demonstrated, for the first time, how one can use in a combined fashion the utility of NMR and thermodynamics to have an insight to the relative binding specificities/affinity between Ca(2+) and Yb(3+). As revealed by the titration experiments, Yb(3+) displaces Ca(2+) from the four metal binding sites present in EhCaBP in a sequential manner. The study provides a structural origin for such a sequential Ca(2+) displacement by Yb(3+) in EhCaBP.

Paramagnetic resonance of biological metal centers

The review deals with recent advances in magnetic resonance spectroscopy (hf EPR and NMR) of paramagnetic metal centers in biological macromolecules. In the first half of our chapter, we present an overview of recent technical developments in the NMR of paramagnetic bio-macromolecules. These are illustrated by a variety of examples deriving mainly from the spectroscopy of metalloproteins and their complexes. The second half focuses on recent developments in high-frequency EPR spectroscopy and the application of the technique to copper, iron, and manganese proteins. Special attention is given to the work on single crystals of copper proteins.

Model-free analysis of protein backbone motion from residual dipolar couplings

On the basis of the measurement of NH residual dipolar couplings (RDCs) in 11 different alignment media, an RDC-based order parameter is derived for each residue in the protein ubiquitin. Dipolar couplings are motionally averaged in the picosecond to millisecond time range and, therefore, reflect motion slower than the inverse overall tumbling correlation time of the protein. It is found that there is considerable motion that is slower than the correlation time and could not be detected with previous NMR methodology. Amplitudes and anisotropies of the motion can be derived from the model-free analysis. The method can be applied provided that at least five sufficiently different alignment media can be found for the biomolecule under investigation.

Paramagnetically induced residual dipolar couplings for solution structure determination of lanthanide binding proteins

Lanthanides may substitute calcium in calcium-binding proteins, such as, for instance, EF-hand proteins. Paramagnetic lanthanides are capable of orienting the protein in high magnetic fields to an extent similar to that obtained by using orienting devices, and each lanthanide orients according to its magnetic susceptibility tensor. Here, Ce(3+), Tb(3+), Dy(3+), Ho(3+), Er(3+), Tm(3+), Yb(3+) in the C-terminal site of calbindin D(9k) have been investigated. Such systems provide (1)H-(15)N residual dipolar couplings (rdc) which can be used for solution structure determinations. Within the frame of optimizing the use of residual dipolar couplings for efficient solution structure determination, it is proposed here to use a number of lanthanides (e.g., >2) to obtain the orientations of the internuclear vectors with respect to an arbitrary reference system. This is facilitated by the independent knowledge of the magnetic susceptibility anisotropy tensor of each metal, obtained from the analysis of the pseudocontact shifts. A further module of the program PARAMAGNETIC-DYANA, called RDCDYANA-ANGLES, is developed to efficiently incorporate such rdc-derived orientations, instead of the rdc themselves, as constraints in the solution structure calculation. This strategy is absolutely general and can be extended to any other pair of dipole-dipole coupled nuclei. The effect of mobility is also assessed. In principle, information on the mobility can be obtained with a number of lanthanide ions >5, or by combining a smaller number of lanthanide ions with a few orienting devices.

Lanthanide modulation of the orientation of macromolecules induced by purple membrane

Addition of Tb3+ to purple membrane (PM) suspensions changes the orientation of the menbrane normal from parallel to perpendicular with respect to the magnetic field. Residual dipolar couplings measured in protein L in the presence of PM are scaled by a factor of -1/2. NMR line broadening and cross-correlation effects induced by the addition of PM are partially reversed by Tb3+ but not by Tm3+ that has no effect on the orientation of PM. This is interpreted as the result of anisotropically restricted motion of protein L in the transiently PM-bound form.

Partial alignment of biomolecules: an aid to NMR characterization

Partial alignment of biomolecules in solution has added a new dimension to structural investigation by high-resolution NMR methods. Applications to proteins, nucleic acids and carbohydrates now abound. Limitations initially associated with compatibility of biomolecules with the liquid-crystal media commonly used to achieve alignment have begun to disappear. This is, in part, a result of the introduction of a wide variety of new media. Future applications to biologically important problems such as the structural organization of multi-domain proteins and multi-protein assemblies look very promising.

Cross correlation between the dipole-dipole interaction and the Curie spin relaxation: the effect of anisotropic magnetic susceptibility

Cross-correlated relaxation caused by the interference of nuclear dipole-dipole interaction and the Curie spin relaxation (DD-CSR cross relaxation) is generalized to treat the case of anisotropic magnetic susceptibility, including the important case where the latter originates from zero-field splitting. It is shown that the phenomenon of DD-CSR cross relaxation is absolutely general and to be expected under any electronic configuration. The results of the generalization are presented for a model system, and the consequences for paramagnetic metalloproteins are illustrated with an example of cerium(III)-substituted calbindin. The effects of the magnetic anisotropy are found to be substantial.

Locating the metal ion in calcium-binding proteins by using cerium(III) as a probe

The detection and assignment of NMR spectroscopic signals of carbon atoms from carbonyl and carboxylate groups in the loop hosting the Ce(III) ion was performed for the cerium-substituted calcium-binding protein calbindin D9k. This provided a tool to characterize in solution the first coordination sphere of the metal ion. Due to the well-documented possibility of replacing calcium with metal ions of the Ln(III) series, this approach turns out to be extremely efficient for characterizing in solution the coordination of calcium ions in proteins, independently of the availability of X-ray crystal structures. The present approach completes the structural characterization of lanthanide-substituted calcium-binding proteins, for which the role of long-range constraints arising from hyperfine interaction and self-orientation has already been assessed.