On the use of pseudocontact shifts in the structure determination of metalloproteins

Contribution of Molecular Dynamics in pNMR for the Structural Determination of AnV and AnVI Complexes in Solution

In this study, we propose to use classical molecular dynamics (MD) coupled with ¹H NMR spectroscopy to study the conformations of different actinyl An^VI (An = U, Np, and Pu) and An^V (An = Np) complexes with tetra-ethyl dyglicolamide (TEDGA) ligands in order to have a better representation of such complexes in solution. Molecular dynamics simulations showed its effectiveness in interpreting the experiments by the calculation of geometric factors needed for the determination of magnetic properties of these complexes. We demonstrated that different conformations of the An^V and An^VI complexes with TEDGA exist in solution with different coordination modes, which is experimentally confirmed by ¹H NMR and EXAFS spectroscopies. Furthermore, MD simulations provide additional insights into the structures of complexes in solution since conformations with fast exchanges, which are not accessible from NMR experiments, have been observed by MD simulations.

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.

What are the methodological and theoretical prospects for paramagnetic NMR in structural biology? A glimpse into the crystal ball

NMR spectroscopy is very sensitive to the presence of unpaired electrons, which perturb the NMR chemical shifts, J splittings and nuclear relaxation rates. These paramagnetic effects have attracted increasing attention over the last decades, and their use is expected to increase further in the future because they can provide structural information not easily achievable with other techniques. In fact, paramagnetic data provide long range structural restraints that can be used to assess the accuracy of crystal structures in solution and to improve them by simultaneous refinements with the X-ray data. They are also precious for obtaining information on the conformational variability of biomolecular systems, possibly in conjunction with SAXS and/or DEER data. We foresee that new tools will be developed in the next years for the simultaneous analysis of the paramagnetic data with data obtained from different techniques, in order to take advantage synergistically of the information content of all of them. Of course, the use of the paramagnetic data for structural purposes requires the knowledge of the relationship between these data and the molecular coordinates. Recently, the equations commonly used, dating back to half a century ago, have been questioned by first principle quantum chemistry calculations. Our prediction is that further theoretical/computational improvements will essentially confirm the validity of the old semi-empirical equations for the analysis of the experimental paramagnetic data.

How Do Nuclei Couple to the Magnetic Moment of a Paramagnetic Center? A New Theory at the Gauntlet of the Experiments

The recent derivation, based on pure quantum chemistry (QC) first-principles, of the pseudocontact shifts (PCSs) caused by a paramagnetic metal center on far away nuclei has cast doubts on the validity of the semiempirical (SE) theory, predicting PCSs to arise from the metal magnetic susceptibility anisotropy. The SE theory has been used and applied countless times, especially in the last 2 decades, to obtain structural information on proteins containing paramagnetic metal ions. We show here that the QC and SE predictions can be directly tested against experiments, provided a suitable macromolecular system is used. The SE approach yields a good prediction of the experimental PCSs while the QC one does not. It appears that the classic theory is able to grasp satisfactorily the underlying physics.

Joint X-ray/NMR structure refinement of multidomain/multisubunit systems

Data integration in structural biology has become a paradigm for the characterization of biomolecular systems, and it is now accepted that combining different techniques can fill the gaps in each other's blind spots. In this frame, one of the combinations, which we have implemented in REFMAC-NMR, is residual dipolar couplings from NMR together with experimental data from X-ray diffraction. The first are exquisitely sensitive to the local details but does not give any information about overall shape, whereas the latter encodes more the information about the overall shape but at the same time tends to miss the local details even at the highest resolutions. Once crystals are obtained, it is often rather easy to obtain a complete X-ray dataset, however it is time-consuming to obtain an exhaustive NMR dataset. Here, we discuss the effect of including a-priori knowledge on the properties of the system to reduce the number of experimental data needed to obtain a more complete picture. We thus introduce a set of new features of REFMAC-NMR that allow for improved handling of RDC data for multidomain proteins and multisubunit biomolecular complexes, and encompasses the use of pseudo-contact shifts as an additional source of NMR-based information. The new feature may either help in improving the refinement, or assist in spotting differences between the crystal and the solution data. We show three different examples where NMR and X-ray data can be reconciled to a unique structural model without invoking mobility.

Towards a protocol for solution structure determination of copper(II) proteins: the case of Cu(II)Zn(II) superoxide dismutase

We have developed an optimized protocol to solve the solution structure of copper(II) proteins. After assignment, proton-proton NOEs are used for the shell where 1H spectra are conveniently observed. In a shell closer to the metal ion, 13C NMR spectra with band-selective homonuclear decoupling provide the assignment of all nuclei except for those of the metal ligands. A convenient method for the measurement of 13C longitudinal-relaxation rates (R1) of carbonyls and carboxylate moieties is proposed. 1H NOEs and 1H and 13C R1 data are sufficient to produce a good/reasonable solution structure, as demonstrated for a monomeric species of superoxide dismutase, a 153-residue protein.

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.

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.

FANTEN: a new web-based interface for the analysis of magnetic anisotropy-induced NMR data

Pseudocontact shifts (PCSs) and residual dipolar couplings (RDCs) arising from the presence of paramagnetic metal ions in proteins as well as RDCs due to partial orientation induced by external orienting media are nowadays routinely measured as a part of the NMR characterization of biologically relevant systems. PCSs and RDCs are becoming more and more popular as restraints (1) to determine and/or refine protein structures in solution, (2) to monitor the extent of conformational heterogeneity in systems composed of rigid domains which can reorient with respect to one another, and (3) to obtain structural information in protein-protein complexes. The use of both PCSs and RDCs proceeds through the determination of the anisotropy tensors which are at the origin of these NMR observables. A new user-friendly web tool, called FANTEN (Finding ANisotropy TENsors), has been developed for the determination of the anisotropy tensors related to PCSs and RDCs and has been made freely available through the WeNMR ( http://fanten-enmr.cerm.unifi.it:8080 ) gateway. The program has many new features not available in other existing programs, among which the possibility of a joint analysis of several sets of PCS and RDC data and the possibility to perform rigid body minimizations.

Simultaneous use of solution NMR and X-ray data in REFMAC5 for joint refinement/detection of structural differences

The program REFMAC5 from CCP4 was modified to allow the simultaneous use of X-ray crystallographic data and paramagnetic NMR data (pseudocontact shifts and self-orientation residual dipolar couplings) and/or diamagnetic residual dipolar couplings. Incorporation of these long-range NMR restraints in REFMAC5 can reveal differences between solid-state and solution conformations of molecules or, in their absence, can be used together with X-ray crystallographic data for structural refinement. Since NMR and X-ray data are complementary, when a single structure is consistent with both sets of data and still maintains reasonably `ideal' geometries, the reliability of the derived atomic model is expected to increase. The program was tested on five different proteins: the catalytic domain of matrix metalloproteinase 1, GB3, ubiquitin, free calmodulin and calmodulin complexed with a peptide. In some cases the joint refinement produced a single model consistent with both sets of observations, while in other cases it indicated, outside the experimental uncertainty, the presence of different protein conformations in solution and in the solid state.

Cd2+ as a Ca2+ surrogate in protein-membrane interactions: isostructural but not isofunctional

Due to its favorable spectroscopic properties, Cd(2+) is frequently used as a probe of Ca(2+) sites in proteins. We investigate the ability of Cd(2+) to act as a structural and functional surrogate of Ca(2+) in protein-membrane interactions. C2 domain from protein kinase Cα (C2α) was chosen as a paradigm for the Ca(2+)-dependent phosphatidylserine-binding peripheral membrane domains. We identified the Cd(2+)-binding sites of C2α using NMR spectroscopy, determined the 1.6 Å crystal structure of Cd(2+)-bound C2α, and characterized metal-ion-dependent interactions between C2α and phospholipid membranes using fluorescence spectroscopy and ultracentrifugation experiments. We show that Cd(2+) forms a tight complex with the membrane-binding loops of C2α but is unable to support its membrane-binding function. This is in sharp contrast with Pb(2+), which is almost as effective as Ca(2+) in driving the C2α-membrane association process. Our results provide the first direct evidence for the specific role of divalent metal ions in mediating protein-membrane interactions, have important implications for metal substitution studies in proteins, and illustrate the potential diversity of functional responses caused by toxic metal ions.

Solution structure of the PAS domain of a thermophilic YybT protein homolog reveals a potential ligand-binding site

The Bacillus subtilis protein YybT (or GdpP) and its homologs were recently established as stress signaling proteins that exert their biological effect by degrading the bacterial messenger cyclic di-AMP. YybT homologs contain a small Per-ARNT-Sim (PAS) domain (~80 amino acids) that can bind b-type heme with 1:1 stoichiometry despite the small size of the domain and the lack of a conserved heme iron-coordinating residue. We determined the solution structure of the PAS domain of GtYybT from Geobacillus thermodenitrificans by NMR spectroscopy to further probe its function. The solution structure confirms that PASGtYybT adopts the characteristic PAS fold composed of a five-stranded antiparallel β sheet and a few short α-helices. One α-helix and three central β-strands of PASGtYybT are noticeably shorter than those of the typical PAS domains. Despite the small size of the protein domain, a hydrophobic pocket is formed by the side chains of nonpolar residues stemming from the β-strands and α-helices. A set of residues in the vicinity of the pocket and in the C-terminal region at the dimeric interface exhibits perturbed NMR parameters in the presence of heme or zinc protoporphyrin. Together, the results unveil a compact PAS domain with a potential ligand-binding pocket and reinforce the view that the PASYybT domains function as regulatory domains in the modulation of cellular cyclic di-AMP concentration.

Solution structure of Pyrococcus furiosus RPP21, a component of the archaeal RNase P holoenzyme, and interactions with its RPP29 protein partner

RNase P is the ubiquitous ribonucleoprotein metalloenzyme responsible for cleaving the 5'-leader sequence of precursor tRNAs during their maturation. While the RNA subunit is catalytically active on its own at high monovalent and divalent ion concentrations, four protein subunits are associated with archaeal RNase P activity in vivo: RPP21, RPP29, RPP30, and POP5. These proteins have been shown to function in pairs: RPP21-RPP29 and POP5-RPP30. We have determined the solution structure of RPP21 from the hyperthermophilic archaeon Pyrococcus furiosus ( Pfu) using conventional and paramagnetic NMR techniques. Pfu RPP21 in solution consists of an unstructured N-terminus, two alpha-helices, a zinc binding motif, and an unstructured C-terminus. Moreover, we have used chemical shift perturbations to characterize the interaction of RPP21 with RPP29. The data show that the primary contact with RPP29 is localized to the two helices of RPP21. This information represents a fundamental step toward understanding structure-function relationships of the archaeal RNase P holoenzyme.

Numbat: an interactive software tool for fitting Deltachi-tensors to molecular coordinates using pseudocontact shifts

Pseudocontact shift (PCS) effects induced by a paramagnetic lanthanide bound to a protein have become increasingly popular in NMR spectroscopy as they yield a complementary set of orientational and long-range structural restraints. PCS are a manifestation of the chi-tensor anisotropy, the Deltachi-tensor, which in turn can be determined from the PCS. Once the Deltachi-tensor has been determined, PCS become powerful long-range restraints for the study of protein structure and protein-ligand complexes. Here we present the newly developed package Numbat (New User-friendly Method Built for Automatic Deltachi-Tensor determination). With a Graphical User Interface (GUI) that allows a high degree of interactivity, Numbat is specifically designed for the computation of the complete set of Deltachi-tensor parameters (including shape, location and orientation with respect to the protein) from a set of experimentally measured PCS and the protein structure coordinates. Use of the program for Linux and Windows operating systems is illustrated by building a model of the complex between the E. coli DNA polymerase III subunits epsilon186 and theta using PCS.

Structure determination of a Galectin-3-carbohydrate complex using paramagnetism-based NMR constraints

The determination of the location and conformation of a natural ligand bound to a protein receptor is often a first step in the rational design of molecules that can modulate receptor function. NMR observables, including NOEs, often provide the basis for these determinations. However, when ligands are carbohydrates, interactions mediated by extensive hydrogen-bonding networks often reduce or eliminate NOEs between ligand and protein protons. In these cases, it is useful to look to other distance- and orientation-dependent observables that can constrain the geometry of ligand-protein complexes. Here we illustrate the use of paramagnetism-based NMR constraints, including pseudo-contact shifts (PCS) and field-induced residual dipolar couplings (RDCs). When a paramagnetic center can be attached to the protein, field-induced RDCs and PCS reflect only bound-state properties of the ligand, even when averages over small fractions of bound states and large fractions of free states are observed. The effects can also be observed over a long range, making it possible to attach a paramagnetic center to a remote part of the protein. The system studied here is a Galectin-3-lactose complex. A lanthanide-binding peptide showing minimal flexibility with respect to the protein was integrated into the C terminus of an expression construct for the Galectin-3-carbohydrate-binding domain. Dysprosium ion, which has a large magnetic susceptibility anisotropy, was complexed to the peptide, making it possible to observe both PCSs and field-induced RDCs for the protein and the ligand. The structure determined from these constraints shows agreement with a crystal structure of a Galectin-3-N-acetyllactosamine complex.