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

Machine Learning-Enhanced Quantum Chemistry-Assisted Refinement of the Active Site Structure of Metalloproteins

Understanding the fine structural details of inhibitor binding at the active site of metalloenzymes can have a profound impact on the rational drug design targeted to this broad class of biomolecules. Structural techniques such as NMR, cryo-EM, and X-ray crystallography can provide bond lengths and angles, but the uncertainties in these measurements can be as large as the range of values that have been observed for these quantities in all the published structures. This uncertainty is far too large to allow for reliable calculations at the quantum chemical (QC) levels for developing precise structure-activity relationships or for improving the energetic considerations in protein-inhibitor studies. Therefore, the need arises to rely upon computational methods to refine the active site structures well beyond the resolution obtained with routine application of structural methods. In a recent paper, we have shown that it is possible to refine the active site of cobalt(II)-substituted MMP12, a metalloprotein that is a relevant drug target, by matching to the experimental pseudocontact shifts (PCS) those calculated using multireference ab initio QC methods. The computational cost of this methodology becomes a significant bottleneck when the starting structure is not sufficiently close to the final one, which is often the case with biomolecular structures. To tackle this problem, we have developed an approach based on a neural network (NN) and a support vector regression (SVR) and applied it to the refinement of the active site structure of oxalate-inhibited human carbonic anhydrase 2 (hCAII), another prototypical metalloprotein target. The refined structure gives a remarkably good agreement between the QC-calculated and the experimental PCS. This study not only contributes to the knowledge of CAII but also demonstrates the utility of combining machine learning (ML) algorithms with QC calculations, offering a promising avenue for investigating other drug targets and complex biological systems in general.

Ligand 1 H NMR Chemical Shifts as Accurate Reporters for Protein-Ligand Binding Interfaces in Solution

The availability of high-resolution 3D structural information is crucial for investigating guest-host systems across a wide range of fields. In the context of drug discovery, the information is routinely used to establish and validate structure-activity relationships, grow initial hits from screening campaigns, and to guide molecular docking. For the generation of protein-ligand complex structural information, X-ray crystallography is the experimental method of choice, however, with limited information on protein flexibility. An experimentally verified structural model of the binding interface in the native solution-state would support medicinal chemists in their molecular design decisions. Here we demonstrate that protein-bound ligand 1 H NMR chemical shifts are highly sensitive and accurate probes for the immediate chemical environment of protein-ligand interfaces. By comparing the experimental ligand 1 H chemical shift values with those computed from the X-ray structure using quantum mechanics methodology, we identify significant disagreements for parts of the ligand between the two experimental techniques. We show that quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) ensembles can be used to refine initial X-ray co-crystal structures resulting in a better agreement with experimental 1 H ligand chemical shift values. Overall, our findings highlight the usefulness of ligand 1 H NMR chemical shift information in combination with a QM/MM MD workflow for generating protein-ligand ensembles that accurately reproduce solution structural data.

Structural Evolution of Paramagnetic Lanthanide Compounds in Solution Compared to Time- and Ensemble-Average Structures

Anisotropy in the magnetic susceptibility strongly influences the paramagnetic shifts seen in nuclear magnetic resonance (NMR) and magnetic resonance imaging (MRI) experiments. A previous study on a series of C3-symmetric prototype MRI contrast agents showed that their magnetic anisotropy was highly sensitive to changes in molecular geometry and concluded that changes in the average angle between the lanthanide-oxygen (Ln-O) bonds and the molecular C3 axis due to solvent interactions had a significant impact on the magnetic anisotropy and, consequently, the paramagnetic shift. However, this study, like many others, was predicated on an idealized C3-symmetric structural model, which may not be representative of the dynamic structure in solution at the single-molecule level. Here, we address this by using ab initio molecular dynamics simulations to simulate how the molecular geometry, in particular the angles between the Ln-O bonds and the pseudo-C3 axis, evolves over time in the solution, mimicking typical experimental conditions. We observe large-amplitude oscillations in the O-Ln-C̃3 angles, and complete active space self-consistent field spin-orbit calculations show that this leads to similarly large oscillations in the pseudocontact (dipolar) paramagnetic NMR shifts. The time-averaged shifts show good agreement with experimental measurements, while the large fluctuations suggest that an idealized structure provides an incomplete description of the solution dynamics. Our observations have significant implications for modeling the electronic and nuclear relaxation times in this and other systems where the magnetic susceptibility is exquisitely sensitive to the molecular structure.

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.

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.

Modeling a unit cell: crystallographic refinement procedure using the biomolecular MD simulation platform Amber

A procedure has been developed for the refinement of crystallographic protein structures based on the biomolecular simulation program Amber. The procedure constructs a model representing a crystal unit cell, which generally contains multiple protein molecules and is fully hydrated with TIP3P water. Periodic boundary conditions are applied to the cell in order to emulate the crystal lattice. The refinement is conducted in the form of a specially designed short molecular-dynamics run controlled by the Amber ff14SB force field and the maximum-likelihood potential that encodes the structure-factor-based restraints. The new Amber-based refinement procedure has been tested on a set of 84 protein structures. In most cases, the new procedure led to appreciably lower R free values compared with those reported in the original PDB depositions or obtained by means of the industry-standard phenix.refine program. In particular, the new method has the edge in refining low-accuracy scrambled models. It has also been successful in refining a number of molecular-replacement models, including one with an r.m.s.d. of 2.15 Å. In addition, Amber-refined structures consistently show superior MolProbity scores. The new approach offers a highly realistic representation of protein-protein interactions in the crystal, as well as of protein-water interactions. It also offers a realistic representation of protein crystal dynamics (akin to ensemble-refinement schemes). Importantly, the method fully utilizes the information from the available diffraction data, while relying on state-of-the-art molecular-dynamics modeling to assist with those elements of the structure that do not diffract well (for example mobile loops or side chains). Finally, it should be noted that the protocol employs no tunable parameters, and the calculations can be conducted in a matter of several hours on desktop computers equipped with graphical processing units or using a designated web service.

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.

Increased usability, algorithmic improvements and incorporation of data mining for structure calculation of proteins with REDCRAFT software package

Background

Traditional approaches to elucidation of protein structures by Nuclear Magnetic Resonance spectroscopy (NMR) rely on distance restraints also known as Nuclear Overhauser effects (NOEs). The use of NOEs as the primary source of structure determination by NMR spectroscopy is time consuming and expensive. Residual Dipolar Couplings (RDCs) have become an alternate approach for structure calculation by NMR spectroscopy. In previous works, the software package REDCRAFT has been presented as a means of harnessing the information containing in RDCs for structure calculation of proteins. However, to meet its full potential, several improvements to REDCRAFT must be made.

Results

In this work, we present improvements to REDCRAFT that include increased usability, better interoperability, and a more robust core algorithm. We have demonstrated the impact of the improved core algorithm in the successful folding of the protein 1A1Z with as high as ±4 Hz of added error. The REDCRAFT computed structure from the highly corrupted data exhibited less than 1.0 Å with respect to the X-ray structure. We have also demonstrated the interoperability of REDCRAFT in a few instances including with PDBMine to reduce the amount of required data in successful folding of proteins to unprecedented levels. Here we have demonstrated the successful folding of the protein 1D3Z (to within 2.4 Å of the X-ray structure) using only N-H RDCs from one alignment medium.

Conclusions

The additional GUI features of REDCRAFT combined with the NEF compliance have significantly increased the flexibility and usability of this software package. The improvements of the core algorithm have substantially improved the robustness of REDCRAFT in utilizing less experimental data both in quality and quantity.

The Photocatalyzed Thiol-ene reaction: A New Tag to Yield Fast, Selective and reversible Paramagnetic Tagging of Proteins

Paramagnetic restraints have been used in biomolecular NMR for the last three decades to elucidate and refine biomolecular structures, but also to characterize protein-ligand interactions. A common technique to generate such restraints in proteins, which do not naturally contain a (paramagnetic) metal, consists in the attachment to the protein of a lanthanide-binding-tag (LBT). In order to design such LBTs, it is important to consider the efficiency and stability of the conjugation, the geometry of the complex (conformational exchanges and coordination) and the chemical inertness of the ligand. Here we describe a photo-catalyzed thiol-ene reaction for the cysteine-selective paramagnetic tagging of proteins. As a model, we designed an LBT with a vinyl-pyridine moiety which was used to attach our tag to the protein GB1 in fast and irreversible fashion. Our tag T1 yields magnetic susceptibility tensors of significant size with different lanthanides and has been characterized using NMR and relaxometry measurements.

On the complementarity of X-ray and NMR data

X-ray crystallography and NMR contain complementary information for the structural characterization of biological macromolecules. X-ray diffraction is primarily sensitive to the overall shape of the molecule, whereas NMR is mostly sensitive to the atomic detail. Their combination can therefore provide a stronger justification for the resulting structure. For their combination we have recently proposed REFMAC-NMR, which relies on primary data from both techniques for joint refinement. This possibility raises the compelling question of how far the complementarity can be extended. In this paper, we describe an integrative approach to the refinement with NMR data of four X-ray structures of hen-egg-white lysozyme, solved at atomic resolution in four different crystal forms, and we demonstrate that the outcome critically depends on the crystal form itself, reflecting the sensitivity of NMR to fine details.

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.

C-terminal Redox Domain of Arabidopsis APR1 is a Non-Canonical Thioredoxin Domain with Glutaredoxin Function

Sulfur is an essential nutrient that can be converted into utilizable metabolic forms to produce sulfur-containing metabolites in plant. Adenosine 5'-phosphosulfate (APS) reductase (APR) plays a vital role in catalyzing the reduction of activated sulfate to sulfite, which requires glutathione. Previous studies have shown that the C-terminal domain of APR acts as a glutathione-dependent reductase. The crystal structure of the C-terminal redox domain of Arabidopsis APR1 (AtAPR1) shows a conserved α/β thioredoxin fold, but not a glutaredoxin fold. Further biochemical studies of the redox domain from AtAPR1 provided evidence to support the structural observation. Collectively, our results provide structural and biochemical information to explain how the thioredoxin fold exerts the glutaredoxin function in APR.

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.

Integrative Approaches in Structural Biology: A More Complete Picture from the Combination of Individual Techniques

With the recent technological and computational advancements, structural biology has begun to tackle more and more difficult questions, including complex biochemical pathways and transient interactions among macromolecules. This has demonstrated that, to approach the complexity of biology, one single technique is largely insufficient and unable to yield thorough answers, whereas integrated approaches have been more and more adopted with successful results. Traditional structural techniques (X-ray crystallography and Nuclear Magnetic Resonance (NMR)) and the emerging ones (cryo-electron microscopy (cryo-EM), Small Angle X-ray Scattering (SAXS)), together with molecular modeling, have pros and cons which very nicely complement one another. In this review, three examples of synergistic approaches chosen from our previous research will be revisited. The first shows how the joint use of both solution and solid-state NMR (SSNMR), X-ray crystallography, and cryo-EM is crucial to elucidate the structure of polyethylene glycol (PEG)ylated asparaginase, which would not be obtainable through any of the techniques taken alone. The second deals with the integrated use of NMR, X-ray crystallography, and SAXS in order to elucidate the catalytic mechanism of an enzyme that is based on the flexibility of the enzyme itself. The third one shows how it is possible to put together experimental data from X-ray crystallography and NMR restraints in order to refine a protein model in order to obtain a structure which simultaneously satisfies both experimental datasets and is therefore closer to the ‘real structure’.

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.

Pseudocontact shifts and paramagnetic susceptibility in semiempirical and quantum chemistry theories

Pseudocontact shifts are traditionally described as a function of the anisotropy of the paramagnetic susceptibility tensor, according to the semiempirical theory mainly developed by Kurland and McGarvey [J. Magn. Reson. 2, 286-301 (1970)]. The paramagnetic susceptibility tensor is required to be symmetric. Applying point-dipole approximation to the quantum chemistry theory of hyperfine shift, pseudocontact shifts are found to scale with a non-symmetric tensor that differs by a factor g^T/g_e from the paramagnetic susceptibility tensor derived within the semiempirical framework. We analyze the foundations of the Kurland-McGarvey pseudocontact shift expression and recall that it is inherently based on the Russell-Saunders (LS) coupling approximation for the spin-orbit coupling. We show that the difference between the semiempirical and quantum chemistry pseudocontact shift expressions arises directly from the different treatment of the orbital contribution to the hyperfine coupling.

West-Life: A Virtual Research Environment for structural biology

•

Data processing and data management services for structural biology.

•

Enhancements to existing web services for structure solution and analysis.

•

New pipelines to link these services into more complex higher-level workflows.

•

New data management facilities.

•

Making the benefits of European e-Infrastructures more accessible to structural biologists.

The West-Life project (https://about.west-life.eu/) is a Horizon 2020 project funded by the European Commission to provide data processing and data management services for the international community of structural biologists, and in particular to support integrative experimental approaches within the field of structural biology. It has developed enhancements to existing web services for structure solution and analysis, created new pipelines to link these services into more complex higher-level workflows, and added new data management facilities. Through this work it has striven to make the benefits of European e-Infrastructures more accessible to life-science researchers in general and structural biologists in particular.

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

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

Paramagnetic NMR as a new tool in structural biology

NMR (nuclear magnetic resonance) investigation through the exploitation of paramagnetic effects is passing from an approach limited to few specialists in the field to a generally applicable method that must be considered, especially for the characterization of systems hardly affordable with other techniques. This is mostly due to the fact that paramagnetic data are long range in nature, thus providing information for the structural and dynamic characterization of complex biomolecular architectures in their native environment. On the other hand, this information usually needs to be complemented by data from other sources. Integration of paramagnetic NMR with other techniques, and the development of protocols for a joint analysis of all available data, is fundamental for achieving a comprehensive characterization of complex biological systems. We describe here a few examples of the new possibilities offered by paramagnetic data used in integrated structural approaches.

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.

Detection of Correlated Protein Backbone and Side-Chain Angle Fluctuations

NMR methods for the characterization of local protein motions have attained a high level of sophistication. Measurement of the synchronization between those motions, however, poses a serious challenge. Such correlated motions are one of the underlying mechanisms for the propagation of local changes to remote sites and as such for information transfer. Here, we demonstrate the experimental detection of the synchronization of motion over an intermediate range. To that purpose, we designed pulse sequences for the measurement of cross-correlated relaxation between the backbone H^N -N and side-chain H^β -C^β dipoles in Ile, Thr, and Val in the protein GB3. These bonds are related through two and three intervening dihedral angles. We show that the correlated motions inherent in a structural ensemble obtained from a large and diverse array of NMR probes are in excellent agreement with our measurements.

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.

Structures of PGAM5 Provide Insight into Active Site Plasticity and Multimeric Assembly

PGAM5 is a mitochondrial membrane protein that functions as an atypical Ser/Thr phosphatase and is a regulator of oxidative stress response, necroptosis, and autophagy. Here we present several crystal structures of PGAM5 including the activating N-terminal regulatory sequences, providing a model for structural plasticity, dimerization of the catalytic domain, and the assembly into an enzymatically active dodecameric form. Oligomeric states observed in structures were supported by hydrogen exchange mass spectrometry, size-exclusion chromatography, and analytical ultracentrifugation experiments in solution. We report that the catalytically important N-terminal WDPNWD motif acts as a structural integrator assembling PGAM5 into a dodecamer, allosterically activating the phosphatase by promoting an ordering of the catalytic loop. Additionally the observed active site plasticity enabled visualization of essential conformational rearrangements of catalytic elements. The comprehensive biophysical characterization offers detailed structural models of this key mitochondrial phosphatase that has been associated with the development of diverse diseases.

•

PGAM5 catalytic domain shares phosphoglycerate mutase fold and forms stable dimer

•

WDPNWD motif allosterically activates the fully active dodecameric form

•

Crystal structures reveal conformational plasticity of the PGAM5 active site

Identification of productive and futile encounters in an electron transfer protein complex

Well-defined, stereospecific states in protein complexes are often in exchange with an ensemble of more dynamic orientations: the encounter states. The structure of the stereospecific complex between cytochrome P450cam and putidaredoxin was solved recently by X-ray diffraction as well as paramagnetic NMR spectroscopy. Other than the stereospecific complex, the NMR data clearly show the presence of additional states in the complex in solution. In these encounter states, populated for a small percentage of the time, putidaredoxin assumes multiple orientations and samples a large part of the surface of cytochrome P450cam. To characterize the nature of the encounter states, an extensive paramagnetic NMR dataset has been analyzed using the Maximum Occurrence of Regions methodology. The analysis reveals the location and maximal spatial extent of the additional states needed to fully explain the NMR data. Under the assumption of sparsity of the size of the conformational ensemble, several minor states can be located quite precisely. The distribution of these minor states correlates with the electrostatic potential map around cytochrome P450cam. Whereas some minor states are on isolated positively charged patches, others are connected to the stereospecific site via positively charged paths. The existence of electrostatically favorable pathways between the stereospecific interaction site and the different minor states or lack thereof suggests a means to discriminate between productive and futile encounter states.

Low Resolution Refinement of Atomic Models Against Crystallographic Data

This review describes some of the problems encountered during low-resolution refinement and map calculation. Refinement is considered as an application of Bayes' theorem, allowing combination of information from various sources including crystallographic experimental data and prior chemical and structural knowledge. The sources of prior knowledge relevant to macromolecules include basic chemical information such as bonds and angles, structural information from reference models of known homologs, knowledge about secondary structures, hydrogen bonding patterns, and similarity of non-crystallographically related copies of a molecule. Additionally, prior information encapsulating local conformational conservation is exploited, keeping local interatomic distances similar to those in the starting atomic model. The importance of designing an accurate likelihood function-the only link between model parameters and observed data-is emphasized. The review also reemphasizes the importance of phases, and describes how the use of raw observed amplitudes could give a better correlation between the calculated and "true" maps. It is shown that very noisy or absent observations can be replaced by calculated structure factors, weighted according to the accuracy of the atomic model. This approach helps to smoothen the map. However, such replacement should be used sparingly, as the bias toward errors in the model could be too much to avoid. It is in general recommended that, whenever a new map is calculated, map quality should be judged by inspection of the parts of the map where there is no atomic model. It is also noted that it is advisable to work with multiple blurred and sharpened maps, as different parts of a crystal may exhibit different degrees of mobility. Doing so can allow accurate building of atomic models, accounting for overall shape as well as finer structural details. Some of the results described in this review have been implemented in the programs REFMAC5, ProSMART and LORESTR, which are available as part of the CCP4 software suite.

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.

In-Cell Protein Structures from 2D NMR Experiments

In-cell NMR spectroscopy provides atomic resolution insights into the structural properties of proteins in cells, but it is rarely used to solve entire protein structures de novo. Here, we introduce a paramagnetic lanthanide-tag to simultaneously measure protein pseudocontact shifts (PCSs) and residual dipolar couplings (RDCs) to be used as input for structure calculation routines within the Rosetta program. We employ this approach to determine the structure of the protein G B1 domain (GB1) in intact Xenopus laevis oocytes from a single set of 2D in-cell NMR experiments. Specifically, we derive well-defined GB1 ensembles from low concentration in-cell NMR samples (∼50 μM) measured at moderate magnetic field strengths (600 MHz), thus offering an easily accessible alternative for determining intracellular protein structures.

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.

Improved Accuracy from Joint X-ray and NMR Refinement of a Protein-RNA Complex Structure

Integrated experimental approaches play an increasingly important role in structural biology, taking advantage of the complementary information provided by different techniques. In particular, the combination of NMR data with X-ray diffraction patterns may provide accurate and precise information about local conformations not available from average-resolution X-ray structures alone. Here, we refined the structure of a ternary protein-protein-RNA complex comprising three domains, Sxl and Unr, bound to a single-stranded region derived in the msl2 mRNA. The joint X-ray and NMR refinement reveals that-despite the poor quality of the fit found for the original structural model-the NMR data can be largely accommodated within the uncertainty in the atom positioning (structural noise) from the primary X-ray data and that the overall domain arrangements and binding interfaces are preserved on passing from the crystalline state to the solution. The refinement highlights local conformational differences, which provide additional information on specific features of the structure. For example, conformational dynamics and heterogeneity observed at the interface between the CSD1 and the Sxl protein components in the ternary complex are revealed by the combination of NMR and crystallographic data. The joint refinement protocol offers unique opportunities to detect structural differences arising from various experimental conditions and reveals static or dynamic differences in the conformation of the biomolecule between the solution and the crystals.

Delicate conformational balance of the redox enzyme cytochrome P450cam

The energy landscapes of proteins are highly complex and can be influenced by changes in physical and chemical conditions under which the protein is studied. The redox enzyme cytochrome P450cam undergoes a multistep catalytic cycle wherein two electrons are transferred to the heme group and the enzyme visits several conformational states. Using paramagnetic NMR spectroscopy with a lanthanoid tag, we show that the enzyme bound to its redox partner, putidaredoxin, is in a closed state at ambient temperature in solution. This result contrasts with recent crystal structures of the complex, which suggest that the enzyme opens up when bound to its partner. The closed state supports a model of catalysis in which the substrate is locked in the active site pocket and the enzyme acts as an insulator for the reactive intermediates of the reaction.

Residue-level global and local ensemble-ensemble comparisons of protein domains

Many methods of protein structure generation such as NMR-based solution structure determination and template-based modeling do not produce a single model, but an ensemble of models consistent with the available information. Current strategies for comparing ensembles lose information because they use only a single representative structure. Here, we describe the ENSEMBLATOR and its novel strategy to directly compare two ensembles containing the same atoms to identify significant global and local backbone differences between them on per-atom and per-residue levels, respectively. The ENSEMBLATOR has four components: eePREP (ee for ensemble-ensemble), which selects atoms common to all models; eeCORE, which identifies atoms belonging to a cutoff-distance dependent common core; eeGLOBAL, which globally superimposes all models using the defined core atoms and calculates for each atom the two intraensemble variations, the interensemble variation, and the closest approach of members of the two ensembles; and eeLOCAL, which performs a local overlay of each dipeptide and, using a novel measure of local backbone similarity, reports the same four variations as eeGLOBAL. The combination of eeGLOBAL and eeLOCAL analyses identifies the most significant differences between ensembles. We illustrate the ENSEMBLATOR's capabilities by showing how using it to analyze NMR ensembles and to compare NMR ensembles with crystal structures provides novel insights compared to published studies. One of these studies leads us to suggest that a "consistency check" of NMR-derived ensembles may be a useful analysis step for NMR-based structure determinations in general. The ENSEMBLATOR 1.0 is available as a first generation tool to carry out ensemble-ensemble comparisons.

Integrative, dynamic structural biology at atomic resolution--it's about time

Biomolecules adopt a dynamic ensemble of conformations, each with the potential to interact with binding partners or perform the chemical reactions required for a multitude of cellular functions. Recent advances in X-ray crystallography, nuclear magnetic resonance (NMR) spectroscopy and other techniques are helping us realize the dream of seeing--in atomic detail--how different parts of biomolecules shift between functional substates using concerted motions. Integrative structural biology has advanced our understanding of the formation of large macromolecular complexes and how their components interact in assemblies by leveraging data from many low-resolution methods. Here, we review the growing opportunities for integrative, dynamic structural biology at the atomic scale, contending there is increasing synergistic potential between X-ray crystallography, NMR and computer simulations to reveal a structural basis for protein conformational dynamics at high resolution.

NMR studies of dynamic biomolecular conformational ensembles

Multidimensional heteronuclear NMR approaches can provide nearly complete sequential signal assignments of isotopically enriched biomolecules. The availability of assignments together with measurements of spin relaxation rates, residual spin interactions, J-couplings and chemical shifts provides information at atomic resolution about internal dynamics on timescales ranging from ps to ms, both in solution and in the solid state. However, due to the complexity of biomolecules, it is not possible to extract a unique atomic-resolution description of biomolecular motions even from extensive NMR data when many conformations are sampled on multiple timescales. For this reason, powerful computational approaches are increasingly applied to large NMR data sets to elucidate conformational ensembles sampled by biomolecules. In the past decade, considerable attention has been directed at an important class of biomolecules that function by binding to a wide variety of target molecules. Questions of current interest are: "Does the free biomolecule sample a conformational ensemble that encompasses the conformations found when it binds to various targets; and if so, on what time scale is the ensemble sampled?" This article reviews recent efforts to answer these questions, with a focus on comparing ensembles obtained for the same biomolecules by different investigators. A detailed comparison of results obtained is provided for three biomolecules: ubiquitin, calmodulin and the HIV-1 trans-activation response RNA.

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.

High accuracy of Karplus equations for relating three-bond J couplings to protein backbone torsion angles

(3) JC'C' and (3) JHNHα couplings are related to the intervening backbone torsion angle ${\varphi }$ by standard Karplus equations. Although these couplings are known to be affected by parameters other than ${\varphi }$, including H-bonding, valence angles and residue type, experimental results and quantum calculations indicate that the impact of these latter parameters is typically very small. The solution NMR structure of protein GB3, newly refined by using extensive sets of residual dipolar couplings, yields 50-60 % better Karplus equation agreement between ${\varphi }$ angles and experimental (3) JC'C' and (3) JHNHα values than does the high-resolution X-ray structure. In intrinsically disordered proteins, (3) JC'C' and (3) JHNHα couplings can be measured at even higher accuracy, and the impact of factors other than the intervening torsion angle on (3) J will be smaller than in folded proteins, making these couplings exceptionally valuable reporters on the ensemble of ${\varphi }$ angles sampled by each residue.

E pluribus unum, no more: from one crystal, many conformations

Several distinct computational approaches have recently been implemented to represent conformational heterogeneity from X-ray crystallography datasets that are averaged in time and space. As these modeling methods mature, newly discovered alternative conformations are being used to derive functional protein mechanisms. Room temperature X-ray data collection is emerging as a key variable for sampling functionally relevant conformations also observed in solution studies. Although concerns about radiation damage are warranted with higher temperature data collection, 'diffract and destroy' strategies on X-ray free electron lasers may permit radiation damage-free data collection. X-ray crystallography need not be confined to 'static unique snapshots'; these experimental and computational advances are revealing how the many conformations populated within a single crystal are used in biological mechanisms.