Visualization of transient ultra-weak protein self-association in solution using paramagnetic relaxation enhancement

Preferential Interactions of a Crowder Protein with the Specific Binding Site of a Native Protein Complex

Nonspecific binding of crowder proteins with functional proteins is likely prevalent in vivo, yet direct quantitative evidence, let alone residue-specific information, is scarce. Here we present nuclear magnetic resonance (NMR) characterization showing that bovine serum albumin weakly but preferentially interacts with the histidine carrier protein (HPr). Notably, the binding interface overlaps with that for HPr's specific partner protein, EIN, leading to competition. The crowder protein thus decreases the EIN-HPr binding affinity and accelerates the dissociation of the native complex. In contrast, Ficoll-70 stabilizes the native complex and slows its dissociation, as one would expect from excluded-volume and microviscosity effects. Our atomistic modeling of macromolecular crowding rationalizes the experimental data and provides quantitative insights into the energetics of protein-crowder interactions. The integrated NMR and modeling study yields benchmarks for the effects of crowded cellular environments on protein-protein specific interactions, with implications for evolution regarding how nonspecific binding can be minimized or exploited.

Probing the Interaction of Huntingtin Exon-1 Polypeptides with the Chaperonin Nanomachine GroEL

Huntington's disease arises from polyQ expansion within the exon-1 region of huntingtin (htt^ex1 ), resulting in an aggregation-prone protein that accumulates in neuronal inclusion bodies. We investigate the interaction of various htt^ex1 constructs with the bacterial analog (GroEL) of the human chaperonin Hsp60. Using fluorescence spectroscopy and electron and atomic force microscopy, we show that GroEL inhibits fibril formation. The binding kinetics of htt^ex1 constructs with intact GroEL and a mini-chaperone comprising the apical domain is characterized by relaxation-based NMR measurements. The lifetimes of the complexes range from 100 to 400 μs with equilibrium dissociation constants (K_D ) of ∼1-2 mM. The binding interface is formed by the N-terminal amphiphilic region of htt^ex1 (which adopts a partially helical conformation) and the H and I helices of the GroEL apical domain. Sequestration of monomeric htt^ex1 by GroEL likely increases the critical concentration required for fibrillization.

Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer's Disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis

Protein misfolding and aggregation is observed in many amyloidogenic diseases affecting either the central nervous system or a variety of peripheral tissues. Structural and dynamic characterization of all species along the pathways from monomers to fibrils is challenging by experimental and computational means because they involve intrinsically disordered proteins in most diseases. Yet understanding how amyloid species become toxic is the challenge in developing a treatment for these diseases. Here we review what computer, in vitro, in vivo, and pharmacological experiments tell us about the accumulation and deposition of the oligomers of the (Aβ, tau), α-synuclein, IAPP, and superoxide dismutase 1 proteins, which have been the mainstream concept underlying Alzheimer's disease (AD), Parkinson's disease (PD), type II diabetes (T2D), and amyotrophic lateral sclerosis (ALS) research, respectively, for many years.

NMR-Based Methods for Protein Analysis

Nuclear magnetic resonance (NMR) spectroscopy is a well-established method for analyzing protein structure, interaction, and dynamics at atomic resolution and in various sample states including solution state, solid state, and membranous environment. Thanks to rapid NMR methodology development, the past decade has witnessed a growing number of protein NMR studies in complex systems ranging from membrane mimetics to living cells, which pushes the research frontier further toward physiological environments and offers unique insights in elucidating protein functional mechanisms. In particular, in-cell NMR has become a method of choice for bridging the huge gap between structural biology and cell biology. Herein, we review the recent developments and applications of NMR methods for protein analysis in close-to-physiological environments, with special emphasis on in-cell protein structural determination and the analysis of protein dynamics, both difficult to be accessed by traditional methods.

PRE-driven protein NMR structures: an alternative approach in highly paramagnetic systems

Metalloproteins play key roles across biology, and knowledge of their structure is essential to understand their physiological role. For those metalloproteins containing paramagnetic states, the enhanced relaxation caused by the unpaired electrons often makes signal detection unfeasible near the metal center, precluding adequate structural characterization right where it is more biochemically relevant. Here, we report a protein structure determination by NMR where two different sets of restraints, one containing Nuclear Overhauser Enhancements (NOEs) and another containing Paramagnetic Relaxation Enhancements (PREs), are used separately and eventually together. The protein PioC from Rhodopseudomonas palustris TIE-1 is a High Potential Iron-Sulfur Protein (HiPIP) where the [4Fe-4S] cluster is paramagnetic in both oxidation states at room temperature providing the source of PREs used as alternative distance restraints. Comparison of the family of structures obtained using NOEs only, PREs only, and the combination of both reveals that the pairwise root-mean-square deviation (RMSD) between them is similar and comparable with the precision within each family. This demonstrates that, under favorable conditions in terms of protein size and paramagnetic effects, PREs can efficiently complement and eventually replace NOEs for the structural characterization of small paramagnetic metalloproteins and de novo-designed metalloproteins by NMR. DATABASES: The 20 conformers with the lowest target function constituting the final family obtained using the full set of NMR restraints were deposited to the Protein Data Bank (PDB ID: 6XYV). The 20 conformers with the lowest target function obtained using NOEs only (PDB ID: 7A58) and PREs only (PDB ID: 7A4L) were also deposited to the Protein Data Bank. The chemical shift assignments were deposited to the BMRB (code 34487).

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

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

Model of a Kinetically Driven Crosstalk between Paralogous Protein Encounter Complexes

Proteins interact with one another across a broad spectrum of affinities. Our understanding of the low end of this spectrum, as characterized by millimolar dissociation constants, relies on a handful of cases in which weak encounters have experimentally been identified. These weak interactions away from the specific target binding site can lead toward a higher-affinity complex. Recently, we detected weak encounters between two paralogous phosphotransferase pathways of Escherichia coli, which regulate various metabolic processes and stress responses. In addition to encounters that are known to occur between cognate proteins, i.e., those that can exchange phosphate groups with each other, surprisingly, encounters involving noncognates were also observed. It is not clear whether these "futile" encounters have a cooperative or competitive role. Using agent-based simulations, we find that the encounter complexes can be cooperative or competitive so as to increase or lower the effective binding affinity of the specific complex under different circumstances. This finding invites further questions into how organisms might exploit such low affinities to connect their signaling components.

Weak Chemical Interactions That Drive Protein Evolution: Crowding, Sticking, and Quinary Structure in Folding and Function

In recent years, better instrumentation and greater computing power have enabled the imaging of elusive biomolecule dynamics in cells, driving many advances in understanding the chemical organization of biological systems. The focus of this Review is on interactions in the cell that affect both biomolecular stability and function and modulate them. The same protein or nucleic acid can behave differently depending on the time in the cell cycle, the location in a specific compartment, or the stresses acting on the cell. We describe in detail the crowding, sticking, and quinary structure in the cell and the current methods to quantify them both in vitro and in vivo. Finally, we discuss protein evolution in the cell in light of current biophysical evidence. We describe the factors that drive protein evolution and shape protein interaction networks. These interactions can significantly affect the free energy, ΔG, of marginally stable and low-population proteins and, due to epistasis, direct the evolutionary pathways in an organism. We finally conclude by providing an outlook on experiments to come and the possibility of collaborative evolutionary biology and biophysical efforts.

Measuring Ultra-Weak Protein Self-Association by Non-ideal Sedimentation Velocity

Ultra-weak self-association can govern the macroscopic solution behavior of concentrated macromolecular solutions ranging from food products to pharmaceutical formulations and the cytosol. For example, it can promote dynamic assembly of multi-protein signaling complexes, lead to intracellular liquid–liquid phase transitions, and seed crystallization or pathological aggregates. Unfortunately, weak self-association is technically extremely difficult to study, as it requires very high protein concentrations where short intermolecular distances cause strongly correlated particle motion. Additionally, protein samples near their solubility limit in vitro frequently show some degree of polydispersity. Here we exploit the strong mass-dependent separation of assemblies in the centrifugal field to study ultra-weak binding, using a sedimentation velocity technique that allows us to determine particle size distributions while accounting for colloidal hydrodynamic interactions and thermodynamic non-ideality (Chaturvedi, S. K.; et al. Nat. Commun.2018, 9, 4415; DOI: 10.1038/s41467-018-06902-x). We show that this approach, applied to self-associating proteins, can reveal a time-average association state for rapidly reversible self-associations from which the free energy of binding can be derived. The method is label-free and allows studying mid-sized proteins at millimolar protein concentrations in a wide range of solution conditions. We examine the performance of this method with hen egg lysozyme as a model system, reproducing its well-known ionic-strength-dependent weak self-association. The application to chicken γS-crystallin reveals weak monomer–dimer self-association with K_D = 24 mM, corresponding to a standard free energy change of approximately −9 kJ/mol, which is a large contribution to the delicate balance of forces ensuring eye lens transparency.

Self-adaptive multiscaling algorithm for efficient simulations of many-protein systems in crowded conditions

A method is described for the efficient simulation of multiprotein systems in crowded environments. It is based on an adaptive, reversible structural coarsening algorithm that preserves relevant physical features of the proteins across scales. Water is treated implicitly whereas all the other components of the aqueous solution, such as ions, cosolutes, or osmolytes, are treated in atomic detail. The focus is on the analytical adaptation of the solvent model to different levels of molecular resolutions, which allows continuous, on-the-fly transitions between scales. This permits the analytical calculation of forces during dynamics and preserves detailed balance in Monte Carlo simulations. A major computational speedup can be achieved in systems containing hundreds of proteins without cutting off the long-range interactions. The method can be combined with a self-adaptive configurational-bias sampling technique described previously, designed to detect strong, weak, or ultra-weak protein associations and shown to improve sampling efficiency and convergence. The implementation aims to simulate early stages of multimeric complexation, aggregation, or self-assembly. The method can be adopted as the basis for a more general algorithm to identify vertices, edges, and hubs in protein interaction networks or to predict critical steps in signal transduction pathways.

Interactions between the Intrinsically Disordered Proteins β-Synuclein and α-Synuclein

Several intrinsically disordered proteins have been implicated in the process of amyloid fibril formation in neurodegenerative disease, and developing approaches to inhibit the aggregation of these intrinsically disordered proteins is critical for establishing effective therapies against disease progression. The aggregation pathway of the intrinsically disordered protein alpha-synuclein, which is implicated in several neurodegenerative diseases known as synucleinopathies, has been extensively characterized. Less attention has been leveraged on beta-synuclein, a homologous intrinsically disordered protein that co-localizes with alpha-synuclein and is known to delay alpha-synuclein fibril formation. In this review, we focus on beta-synuclein and the molecular-level interactions between alpha-synuclein and beta-synuclein that underlie the delay of fibril formation. We highlight studies that begin to define alpha-synuclein and beta-synuclein interactions at the monomer, oligomer, and surface levels, and suggest that beta-synuclein plays a role in regulation of inhibition at many different stages of alpha-synuclein aggregation.

Investigating energy-based pool structure selection in the structure ensemble modeling with experimental distance constraints: The example from a multidomain protein Pub1

The structural variations of multidomain proteins with flexible parts mediate many biological processes, and a structure ensemble can be determined by selecting a weighted combination of representative structures from a simulated structure pool, producing the best fit to experimental constraints such as interatomic distance. In this study, a hybrid structure-based and physics-based atomistic force field with an efficient sampling strategy is adopted to simulate a model di-domain protein against experimental paramagnetic relaxation enhancement (PRE) data that correspond to distance constraints. The molecular dynamics simulations produce a wide range of conformations depicted on a protein energy landscape. Subsequently, a conformational ensemble recovered with low-energy structures and the minimum-size restraint is identified in good agreement with experimental PRE rates, and the result is also supported by chemical shift perturbations and small-angle X-ray scattering data. It is illustrated that the regularizations of energy and ensemble-size prevent an arbitrary interpretation of protein conformations. Moreover, energy is found to serve as a critical control to refine the structure pool and prevent data overfitting, because the absence of energy regularization exposes ensemble construction to the noise from high-energy structures and causes a more ambiguous representation of protein conformations. Finally, we perform structure-ensemble optimizations with a topology-based structure pool, to enhance the understanding on the ensemble results from different sources of pool candidates.

Application of Lysine-specific Labeling to Detect Transient Interactions Present During Human Lysozyme Amyloid Fibril Formation

Populating transient and partially unfolded species is a crucial step in the formation and accumulation of amyloid fibrils formed from pathogenic variants of human lysozyme linked with a rare but fatal hereditary systemic amyloidosis. The partially unfolded species possess an unstructured β-domain and C-helix with the rest of the α-domain remaining native-like. Here we use paramagnetic relaxation enhancement (PRE) measured by NMR spectroscopy to study the transient intermolecular interactions between such intermediate species. Nitroxide spin labels, introduced specifically at three individual lysine residues, generate distinct PRE profiles, indicating the presence of intermolecular interactions between residues within the unfolded β-domain. This study describes the applicability to PRE NMR measurements of selective lysine labeling, at different sites within a protein, as an alternative to the introduction of spin labels via engineered cysteine residues. These results reveal the importance of the β-sheet region of lysozyme for initiating self-assembly into amyloid fibrils.

Impact of spin label rigidity on extent and accuracy of distance information from PRE data

Paramagnetic relaxation enhancement (PRE) is a versatile tool for NMR spectroscopic structural and kinetic studies in biological macromolecules. Here, we compare the quality of PRE data derived from two spin labels with markedly different dynamic properties for large RNAs using the I-A riboswitch aptamer domain (78 nt) from Mesoplamsa florum as model system. We designed two I-A aptamer constructs that were spin-labeled by noncovalent hybridization of short spin-labeled oligomer fragments. As an example of a flexible spin label, ^UreidoU-TEMPO was incorporated into the 3' terminal end of helix P1 while, the recently developed rigid spin-label Çm was incorporated in the 5' terminal end of helix P1. We determined PRE rates obtained from aromatic ¹³C bound proton intensities and compared these rates to PREs derived from imino proton intensities in this sizeable RNA (~78 nt). PRE restraints derived from both imino and aromatic protons yielded similar data quality, and hence can both be reliably used for PRE determination. For NMR, the data quality derived from the rigid spin label Çm is slightly better than the data quality for the flexible ^UreidoTEMPO as judged by comparison of the structural agreement with the I-A aptamer crystal structure (3SKI).

Change in the Molecular Dimension of a RAGE-Ligand Complex Triggers RAGE Signaling

The weak oligomerization exhibited by many transmembrane receptors has a profound effect on signal transduction. The phenomenon is difficult to characterize structurally due to the large sizes of and transient interactions between monomers. The receptor for advanced glycation end products (RAGE), a signaling molecule central to the induction and perpetuation of inflammatory responses, is a weak constitutive oligomer. The RAGE domain interaction surfaces that mediate homo-dimerization were identified by combining segmental isotopic labeling of extracellular soluble RAGE (sRAGE) and nuclear magnetic resonance spectroscopy with chemical cross-linking and mass spectrometry. Molecular modeling suggests that two sRAGE monomers orient head to head forming an asymmetric dimer with the C termini directed toward the cell membrane. Ligand-induced association of RAGE homo-dimers on the cell surface increases the molecular dimension of the receptor, recruiting Diaphanous 1 (DIAPH1) and activating signaling pathways.

Computational study of the inhibitory mechanism of the kinase CDK5 hyperactivity by peptide p5 and derivation of a pharmacophore

The hyperactivity of the cyclic dependent kinase 5 (CDK5) induced by the activator protein p25 has been linked to a number of pathologies of the brain. The CDK5-p25 complex has thus emerged as a major therapeutic target for Alzheimer's disease (AD) and other neurodegenerative conditions. Experiments have shown that the peptide p5 reduces the CDK5-p25 activity without affecting the endogenous CDK5-p35 activity, whereas the peptide TFP5, obtained from p5, elicits similar inhibition, crosses the blood-brain barrier, and exhibits behavioral rescue of AD mice models with no toxic side effects. The molecular basis of the kinase inhibition is not currently known, and is here investigated by computer simulations. It is shown that p5 binds the kinase at the same CDK5/p25 and CDK5/p35 interfaces, and is thus a non-selective competitor of both activators, in agreement with available experimental data in vitro. Binding of p5 is enthalpically driven with an affinity estimated in the low µM range. A quantitative description of the binding site and pharmacophore is presented, and options are discussed to increase the binding affinity and selectivity in the design of drug-like compounds against AD.

Intermolecular Paramagnetic Relaxation Enhancement (PRE) Studies of Transient Complexes in Intrinsically Disordered Proteins

NMR interchain paramagnetic relaxation enhancement (PRE) techniques are a very powerful approach for detecting transient interchain interactions between intrinsically disordered proteins. These experiments, requiring a mixed sample containing a 1:1 ratio of isotope-labeled (15)N protein and natural abundance (14)N protein with a paramagnetic spin label, provide data that is limited to interchain interactions only. Application of these experiments to weakly associated transient species such as those that are present in the very early stages of self-assembly processes will aid our understanding of protein aggregation or fibril formation processes.

Weak self-association of cytochrome c peroxidase molecules observed by paramagnetic NMR

There is growing experimental evidence that many proteins exhibit a tendency for (ultra)weak homo- or hetero- oligomerization interactions. With the development of paramagnetic relaxation enhancement NMR spectroscopy it has become possible to characterize weak complexes experimentally and even detect complexes with affinities in the 1-25 mM range. We present evidence for a weak complex between cytochrome c peroxidase (CcP) molecules. In a previous study, we attached nitroxide based spin labels at three positions on CcP with the intent of observing intramolecular PRE effects. However, several intermolecular PRE effects were also observed suggesting a weak self-association between CcP molecules. The CcP-CcP complex was characterized using paramagnetic NMR and protein docking. The interaction occurs between the surface that is also part of the stereo-specific binding site for its physiological partner, cytochrome c (Cc), and several small, positively charged patches on the "back" of CcP. The CcP-CcP complex is not a stereo-specific complex. It is a dynamic ensemble of orientations, characteristic of an encounter state. The contact areas resemble those observed for CcP molecules in crystals. The CcP-CcP complex formation competes with that of the CcP-Cc complex. However, the affinity for Cc is much larger and thus it is expected that, under physiological conditions, auto-inhibition will be limited. A weak self-association between cytochrome c peroxidase molecules was characterized using paramagnetic NMR.

Advanced computational modeling for in vitro nanomaterial dosimetry

Background

Accurate and meaningful dose metrics are a basic requirement for in vitro screening to assess potential health risks of engineered nanomaterials (ENMs). Correctly and consistently quantifying what cells “see,” during an in vitro exposure requires standardized preparation of stable ENM suspensions, accurate characterizatoin of agglomerate sizes and effective densities, and predictive modeling of mass transport. Earlier transport models provided a marked improvement over administered concentration or total mass, but included assumptions that could produce sizable inaccuracies, most notably that all particles at the bottom of the well are adsorbed or taken up by cells, which would drive transport downward, resulting in overestimation of deposition.

Methods

Here we present development, validation and results of two robust computational transport models. Both three-dimensional computational fluid dynamics (CFD) and a newly-developed one-dimensional Distorted Grid (DG) model were used to estimate delivered dose metrics for industry-relevant metal oxide ENMs suspended in culture media. Both models allow simultaneous modeling of full size distributions for polydisperse ENM suspensions, and provide deposition metrics as well as concentration metrics over the extent of the well. The DG model also emulates the biokinetics at the particle-cell interface using a Langmuir isotherm, governed by a user-defined dissociation constant, K_D, and allows modeling of ENM dissolution over time.

Results

Dose metrics predicted by the two models were in remarkably close agreement. The DG model was also validated by quantitative analysis of flash-frozen, cryosectioned columns of ENM suspensions. Results of simulations based on agglomerate size distributions differed substantially from those obtained using mean sizes. The effect of cellular adsorption on delivered dose was negligible for K_D values consistent with non-specific binding (> 1 nM), whereas smaller values (≤ 1 nM) typical of specific high-affinity binding resulted in faster and eventual complete deposition of material.

Conclusions

The advanced models presented provide practical and robust tools for obtaining accurate dose metrics and concentration profiles across the well, for high-throughput screening of ENMs. The DG model allows rapid modeling that accommodates polydispersity, dissolution, and adsorption. Result of adsorption studies suggest that a reflective lower boundary condition is appropriate for modeling most in vitro ENM exposures.

Electronic supplementary material

The online version of this article (doi:10.1186/s12989-015-0109-1) contains supplementary material, which is available to authorized users.

Unveiling transient protein-protein interactions that modulate inhibition of alpha-synuclein aggregation by beta-synuclein, a pre-synaptic protein that co-localizes with alpha-synuclein

Pathology in Parkinson’s disease is linked to self-association of α-Synuclein (αS) into pathogenic oligomeric species and highly ordered amyloid fibrils. Developing effective therapeutic strategies against this debilitating disease is critical and βS, a pre-synaptic protein that co-localizes with αS, can act as an inhibitor of αS assembly. Despite the potential importance of βS as an inhibitor of αS, the nature, location and specificity of the molecular interactions between these two proteins is unknown. Here we use NMR paramagnetic relaxation enhancement experiments, to demonstrate that βS interacts directly with αS in a transient dimer complex with high specificity and weak affinity. Inhibition of αS by βS arises from transient αS/βS heterodimer species that exist primarily in head- to- tail configurations while αS aggregation arises from a more heterogeneous and weaker range of transient interactions that include both head-to-head and head-to-tail configurations. Our results highlight that intrinsically disordered proteins can interact directly with one another at low affinity and that the transient interactions that drive inhibition versus aggregation are distinct by virtue of their plasticity and specificity.

Detection and characterization of nonspecific, sparsely populated binding modes in the early stages of complexation

A method is proposed to study protein-ligand binding in a system governed by specific and nonspecific interactions. Strong associations lead to narrow distributions in the proteins configuration space; weak and ultraweak associations lead instead to broader distributions, a manifestation of nonspecific, sparsely populated binding modes with multiple interfaces. The method is based on the notion that a discrete set of preferential first-encounter modes are metastable states from which stable (prerelaxation) complexes at equilibrium evolve. The method can be used to explore alternative pathways of complexation with statistical significance and can be integrated into a general algorithm to study protein interaction networks. The method is applied to a peptide-protein complex. The peptide adopts several low-population conformers and binds in a variety of modes with a broad range of affinities. The system is thus well suited to analyze general features of binding, including conformational selection, multiplicity of binding modes, and nonspecific interactions, and to illustrate how the method can be applied to study these problems systematically. The equilibrium distributions can be used to generate biasing functions for simulations of multiprotein systems from which bulk thermodynamic quantities can be calculated.

The low-affinity complex of cytochrome c and its peroxidase

The complex of yeast cytochrome c peroxidase and cytochrome c is a paradigm of the biological electron transfer (ET). Building on seven decades of research, two different models have been proposed to explain its functional redox activity. One postulates that the intermolecular ET occurs only in the dominant, high-affinity protein–protein orientation, while the other posits formation of an additional, low-affinity complex, which is much more active than the dominant one. Unlike the high-affinity interaction—extensively studied by X-ray crystallography and NMR spectroscopy—until now the binding of cytochrome c to the low-affinity site has not been observed directly, but inferred mainly from kinetics experiments. Here we report the structure of this elusive, weak protein complex and show that it consists of a dominant, inactive bound species and an ensemble of minor, ET-competent protein–protein orientations, which summarily account for the experimentally determined value of the ET rate constant.

The redox activity of cytochrome c in complex with its peroxidase has been rationalized by two possible models; a readily observable high-affinity complex and a more elusive but potentially more active low-affinity complex. Here, the authors provide an NMR-based structural mapping of this low-affinity complex.

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.

Specific and nonspecific interactions in ultraweak protein-protein associations revealed by solvent paramagnetic relaxation enhancements

Weak and transient protein–protein interactions underlie numerous biological processes. However, the location of the interaction sites of the specific complexes and the effect of transient, nonspecific protein–protein interactions often remain elusive. We have investigated the weak self-association of human growth hormone (hGH, K_D = 0.90 ± 0.03 mM) at neutral pH by the paramagnetic relaxation enhancement (PRE) of the amide protons induced by the soluble paramagnetic relaxation agent, gadodiamide (Gd(DTPA-BMA)). Primarily, it was found that the PREs are in agreement with the general Hwang-Freed model for relaxation by translational diffusion (J. Chem. Phys.1975, 63, 4017–4025), only if crowding effects on the diffusion in the protein solution are taken into account. Second, by measuring the PREs of the amide protons at increasing hGH concentrations and a constant concentration of the relaxation agent, it is shown that a distinction can be made between residues that are affected only by transient, nonspecific protein–protein interactions and residues that are involved in specific protein–protein associations. Thus, the PREs of the former residues increase linearly with the hGH concentration in the entire concentration range because of a reduction of the diffusion caused by the transient, nonspecific protein–protein interactions, while the PREs of the latter residues increase only at the lower hGH concentrations but decrease at the higher concentrations because of specific protein–protein associations that impede the access of gadodiamide to the residues of the interaction surface. Finally, it is found that the ultraweak aggregation of hGH involves several interaction sites that are located in patches covering a large part of the protein surface.

Visualization of transient protein-protein interactions that promote or inhibit amyloid assembly

In the early stages of amyloid formation, heterogeneous populations of oligomeric species are generated, the affinity, specificity, and nature of which may promote, inhibit, or define the course of assembly. Despite the importance of the intermolecular interactions that initiate amyloid assembly, our understanding of these events remains poor. Here, using amyloidogenic and nonamyloidogenic variants of β₂-microglobulin, we identify the interactions that inhibit or promote fibril formation in atomic detail. The results reveal that different outcomes of assembly result from biomolecular interactions involving similar surfaces. Specifically, inhibition occurs via rigid body docking of monomers in a head-to-head orientation to form kinetically trapped dimers. By contrast, the promotion of fibrillation involves relatively weak protein association in a similar orientation, which results in conformational changes in the initially nonfibrillogenic partner. The results highlight the complexity of interactions early in amyloid assembly and reveal atomic-level information about species barriers in amyloid formation.

•

Dissection of protein-protein interactions in the early stages of amyloid assembly

•

Rare biomolecular collisions and the course of amyloid assembly

•

Interaction surfaces and different outcomes of amyloid assembly

•

Molecular description of species barriers in amyloid assembly

Noncovalent spin labeling of riboswitch RNAs to obtain long-range structural NMR restraints

Paramagnetic relaxation enhancement (PRE) NMR is a powerful method to study structure, dynamics and function of proteins. Up to now, the application of PRE NMR on RNAs is a significant challenge due to the limited size of chemically synthesized RNA. Here, we present a noncovalent spin labeling strategy to spin label RNAs in high yields required for NMR studies. The approach requires the presence of a helix segment composed of about 10 nucleotides (nt) but is not restricted by the size of the RNA. We show successful application of this strategy on the 2'dG sensing aptamer domain of Mesoplasma florum (78 nt). The aptamer domain was prepared in two fragments. A larger fragment was obtained by biochemical means, while a short spin labeled fragment was prepared by chemical solid-phase synthesis. The two fragments were annealed noncovalently by hybridization. We performed NMR, cw-EPR experiments and gel shift assays to investigate the stability of the two-fragment complex. NMR analysis in (15)N-TROSY and (1)H,(1)H-NOESY spectra of both unmodified and spin labeled aptamer domain show that the fragmented system forms a stable hybridization product, is in structural agreement with the full length aptamer domain and maintains its function. Together with structure modeling, experimentally determined (1)H-Γ2 rates are in agreement with reported crystal structure data and show that distance restraints up to 25 Å can be obtained from NMR PRE data of RNA.

Encounter complexes and dimensionality reduction in protein-protein association

An outstanding challenge has been to understand the mechanism whereby proteins associate. We report here the results of exhaustively sampling the conformational space in protein–protein association using a physics-based energy function. The agreement between experimental intermolecular paramagnetic relaxation enhancement (PRE) data and the PRE profiles calculated from the docked structures shows that the method captures both specific and non-specific encounter complexes. To explore the energy landscape in the vicinity of the native structure, the nonlinear manifold describing the relative orientation of two solid bodies is projected onto a Euclidean space in which the shape of low energy regions is studied by principal component analysis. Results show that the energy surface is canyon-like, with a smooth funnel within a two dimensional subspace capturing over 75% of the total motion. Thus, proteins tend to associate along preferred pathways, similar to sliding of a protein along DNA in the process of protein-DNA recognition.

DOI:http://dx.doi.org/10.7554/eLife.01370.001

Proteins rarely act alone. Instead, they tend to bind to other proteins to form structures known as complexes. When two proteins come together to form a complex, they twist and turn through a series of intermediate states before they form the actual complex. These intermediate states are difficult to study because they don’t last for very long, which means that our knowledge of how complexes are formed remains incomplete.

One promising approach for studying the formation of complexes is called paramagnetic relaxation enhancement. In this technique certain areas in one of the proteins are labelled with magnetic particles, which produce signals when the two proteins are close to each other. Repeating the measurement several times with the magnetic particles in different positions provides information about the overall structure of the complex. Computational modelling can then be used to work out the fine details of the structure, including the shapes of the intermediate structures made by the proteins as they interact.

A computer method called docking can be used to predict the most favourable positions that the proteins can take, relative to one another, in a complex. This involves calculating the energy contained in the system, with the correct structure having the lowest energy. Docking methods also predict protein models with slightly higher energies, but with structures that are radically different. Modellers usually ignore these structures, but comparing the docking results to paramagnetic relaxation enhancement data, Kozakov et al. found that these structures actually represent the intermediate states.

Analysing the structure of the intermediate states revealed that the movement of the two proteins relative to one another is severely restricted as they form the final complex. Kozakov et al. found that proteins associate along preferred pathways, similar to the way a protein slides along DNA in the process of protein-DNA recognition. Knowing that the movement of the proteins is restricted in this way will enable researchers to improve the efficiency of docking calculations.

DOI:http://dx.doi.org/10.7554/eLife.01370.002

Specific and non-specific protein association in solution: computation of solvent effects and prediction of first-encounter modes for efficient configurational bias Monte Carlo simulations

Weak and ultraweak protein-protein association play a role in molecular recognition and can drive spontaneous self-assembly and aggregation. Such interactions are difficult to detect experimentally, and are a challenge to the force field and sampling technique. A method is proposed to identify low-population protein-protein binding modes in aqueous solution. The method is designed to identify preferential first-encounter complexes from which the final complex(es) at equilibrium evolve. A continuum model is used to represent the effects of the solvent, which accounts for short- and long-range effects of water exclusion and for liquid-structure forces at protein/liquid interfaces. These effects control the behavior of proteins in close proximity and are optimized on the basis of binding enthalpy data and simulations. An algorithm is described to construct a biasing function for self-adaptive configurational-bias Monte Carlo of a set of interacting proteins. The function allows mixing large and local changes in the spatial distribution of proteins, thereby enhancing sampling of relevant microstates. The method is applied to three binary systems. Generalization to multiprotein complexes is discussed.

NMR-based conformational ensembles explain pH-gated opening and closing of OmpG channel

The outer membrane protein G (OmpG) is a monomeric 33 kDa 14-stranded β-barrel membrane protein functioning as a nonspecific porin for the uptake of oligosaccharides in Escherichia coli. Two different crystal structures of OmpG obtained at different values of pH suggest a pH-gated pore opening mechanism. In these structures, extracellular loop 6 extends away from the barrel wall at neutral pH but is folded back into the pore lumen at low pH, blocking transport through the pore. Loop 6 was invisible in a previously published solution NMR structure of OmpG in n-dodecylphosphocholine micelles, presumably due to conformational exchange on an intermediate NMR time scale. Here we present an NMR paramagnetic relaxation enhancement (PRE)-based approach to visualize the conformational dynamics of loop 6 and to calculate conformational ensembles that explain the pH-gated opening and closing of the OmpG channel. The different loop conformers detected by the PRE ensemble calculations were validated by disulfide cross-linking of strategically engineered cysteines and electrophysiological single channel recordings. The results indicate a more dynamically regulated channel opening and closing than previously thought and reveal additional membrane-associated conformational ensembles at pH 6.3 and 7.0. We anticipate this approach to be generally applicable to detect and characterize functionally important conformational ensembles of membrane proteins.

Measurement of rate constants for homodimer subunit exchange using double electron-electron resonance and paramagnetic relaxation enhancements

Here, we report novel methods to measure rate constants for homodimer subunit exchange using double electron-electron resonance (DEER) electron paramagnetic resonance spectroscopy measurements and nuclear magnetic resonance spectroscopy based paramagnetic relaxation enhancement (PRE) measurements. The techniques were demonstrated using the homodimeric protein Dsy0195 from the strictly anaerobic bacterium Desulfitobacterium hafniense Y51. At specific times following mixing site-specific MTSL-labeled Dsy0195 with uniformly (15)N-labeled Dsy0195, the extent of exchange was determined either by monitoring the decrease of MTSL-labeled homodimer from the decay of the DEER modulation depth or by quantifying the increase of MTSL-labeled/(15)N-labeled heterodimer using PREs. Repeated measurements at several time points following mixing enabled determination of the homodimer subunit dissociation rate constant, k (-1), which was 0.037 ± 0.005 min(-1) derived from DEER experiments with a corresponding half-life time of 18.7 min. These numbers agreed with independent measurements obtained from PRE experiments. These methods can be broadly applied to protein-protein and protein-DNA complex studies.

Jahn-Teller dynamics in a series of high-symmetry Co(II) chelates determine paramagnetic relaxation enhancements

NMR paramagnetic relaxation enhancements (PREs) of a series of structurally characterized, trigonal bis-trispyrazolylborate (Tp) chelates of high-spin Co(II), spanning 100-850 MHz in field, are reported. Prior knowledge of the metal-nucleus distances allows numerical extraction of position-dependent electron spin relaxation rates (τ(c)(-1)) from direct measurement of the individual PREs of the four symmetry distinct protons in Co(Tp)(2), using available closed-form expressions. The data for this electronically complex system where spin-orbit coupling defines the ground state electronic structure are analyzed in terms of the Solomon-Bloembergen-Morgan (SBM) relations, as well as available zero-field splitting limit theories. A simple angular correction is shown to be sufficient to reconcile the individual τ(c)(T) data for the four classes of protons. The data identify a previously unrecognized dynamic Jahn-Teller effect in these historically important complexes, with a barrier of ~230 cm(-1), pointing to a level of dynamics in trispyrazolylborate chemistry that has not been described before, and further show that it is the Jahn-Teller that is responsible for the PREs in fluid solution. A field-dependent component is also identified for the two protons nearest g(//), which is suggested to arise due to Zeeman mixing of excited state character into the ground level.

Spatial distribution of cytoplasmic domains of the Mg(2+)-transporter MgtE, in a solution lacking Mg(2+), revealed by paramagnetic relaxation enhancement

MgtE is a prokaryotic Mg(2+) transporter that controls cellular Mg(2+) concentrations. We previously reported crystal structures of the cytoplasmic region of MgtE, consisting of 2 domains, that is, N and CBS, in the Mg(2+)-free and Mg(2+)-bound forms. The Mg(2+)-binding sites lay at the interface of the 2 domains, making the Mg(2+)-bound form compact and globular. In the Mg(2+)-free structure, however, the domains are far apart, and the Mg(2+)-binding sites are destroyed. Therefore, it is unclear how Mg(2+)-free MgtE changes its conformation to accommodate Mg(2+) ions. Here, we used paramagnetic relaxation enhancement (PRE) to characterize the relative orientation of the N and CBS domains in the absence of Mg(2+) in solution. When the residues on the surface of the CBS domain were labeled with nitroxide tags, significant PRE effects were observed for the residues in the N domain. No single structure satisfied the PRE profiles, suggesting that the N and CBS domains are not fixed in a particular orientation in solution. We then conducted ensemble simulated annealing calculations in order to obtain the atomic probability density and visualize the spatial distribution of the N domain in solution. The results indicate that the N domain tends to occupy the space near its position in the Mg(2+)-bound crystal structure, facilitating efficient capture of Mg(2+) with increased intracellular Mg(2+) concentration, which is necessary to close the gate.

Fast mapping of biomolecular interfaces by Random Spin Labeling (RSL)

Random spin labeling (RSL) is a method for rapid mapping of biomolecular interaction surfaces using an interaction partner with SL and an interaction partner enriched in (13)C or (15)N nuclei for paramagnetic relaxation enhanced NMR-based detection. The SL reaction is conducted in a manner resulting in a heterogeneous reaction product consisting of different populations of the protein carrying a varying number of spin labels at different positions. Preparation of the paramagnetic probe is complete within a few hours and hence much faster than site selective SL. RSL is applicable to tightly interacting systems but shows its particular strength when applied to systems involving weak or transient contacts.

Noncovalent dimerization of ubiquitin

Another kind of dynamics: Ubiquitin noncovalently dimerizes with a dissociation constant of approximately 5 mM. The two subunits adopt an array of relative orientations, utilizing an interface also for binding to other proteins (see picture). Quaternary fluctuation among members of the dimer ensemble constitutes a different kind of dynamics that complements the tertiary dynamics of each ubiquitin subunit.

Transient weak protein-protein complexes transfer heme across the cell wall of Staphylococcus aureus

Iron is an essential nutrient for the bacterial pathogen Staphylococcus aureus . Heme in hemoglobin (Hb) is the most abundant source of iron in the human body and during infections is captured by S. aureus using iron-regulated surface determinant (Isd) proteins. A central step in this process is the transfer of heme between the cell wall associated IsdA and IsdC hemoproteins. Biochemical evidence indicates that heme is transferred via an activated IsdA:heme:IsdC heme complex. Transfer is rapid and occurs up to 70,000 times faster than indirect mechanisms in which heme is released into the solvent. To gain insight into the mechanism of transfer, we modeled the structure of the complex using NMR paramagnetic relaxation enhancement (PRE) methods. Our results indicate that IsdA and IsdC transfer heme via an ultraweak affinity "handclasp" complex that juxtaposes their respective 3(10) helices and β7/β8 loops. Interestingly, PRE also identified a set of transient complexes that could represent high-energy pre-equilibrium encounter species that form prior to the stereospecific handclasp complex. Targeted amino acid mutagenesis and stopped-flow measurements substantiate the functional relevance of a PRE-derived model, as mutation of interfacial side chains significantly slows the rate of transfer. IsdA and IsdC bind heme using NEAr Transporter (NEAT) domains that are conserved in many species of pathogenic Gram-positive bacteria. Heme transfer in these microbes may also occur through structurally similar transient stereospecific complexes.

Exploring sparsely populated states of macromolecules by diamagnetic and paramagnetic NMR relaxation

Sparsely populated states of macromolecules, characterized by short lifetimes and high free-energies relative to the predominant ground state, often play a key role in many biological, chemical, and biophysical processes. In this review, we briefly summarize various new developments in NMR spectroscopy that permit these heretofore invisible, sparsely populated states to be detected, characterized, and in some instances visualized. Relaxation dispersion spectroscopy yields detailed kinetic information on processes involving species characterized by distinct chemical shifts with lifetimes in the ∼50 μs-10 ms range and populations as low as 0.5%. In the fast exchange regime (time scale less than ∼250-500 μs), the footprint of sparsely populated states can be observed on paramagnetic relaxation enhancement profiles measured on the resonances of the major species, thereby yielding structural information that is directly related to paramagnetic center-nuclei distances from which it is possible, under suitable circumstances, to compute a structure or ensemble of structures for the minor species. Finally, differential transverse relaxation measurements can be used to detect lifetime broadening effects that directly reflect the unidirectional rates for the conversion of NMR-visible into high-molecular weight NMR-invisible species. Examples of these various approaches are presented.

Characterization of the interaction between the Salmonella type III secretion system tip protein SipD and the needle protein PrgI by paramagnetic relaxation enhancement

Many Gram-negative bacteria that cause major diseases and mortality worldwide require the type III secretion system (T3SS) to inject virulence proteins into their hosts and cause infections. A structural component of the T3SS is the needle apparatus, which consists of a base, an external needle, and a tip complex. In Salmonella typhimurium, the external needle is assembled by the polymerization of the needle protein PrgI. On top of this needle sits a tip complex, which is partly formed by the tip protein SipD. How SipD interacts with PrgI during the assembly of the T3SS needle apparatus remains unknown. The central region of PrgI forms an α-helical hairpin, whereas SipD has a long central coiled-coil, which is a defining structural feature of other T3SS tip proteins as well. Using NMR paramagnetic relaxation enhancement, we have identified a specific region on the SipD coiled-coil that interacts directly with PrgI. We present a model of how SipD might dock at the tip of the needle based on our paramagnetic relaxation enhancement results, thus offering new insight about the mechanism of assembly of the T3SS needle apparatus.

Mapping the encounter state of a transient protein complex by PRE NMR spectroscopy

Many biomolecular interactions proceed via a short-lived encounter state, consisting of multiple, lowly-populated species invisible to most experimental techniques. Recent development of paramagnetic relaxation enhancement (PRE) nuclear magnetic resonance (NMR) spectroscopy has allowed to directly visualize such transient intermediates in a number of protein-protein and protein-DNA complexes. Here we present an analysis of the recently published PRE NMR data for a protein complex of yeast cytochrome c (Cc) and cytochrome c peroxidase (CcP). First, we describe a simple, general method to map out the spatial and temporal distributions of binding geometries constituting the Cc-CcP encounter state. We show that the spatiotemporal mapping provides a reliable estimate of the experimental coverage and, at higher coverage levels, allows to delineate the conformational space sampled by the minor species. To further refine the encounter state, we performed PRE-based ensemble simulations. The generated solutions reproduce well the experimental data and lie within the allowed regions of the encounter maps, confirming the validity of the mapping approach. The refined encounter ensembles are distributed predominantly in a region encompassing the dominant form of the complex, providing experimental proof for the results of classical theoretical simulations.