Paramagnetism-based restraints for Xplor-NIH

The HADDOCK2.4 web server for integrative modeling of biomolecular complexes

Interactions between macromolecules, such as proteins and nucleic acids, are essential for cellular functions. Experimental methods can fail to provide all the information required to fully model biomolecular complexes at atomic resolution, particularly for large and heterogeneous assemblies. Integrative computational approaches have, therefore, gained popularity, complementing traditional experimental methods in structural biology. Here, we introduce HADDOCK2.4, an integrative modeling platform, and its updated web interface ( https://wenmr.science.uu.nl/haddock2.4 ). The platform seamlessly integrates diverse experimental and theoretical data to generate high-quality models of macromolecular complexes. The user-friendly web server offers automated parameter settings, access to distributed computing resources, and pre- and post-processing steps that enhance the user experience. To present the web server's various interfaces and features, we demonstrate two different applications: (i) we predict the structure of an antibody-antigen complex by using NMR data for the antigen and knowledge of the hypervariable loops for the antibody, and (ii) we perform coarse-grained modeling of PRC1 with a nucleosome particle guided by mutagenesis and functional data. The described protocols require some basic familiarity with molecular modeling and the Linux command shell. This new version of our widely used HADDOCK web server allows structural biologists and non-experts to explore intricate macromolecular assemblies encompassing various molecule types.

The evolution of paramagnetic NMR as a tool in structural biology

Paramagnetic NMR data contain extremely accurate long-range information on metalloprotein structures and, when used in the frame of integrative structural biology approaches, they allow for the retrieval of structural details to a resolution that is not achievable using other techniques. Paramagnetic data thus represent an extremely powerful tool to refine protein models in solution, especially when coupled to X-ray or cryoelectron microscopy data, to monitor the formation of complexes and determine the relative arrangements of their components, and to highlight the presence of conformational heterogeneity. More recently, theoretical and computational advancements in quantum chemical calculations of paramagnetic NMR observables are progressively opening new routes in structural biology, because they allow for the determination of the structure within the coordination sphere of the metal center, thus acting as a loupe on sites that are difficult to observe but very important for protein function.

Incorporation of residual chemical shift anisotropy into the treatment of 15N pseudocontact shifts for structural refinement

Paramagnetic NMR experiments, including the pseudocontact shift experiment, have seen increasing use due to recently developed probes and labeling strategies. The pseudocontact shift experiment can provide valuable intra- or inter-molecular distance and orientation information. However, the use of ¹H/¹³C or ¹H/¹⁵N PCS data in structure calculations is currently complicated by the contribution of residual chemical shift anisotropy to the ¹³C or ¹⁵N datasets. Here, we present a corrected PCS energy term for the software package Xplor-NIH with the appropriate residual chemical shift anisotropy correction and show its suitability for model refinements of ubiquitin labeled at residue 57 with a Tm-M8-SPy tag. For data taken at 800 MHz, the improvement with the corrected energy term is sufficient to make the quality of the fit for the ¹⁵N dataset comparable to that of the ¹H dataset, for which no correction is needed. The corrected energy term is expected to become more relevant with increased use of higher field instruments and as new paramagnetic probes with larger magnetic susceptibility tensors continue to be developed.

Paramagnetic Chemical Probes for Studying Biological Macromolecules

Paramagnetic chemical probes have been used in electron paramagnetic resonance (EPR) and nuclear magnetic resonance (NMR) spectroscopy for more than four decades. Recent years witnessed a great increase in the variety of probes for the study of biological macromolecules (proteins, nucleic acids, and oligosaccharides). This Review aims to provide a comprehensive overview of the existing paramagnetic chemical probes, including chemical synthetic approaches, functional properties, and selected applications. Recent developments have seen, in particular, a rapid expansion of the range of lanthanoid probes with anisotropic magnetic susceptibilities for the generation of structural restraints based on residual dipolar couplings and pseudocontact shifts in solution and solid state NMR spectroscopy, mostly for protein studies. Also many new isotropic paramagnetic probes, suitable for NMR measurements of paramagnetic relaxation enhancements, as well as EPR spectroscopic studies (in particular double resonance techniques) have been developed and employed to investigate biological macromolecules. Notwithstanding the large number of reported probes, only few have found broad application and further development of probes for dedicated applications is foreseen.

Characterizing conformational ensembles of multi-domain proteins using anisotropic paramagnetic NMR restraints

It has been over two decades since paramagnetic NMR started to form part of the essential techniques for structural analysis of proteins under physiological conditions. Paramagnetic NMR has significantly expanded our understanding of the inherent flexibility of proteins, in particular, those that are formed by combinations of two or more domains. Here, we present a brief overview of techniques to characterize conformational ensembles of such multi-domain proteins using paramagnetic NMR restraints produced through anisotropic metals, with a focus on the basics of anisotropic paramagnetic effects, the general procedures of conformational ensemble reconstruction, and some representative reweighting approaches.

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.

Lanthanide complexes with polyaminopolycarboxylates as prospective NMR/MRI diagnostic probes: peculiarities of molecular structure, dynamics and paramagnetic properties

The paramagnetic lanthanide complexes with polyaminopolycarboxylate (PAPC) ligands attract considerable attention from the standpoint of potential applications thereof as relaxation agents used in medical magnetic resonance imaging (MRI) and in luminescent materials, as well as owing to promising use thereof as paramagnetic labels for studying the properties of biopolymers since they exhibit thermodynamic stability, good solubility in aqueous media and moderate toxicity. For the last decades, the NMR methods have been used to determine the physical and chemical properties of paramagnetic Ln compounds. The studies concerning paramagnetic NMR lanthanide-induced shifts (LISs) in dissolved Ln complexes, as well as the analysis of band shape as a function of temperature make it possible to obtain valuable information on the structure, intra- and intermolecular dynamics and paramagnetic properties thereof. This review is devoted solely to the following features: firstly, the processes of intramolecular dynamics of lanthanide complexes with polyamino-polycarboxylate ligands such as DOTA, EDTA and DTPA and their derivatives studied by NMR; secondly, the LISs of lanthanide complexes with EDTA, DOTA, DTPA and some of their derivatives depending on temperature and pH. Moreover, in this review, for the first time, the dependence of the activation energy of molecular dynamics in complexes with polydentate ligands on the atomic number of the lanthanide cation is analyzed and a monotonic change in energy is detected, which is due to the effect of lanthanide contraction. It should be noted that this phenomenon is quite general and may also appear in the future in many other series of lanthanide complexes with both other multidentate ligands and with bidentate and monodentate ligands. In the future, it is possible to predict the dependence of the properties of certain lanthanide complexes on the ionic radius of the lanthanide cation based on the approaches presented in the review. In this review, we have also presented the dynamic NMR as the main research method widely used to analyze the processes of molecular dynamics, and the structural studies based on the NMR relaxation spectroscopy and LIS analysis.

Conformational and Structural characterization of carbohydrates and their interactions studied by NMR

Carbohydrates, either free or as glycans conjugated with other biomolecules, participate in many essential biological processes. Their apparent simplicity in terms of chemical functionality hides an extraordinary diversity and structural complexity. Deeply deciphering at the atomic level their structures is essential to understand their biological function and activities, but it is still a challenging task in need of complementary approaches and no generalized procedures are available to address the study of such complex, natural glycans. The versatility of Nuclear Magnetic Resonance spectroscopy (NMR) often makes it the preferred choice to study glycans and carbohydrates in solution media. The most basic NMR parameters, namely chemical shifts, coupling constants and nuclear Overhauser effects, allow defining short or repetitive chain sequences and characterize their structures and local geometries either in the free state or when interacting with other biomolecules, rendering additional information on the molecular recognition processes. The increased accessibility to carbohydrate molecules extensively or selectively labeled with 13C boosts the resolution and detail that analyzed glycan structures can reach. In turn, structural information derived from NMR, complemented with molecular modeling and theoretical calculations can also provide dynamic information on the conformational flexibility of carbohydrate structures. Furthermore, using partially oriented media or paramagnetic perturbations, it has been possible to introduce additional long-range observables rendering structural information on longer and branched glycan chains. In this review, we provide examples of these studies and an overview of the recent and most relevant NMR applications in the glycobiology field.

A High-Resolution View of the Coordination Environment in a Paramagnetic Metalloprotein from its Magnetic Properties

Metalloproteins constitute a significant fraction of the proteome of all organisms and their characterization is critical for both basic sciences and biomedical applications. A large portion of metalloproteins bind paramagnetic metal ions, and paramagnetic NMR spectroscopy has been widely used in their structural characterization. However, the signals of nuclei in the immediate vicinity of the metal center are often broadened beyond detection. In this work, we show that it is possible to determine the coordination environment of the paramagnetic metal in the protein at a resolution inaccessible to other techniques. Taking the structure of a diamagnetic analogue as a starting point, a geometry optimization is carried out by fitting the pseudocontact shifts obtained from first principles quantum chemical calculations to the experimental ones.

Characterization of lanthanoid-binding proteins using NMR spectroscopy

The variety of magnetic properties exhibited by paramagnetic lanthanoids provides outstanding information in NMR-based structural biology and therefore can be a very useful tool for characterizing lanthanoid-binding proteins. Because of their dependence on the relative positions of the protein nuclei and of the lanthanoid ion, the paramagnetic restraints (PCS, PRDC and PRE) provide information on structure and dynamics of proteins. In this Chapter, we cover the use of lanthanoids in structural biology including protein sample preparation, NMR experiments and data interpretation.

Paramagnetic NMR Spectroscopy Is a Tool to Address Reactivity, Structure, and Protein–Protein Interactions of Metalloproteins: The Case of Iron–Sulfur Proteins

The study of cellular machineries responsible for the iron–sulfur (Fe–S) cluster biogenesis has led to the identification of a large number of proteins, whose importance for life is documented by an increasing number of diseases linked to them. The labile nature of Fe–S clusters and the transient protein–protein interactions, occurring during the various steps of the maturation process, make their structural characterization in solution particularly difficult. Paramagnetic nuclear magnetic resonance (NMR) has been used for decades to characterize chemical composition, magnetic coupling, and the electronic structure of Fe–S clusters in proteins; it represents, therefore, a powerful tool to study the protein–protein interaction networks of proteins involving into iron–sulfur cluster biogenesis. The optimization of the various NMR experiments with respect to the hyperfine interaction will be summarized here in the form of a protocol; recently developed experiments for measuring longitudinal and transverse nuclear relaxation rates in highly paramagnetic systems will be also reviewed. Finally, we will address the use of extrinsic paramagnetic centers covalently bound to diamagnetic proteins, which contributed over the last twenty years to promote the applications of paramagnetic NMR well beyond the structural biology of metalloproteins.

Paramagnetic NMR in drug discovery

The presence of an unpaired electron in paramagnetic molecules generates significant effects in NMR spectra, which can be exploited to provide restraints complementary to those used in standard structure-calculation protocols. NMR already occupies a central position in drug discovery for its use in fragment screening, structural biology and validation of ligand-target interactions. Paramagnetic restraints provide unique opportunities, for example, for more sensitive screening to identify weaker-binding fragments. A key application of paramagnetic NMR in drug discovery, however, is to provide new structural restraints in cases where crystallography proves intractable. This is particularly important at early stages in drug-discovery programs where crystal structures of weakly-binding fragments are difficult to obtain and crystallization artefacts are probable, but structural information about ligand poses is crucial to guide medicinal chemistry. Numerous applications show the value of paramagnetic restraints to filter computational docking poses and to generate interaction models. Paramagnetic relaxation enhancements (PREs) generate a distance-dependent effect, while pseudo-contact shift (PCS) restraints provide both distance and angular information. Here, we review strategies for introducing paramagnetic centers and discuss examples that illustrate the utility of paramagnetic restraints in drug discovery. Combined with standard approaches, such as chemical shift perturbation and NOE-derived distance information, paramagnetic NMR promises a valuable source of information for many challenging drug-discovery programs.

Paramagpy: software for fitting magnetic susceptibility tensors using paramagnetic effects measured in NMR spectra

Paramagnetic metal ions with fast-relaxing electrons generate pseudocontact shifts (PCSs), residual dipolar couplings (RDCs), paramagnetic relaxation enhancements (PREs) and cross-correlated relaxation (CCR) in the nuclear magnetic resonance (NMR) spectra of the molecules they bind to. These effects offer long-range structural information in molecules equipped with binding sites for such metal ions. Here we present the new open-source software Paramagpy, which has been written in Python 3 with a graphic user interface. Paramagpy combines the functionalities of different currently available programs to support the fitting of magnetic susceptibility tensors using PCS, RDC, PRE and CCR data and molecular coordinates in Protein Data Bank (PDB) format, including a convenient graphical user interface. Paramagpy uses efficient fitting algorithms to avoid local minima and supports corrections to back-calculated PCS and PRE data arising from cross-correlation effects with chemical shift tensors. The source code is available from 10.5281/zenodo.3594568 .

Assessing Structural Preferences of Unstructured Protein Regions by NMR

Biomacromolecules, such as proteins, often exhibit significant motions intimately associated with their function. Intrinsically disordered proteins and proteins with intrinsically disordered regions, although extremely important for a plethora of cellular functions, are difficult to structurally characterize at the atomic level because the experimental parameters report on ensemble and time averages. Here, we demonstrate for the C-terminal domain of the human immunodeficiency virus type 1 capsid protein that NMR and, in particular, residual dipolar couplings (RDCs) measured for the folded portion of the protein can inform on the structural preferences of the unstructured portion using RDC-prediction tools and the maximum occurrence approach.

Resonance Assignments of Lowly Populated and Unstable Enzyme Intermediate Complex under Real-Time Conditions

Unstable and low-abundance protein complexes represent a large family of transient protein complexes that are difficult to characterize, even by means of high-resolution NMR spectroscopy. A method to assign the NMR signals of these unstable complexes through a combination of selective isotope labeling of amino acids in a protein and site-specific labeling the protein with a paramagnetic tag is presented herein. By using this method, the resonances of unstable thioester intermediate complex (lifetime <5 h and highest concentration ≈20 μm) generated by Staphylococcus aureus sortase A and its peptide substrate under a real-time reaction have been assigned.

Capturing a dynamically interacting inhibitor by paramagnetic NMR spectroscopy

Transient and fuzzy intermolecular interactions are fundamental to many biological processes. Despite their importance, they are notoriously challenging to characterize. Effects induced by paramagnetic ligands in the NMR spectra of interacting biomolecules provide an opportunity to amplify subtle manifestations of weak intermolecular interactions observed for diamagnetic ligands. Here, we present an approach to characterizing dynamic interactions between a partially flexible dimeric protein, HIV-1 protease, and a metallacarborane-based ligand, a system for which data obtained by standard NMR approaches do not enable detailed structural interpretation. We show that for the case where the experimental data are significantly averaged to values close to zero the standard fitting of pseudocontact shifts cannot provide reliable structural information. We based our approach on generating a large ensemble of full atomic models, for which the experimental data can be predicted, ensemble averaged and finally compared to the experiment. We demonstrate that a combination of paramagnetic NMR experiments, quantum chemical calculations, and molecular dynamics simulations offers a route towards structural characterization of dynamic protein-ligand complexes.

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.

Determination of the magnetic properties and orientation of the heme axial ligands of PpcA from Geobacter metallireducens by paramagnetic NMR

The rising interest in the use of Geobacter bacteria for biotechnological applications demands a deep understanding of how these bacteria are able to thrive in a variety of environments and perform extracellular electron transfer. The Geobacter metallireducens bacterium can couple the oxidation of a wide range of compounds to the reduction of several extracellular acceptors, including heavy metals, toxic organic compounds or electrode surfaces. The periplasmic c-type cytochrome PpcA from this bacterium is a member of a family composed of five periplasmic triheme cytochromes, which are important to bridge the electron transfer between the cytoplasm and the extracellular environment. To better understand the functional mechanism of PpcA it is essential to obtain structural data for this cytochrome. In this work, the geometry of the heme axial ligands, as well as the magnetic properties of the hemes were determined for the oxidized form of the cytochrome, using the ¹³C NMR chemical shifts of the heme α-substituents. The results were further compared with those previously obtained for the homologous cytochrome from Geobacter sulfurreducens. The orientations of the axial histidine planes and the magnetic properties of the hemes are conserved in both proteins. Overall, the results obtained allowed the definition of the orientation of the magnetic axes of PpcA from G. metallireducens, which will be used as constraints to assist the solution structure determination of the cytochrome in the oxidized form.

Assessing protein conformational landscapes: integration of DEER data in Maximum Occurrence analysis

The properties of the conformational landscape of a biomolecule are of capital importance to understand its function. It is widely accepted that a statistical ensemble is far more representative than a single structure, especially for proteins with disordered regions. While experimental data provide the most important handle on the conformational variability that the system is experiencing, they usually report on either time or ensemble averages. Since the available conformations largely outnumber the (independent) available experimental data, the latter can be equally well reproduced by a variety of ensembles. We have proposed the Maximum Occurrence (MaxOcc) approach to provide an upper bound of the statistical weight of each conformation. This method is expected to converge towards the true statistical weights by increasing the number of independent experimental datasets. In this paper we explore the ability of DEER (Double Electron Electron Resonance) data, which report on the distance distribution between two spin labels attached to a biomolecule, to restrain the MaxOcc values and its complementarity to previously introduced experimental techniques such as NMR and Small-Angle X-ray Scattering. We here present the case of Ca2+ bound calmodulin (CaM) as a test case and show that DEER data impose a sizeable reduction of the conformational space described by high MaxOcc conformations.

Simultaneous detection of intra- and inter-molecular paramagnetic relaxation enhancements in protein complexes

Paramagnetic relaxation enhancement (PRE) measurements constitute a powerful approach for detecting both permanent and transient protein-protein interactions. Typical PRE experiments require an intrinsic or engineered paramagnetic site on one of the two interacting partners; while a second, diamagnetic binding partner is labeled with stable isotopes (¹⁵N or ¹³C). Multiple paramagnetic labeled centers or reversed labeling schemes are often necessary to obtain sufficient distance restraints to model protein-protein complexes, making this approach time consuming and expensive. Here, we show a new strategy that combines a modified pulse sequence (¹H_N-Γ₂-CCLS) with an asymmetric labeling scheme to enable the detection of both intra- and inter-molecular PREs simultaneously using only one sample preparation. We applied this strategy to the non-covalent dimer of ubiquitin. Our method confirmed the previously identified binding interface for the transient di-ubiquitin complex, and at the same time, unveiled the internal structural dynamics rearrangements of ubiquitin upon interaction. In addition to reducing the cost of sample preparation and speed up PRE measurements, by detecting the intra-molecular PRE this new strategy will make it possible to measure and calibrate inter-molecular distances more accurately for both symmetric and asymmetric protein-protein complexes.

Long-range paramagnetic NMR data can provide a closer look on metal coordination in metalloproteins

Paramagnetic NMR data can be profitably incorporated in structural refinement protocols of metalloproteins or metal-substituted proteins, mostly as distance or angle restraints. However, they could in principle provide much more information, because the magnetic susceptibility of a paramagnetic metal ion is largely determined by its coordination sphere. This information can in turn be used to evaluate changes occurring in the coordination sphere of the metal when ligands (e.g.: inhibitors) are bound to the protein. This gives an experimental handle on the molecular structure in the vicinity of the metal which falls in the so-called blind sphere. The magnetic susceptibility anisotropy tensors of cobalt(II) and nickel(II) ions bound to human carbonic anhydrase II in free and inhibited forms have been determined. The change of the magnetic susceptibility anisotropy is directly linked to the binding mode of different ligands in the active site of the enzyme. Indication about the metal coordination sphere in the presence of an inhibitor in pharmaceutically relevant proteins could be important in the design of selective drugs with a structure-based approach.

Perspectives on paramagnetic NMR from a life sciences infrastructure

The effects arising in NMR spectroscopy because of the presence of unpaired electrons, collectively referred to as "paramagnetic NMR" have attracted increasing attention over the last decades. From the standpoint of the structural and mechanistic biology, paramagnetic NMR provides long range restraints that can be used to assess the accuracy of crystal structures in solution and to improve them by simultaneous refinements through NMR and X-ray data. These restraints also provide information on structure rearrangements and conformational variability in biomolecular systems. Theoretical improvements in quantum chemistry calculations can nowadays allow for accurate calculations of the paramagnetic data from a molecular structural model, thus providing a tool to refine the metal coordination environment by matching the paramagnetic effects observed far away from the metal. Furthermore, the availability of an improved technology (higher fields and faster magic angle spinning) has promoted paramagnetic NMR applications in the fast-growing area of biomolecular solid-state NMR. Major improvements in dynamic nuclear polarization have been recently achieved, especially through the exploitation of the Overhauser effect occurring through the contact-driven relaxation mechanism: the very large enhancement of the ¹³C signal observed in a variety of liquid organic compounds at high fields is expected to open up new perspectives for applications of solution NMR.

Molecular modeling of biomolecules by paramagnetic NMR and computational hybrid methods

The 3D atomic structures of biomolecules and their complexes are key to our understanding of biomolecular function, recognition, and mechanism. However, it is often difficult to obtain structures, particularly for systems that are complex, dynamic, disordered, or exist in environments like cell membranes. In such cases sparse data from a variety of paramagnetic NMR experiments offers one possible source of structural information. These restraints can be incorporated in computer modeling algorithms that can accurately translate the sparse experimental data into full 3D atomic structures. In this review, we discuss various types of paramagnetic NMR/computational hybrid modeling techniques that can be applied to successful modeling of not only the atomic structure of proteins but also their interacting partners. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.

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.

Characterizing the magnetic susceptibility tensor of lanthanide-containing polymethylated-DOTA complexes

Lanthanide complexes based on the DOTA (1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid) cage are commonly used as phase contrast agents in magnetic resonance imaging, but can also be utilized in structural NMR applications due to their ability to induce either paramagnetic relaxation enhancement or a pseudocontact shift (PCS) depending on the choice of the lanthanide. The size and sign of the PCS for any given atom is determined by its coordinates relative to the metal center, and the characteristics of the lanthanide's magnetic susceptibility tensor. Using a polymethylated DOTA tag (Ln-M8-SPy) conjugated to ubiquitin, we calculated the position of the metal center and characterized the susceptibility tensor for a number of lanthanides (dysprosium, thulium, and ytterbium) under a range of pH and temperature conditions. We found that there was a difference in temperature sensitivity for each of the complexes studied, which depended on the size of the lanthanide ion as well as the isomeric state of the cage. Using ¹⁷O-NMR, we confirmed that the temperature sensitivity of the compounds was enhanced by the presence of an apically bound water molecule. Since amide-containing lanthanide complexes are known to be pH sensitive and can be used as probes of physiological pH, we also investigated the effect of pH on the Ln-M8-SPy susceptibility tensor, but we found that the changes in this pH range (5.0-7.4) were not significant.

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.

Structure of the transmembrane domain of human nicastrin-a component of γ-secretase

Nicastrin is the largest component of γ-secretase that is an intramembrane protease important in the development of Alzheimer's disease. Nicastrin contains a large extracellular domain, a single transmembrane (TM) domain, and a short C-terminus. Its TM domain is important for the γ-secretase complex formation. Here we report nuclear magnetic resonance (NMR) studies of the TM and C-terminal regions of human nicastrin in both sodium dodecyl sulfate (SDS) and dodecylphosphocholine (DPC) micelles. Structural study and dynamic analysis reveal that the TM domain is largely helical and stable under both SDS and DPC micelles with its N-terminal region undergoing intermediate time scale motion. The TM helix contains a hydrophilic patch that is important for TM-TM interactions. The short C-terminus is not structured in solution and a region formed by residues V697-A702 interacts with the membrane, suggesting that these residues may play a role in the γ-secretase complex formation. Our study provides structural insight into the function of the nicastrin TM domain and the C-terminus in γ-secretase complex.

Defining the Apoptotic Trigger: THE INTERACTION OF CYTOCHROME c AND CARDIOLIPIN

The interaction between cytochrome c and the anionic lipid cardiolipin has been proposed as a primary event in the apoptotic signaling cascade. Numerous studies that have examined the interaction of cytochrome c with cardiolipin embedded in a variety of model phospholipid membranes have suggested that partial unfolding of the protein is a precursor to the apoptotic response. However, these studies lacked site resolution and used model systems with negligible or a positive membrane curvature, which is distinct from the large negative curvature of the invaginations of the inner mitochondrial membrane where cytochrome c resides. We have used reverse micelle encapsulation to mimic the potential effects of confinement on the interaction of cytochrome c with cardiolipin. Encapsulation of oxidized horse cytochrome c in 1-decanoyl-rac-glycerol/lauryldimethylamine-N-oxide/hexanol reverse micelles prepared in pentane yields NMR spectra essentially identical to the protein in free aqueous solution. The structure of encapsulated ferricytochrome c was determined to high precision (bb ∼ 0.23 Å) using NMR-based methods and is closely similar to the cryogenic crystal structure (bb ∼ 1.2 Å). Incorporation of cardiolipin into the reverse micelle surfactant shell causes localized chemical shift perturbations of the encapsulated protein, providing the first view of the cardiolipin/cytochrome c interaction interface at atomic resolution. Three distinct sites of interaction are detected: the so-called A- and L-sites, plus a previously undocumented interaction centered on residues Phe-36, Gly-37, Thr-58, Trp-59, and Lys-60. Importantly, in distinct contrast to earlier studies of this interaction, the protein is not significantly disturbed by the binding of cardiolipin in the context of the reverse micelle.

Nitrogen detected TROSY at high field yields high resolution and sensitivity for protein NMR

Detection of (15)N in multidimensional NMR experiments of proteins has sparsely been utilized because of the low gyromagnetic ratio (γ) of nitrogen and the presumed low sensitivity of such experiments. Here we show that selecting the TROSY components of proton-attached (15)N nuclei (TROSY (15)NH) yields high quality spectra in high field magnets (>600 MHz) by taking advantage of the slow (15)N transverse relaxation and compensating for the inherently low (15)N sensitivity. The (15)N TROSY transverse relaxation rates increase modestly with molecular weight but the TROSY gain in peak heights depends strongly on the magnetic field strength. Theoretical simulations predict that the narrowest line width for the TROSY (15)NH component can be obtained at 900 MHz, but sensitivity reaches its maximum around 1.2 GHz. Based on these considerations, a (15)N-detected 2D (1)H-(15)N TROSY-HSQC ((15)N-detected TROSY-HSQC) experiment was developed and high-quality 2D spectra were recorded at 800 MHz in 2 h for 1 mM maltose-binding protein at 278 K (τc ~ 40 ns). Unlike for (1)H detected TROSY, deuteration is not mandatory to benefit (15)N detected TROSY due to reduced dipolar broadening, which facilitates studies of proteins that cannot be deuterated, especially in cases where production requires eukaryotic expression systems. The option of recording (15)N TROSY of proteins expressed in H2O media also alleviates the problem of incomplete amide proton back exchange, which often hampers the detection of amide groups in the core of large molecular weight proteins that are expressed in D2O culture media and cannot be refolded for amide back exchange. These results illustrate the potential of (15)NH-detected TROSY experiments as a means to exploit the high resolution offered by high field magnets near and above 1 GHz.

Ligand-driven conformational changes of MurD visualized by paramagnetic NMR

Proteins, especially multi-domain proteins, often undergo drastic conformational changes upon binding to ligands or by post-translational modifications, which is a key step to regulate their function. However, the detailed mechanisms of such dynamic regulation of the functional processes are poorly understood because of the lack of an efficient tool. We here demonstrate detailed characterization of conformational changes of MurD, a 47 kDa protein enzyme consisting of three domains, by the use of solution NMR equipped with paramagnetic lanthanide probe. Quantitative analysis of pseudocontact shifts has identified a novel conformational state of MurD, named semi-closed conformation, which is found to be the key to understand how MurD regulates the binding of the ligands. The modulation of the affinity coupled with conformational changes accentuates the importance of conformational state to be evaluated in drug design.

Solution structure of the transmembrane domain of the mouse erythropoietin receptor in detergent micelles

Erythropoiesis is regulated by the erythropoietin receptor (EpoR) binding to its ligand. The transmembrane domain (TMD) and the juxtamembrane (JM) regions of the EpoR are important for signal transduction across the cell membrane. We report a solution NMR study of the mouse erythropoietin receptor (mEpoR) comprising the TMD and the JM regions reconstituted in dodecylphosphocholine (DPC) micelles. The TMD and the C-terminal JM region of the mEpoR are mainly α-helical, adopting a similar structure to those of the human EpoR. Residues from S216 to T219 in mEpoR form a short helix. Relaxation study demonstrates that the TMD of the mEpoR is rigid whilst the N-terminal region preceding the TMD is flexible. Fluorescence spectroscopy and sequence analysis indicate that the C-terminal JM region is exposed to the solvent. Helix wheel result shows that there is hydrophilic patch in the TMD of the mEpoR formed by residues S231, S238 and T242, and these residues might be important for the receptor dimerization.

Structural insight into the transmembrane domain and the juxtamembrane region of the erythropoietin receptor in micelles

Erythropoietin receptor (EpoR) dimerization is an important step in erythrocyte formation. Its transmembrane domain (TMD) and juxtamembrane (JM) region are essential for signal transduction across the membrane. A construct compassing residues S212-P259 and containing the TMD and JM region of the human EpoR was purified and reconstituted in detergent micelles. The solution structure of the construct was determined in dodecylphosphocholine (DPC) micelles by solution NMR spectroscopy. Structural and dynamic studies demonstrated that the TMD and JM region are an ?-helix in DPC micelles, whereas residues S212-D224 at the N-terminus of the construct are not structured. The JM region is a helix that contains a hydrophobic patch formed by conserved hydrophobic residues (L253, I257, and W258). Nuclear Overhauser effect analysis, fluorescence spectroscopy, and paramagnetic relaxation enhancement experiments suggested that the JM region is exposed to the solvent. The structures of the TMD and JM region of the mouse EpoR were similar to those of the human EpoR.

Integral membrane protein structure determination using pseudocontact shifts

Obtaining enough experimental restraints can be a limiting factor in the NMR structure determination of larger proteins. This is particularly the case for large assemblies such as membrane proteins that have been solubilized in a membrane-mimicking environment. Whilst in such cases extensive deuteration strategies are regularly utilised with the aim to improve the spectral quality, these schemes often limit the number of NOEs obtainable, making complementary strategies highly beneficial for successful structure elucidation. Recently, lanthanide-induced pseudocontact shifts (PCSs) have been established as a structural tool for globular proteins. Here, we demonstrate that a PCS-based approach can be successfully applied for the structure determination of integral membrane proteins. Using the 7TM α-helical microbial receptor pSRII, we show that PCS-derived restraints from lanthanide binding tags attached to four different positions of the protein facilitate the backbone structure determination when combined with a limited set of NOEs. In contrast, the same set of NOEs fails to determine the correct 3D fold. The latter situation is frequently encountered in polytopical α-helical membrane proteins and a PCS approach is thus suitable even for this particularly challenging class of membrane proteins. The ease of measuring PCSs makes this an attractive route for structure determination of large membrane proteins in general.

NMR-based modeling and refinement of protein 3D structures

NMR is a well-established method to characterize the structure and dynamics of biomolecules in solution. High-quality structures can now be produced thanks to both experimental advances and computational developments that incorporate new NMR parameters and improved protocols and force fields in the structure calculation and refinement process. In this chapter, we give a short overview of the various types of NMR data that can provide structural information, and then focus on the structure calculation methodology itself. We discuss and illustrate with tutorial examples "classical" structure calculation, refinement, and structure validation approaches.

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.

G-triplex structure and formation propensity

The occurrence of a G-triplex folding intermediate of thrombin binding aptamer (TBA) has been recently predicted by metadynamics calculations, and experimentally supported by Nuclear Magnetic Resonance (NMR), Circular Dichroism (CD) and Differential Scanning Calorimetry (DSC) data collected on a 3' end TBA-truncated 11-mer oligonucleotide (11-mer-3'-t-TBA). Here we present the solution structure of 11-mer-3'-t-TBA in the presence of potassium ions. This structure is the first experimental example of a G-triplex folding, where a network of Hoogsteen-like hydrogen bonds stabilizes six guanines to form two G:G:G triad planes. The G-triplex folding of 11-mer-3'-t-TBA is stabilized by the potassium ion and destabilized by increasing the temperature. The superimposition of the experimental structure with that predicted by metadynamics shows a great similarity, with only significant differences involving two loops. These new structural data show that 11-mer-3'-t-TBA assumes a G-triplex DNA conformation as its stable form, reinforcing the idea that G-triplex folding intermediates may occur in vivo in human guanine-rich sequences. NMR and CD screening of eight different constructs obtained by removing from one to four bases at either the 3' and the 5' ends show that only the 11-mer-3'-t-TBA yields a relatively stable G-triplex.

Structural insight into the transmembrane segments 3 and 4 of the hERG potassium channel

The hERG (human ether-a-go-go related gene) potassium channel is a voltage-gated potassium channel containing an N-terminal domain, a voltage-sensor domain, a pore domain and a C-terminal domain. The transmembrane segment 4 (S4) is important for sensing changes of membrane potentials through positively charge residues. A construct containing partial S2-S3 linker, S3, S4 and the S4-S5 linker of the hERG channel was purified into detergent micelles. This construct exhibits good quality NMR spectrum when it was purified in lyso-myristoyl phosphatidylglycerol (LMPG) micelles. Structural study showed that S3 contains two short helices with a negatively charged surface. The S4 and S4-S5 linker adopt helical structures. The six positively charged residues in S4 localize at different sides, suggesting that they may have different functions in channel gating. Relaxation studies indicated that S3 is more flexible than S4. The boundaries of S3-S4 and S4-S4-S5 linker were identified. Our results provided structural information of the S3 and S4, which will be helpful to understand their roles in channel gating.

NMR approaches for structural analysis of multidomain proteins and complexes in solution

NMR spectroscopy is a key method for studying the structure and dynamics of (large) multidomain proteins and complexes in solution. It plays a unique role in integrated structural biology approaches as especially information about conformational dynamics can be readily obtained at residue resolution. Here, we review NMR techniques for such studies focusing on state-of-the-art tools and practical aspects. An efficient approach for determining the quaternary structure of multidomain complexes starts from the structures of individual domains or subunits. The arrangement of the domains/subunits within the complex is then defined based on NMR measurements that provide information about the domain interfaces combined with (long-range) distance and orientational restraints. Aspects discussed include sample preparation, specific isotope labeling and spin labeling; determination of binding interfaces and domain/subunit arrangements from chemical shift perturbations (CSP), nuclear Overhauser effects (NOEs), isotope editing/filtering, cross-saturation, and differential line broadening; and based on paramagnetic relaxation enhancements (PRE) using covalent and soluble spin labels. Finally, the utility of complementary methods such as small-angle X-ray or neutron scattering (SAXS, SANS), electron paramagnetic resonance (EPR) or fluorescence spectroscopy techniques is discussed. The applications of NMR techniques are illustrated with studies of challenging (high molecular weight) protein complexes.