Paramagnetism-based NMR restraints provide maximum allowed probabilities for the different conformations of partially independent protein domains

Are Protein Conformational Ensembles in Agreement with Experimental Data? A Geometrical Interpretation of the Problem

The conformational variability of biological macromolecules can play an important role in their biological function. Therefore, understanding conformational variability is expected to be key for predicting the behavior of a particular molecule in the context of organism-wide studies. Several experimental methods have been developed and deployed for accessing this information, and computational methods are continuously updated for the profitable integration of different experimental sources. The outcome of this endeavor is conformational ensembles, which may vary significantly in properties and composition when different ensemble reconstruction methods are used, and this raises the issue of comparing the predicted ensembles against experimental data. In this article, we discuss a geometrical formulation to provide a framework for understanding the agreement of an ensemble prediction to the experimental observations.

Exploring the impact of alignment media on RDC analysis of phosphorus-containing compounds: a molecular docking approach

Residual dipolar couplings (RDCs) are employed in NMR analysis when conventional methods, such as J-couplings and nuclear Overhauser effects (NOEs) fail. Low-energy (optimized) conformers are often used as input structures in RDC analysis programs. However, these low-energy structures do not necessarily resemble conformations found in anisotropic environments due to interactions with the alignment medium, especially if the analyte molecules are flexible. Considering interactions with alignment media in RDC analysis, we developed and evaluated a molecular docking-based approach to generate more accurate conformer ensembles for compounds in the presence of the poly-γ-benzyl-L-glutamate alignment medium. We designed chiral phosphorus-containing compounds that enabled us to utilize 31P NMR parameters for the stereochemical analysis. Using P3D/PALES software to evaluate diastereomer discrimination, we found that our conformer ensembles outperform moderately the standard, low-energy conformers in RDC analysis. To further improve our results, we (i) averaged the experimental values of the molecular docking-based conformers by applying the Boltzmann distribution and (ii) optimized the structures through normal mode relaxation, thereby enhancing the Pearson correlation factor R and even diastereomer discrimination in some cases. Nevertheless, we presume that significant differences between J-couplings in isotropic and in anisotropic environments may preclude RDC measurements for flexible molecules. Therefore, generating conformer ensembles based on molecular docking enhances RDC analysis for mildly flexible systems while flexible molecules may require applying more advanced approaches, in particular approaches including dynamical effects.

Functional and Pathological Effects of α-Synuclein on Synaptic SNARE Complexes

α-Synuclein is an abundant protein at the neuronal synapse that has been implicated in Parkinson's disease for over 25 years and characterizes the hallmark pathology of a group of neurodegenerative diseases now known as the synucleinopathies. Physiologically, α-synuclein exists in an equilibrium between a synaptic vesicle membrane-bound α-helical multimer and a cytosolic largely unstructured monomer. Through its membrane-bound state, α-synuclein functions in neurotransmitter release by modulating several steps in the synaptic vesicle cycle, including synaptic vesicle clustering and docking, SNARE complex assembly, and homeostasis of synaptic vesicle pools. These functions have been ascribed to α-synuclein's interactions with the synaptic vesicle SNARE protein VAMP2/synaptobrevin-2, the synaptic vesicle-attached synapsins, and the synaptic vesicle membrane itself. How α-synuclein affects these processes, and whether disease is due to loss-of-function or gain-of-toxic-function of α-synuclein remains unclear. In this review, we provide an in-depth summary of the existing literature, discuss possible reasons for the discrepancies in the field, and propose a working model that reconciles the findings in the literature.

Quantitative Agreement between Conformational Substates of Holo Calcium-Loaded Calmodulin Detected by Double Electron-Electron Resonance EPR and Predicted by Molecular Dynamics Simulations

Calcium-loaded calmodulin (CaM/4Ca²⁺) comprises two domains that undergo rigid body reorientation from a predominantly extended conformation to a compact one upon binding target peptides. A recent replica-exchange molecular dynamics (MD) simulation on holo CaM/4Ca²⁺ suggested the existence of distinct structural clusters (substates) along the path from extended to compact conformers in the absence of substrates. Here, we experimentally demonstrate the existence of CaM/4Ca²⁺ substates trapped in local minima by three freezing/annealing regimes (slow, 40 s; intermediate, 1.5 s; fast, 0.5 ms) using pulsed Q-band double electron-electron resonance (DEER) EPR spectroscopy to measure interdomain distances between nitroxide spin-labels positioned at A17C and A128C in the N- and C-terminal domains, respectively. The DEER echo curves were directly fit to population-optimized P(r) pairwise distance distributions calculated from the coordinates of the MD clusters and compact crystal structure. DEER data on fully deuterated CaM/4Ca²⁺ were acquired at multiple values of the second echo period (10-35 μs) and analyzed globally to eliminate instrumental and overfitting artifacts and ensure accurate populations, peak positions, and widths. The DEER data for all three freezing regimes are quantitatively accounted for within experimental error by 5-6 distinct conformers comprising a predominantly populated extended form (60-75%) and progressively more compact states whose populations decrease as the degree of compactness increases. The shortest interdomain separation is found in the compact crystal structure, which has an occupancy of 4-6%. Thus, CaM/4Ca²⁺ samples high energy local minima comprising a few discrete substates of increasing compactness in a rugged energy landscape.

Pathogenic Impact of α-Synuclein Phosphorylation and Its Kinases in α-Synucleinopathies

α-Synuclein is a protein with a molecular weight of 14.5 kDa and consists of 140 amino acids encoded by the SNCA gene. Missense mutations and gene duplications in the SNCA gene cause hereditary Parkinson’s disease. Highly phosphorylated and abnormally aggregated α-synuclein is a major component of Lewy bodies found in neuronal cells of patients with sporadic Parkinson’s disease, dementia with Lewy bodies, and glial cytoplasmic inclusion bodies in oligodendrocytes with multiple system atrophy. Aggregated α-synuclein is cytotoxic and plays a central role in the pathogenesis of the above-mentioned synucleinopathies. In a healthy brain, most α-synuclein is unphosphorylated; however, more than 90% of abnormally aggregated α-synuclein in Lewy bodies of patients with Parkinson’s disease is phosphorylated at Ser129, which is presumed to be of pathological significance. Several kinases catalyze Ser129 phosphorylation, but the role of phosphorylation enzymes in disease pathogenesis and their relationship to cellular toxicity from phosphorylation are not fully understood in α-synucleinopathy. Consequently, this review focuses on the pathogenic impact of α-synuclein phosphorylation and its kinases during the neurodegeneration process in α-synucleinopathy.

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.

On the use of residual dipolar couplings in multi-state structure calculation of two-domain proteins

Residual dipolar couplings (RDCs) are powerful nuclear magnetic resonance (NMR) probes for the structure calculation of biomacromolecules. Typically, an alignment tensor that defines the orientation of the entire molecule relative to the magnetic field is determined either before refinement of individual bond vectors or simultaneously with this refinement. For single-domain proteins this approach works well since all bond vectors can be described within the same coordinate frame, which is given by the alignment tensor. However, novel approaches are sought after for systems where no universal alignment tensor can be used. Here, we present an approach that can be applied to two-domain proteins that enables the calculation of multiple states within each domain as well as with respect to the relative positions of the two domains.

DEER experiments reveal fundamental differences between calmodulin complexes with IQ and MARCKS peptides in solution

Calmodulin (CaM) is a calcium-binding protein that regulates the function of many proteins by indirectly conferring Ca²⁺ sensitivity, and it undergoes a large conformational change on partners' binding. We compared the solution binding mode of the target peptides MARCKS and IQ by double electron-electron resonance (DEER) distance measurements and paramagnetic NMR. We combined nitroxide and Gd(III) spin labels, including specific substitution of one of the Ca²⁺ ions in the CaM mutant N60D by a Gd(III) ion. The binding of MARCKS to holo-CaM resulted neither in a closed conformation nor in a unique relative orientation between the two CaM domains, in contrast with the crystal structure. Binding of IQ to holo-CaM did generate a closed conformation. Using elastic network modeling and 12 distance restraints obtained from multiple holo-CaM/IQ DEER data, we derived a model of the solution structure, which is in reasonable agreement with the crystal structure.

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.

Reconstruction of Coupled Intra- and Interdomain Protein Motion from Nuclear and Electron Magnetic Resonance

Proteins composed of multiple domains allow for structural heterogeneity and interdomain dynamics that may be vital for function. Intradomain structures and dynamics can influence interdomain conformations and vice versa. However, no established structure determination method is currently available that can probe the coupling of these motions. The protein Pin1 contains separate regulatory and catalytic domains that sample "extended" and "compact" states, and ligand binding changes this equilibrium. Ligand binding and interdomain distance have been shown to impact the activity of Pin1, suggesting interdomain allostery. In order to characterize the conformational equilibrium of Pin1, we describe a novel method to model the coupling between intra- and interdomain dynamics at atomic resolution using multistate ensembles. The method uses time-averaged nuclear magnetic resonance (NMR) restraints and double electron-electron resonance (DEER) data that resolve distance distributions. While the intradomain calculation is primarily driven by exact nuclear Overhauser enhancements (eNOEs), J couplings, and residual dipolar couplings (RDCs), the relative domain distribution is driven by paramagnetic relaxation enhancement (PREs), RDCs, interdomain NOEs, and DEER. Our data support a 70:30 population of the compact and extended states in apo Pin1. A multistate ensemble describes these conformations simultaneously, with distinct conformational differences located in the interdomain interface stabilizing the compact or extended states. We also describe correlated conformations between the catalytic site and interdomain interface that may explain allostery driven by interdomain contact.

NMR for Single Ion Magnets

Nuclear Magnetic Resonance is particularly sensitive to the electronic structure of matter and is thus a powerful tool to characterize in-depth the magnetic properties of a system. NMR is indeed increasingly recognized as an ideal tool to add precious structural information for the development of Single Ion Magnets, small complexes that are recently gaining much popularity due to their quantum computing and spintronics applications. In this review, we recall the theoretical principles of paramagnetic NMR, with particular attention to lanthanoids, and we give an overview of the recent advances in this field.

De novo designed coiled coils as scaffolds for lanthanides, including novel imaging agents with a twist

The design of artificial miniature lanthanide proteins, provide an opportunity to access new functional metalloproteins as well as insight into native lanthanide biochemistry.

Mechanism and Inhibition of Matrix Metalloproteinases

Matrix metalloproteinases hydrolyze proteins and glycoproteins forming the extracellular matrix, cytokines and growth factors released in the extracellular space, and membrane-bound receptors on the outer cell membrane. The pathological relevance of MMPs has prompted the structural and functional characterization of these enzymes and the development of synthetic inhibitors as possible drug candidates. Recent studies have provided a better understanding of the substrate preference of the different members of the family, and structural data on the mechanism by which these enzymes hydrolyze the substrates. Here, we report the recent advancements in the understanding of the mechanism of collagenolysis and elastolysis, and we discuss the perspectives of new therapeutic strategies for targeting MMPs.

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’.

α-synuclein in the pathophysiology of Alzheimer's disease

The Alzheimer’s disease (AD) afflicted brain is neuropathologically defined by extracellular amyloid-β (Aβ) plaques and intraneuronal neurofibrillary tangles composed of hyperphosphorylated tau protein. However, accumulating evidence suggests that the presynaptic protein α-synuclein (αSyn), mainly associated with synucleinopathies like Parkinson’s disease (PD), dementia with Lewy bodies (DLB) and multiple system atrophy (MSA), is involved in the pathophysiology of AD. Lewy-related pathology (LRP), primarily comprised of αSyn, is present in a majority of autopsied AD brains, and higher levels of αSyn in the cerebrospinal fluid (CSF) of patients with mild cognitive impairment (MCI) and AD have been linked to cognitive decline. Recent studies also suggest that the asymptomatic accumulation of Aβ plaques is associated with higher CSF αSyn levels in subjects at risk of sporadic AD and in individuals carrying autosomal dominant AD mutations. Experimental evidence has further linked αSyn mainly to tau hyperphosphorylation, but also to the pathological actions of Aβ and the APOEε4 allele, the latter being a major genetic risk factor for both AD and DLB. In this review, we provide a summary of the current evidence proposing an involvement of αSyn either as an active or passive player in the pathophysiological ensemble of AD, and furthermore describe in detail the current knowledge of αSyn structure and inferred function.

13C NMR-Based Approaches for Solving Challenging Stereochemical Problems

Determining the configuration of proton-deficient molecules is challenging using conventional NMR methods including nuclear Overhauser effect (NOE) and the proton-dependent J-based configuration analysis (JBCA). The problem is exacerbated when only one stereoisomer is available. Alternative methods based on the utilization of ¹³C NMR chemical shifts, ¹³C-¹³C homonuclear couplings measured at natural abundance, and residual chemical shift anisotropy measurements in conjunction with density functional theory calculations are illustrated with a proton-deficient model compound.

The theory of frame ordering: observing motions in calmodulin complexes

Large scale functional motions of molecules are studied experimentally using numerous molecular and biophysics techniques, the data from which are subsequently interpreted using diverse models of Brownian molecular dynamics. To unify all rotational physics techniques and motional models, the frame order tensor - a universal statistical mechanics theory based on the rotational ordering of rigid body frames - is herein formulated. The frame ordering is the fundamental physics that governs how motions modulate rotational molecular physics and it defines the properties and maximum information content encoded in the observable physics. Using the tensor to link residual dipolar couplings and pseudo-contact shifts, two distinct information-rich and atomic-level biophysical measurements from the field of nuclear magnetic resonance spectroscopy, to a number of basic mechanical joint models, a highly dynamic state of calmodulin (CaM) bound to a target peptide in a tightly closed conformation was observed. Intra- and inter-domain motions reveal the CaM complex to be entropically primed for peptide release.

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.

Treating Biomacromolecular Conformational Variability

The function of a biomacromolecule is related not only to its structure but also to the different conformations that its structural elements can sample. It is therefore important to determine the extent of the structural fluctuations and to identify the states that are actually populated as a result of the rearrangement. However, this accomplishment is undermined by an intrinsic limitation: the amount of experimental data is by and large inferior to the number of the states that a biomacromolecule can actually sample. This means that additional, a priori information must be applied in order to derive the most from the available experimental data but not to run into overinterpretation. In this chapter we will give a summary of the experimental observables that can be used towards the reconstruction of structural ensembles, how the data can be profitably combined and to what extent the data are affected by error; finally we will give an overview of the computational methods that have been developed to model structural ensembles, highlighting their difference and similarities, advantages and disadvantages.

Ligand Proton Pseudocontact Shifts Determined from Paramagnetic Relaxation Dispersion in the Limit of NMR Intermediate Exchange

Delineation of protein-ligand interaction modes is key for rational drug discovery. The availability of complex crystal structures is often limited by the aqueous solubility of the compounds, while lead-like compounds with micromolar affinities normally fall into the NMR intermediate exchange regime, in which severe line broadening to beyond the detection of interfacial resonances limits NMR applications. Here, we developed a new method to retrieve low-populated bound-state ¹H pseudocontact shifts (PCSs) using paramagnetic relaxation dispersion (RD). We evaluated using a ¹H PCS-RD approach in a BRM bromodomain lead-like inhibitor to filter molecular docking poses using multiple intermolecular structural restraints. Considering the universal presence of proton atoms in druglike compounds, our work will have wide application in structure-guided drug discovery even under an extreme condition of NMR intermediate exchange and low aqueous solubility of ligands.

Continuous Interdomain Orientation Distributions Reveal Components of Binding Thermodynamics

The flexibility of biological macromolecules is an important structural determinant of function. Unfortunately, the correlations between different motional modes are poorly captured by discrete ensemble representations. Here, we present new ways to both represent and visualize correlated interdomain motions. Interdomain motions are determined directly from residual dipolar couplings, represented as a continuous conformational distribution, and visualized using the disk-on-sphere representation. Using the disk-on-sphere representation, features of interdomain motions, including correlations, are intuitively visualized. The representation works especially well for multidomain systems with broad conformational distributions.This analysis also can be extended to multiple probability density modes, using a Bingham mixture model. We use this new paradigm to study the interdomain motions of staphylococcal protein A, which is a key virulence factor contributing to the pathogenicity of Staphylococcus aureus. We capture the smooth transitions between important states and demonstrate the utility of continuous distribution functions for computing the reorientational components of binding thermodynamics. Such insights allow for the dissection of the dynamic structural components of functionally important intermolecular interactions.

Modulation of Correlated Segment Fluctuations in IDPs upon Complex Formation as an Allosteric Regulatory Mechanism

Molecular recognition of and by intrinsically disordered proteins (IDPs) is an intriguing and still largely elusive phenomenon. Typically, protein recognition involving IDPs requires either folding upon binding or, alternatively, the formation of "fuzzy complexes." Here we show via correlation analyses of paramagnetic relaxation enhancement data unprecedented and striking alterations of the concerted fluctuations within the conformational ensemble of IDPs upon ligand binding. We study the binding of α-synuclein to calmodulin, a ubiquitous calcium-binding protein, and the binding of the extracellular matrix IDP osteopontin to heparin, a mimic of the extracellular matrix ligand hyaluronic acid. In both cases, binding leads to reduction of correlated long-range motions in these two IDPs and thus indicates a loosening of structural compaction upon binding. Most importantly, however, the simultaneous presence of correlated and anti-correlated fluctuations in IDPs suggests the prevalence of "energetic frustration" and provides an explanation for the puzzling observation of disordered allostery in IDPs.

Cell Biology and Pathophysiology of α-Synuclein

α-Synuclein is an abundant neuronal protein that is highly enriched in presynaptic nerve terminals. Genetics and neuropathology studies link α-synuclein to Parkinson's disease (PD) and other neurodegenerative disorders. Accumulation of misfolded oligomers and larger aggregates of α-synuclein defines multiple neurodegenerative diseases called synucleinopathies, but the mechanisms by which α-synuclein acts in neurodegeneration are unknown. Moreover, the normal cellular function of α-synuclein remains debated. In this perspective, we review the structural characteristics of α-synuclein, its developmental expression pattern, its cellular and subcellular localization, and its function in neurons. We also discuss recent progress on secretion of α-synuclein, which may contribute to its interneuronal spread in a prion-like fashion, and describe the neurotoxic effects of α-synuclein that are thought to be responsible for its role in neurodegeneration.

Pseudocontact shifts from mobile spin labels

This paper presents a detailed analysis of the pseudocontact shift (PCS) field induced by a mobile spin label that is viewed as a probability density distribution with an associated effective magnetic susceptibility anisotropy. It is demonstrated that non-spherically symmetric density can lead to significant deviations from the commonly used point dipole approximation for the PCS. Analytical and numerical solutions are presented for the general partial differential equation that describes the non-point case. It is also demonstrated that it is possible, with some reasonable approximations, to reconstruct paramagnetic centre probability distributions from the experimental PCS data.

Conformational landscape mapping the difference between N-lobes and C-lobes of calmodulin

Calmodulin is a calcium binding protein that consists of four EF-hand domains. The two EF-lobes of calmodulin, called the N-lobe and the C-lobe, arose from duplication and fusion of a precursor EF-hand. The amino acid sequences and the structures of the N-lobe and of the C-lobe are quite similar to each other. The N-lobe and the C-lobe, however, have subtle differences in structure and function. We analyzed the helix positions of calmodulin lobes by the alignment with the pseudo-two fold axis of the EF-lobe. We made a map of conformational landscape of helix positions. The four states of the EF-lobe appeared on two lines in the landscape; these two lines show the trajectory of opening and closing of the EF-lobe. For the N-lobe of calmodulin, the calcium bound form and the apo-forms are on the lower line. The two apo-forms of the C-lobe of calmodulin, with target and without target, are on the upper line. The calcium bound form of the C-lobe is on the lower line. The rearrangement of helix interaction between two the EF-hands is necessary for calcium binding in the C-lobe. The hydrophobic packing in the apo-form of the N-lobe is similar to the packing of the N- and C-lobes of the calcium bound form. However, the packing of C-lobe side chains in the apo-form is different from these other three structures. Our detailed analysis should serve as an example that can be applied to other proteins that undergo changes in conformation upon binding effectors.

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.

Recent excitements in protein NMR: Large proteins and biologically relevant dynamics

The advent of Transverse Relaxation Optimized SpectroscopY (TROSY) and perdeuteration allowed biomolecular NMR spectroscopists to overcome the size limitation barrier (approx. 20 kDa) in de novo structure determination of proteins. The utility of these techniques was immediately demonstrated on large proteins and protein complexes (e.g. GroELGroES, ClpP protease, Hsp90-p53, 20S proteasome, etc.). Further, recent methodological developments such as Residual Dipolar Couplings and Paramagnetic Relaxation Enhancement allowed accurate measurement of long-range structural restraints. Additionally, Carr-Purcell-Meiboom-Gill (CPMG), rotating frame relaxation experiments (R1(rho)) and saturation transfer experiments (CEST and DEST) created never-before accessibility to the (mu)s-ms timescale dynamic parameters that led to the deeper understanding of biological processes. Meanwhile, the excitement in the field continued with a series of developments in the fast data acquisition methods allowing rapid structural studies on less stable proteins. This review aims to discuss important developments in the field of biomolecular NMR spectroscopy in the recent past, i.e., in the post TROSY era. These developments not only gave access to the structural studies of large protein assemblies, but also revolutionized tools in the arsenal of today's biomolecular NMR and point to a bright future of biomolecular NMR spectroscopy.

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.

Implication of Alpha-Synuclein Phosphorylation at S129 in Synucleinopathies: What Have We Learned in the Last Decade?

Abnormal accumulation of proteinaceous intraneuronal inclusions called Lewy bodies (LBs) is the neurpathological hallmark of Parkinson’s disease (PD) and related synucleinopathies. These inclusions are mainly constituted of a presynaptic protein, α-synuclein (α-syn). Over the past decade, growing amounts of studies reported an aberrant accumulation of phosphorylated α-syn at the residue S129 (pS129) in the brain of patients suffering from PD, as well as in transgenic animal models of synucleinopathies. Whereas only a small fraction of α-syn (<4%) is phosphorylated in healthy brains, a dramatic accumulation of pS129 (>90%) has been observed within LBs, suggesting that this post-translational modification may play an important role in the regulation of α-syn aggregation, LBs formation and neuronal degeneration. However, whether phosphorylation at S129 suppresses or enhances α-syn aggregation and toxicity in vivo remains a subject of active debate. The answer to this question has important implications for understanding the role of phosphorylation in the pathogenesis of synucleinopathies and determining if targeting kinases or phosphatases could be a viable therapeutic strategy for the treatment of these devastating neurological disorders. In the present review, we explore recent findings from in vitro, cell-based assays and in vivo studies describing the potential implications of pS129 in the regulation of α-syn physiological functions, as well as its implication in synucleinopathies pathogenesis and diagnosis.

Inter-helical conformational preferences of HIV-1 TAR-RNA from maximum occurrence analysis of NMR data and molecular dynamics simulations

Detecting conformational heterogeneity in biological macromolecules is a key for the understanding of their biological function. We here provide a comparison between two independent approaches to assess conformational heterogeneity: molecular dynamics simulations, performed without inclusion of any experimental data, and maximum occurrence (MaxOcc) distribution over the topologically available conformational space. The latter only reflects the extent of the averaging and identifies regions which are most compliant with the experimentally measured NMR Residual Dipolar Couplings (RDCs). The analysis was performed for the HIV-1 TAR RNA, consisting of two helical domains connected by a flexible bulge junction, for which four sets of RDCs were available as well as an 8.2 μs all-atom molecular dynamics simulation. A sample and select approach was previously applied to extract from the molecular dynamics trajectory conformational ensembles in agreement with the four sets of RDCs. The MaxOcc analysis performed here identifies the most likely sampled region in the conformational space of the system which, strikingly, overlaps well with the structures independently sampled in the molecular dynamics calculations and even better with the RDC selected ensemble.

A critical assessment of methods to recover information from averaged data

Conformational heterogeneity is key to the function of many biomacromolecules, but only a few groups have tried to characterize it until recently. Now, thanks to the increased throughput of experimental data and the increased computational power, the problem of the characterization of protein structural variability has become more and more popular. Several groups have devoted their efforts in trying to create quantitative, reliable and accurate protocols for extracting such information from averaged data. We analyze here different approaches, discussing strengths and weaknesses of each. All approaches can roughly be clustered into two groups: those satisfying the maximum entropy principle and those recovering ensembles composed of a restricted number of molecular conformations. In the first case, the solution focuses on the features that are common to all the infinite solutions satisfying the experimental data; in the second case, the reconstructed ensemble shows the conformational regions where a large probability can be placed. The upper limits for conformational probabilities (MaxOcc) can also be calculated. We also give an overview of the mainstream experimental observables, with considerations on the assumptions underlying their usage.

Investigating the Role of Large-Scale Domain Dynamics in Protein-Protein Interactions

Intrinsically disordered linkers provide multi-domain proteins with degrees of conformational freedom that are often essential for function. These highly dynamic assemblies represent a significant fraction of all proteomes, and deciphering the physical basis of their interactions represents a considerable challenge. Here we describe the difficulties associated with mapping the large-scale domain dynamics and describe two recent examples where solution state methods, in particular NMR spectroscopy, are used to investigate conformational exchange on very different timescales.

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.

Investigating protein conformational energy landscapes and atomic resolution dynamics from NMR dipolar couplings: a review

Nuclear magnetic resonance spectroscopy is exquisitely sensitive to protein dynamics. In particular inter-nuclear dipolar couplings, that become measurable in solution when the protein is dissolved in a dilute liquid crystalline solution, report on all conformations sampled up to millisecond timescales. As such they provide the opportunity to describe the Boltzmann distribution present in solution at atomic resolution, and thereby to map the conformational energy landscape in unprecedented detail. The development of analytical methods and approaches based on numerical simulation and their application to numerous biologically important systems is presented.

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

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

Compact, hydrophilic, lanthanide-binding tags for paramagnetic NMR spectroscopy

The design, synthesis and evaluation of four novel lanthanide-binding tags for paramagnetic NMR spectroscopy are reported.

The design, synthesis and evaluation of four novel lanthanide-binding tags for paramagnetic NMR spectroscopy are reported. Each tag is based on the ((2S,2′S,2′′S,2′′′S)-1,1′,1′′,1′′′-(1,4,7,10-tetraazacyclododecane-1,4,7,10-tetrayl)tetrakis(propan-2-ol)) scaffold, featuring small chiral alcohol coordinating pendants to minimise the size and hydrophobic character of each tag. The tags feature different linkers of variable length for conjugation to protein via a single cysteine residue. Each tag's ability to induce pseudocontact shifts (PCS) was assessed on a ubiquitin A28C mutant. Two enantiomeric tags of particular note, C7 and C8, produced significantly larger Δχ-tensors compared to a previously developed tag, C1, attributed to the extremely short linker utilised, limiting the mobility of the bound lanthanide ion. The C7 and C8 tags' capacity to induce PCSs was further demonstrated on GB1 Q32C and 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) S112C/C80A mutants. Whilst factors such as the choice of lanthanide ion, pH and site of conjugation influence the size of the PCSs obtained, the tags represent a significant advance in the field.

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.

Visualizing transient dark states by NMR spectroscopy

Myriad biological processes proceed through states that defy characterization by conventional atomic-resolution structural biological methods. The invisibility of these 'dark' states can arise from their transient nature, low equilibrium population, large molecular weight, and/or heterogeneity. Although they are invisible, these dark states underlie a range of processes, acting as encounter complexes between proteins and as intermediates in protein folding and aggregation. New methods have made these states accessible to high-resolution analysis by nuclear magnetic resonance (NMR) spectroscopy, as long as the dark state is in dynamic equilibrium with an NMR-visible species. These methods - paramagnetic NMR, relaxation dispersion, saturation transfer, lifetime line broadening, and hydrogen exchange - allow the exploration of otherwise invisible states in exchange with a visible species over a range of timescales, each taking advantage of some unique property of the dark state to amplify its effect on a particular NMR observable. In this review, we introduce these methods and explore two specific techniques - paramagnetic relaxation enhancement and dark state exchange saturation transfer - in greater detail.

Ensemble Calculation for Intrinsically Disordered Proteins Using NMR Parameters

Intrinsically disordered proteins (IDPs) perform their function despite their lack of well-defined tertiary structure. Residual structure has been observed in IDPs, commonly described as transient/dynamic or expressed in terms of fractional populations. In order to understand how the protein primary sequence dictates the dynamic and structural properties of IDPs and in general to understand how IDPs function, atomic-level descriptions are needed. Nuclear magnetic resonance spectroscopy provides information about local and long-range structure in IDPs at amino acid specific resolution and can be used in combination with ensemble descriptions to represent the dynamic nature of IDPs. In this chapter we describe sample-and-select approaches for ensemble modelling of local structural propensities in IDPs with specific emphasis on validation of these ensembles.

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.

Alpha-synuclein function and dysfunction on cellular membranes

Alpha-synuclein is a small neuronal protein that is closely associated with the etiology of Parkinson's disease. Mutations in and alterations in expression levels of alpha-synuclein cause autosomal dominant early onset heredity forms of Parkinson's disease, and sporadic Parkinson's disease is defined in part by the presence of Lewy bodies and Lewy neurites that are composed primarily of alpha-synuclein deposited in an aggregated amyloid fibril state. The normal function of alpha-synuclein is poorly understood, and the precise mechanisms by which it leads to toxicity and cell death are also unclear. Although alpha-synuclein is a highly soluble, cytoplasmic protein, it binds to a variety of cellular membranes of different properties and compositions. These interactions are considered critical for at least some normal functions of alpha-synuclein, and may well play critical roles in both the aggregation of the protein and its mechanisms of toxicity. Here we review the known features of alpha-synuclein membrane interactions in the context of both the putative functions of the protein and of its pathological roles in disease.