Conformational Ensemble of Disordered Proteins Probed by Solvent Paramagnetic Relaxation Enhancement (sPRE)

Characterization of the conformational ensemble of disordered proteins is highly important for understanding protein folding and aggregation mechanisms, but remains a computational and experimental challenge owing to the dynamic nature of these proteins. New observables that can provide unique insights into transient residual structures in disordered proteins are needed. Here using denatured ubiquitin as a model system, NMR solvent paramagnetic relaxation enhancement (sPRE) measurements provide an accurate and highly sensitive probe for detecting low populations of residual structure in a disordered protein. Furthermore, a new ensemble calculation approach based on sPRE restraints in conjunction with residual dipolar couplings (RDCs) and small‐angle X‐ray scattering (SAXS) is used to define the conformational ensemble of disordered proteins at atomic resolution. The approach presented should be applicable to a wide range of dynamic macromolecules.     

Solid-state NMR on a large multidomain integral membrane protein: the outer membrane protein assembly factor BamA

Multidomain proteins constitute a large part of prokaryotic and eukaryotic proteomes and play fundamental roles in various physiological processes. However, their structural characterization is challenging because of their large size and intrinsic flexibility. We show here that motional-filtered high-resolution solid-state NMR (ssNMR) experiments allow for the observation and structural analysis of very large multidomain membrane proteins that are characterized by different motional time scales. This approach was used to probe the folding of the 790-residue membrane protein BamA, which is the core component of the Escherichia coli outer membrane protein assembly machinery. A combination of dipolar- and scalar-based two-dimensional ssNMR experiments applied to two uniformly (13)C,(15)N-labeled BamA variants revealed characteristic secondary structure elements and distinct dynamics within the BamA transmembrane protein segment and the periplasmic POTRA domains. This approach hence provides a general strategy for collecting atomic-scale structural information on multidomain (membrane) proteins in a native-like environment.     

Residual dipolar coupling measurements of transmembrane proteins using aligned low-q bicelles and high-resolution magic angle spinning NMR spectroscopy

Bicelles are a major medium form to produce weak alignment of soluble proteins for residual dipolar coupling (RDC) measurements. The obstacle to using the same type of bicelles for transmembrane proteins with solution-state NMR spectroscopy is the loss of signals due to the adhesion or penetration of the proteins into large bicelles, resulting in slow protein tumbling. In this study, weak alignment of the second and third transmembrane domains (TM23) of the human glycine receptor (GlyR) was achieved in low-q bicelles (q = DMPC/DHPC). Although protein-free bicelles with such low q would likely show isotropic properties, the insertion of TM23 induced weakly preferred orientations so that the RDC of the embedded protein can be measured. The extent of the alignment increased but the TM23 signal intensity decreased when q was varied from 0.19 to 0.60. A q of 0.50 was found to be an optimal compromise between alignment and the signal-to-noise ratio. In each pair of NMR experiments for RDC measurements, the same sample and pulse sequence were used, with one being performed at high-resolution magic-angle spinning to obtain pure J-couplings without RDC. A meaningful structure refinement in bicelles was possible by iteratively fitting the experimental RDCs to the back-calculated RDCs using the high-resolution NMR structure of GlyR TM23 in trifluoroethanol as the starting template. Combination of this method with the conventional high-resolution NMR in membrane mimicking mixtures of water and organic solvents offers an attractive way to derive structural information for membrane proteins in their native environment.     

Site-directed parallel spin-labeling and paramagnetic relaxation enhancement in structure determination of membrane proteins by solution NMR spectroscopy

A major challenge for the structure determination of integral membrane proteins by solution NMR spectroscopy is the limited number of NOE restraints in these systems stemming from extensive deuteration. Paramagnetic relaxation enhancement (PRE) by means of nitroxide spin-labels can provide valuable long-range distance information but, in practice, has limits in its application to membrane proteins because spin-labels are often incompletely reduced in highly apolar environments. Using the integral membrane protein OmpA as a model system, we introduce a method of parallel spin-labeling with paramagnetic and diamagnetic labels and show that distances in the range 15-24 Angstroms can be readily determined. The protein was labeled at 11 water-exposed and lipid-covered sites, and 320 PRE distance restraints were measured. The addition of these restraints resulted in significant improvement of the calculated backbone structure of OmpA. Structures of reasonable quality can even be calculated with PRE distance restraints only, i.e., in the absence of NOE distance restraints.     

Bioconjugation of Proteins with a Paramagnetic NMR and Fluorescent Tag

Site‐specific labeling of proteins with lanthanide ions offers great opportunities for investigating the structure, function, and dynamics of proteins by virtue of the unique properties of lanthanides. Lanthanide‐tagged proteins can be studied by NMR, X‐ray, fluorescence, and EPR spectroscopy. However, the rigidity of a lanthanide tag in labeling of proteins plays a key role in the determination of protein structures and interactions. Pseudocontact shift (PCS) and paramagnetic relaxation enhancement (PRE) are valuable long‐range structure restraints in structural‐biology NMR spectroscopy. Generation of these paramagnetic restraints generally relies on site‐specific tagging of the target proteins with paramagnetic species. To avoid nonspecific interaction between the target protein and paramagnetic tag and achieve reliable paramagnetic effects, the rigidity, stability, and size of lanthanide tag is highly important in paramagnetic labeling of proteins. Here 4′‐mercapto‐2,2′: 6′,2′′‐terpyridine‐6,6′′‐dicarboxylic acid (4MTDA) is introduced as a a rigid paramagnetic and fluorescent tag which can be site‐specifically attached to a protein by formation of a disulfide bond. 4MTDA can be readily immobilized by coordination of the protein side chain to the lanthanide ion. Large PCSs and RDCs were observed for 4MTDA‐tagged proteins in complexes with paramagnetic lanthanide ions. At an excitation wavelength of 340 nm, the complex formed by protein–4MTDA and Tb³⁺ produces high fluorescence with the main emission at 545 nm. These interesting features of 4MTDA make it a very promising tag that can be exploited in NMR, fluorescence, and EPR spectroscopic studies on protein structure, interaction, and dynamics.     

Mapping the potential energy landscape of intrinsically disordered proteins at amino Acid resolution

Intrinsically disordered regions are predicted to exist in a significant fraction of proteins encoded in eukaryotic genomes. The high levels of conformational plasticity of this class of proteins endows them with unique capacities to act in functional modes not achievable by folded proteins, but also places their molecular characterization beyond the reach of classical structural biology. New techniques are therefore required to understand the relationship between primary sequence and biological function in this class of proteins. Although dependences of some NMR parameters such as chemical shifts (CSs) or residual dipolar couplings (RDCs) on structural propensity are known, so that sampling regimes are often inferred from experimental observation, there is currently no framework that allows for a statistical mapping of the available Ramachandran space of each amino acid in terms of conformational propensity. In this study we develop such an approach, combining highly efficient conformational sampling with ensemble selection to map the backbone conformational sampling of IDPs on a residue specific level. By systematically analyzing the ability of NMR data to map the conformational landscape of disordered proteins, we identify combinations of RDCs and CSs that can be used to raise conformational degeneracies inherent to different data types, and apply these approaches to characterize the conformational behavior of two intrinsically disordered proteins, the K18 domain from Tau protein and N(TAIL) from measles virus nucleoprotein. In both cases, we identify the enhanced populations of turn and helical regions in key regions of the proteins, as well as contiguous strands that show clear and enhanced polyproline II sampling.     

Engineering encodable lanthanide-binding tags into loop regions of proteins

Lanthanide-binding tags (LBTs) are valuable tools for investigation of protein structure, function, and dynamics by NMR spectroscopy, X-ray crystallography, and luminescence studies. We have inserted LBTs into three different loop positions (denoted L, R, and S) of the model protein interleukin-1β (IL1β) and varied the length of the spacer between the LBT and the protein (denoted 1?3). Luminescence studies demonstrate that all nine constructs bind Tb3+ tightly in the low nanomolar range. No significant change in the fusion protein occurs from insertion of the LBT, as shown by two X-ray crystallographic structures of the IL1β-S1 and IL1β-L3 constructs and for the remaining constructs by comparing the 1H?15N heteronuclear single-quantum coherence NMR spectra with that of the wild-type IL1β. Additionally, binding of LBT-loop IL1β proteins to their native binding partner in vitro remains unaltered. X-ray crystallographic phasing was successful using only the signal from the bound lanthanide. Large residual dipolar couplings (RDCs) could be determined by NMR spectroscopy for all LBT-loop constructs and revealed that the LBT-2 series were rigidly incorporated into the interleukin-1β structure. The paramagnetic NMR spectra of loop-LBT mutant IL1β-R2 were assigned and the Δχ tensor components were calculated on the basis of RDCs and pseudocontact shifts. A structural model of the IL1β-R2 construct was calculated using the paramagnetic restraints. The current data provide support that encodable LBTs serve as versatile biophysical tags when inserted into loop regions of proteins of known structure or predicted via homology modeling.     

Searching and optimizing structure ensembles for complex flexible sugars

NMR restrictions are suitable to specify the geometry of a molecule when a single well-defined global free energy minimum exists that is significantly lower than other local minima. Carbohydrates are quite flexible, and therefore, NMR observables do not always correlate with a single conformer but instead with an ensemble of low free energy conformers that can be accessed by thermal fluctuations. In this communication, we describe a novel procedure to identify and weight the contribution to the ensemble of local minima conformers based on comparison to residual dipolar couplings (RDCs) or other NMR observables, such as scalar couplings. A genetic algorithm is implemented to globally minimize the R factor comparing calculated RDCs to experiment. This is done by optimizing the weights of different conformers derived from the exhaustive local minima conformational search program, fast sugar structure prediction software (FSPS). We apply this framework to six human milk sugars, LND-1, LNF-1, LNF-2, LNF-3, LNnT, and LNT, and are able to determine corresponding population weights for the ensemble of conformers. Interestingly, our results indicate that in all cases the RDCs can be well represented by only a few most important conformers. This confirms that several, but not all of the glycosidic linkages in histo-blood group "epitopes" are quite rigid.     

Refinement of NMR structures using implicit solvent and advanced sampling techniques

NMR biomolecular structure calculations exploit simulated annealing methods for conformational sampling and require a relatively high level of redundancy in the experimental restraints to determine quality three-dimensional structures. Recent advances in generalized Born (GB) implicit solvent models should make it possible to combine information from both experimental measurements and accurate empirical force fields to improve the quality of NMR-derived structures. In this paper, we study the influence of implicit solvent on the refinement of protein NMR structures and identify an optimal protocol of utilizing these improved force fields. To do so, we carry out structure refinement experiments for model proteins with published NMR structures using full NMR restraints and subsets of them. We also investigate the application of advanced sampling techniques to NMR structure refinement. Similar to the observations of Xia et al. (J.Biomol. NMR 2002, 22, 317-331), we find that the impact of implicit solvent is rather small when there is a sufficient number of experimental restraints (such as in the final stage of NMR structure determination), whether implicit solvent is used throughout the calculation or only in the final refinement step. The application of advanced sampling techniques also seems to have minimal impact in this case. However, when the experimental data are limited, we demonstrate that refinement with implicit solvent can substantially improve the quality of the structures. In particular, when combined with an advanced sampling technique, the replica exchange (REX) method, near-native structures can be rapidly moved toward the native basin. The REX method provides both enhanced sampling and automatic selection of the most native-like (lowest energy) structures. An optimal protocol based on our studies first generates an ensemble of initial structures that maximally satisfy the available experimental data with conventional NMR software using a simplified force field and then refines these structures with implicit solvent using the REX method. We systematically examine the reliability and efficacy of this protocol using four proteins of various sizes ranging from the 56-residue B1 domain of Streptococcal protein G to the 370-residue Maltose-binding protein. Significant improvement in the structures was observed in all cases when refinement was based on low-redundancy restraint data. The proposed protocol is anticipated to be particularly useful in early stages of NMR structure determination where a reliable estimate of the native fold from limited data can significantly expedite the overall process. This refinement procedure is also expected to be useful when redundant experimental data are not readily available, such as for large multidomain biomolecules and in solid-state NMR structure determination.     

4D solid-state NMR for protein structure determination

Solid-state NMR offers the chance to extend structural studies to proteins that are otherwise difficult to study at atomic resolution, such as protein fibrils, membrane proteins or poorly diffracting crystals. As two-dimensional spatial correlation NMR spectra of proteins suffer from severe resonance overlap, we analyze in this perspective article the potential of higher-dimensional (3D and 4D) proton-detected experiments, which have an increased number of identifiable and assignable distance restraints for solid-state structural studies. We discuss practical considerations for the NMR measurements and the preparation of suitable protein samples and show results of structure calculations from 4D solid-state NMR spectra.     

Anisotropic nuclear spin interactions for the morphology analysis of proteins in solution by NMR spectroscopy.

Determining the relative orientation of domains within a protein is an important problem in structural biology, which has been difficult to address by either X-ray crystallography or NMR. The structure of a multidomain protein in a crystal lattice can be altered by crystal packing forces, resulting in different domain arrangements from those in solution. On the other hand, conventional NMR primarily provides short-range structural information, including proton-proton distances derived from nuclear Overhauser effects (NOEs) and torsion angles through vicinal spin couplings. Thus, NMR cannot always determine the precise interdomain arrangements due to the sparsely observed spin interactions between domains. However, the weak alignment of proteins in solution has enabled a new NMR technique to determine the domain arrangement based on the different structural information, which defines the orientation of a structural unit in protein against the magnetic field. This technique relies on the anisotropic nuclear spin interactions that only occur for a molecule in a weakly aligned state. In this review, the basics of the new NMR approach are described with focusing on its application to domain orientation analysis. We also describe our recently established NMR approach using the same spin interactions, which expands the domain arrangement analysis to higher-molecular weight proteins over 100 kDa.     

Combined Use of Residual Dipolar Couplings and Solution X-ray Scattering To Rapidly Probe Rigid-Body Conformational Transitions in a Non-phosphorylatable Active-Site Mutant of the 128 kDa Enzyme I Dimer

The first component of the bacterial phosphotransferase system, enzyme I (EI), is a multidomain 128 kDa dimer that undergoes large rigid-body conformational transitions during the course of its catalytic cycle. Here we investigate the solution structure of a non-phosphorylatable active-site mutant in which the active-site histidine is substituted by glutamine. We show that perturbations in the relative orientations and positions of the domains and subdomains can be rapidly and reliably determined by conjoined rigid-body/torsion angle/Cartesian simulated annealing calculations driven by orientational restraints from residual dipolar couplings and shape and translation information afforded by small- and wide-angle X-ray scattering. Although histidine and glutamine are isosteric, the conformational space available to a Gln side chain is larger than that for the imidazole ring of His. An additional hydrogen bond between the side chain of Gln189 located on the EIN(α/β) subdomain and an aspartate (Asp129) on the EIN(α) subdomain results in a small (～9°) reorientation of the EIN(α) and EIN(α/β) subdomains that is in turn propagated to a larger reorientation (～26°) of the EIN domain relative to the EIC dimerization domain, illustrating the positional sensitivity of the EIN domain and its constituent subdomains to small structural perturbations.     

Characterization of the conformational equilibrium between the two major substates of RNase A using NMR chemical shifts

Following the recognition that NMR chemical shifts can be used for protein structure determination, rapid advances have recently been made in methods for extending this strategy for proteins and protein complexes of increasing size and complexity. A remaining major challenge is to develop approaches to exploit the information contained in the chemical shifts about conformational fluctuations in native states of proteins. In this work we show that it is possible to determine an ensemble of conformations representing the free energy surface of RNase A using chemical shifts as replica-averaged restraints in molecular dynamics simulations. Analysis of this surface indicates that chemical shifts can be used to characterize the conformational equilibrium between the two major substates of this protein.     

Determination of the structures of symmetric protein oligomers from NMR chemical shifts and residual dipolar couplings

Symmetric protein dimers, trimers, and higher-order cyclic oligomers play key roles in many biological processes. However, structural studies of oligomeric systems by solution NMR can be difficult due to slow tumbling of the system and the difficulty in identifying NOE interactions across protein interfaces. Here, we present an automated method (RosettaOligomers) for determining the solution structures of oligomeric systems using only chemical shifts, sparse NOEs, and domain orientation restraints from residual dipolar couplings (RDCs) without a need for a previously determined structure of the monomeric subunit. The method integrates previously developed Rosetta protocols for solving the structures of monomeric proteins using sparse NMR data and for predicting the structures of both nonintertwined and intertwined symmetric oligomers. We illustrated the performance of the method using a benchmark set of nine protein dimers, one trimer, and one tetramer with available experimental data and various interface topologies. The final converged structures are found to be in good agreement with both experimental data and previously published high-resolution structures. The new approach is more readily applicable to large oligomeric systems than conventional structure-determination protocols, which often require a large number of NOEs, and will likely become increasingly relevant as more high-molecular weight systems are studied by NMR.     

Calculation of residual dipolar couplings from disordered state ensembles using local alignment

Residual dipolar couplings (RDCs) have been observed in disordered states of several proteins. While their nonuniform values were initially surprising, it has been shown that reasonable approximation of experimental RDCs can be obtained using simple statistical coil models and assuming global alignment of each structure, provided that many thousands of conformers are averaged. Here we show that, by using short local alignment tensors, we can achieve good agreement between experimental and simulated RDCs with far fewer structures than required when using global alignment. This makes the possibility of using RDCs as direct restraints in structural calculations of disordered proteins much more feasible. In addition, it provides insight into the nature of RDCs in disordered states, suggesting that they are primarily reporting on local structure.     

Refinement of multidomain protein structures by combination of solution small-angle X-ray scattering and NMR data

Determination of the 3D structures of multidomain proteins by solution NMR methods presents a number of unique challenges related to their larger molecular size and the usual scarcity of constraints at the interdomain interface, often resulting in a decrease in structural accuracy. In this respect, experimental information from small-angle scattering of X-ray radiation in solution (SAXS) presents a suitable complement to the NMR data, as it provides an independent constraint on the overall molecular shape. A computational procedure is described that allows incorporation of such SAXS data into the mainstream high-resolution macromolecular structure refinement. The method is illustrated for a two-domain 177-amino-acid protein, gammaS crystallin, using an experimental SAXS data set fitted at resolutions from approximately 200 A to approximately 30 A. Inclusion of these data during structure refinement decreases the backbone coordinate root-mean-square difference between the derived model and the high-resolution crystal structure of a 54% homologous gammaB crystallin from 1.96 +/- 0.07 A to 1.31 +/- 0.04 A. Combining SAXS data with NMR restraints can be accomplished at a moderate computational expense and is expected to become useful for multidomain proteins, multimeric assemblies, and tight macromolecular complexes.     

Protein Structural Ensembles Visualized by Solvent Paramagnetic Relaxation Enhancement

A protein can be in different conformations when fulfilling its function. Yet depiction of protein structural ensembles remains difficult. Here we show that the accurate measurement of solvent paramagnetic relaxation enhancement (sPRE) in the presence of an inert paramagnetic cosolute allows the assessment of protein dynamics. Demonstrated with two multi‐domain proteins, we present a method to characterize protein microsecond–millisecond dynamics based on the analysis of the sPRE. Provided with the known structures of a protein, our method uncovers an ensemble of structures that fully accounts for the observed sPRE. In conjunction with molecular dynamics simulations, our method can identify protein alternative conformation that has only been theorized before. Together, our method expands the application of sPRE beyond structural characterization of rigid proteins and complements the established PRE NMR technique.