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

Efficient determination of the accessible conformation space of multi-domain complexes based on EPR PELDOR data

Many proteins can adopt multiple conformations which are important for their function. This is also true for proteins and domains that are covalently linked to each other. One important example is ubiquitin, which can form chains of different conformations depending on which of its lysine side chains is used to form an isopeptide bond with the C-terminus of another ubiquitin molecule. Similarly, ubiquitin gets covalently attached to active-site residues of E2 ubiquitin-conjugating enzymes. Due to weak interactions between ubiquitin and its interaction partners, these covalent complexes adopt multiple conformations. Understanding the function of these complexes requires the characterization of the entire accessible conformation space and its modulation by interaction partners. Long-range (1.8-10 nm) distance restraints obtained by EPR spectroscopy in the form of probability distributions are ideally suited for this task as not only the mean distance but also information about the conformation dynamics is encoded in the experimental data. Here we describe a computational method that we have developed based on well-established structure determination software using NMR restraints to calculate the accessible conformation space using PELDOR/DEER data.

Paramagnetic NMR restraints for the characterization of protein structural rearrangements

Mobility is a common feature of biomacromolecules, often fundamental for their function. Thus, in many cases, biomacromolecules cannot be described by a single conformation, but rather by a conformational ensemble. NMR paramagnetic data demonstrated quite informative to monitor this conformational variability, especially when used in conjunction with data from different sources. Due to their long-range nature, paramagnetic data can, for instance, i) clearly demonstrate the occurrence of conformational rearrangements, ii) reveal the presence of minor conformational states, sampled only for a short time, iii) indicate the most representative conformations within the conformational ensemble sampled by the molecule, iv) provide an upper limit to the weight of each conformation.

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.

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.

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.

Mechanism of activation and the rewired network: New drug design concepts

Precision oncology benefits from effective early phase drug discovery decisions. Recently, drugging inactive protein conformations has shown impressive successes, raising the cardinal questions of which targets can profit and what are the principles of the active/inactive protein pharmacology. Cancer driver mutations have been established to mimic the protein activation mechanism. We suggest that the decision whether to target an inactive (or active) conformation should largely rest on the protein mechanism of activation. We next discuss the recent identification of double (multiple) same-allele driver mutations and their impact on cell proliferation and suggest that like single driver mutations, double drivers also mimic the mechanism of activation. We further suggest that the structural perturbations of double (multiple) in cis mutations may reveal new surfaces/pockets for drug design. Finally, we underscore the preeminent role of the cellular network which is deregulated in cancer. Our structure-based review and outlook updates the traditional Mechanism of Action, informs decisions, and calls attention to the intrinsic activation mechanism of the target protein and the rewired tumor-specific network, ushering innovative considerations in precision medicine.

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.

Conformational Space Sampled by Domain Reorientation of Linear Diubiquitin Reflected in Its Binding Mode for Target Proteins

Linear polyubiquitin chains regulate diverse signaling proteins, in which the chains adopt various conformations to recognize different target proteins. Thus, the structural plasticity of the chains plays an important role in controlling the binding events. Herein, paramagnetic NMR spectroscopy is employed to explore the conformational space sampled by linear diubiquitin, a minimal unit of linear polyubiquitin, in its free state. Rigorous analysis of the data suggests that, regarding the relative positions of the ubiquitin units, particular regions of conformational space are preferentially sampled by the molecule. By combining these results with further data collected for charge-reversal derivatives of linear diubiquitin, structural insights into the factors underlying the binding events of linear diubiquitin are obtained.

Characterization of Weak Protein Domain Structure by Spin-Label Distance Distributions

Function of intrinsically disordered proteins may depend on deviation of their conformational ensemble from that of a random coil. Such deviation may be hard to characterize and quantify, if it is weak. We explored the potential of distance distributions between spin labels, as they can be measured by electron paramagnetic resonance techniques, for aiding such characterization. On the example of the intrinsically disordered N-terminal domain 1-267 of fused in sarcoma (FUS) we examined what such distance distributions can and cannot reveal on the random-coil reference state. On the example of the glycine-rich domain 188-320 of heterogeneous nuclear ribonucleoprotein A1 (hnRNP A1) we studied whether deviation from a random-coil ensemble can be robustly detected with 19 distance distribution restraints. We discuss limitations imposed by ill-posedness of the conversion of primary data to distance distributions and propose overlap of distance distributions as a fit criterion that can tackle this problem. For testing consistency and size sufficiency of the restraint set, we propose jack-knife resampling. At current desktop computers, our approach is expected to be viable for domains up to 150 residues and for between 10 and 50 distance distribution restraints.

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.

Magnetic susceptibility and paramagnetism-based NMR

The magnetic interactions between the nuclear magnetic moment and the magnetic moment of unpaired electron(s) depend on the structure and dynamics of the molecules where the paramagnetic center is located and of their partners. The long-range nature of the magnetic interactions is thus a reporter of invaluable information for structural biology studies, when other techniques often do not provide enough data for the atomic-level characterization of the system. This precious information explains the flourishing of paramagnetism-assisted NMR studies in recent years. Many paramagnetic effects are related to the magnetic susceptibility of the paramagnetic metal. Although these effects have been known for more than half a century, different theoretical models and new approaches have been proposed in the last decade. In this review, we have summarized the consequences for NMR spectroscopy of magnetic interactions between nuclear and electron magnetic moments, and thus of the presence of a magnetic susceptibility due to metals, and we do so using a unified notation.

Tracking Conformational Changes in Calmodulin in vitro, in Cell Extract, and in Cells by Electron Paramagnetic Resonance Distance Measurements

It is an open question whether the conformations of proteins sampled in dilute solutions are the same as in the cellular environment. Here we address this question by double electron-electron resonance (DEER) distance measurements with Gd(III) spin labels to probe the conformations of calmodulin (CaM) in vitro, in cell extract, and in human HeLa cells. Using the CaM mutants N53C/T110C and T34C/T117C labeled with maleimide-DOTA-Gd(III) in the N- and C-terminal domains, we observed broad and varied interdomain distance distributions. The in vitro distance distributions of apo-CaM and holo-CaM in the presence and absence of the IQ target peptide can be described by combinations of closed, open, and collapsed conformations. In cell extract, apo- and holo-CaM bind to target proteins in a similar way as apo- and holo-CaM bind to IQ peptide in vitro. In HeLa cells, however, in the presence or absence of elevated in-cell Ca²⁺ levels CaM unexpectedly produced more open conformations and very broad distance distributions indicative of many different interactions with in-cell components. These results show-case the importance of in-cell analyses of protein structures.

Site-Specific Tagging of Proteins with Paramagnetic Ions for Determination of Protein Structures in Solution and in Cells

High-resolution NMR spectroscopy is sensitive to local structural variations and subtle dynamics of biomolecules and is an important technique for studying the structures, dynamics, and interactions of these molecules. Small-molecule probes, including paramagnetic tags, have been developed for this purpose. Paramagnetic effects manifested in magnetic resonance spectra have long been recognized as valuable tools for chemical analysis of small molecules, and these effects were later applied in the fields of chemical biology and structural biology. However, such applications require the installation of a paramagnetic center in the biomolecules of interest. Paramagnetic metal ions and stable free radicals are the most widely used paramagnetic probes for biological magnetic resonance spectroscopy, and therefore mild, high-yielding approaches for chemically attaching paramagnetic tags to biomolecules are in high demand. In this Account, we begin by discussing paramagnetic species, especially transition metal ions and lanthanide ions, that are suitable for NMR and EPR studies, particularly for in-cell applications. Thereafter, we describe approaches for site-specific tagging of proteins with paramagnetic ions and discuss considerations involved in designing high-quality paramagnetic tags, including the strength of the binding between the metal-chelating moiety and the paramagnetic ion, the chemical stability, and the flexibility of the tether between the paramagnetic tag and the target protein. The flexibility of a tag correlates strongly with the averaging of paramagnetic effects observed in NMR spectra, and we describe methods for increasing tag rigidity and applications of such tags in biological systems. We also describe specific applications of established site-specific tagging approaches and newly developed paramagnetic tags for the elucidation of protein structures and dynamics at atomic resolution both in solution and in cells. First, we describe the determination of the 3D structure of a short-lived, low-abundance enzyme intermediate complex in real time by using pseudocontact shifts as structural restraints. Second, we demonstrate the utility of stable paramagnetic tags for determining 3D structures of proteins in live cells, and pseudocontact shifts are shown to be valuable structural restraints for in-cell protein analysis. Third, we show that a NMR optimized paramagnetic tag allows one to determine distance restraints on proteins by double electron-electron resonance (DEER) measurements with high spatial resolution both in vitro and in cells. Finally, we summarize recent advances in site-specific tagging of proteins to achieve atomic-resolution information about structural changes of proteins, and the advantages and challenges of magnetic resonance spectroscopy in biological systems.

Hypothesis: entaticversusecstatic states in metalloproteins

Biological metal sites are optimized for function by conformational properties of the protein macroligand.