Identification and removal of nitroxide spin label contaminant: impact on PRE studies of α-helical membrane proteins in detergent

NMR paramagnetic relaxation enhancement (PRE) provides long‐range distance constraints (～15–25 Å) that can be critical to determining overall protein topology, especially where long‐range NOE information is unavailable such as in the case of larger proteins that require deuteration. However, several challenges currently limit the use of NMR PRE for α‐helical membrane proteins. One challenge is the nonspecific association of the nitroxide spin label to the protein‐detergent complex that can result in spurious PRE derived distance restraints. The effect of the nitroxide spin label contaminant is evaluated and quantified and a robust method for the removal of the contaminant is provided to advance the application of PRE restraints to membrane protein NMR structure determination.     

Structure of hen egg-white lysozyme solvated in TFE/water: a molecular dynamics simulation study based on NMR data

Various experimental studies of hen egg white lysozyme (HEWL) in water and TFE/water clearly indicate structural differences between the native state and TFE state of HEWL, e.g. the helical content of the protein in the TFE state is much higher than in the native state. However, the available detailed NMR studies were not sufficient to determine fully a structure of HEWL in the TFE state. Different molecular dynamics (MD) simulations, i.e. at room temperature, at increased temperature and using proton–proton distance restraints derived from NMR NOE data, have been used to generate configurational ensembles corresponding to the TFE state of HEWL. The configurational ensemble obtained at room temperature using atom-atom distance restraints measured for HEWL in TFE/water solution satisfies the experimental data and has the lowest protein energy. In this ensemble residues 50–58, which are part of the β-sheet in native HEWL, adopt fluctuating α-helical secondary structure.     

The impact of window functions on NMR-based paramagnetic relaxation enhancement measurements in membrane proteins

Though challenging, solution NMR spectroscopy allows fundamental interrogation of the structure and dynamics of membrane proteins. One major technical hurdle in studies of helical membrane proteins by NMR is the difficulty of obtaining sufficient long range NOEs to determine tertiary structure. For this reason, long range distance information is sometimes sought through measurement of paramagnetic relaxation enhancements (PRE) of NMR nuclei as a function of distance from an introduced paramagnetic probe. Current PRE interpretation is based on the assumption of Lorentzian resonance lineshapes. However, in order to optimize spectral resolution, modern multidimensional NMR spectra are almost always subjected to resolution-enhancement, leading to distortions in the Lorentizian peak shape. Here it is shown that when PREs are derived using peak intensities (i.e., peak height) and linewidths from both real and simulated spectra that were produced using a wide range of apodization/window functions, that there is little variation in the distances determined (< 1 Å at the extremes). This indicates that the high degree of resolution enhancement required to obtain well-resolved spectra from helical membrane proteins is compatible with the use of PRE data as a source of distance restraints. While these conclusions are particularly important for helical membrane proteins, they are generally applicable to all PRE measurements made using resolution-enhanced data.     

A rigid disulfide-linked nitroxide side chain simplifies the quantitative analysis of PRE data

The measurement of (1)H transverse paramagnetic relaxation enhancement (PRE) has been used in biomolecular systems to determine long-range distance restraints and to visualize sparsely-populated transient states. The intrinsic flexibility of most nitroxide and metal-chelating paramagnetic spin-labels, however, complicates the quantitative interpretation of PREs due to delocalization of the paramagnetic center. Here, we present a novel, disulfide-linked nitroxide spin label, R1p, as an alternative to these flexible labels for PRE studies. When introduced at solvent-exposed α-helical positions in two model proteins, calmodulin (CaM) and T4 lysozyme (T4L), EPR measurements show that the R1p side chain exhibits dramatically reduced internal motion compared to the commonly used R1 spin label (generated by reacting cysteine with the spin labeling compound often referred to as MTSL). Further, only a single nitroxide position is necessary to account for the PREs arising from CaM S17R1p, while an ensemble comprising multiple conformations is necessary for those observed for CaM S17R1. Together, these observations suggest that the nitroxide adopts a single, fixed position when R1p is placed at solvent-exposed α-helical positions, greatly simplifying the interpretation of PRE data by removing the need to account for the intrinsic flexibility of the spin label.     

Algorithm for selection of optimized EPR distance restraints for de novo protein structure determination

A hybrid protein structure determination approach combining sparse Electron Paramagnetic Resonance (EPR) distance restraints and Rosetta de novo protein folding has been previously demonstrated to yield high quality models (Alexander et al. (2008)). However, widespread application of this methodology to proteins of unknown structures is hindered by the lack of a general strategy to place spin label pairs in the primary sequence. In this work, we report the development of an algorithm that optimally selects spin labeling positions for the purpose of distance measurements by EPR. For the α-helical subdomain of T4 lysozyme (T4L), simulated restraints that maximize sequence separation between the two spin labels while simultaneously ensuring pairwise connectivity of secondary structure elements yielded vastly improved models by Rosetta folding. 54% of all these models have the correct fold compared to only 21% and 8% correctly folded models when randomly placed restraints or no restraints are used, respectively. Moreover, the improvements in model quality require a limited number of optimized restraints, which is determined by the pairwise connectivities of T4L α-helices. The predicted improvement in Rosetta model quality was verified by experimental determination of distances between spin labels pairs selected by the algorithm. Overall, our results reinforce the rationale for the combined use of sparse EPR distance restraints and de novo folding. By alleviating the experimental bottleneck associated with restraint selection, this algorithm sets the stage for extending computational structure determination to larger, traditionally elusive protein topologies of critical structural and biochemical importance.     

Intrinsic order and disorder in the bcl-2 member harakiri: insights into its proapoptotic activity

Harakiri is a BH3-only member of the Bcl-2 family that localizes in membranes and induces cell death by binding to prosurvival Bcl-x(L) and Bcl-2. The cytosolic domain of Harakiri is largely disorder with residual α-helical conformation according to previous structural studies. As these helical structures could play an important role in Harakiri's function, we have used NMR and circular dichroism to fully characterize them at the residue-atomic level. In addition, we report structural studies on a peptide fragment spanning Harakiri's C-terminal hydrophobic sequence, which potentially operates as a transmembrane domain. We initially checked by enzyme immunoassays and NMR that peptides encompassing different lengths of the cytosolic domain are functional as they bind Bcl-x(L) and Bcl-2. The structural data in water indicate that the α-helical conformation is restricted to a 25-residue segment comprising the BH3 domain. However, structure calculation was precluded because of insufficient NMR restraints. To bypass this problem we used alcohol-water mixture to increase structure population and confirmed by NMR that the conformation in both milieus is equivalent. The resulting three-dimensional structure closely resembles that of peptides encompassing the BH3 domain of BH3-only members in complex with their prosurvival partners, suggesting that preformed structural elements in the disordered protein are central to binding. In contrast, the transmembrane domain forms in micelles a monomeric α-helix with a population close to 100%. Its three-dimensional structure here reported reveals features that explain its function as membrane anchor. Altogether these results are used to propose a tentative structural model of how Harakiri works.     

NMR Model of PrgI–SipD Interaction and Its Implications in the Needle-Tip Assembly of the Salmonella Type III Secretion System

Salmonella and other pathogenic bacteria use the type III secretion system (T3SS) to inject virulence proteins into human cells to initiate infections. The structural component of the T3SS contains a needle and a needle tip. The needle is assembled from PrgI needle protomers and the needle tip is capped with several copies of the SipD tip protein. How a tip protein docks on the needle is unclear. A crystal structure of a PrgI–SipD fusion protein docked on the PrgI needle results in steric clash of SipD at the needle tip when modeled on the recent atomic structure of the needle. Thus, there is currently no good model of how SipD is docked on the PrgI needle tip. Previously, we showed by NMR paramagnetic relaxation enhancement (PRE) methods that a specific region in the SipD coiled coil is the binding site for PrgI. Others have hypothesized that a domain of the tip protein—the N-terminal α-helical hairpin—has to swing away during the assembly of the needle apparatus. Here, we show by PRE methods that a truncated form of SipD lacking the α-helical hairpin domain binds more tightly to PrgI. Further, PRE-based structure calculations revealed multiple PrgI binding sites on the SipD coiled coil. Our PRE results together with the recent NMR-derived atomic structure of the Salmonella needle suggest a possible model of how SipD might dock at the PrgI needle tip. SipD and PrgI are conserved in other bacterial T3SSs; thus, our results have wider implication in understanding other needle-tip complexes.     

How Amphipols Embed Membrane Proteins: Global Solvent Accessibility and Interaction with a Flexible Protein Terminus

Amphipathic polymers called amphipols provide a valuable alternative to detergents for keeping integral membrane proteins soluble in aqueous buffers. Here, we characterize spatial contacts of amphipol A8-35 with membrane proteins from two architectural classes: The 8-stranded β-barrel outer membrane protein OmpX and the α-helical protein bacteriorhodopsin. OmpX is well structured in A8-35, with its barrel adopting a fold closely similar to that in dihexanoylphosphocholine micelles. The accessibility of A8-35-trapped OmpX by a water-soluble paramagnetic molecule is highly similar to that in detergent micelles and resembles the accessibility in the natural membrane. For the α-helical protein bacteriorhodopsin, previously shown to keep its fold and function in amphipols, NMR data show that the imidazole protons of a polyhistidine tag at the N-terminus of the protein are exchange protected in the presence of detergent and lipid bilayer nanodiscs, but not in amphipols, indicating the absence of an interaction in the latter case. Overall, A8-35 exhibits protein interaction properties somewhat different from detergents and lipid bilayer nanodiscs, while maintaining the structure of solubilized integral membrane proteins.     

Solution NMR studies of polytopic α-helical membrane proteins

NMR spectroscopy has established itself as one of the main techniques for the structural study of integral membrane proteins. Remarkably, over the last few years, substantial progress has been achieved in the structure determination of increasingly complex polytopical α-helical membrane proteins, with their size approaching ～100kDa. Such advances are the result of significant improvements in NMR methodology, sample preparation and powerful selective isotope labelling schemes. We review the requirements facilitating such work based on the more recent solution NMR studies of α-helical proteins. While the majority of such studies still use detergent-solubilized proteins, alternative more native-like lipid-based media are emerging. Recent interaction, dynamics and conformational studies are discussed that cast a promising light on the future role of NMR in this important and exciting area.     

NMR structures of polytopic integral membrane proteins

Membrane proteins represent up to 30% of the proteins in all organisms, they are involved in many biological processes and are the molecular targets for around 50% of validated drugs. Despite this, membrane proteins represent less than 1% of all high-resolution protein structures due to various challenges associated with applying the main biophysical techniques used for protein structure determination. Recent years have seen an explosion in the number of high-resolution structures of membrane proteins determined by NMR spectroscopy, especially for those with multiple transmembrane-spanning segments. This is a review of the structures of polytopic integral membrane proteins determined by NMR spectroscopy up to the end of the year 2010, which includes both β-barrel and α-helical proteins from a number of different organisms and with a range in types of function. It also considers the challenges associated with performing structural studies by NMR spectroscopy on membrane proteins and how some of these have been overcome, along with its exciting potential for contributing new knowledge about the molecular mechanisms of membrane proteins, their roles in human disease, and for assisting drug design.     

Backbone structure of a small helical integral membrane protein: A unique structural characterization

The structural characterization of small integral membrane proteins pose a significant challenge for structural biology because of the multitude of molecular interactions between the protein and its heterogeneous environment. Here, the three‐dimensional backbone structure of Rv1761c from Mycobacterium tuberculosis has been characterized using solution NMR spectroscopy and dodecylphosphocholine (DPC) micelles as a membrane mimetic environment. This 127 residue single transmembrane helix protein has a significant (10 kDa) C‐terminal extramembranous domain. Five hundred and ninety distance, backbone dihedral, and orientational restraints were employed resulting in a 1.16 Å rmsd backbone structure with a transmembrane domain defined at 0.40 Å. The structure determination approach utilized residual dipolar coupling orientation data from partially aligned samples, long‐range paramagnetic relaxation enhancement derived distances, and dihedral restraints from chemical shift indices to determine the global fold. This structural model of Rv1761c displays some influences by the membrane mimetic illustrating that the structure of these membrane proteins is dictated by a combination of the amino acid sequence and the protein's environment. These results demonstrate both the efficacy of the structural approach and the necessity to consider the biophysical properties of membrane mimetics when interpreting structural data of integral membrane proteins and, in particular, small integral membrane proteins.     

NMR observable-based structure refinement of DAP12-NKG2C activating immunoreceptor complex in explicit membranes

NMR observables, such as NOE-based distance measurements, are increasingly being used to characterize membrane protein structures. However, challenges in membrane protein NMR studies often yield a relatively small number of such restraints that can create ambiguities in defining critical side chain-side chain interactions. In the recent solution NMR structure of the DAP12-NKG2C immunoreceptor transmembrane helix complex, five functionally required interfacial residues (two Asps and two Thrs in the DAP12 dimer and one Lys in NKG2C) display a surprising arrangement in which one Asp side chain faces the membrane hydrophobic core. To explore whether these side-chain interactions are energetically optimal, we used the published distance restraints for molecular dynamics simulations in explicit micelles and bilayers. The structures refined by this protocol are globally similar to the published structure, but the side chains of those five residues form a stable network of salt bridges and hydrogen bonds, leaving the Asp side chain shielded from the hydrophobic core, which is also consistent with available experimental observations. Moreover, the simulations show similar short-range interactions between the transmembrane complex and lipid/detergent molecules in micelles and bilayers, respectively. This study illustrates the efficacy of NMR membrane protein structure refinements in explicit membrane systems.     

High-resolution membrane protein structure by joint calculations with solid-state NMR and X-ray experimental data

X-ray diffraction and nuclear magnetic resonance spectroscopy (NMR) are the staple methods for revealing atomic structures of proteins. Since crystals of biomolecular assemblies and membrane proteins often diffract weakly and such large systems encroach upon the molecular tumbling limit of solution NMR, new methods are essential to extend structures of such systems to high resolution. Here we present a method that incorporates solid-state NMR restraints alongside of X-ray reflections to the conventional model building and refinement steps of structure calculations. Using the 3.7 Å crystal structure of the integral membrane protein complex DsbB-DsbA as a test case yielded a significantly improved backbone precision of 0.92 Å in the transmembrane region, a 58% enhancement from using X-ray reflections alone. Furthermore, addition of solid-state NMR restraints greatly improved the overall quality of the structure by promoting 22% of DsbB transmembrane residues into the most favored regions of Ramachandran space in comparison to the crystal structure. This method is widely applicable to any protein system where X-ray data are available, and is particularly useful for the study of weakly diffracting crystals.     

Solid-state NMR structures of integral membrane proteins

Abstract

Solid-state NMR is unique for its ability to obtain three-dimensional structures and to measure atomic-resolution structural and dynamic information for membrane proteins in native lipid bilayers. An increasing number and complexity of integral membrane protein structures have been determined by solid-state NMR using two main methods. Oriented sample solid-state NMR uses macroscopically aligned lipid bilayers to obtain orientational restraints that define secondary structure and global fold of embedded peptides and proteins and their orientation and topology in lipid bilayers. Magic angle spinning (MAS) solid-state NMR uses unoriented rapidly spinning samples to obtain distance and torsion angle restraints that define tertiary structure and helix packing arrangements. Details of all current protein structures are described, highlighting developments in experimental strategy and other technological advancements. Some structures originate from combining solid- and solution-state NMR information and some have used solid-state NMR to refine X-ray crystal structures. Solid-state NMR has also validated the structures of proteins determined in different membrane mimetics by solution-state NMR and X-ray crystallography and is therefore complementary to other structural biology techniques. By continuing efforts in identifying membrane protein targets and developing expression, isotope labelling and sample preparation strategies, probe technology, NMR experiments, calculation and modelling methods and combination with other techniques, it should be feasible to determine the structures of many more membrane proteins of biological and biomedical importance using solid-state NMR. This will provide three-dimensional structures and atomic-resolution structural information for characterising ligand and drug interactions, dynamics and molecular mechanisms of membrane proteins under physiological lipid bilayer conditions.     

NMR structures of polytopic integral membrane proteins

Abstract

Membrane proteins represent up to 30% of the proteins in all organisms, they are involved in many biological processes and are the molecular targets for around 50% of validated drugs. Despite this, membrane proteins represent less than 1% of all high-resolution protein structures due to various challenges associated with applying the main biophysical techniques used for protein structure determination. Recent years have seen an explosion in the number of high-resolution structures of membrane proteins determined by NMR spectroscopy, especially for those with multiple transmembrane-spanning segments. This is a review of the structures of polytopic integral membrane proteins determined by NMR spectroscopy up to the end of the year 2010, which includes both β-barrel and α-helical proteins from a number of different organisms and with a range in types of function. It also considers the challenges associated with performing structural studies by NMR spectroscopy on membrane proteins and how some of these have been overcome, along with its exciting potential for contributing new knowledge about the molecular mechanisms of membrane proteins, their roles in human disease, and for assisting drug design.     

Folding determinants of disulfide bond forming protein B explored by solution nuclear magnetic resonance spectroscopy

The two-stage model for membrane protein folding postulates that individual helices form first and are subsequently packed against each other. To probe the two-stage model, the structures of peptides representing individual transmembrane helices of the disulfide bond forming protein B have been studied in trifluoroethanol solution as well as in detergent micelles using nuclear magnetic resonance (NMR) and circular dichroism spectroscopy. In TFE solution, peptides showed well-defined α-helical structures. Peptide structures in TFE were compared to the structures of full-length protein obtained by X-ray crystallography and NMR. The extent of α-helical secondary structure coincided well, lending support for the two-stage model for membrane protein folding. However, the conformation of some amino acid side chains differs between the structures of peptide and full-length protein. In micellar solution, the peptides also adopted a helical structure, albeit of reduced definition. Using measurements of paramagnetic relaxation enhancement, peptides were confirmed to be embedded in micelles. These observations may indicate that in the native protein, tertiary interactions additionally stabilize the secondary structure of the individual transmembrane helices.     

Solid-State NMR-Restrained Ensemble Dynamics of a Membrane Protein in Explicit Membranes

Solid-state NMR has been used to determine the structures of membrane proteins in native-like lipid bilayer environments. Most structure calculations based on solid-state NMR observables are performed using simulated annealing with restrained molecular dynamics and an energy function, where all nonbonded interactions are represented by a single, purely repulsive term with no contributions from van der Waals attractive, electrostatic, or solvation energy. To our knowledge, this is the first application of an ensemble dynamics technique performed in explicit membranes that uses experimental solid-state NMR observables to obtain the refined structure of a membrane protein together with information about its dynamics and its interactions with lipids. Using the membrane-bound form of the fd coat protein as a model membrane protein and its experimental solid-state NMR data, we performed restrained ensemble dynamics simulations with different ensemble sizes in explicit membranes. For comparison, a molecular dynamics simulation of fd coat protein was also performed without any restraints. The average orientation of each protein helix is similar to a structure determined by traditional single-conformer approaches. However, their variations are limited in the resulting ensemble of structures with one or two replicas, as they are under the strong influence of solid-state NMR restraints. Although highly consistent with all solid-state NMR observables, the ensembles of more than two replicas show larger orientational variations similar to those observed in the molecular dynamics simulation without restraints. In particular, in these explicit membrane simulations, Lys⁴⁰, residing at the C-terminal side of the transmembrane helix, is observed to cause local membrane curvature. Therefore, compared to traditional single-conformer approaches in implicit environments, solid-state NMR restrained ensemble simulations in explicit membranes readily characterize not only protein dynamics but also protein-lipid interactions in detail.     

Exhaustive enumeration of protein conformations using experimental restraints.

We present an efficient new algorithm that enumerates all possible conformations of a protein that satisfy a given set of distance restraints. Rapid growth of all possible self-avoiding conformations on the diamond lattice provides construction of alpha-carbon representations of a protein fold. We investigated the dependence of the number of conformations on pairwise distance restraints for the proteins crambin, pancreatic trypsin inhibitor, and ubiquitin. Knowledge of between one and two contacts per monomer is shown to be sufficient to restrict the number of candidate structures to approximately 1,000 conformations. Pairwise RMS deviations of atomic position comparisons between pairs of these 1,000 structures revealed that these conformations can be grouped into about 25 families of structures. These results suggest a new approach to assessing alternative protein folds given a very limited number of distance restraints. Such restraints are available from several experimental techniques such as NMR, NOESY, energy transfer fluorescence spectroscopy, and crosslinking experiments. This work focuses on exhaustive enumeration of protein structures with emphasis on the possible use of NOESY-determined distance restraints.     

A computational method for NMR-constrained protein threading.

Protein threading provides an effective method for fold recognition and backbone structure prediction. But its application is currently limited due to its level of prediction accuracy and scope of applicability. One way to significantly improve its usefulness is through the incorporation of underconstrained (or partial) NMR data. It is well known that the NMR method for protein structure determination applies only to small proteins and that its effectiveness decreases rapidly as the protein mass increases beyond about 30 kD. We present, in this paper, a computational framework for applying underconstrained NMR data (that alone are insufficient for structure determination) as constraints in protein threading and also in all-atom model construction. In this study, we consider both secondary structure assignments from chemical shifts and NOE distance restraints. Our results have shown that both secondary structure assignments and a small number of long-range NOEs can significantly improve the threading quality in both fold recognition and threading-alignment accuracy, and can possibly extend threading's scope of applicability from homologs to analogs. An accurate backbone structure generated by NMR-constrained threading can then provide a great amount of structural information, equivalent to that provided by many NMR data; and hence can help reduce the number of NMR data typically required for an accurate structure determination. This new technique can potentially accelerate current NMR structure determination processes and possibly expand NMR's capability to larger proteins.     

Structural Interpretation of Paramagnetic Relaxation Enhancement-Derived Distances for Disordered Protein States

Paramagnetic relaxation enhancement (PRE) is a powerful technique for studying transient tertiary organizations of unfolded and partially folded proteins. The heterogeneous and dynamic nature of disordered protein states, together with the r⁻⁶ dependence of PRE, presents significant challenges for reliable structural interpretation of PRE-derived distances. Without additional knowledge of accessible conformational substates, ensemble-simulation-based protocols have been used to calculate structure ensembles that appear to be consistent with the PRE distance restraints imposed on the ensemble level with the proper r⁻⁶ weighting. However, rigorous assessment of the reliability of such protocols has been difficult without intimate knowledge of the true nature of disordered protein states. Here we utilize sets of theoretical PRE distances derived from simulated structure ensembles that represent the folded, partially folded and unfolded states of a small protein to investigate the efficacy of ensemble-simulation-based structural interpretation of PRE distances. The results confirm a critical limitation that, due to r⁻⁶ weighting, only one or a few members need to satisfy the distance restraints and the rest of the ensemble are essentially unrestrained. Consequently, calculated structure ensembles will appear artificially heterogeneous no matter whether the PRE distances are derived from the folded, partially unfolded or unfolded state. Furthermore, the nature of the heterogeneous ensembles is largely determined by the protein model employed in structure calculation and reflects little on the true nature of the underlying disordered state. These findings suggest that PRE measurements on disordered protein states alone generally do not contain enough information for a reliable structural interpretation and that the latter will require additional knowledge of accessible conformational substates. Interestingly, when a very large number of PRE measurements is available, faithful structural interpretation might be possible with intermediate ensemble sizes under ideal conditions.