Use of chemical shifts in macromolecular structure determination

Using chemical shifts to assess transient secondary structure and generate ensemble structures of intrinsically disordered proteins

The chemical shifts of backbone atoms in polypeptides are sensitive to the dihedral angles phi and psi and can be used to estimate transient secondary structure and to generate structural ensembles of intrinsically disordered proteins (IDPs). In this chapter, several of the random coil reference databases used to estimate transient secondary structure are described, and the procedure is outlined for using these databases to estimate transient secondary structure. A new protocol is also presented for generating a diverse ensemble of structures for an IDP and reweighting these structures to optimize the fit between simulated and experimental chemical shift values.

Mechanisms of allosteric gene regulation by NMR quantification of microsecond-millisecond protein dynamics

The trp RNA-binding attenuation protein (TRAP) is a paradigmatic allosteric protein that regulates the tryptophan biosynthetic genes associated with the trp operon in bacilli. The ring-shaped 11-mer TRAP is activated for recognition of a specific trp-mRNA target by binding up to 11 tryptophan molecules. To characterize the mechanisms of tryptophan-induced TRAP activation, we have performed methyl relaxation dispersion (MRD) nuclear magnetic resonance (NMR) experiments that probe the time-dependent structure of TRAP in the microsecond-to-millisecond "chemical exchange" time window. We find significant side chain flexibility localized to the RNA and tryptophan binding sites of the apo protein and that these dynamics are dramatically reduced upon ligand binding. Analysis of the MRD NMR data provides insights into the structural nature of transiently populated conformations sampled in solution by apo TRAP. The MRD data are inconsistent with global two-state exchange, indicating that conformational sampling in apo TRAP is asynchronous. These findings imply a temporally heterogeneous population of structures that are incompatible with RNA binding and substantiate the study of TRAP as a paradigm for probing and understanding essential dynamics in allosteric, regulatory proteins.

Assembly of the complex between archaeal RNase P proteins RPP30 and Pop5

RNase P is a highly conserved ribonucleoprotein enzyme that represents a model complex for understanding macromolecular RNA-protein interactions. Archaeal RNase P consists of one RNA and up to five proteins (Pop5, RPP30, RPP21, RPP29, and RPP38/L7Ae). Four of these proteins function in pairs (Pop5-RPP30 and RPP21–RPP29). We have used nuclear magnetic resonance (NMR) spectroscopy and isothermal titration calorimetry (ITC) to characterize the interaction between Pop5 and RPP30 from the hyperthermophilic archaeon Pyrococcus furiosus (Pfu). NMR backbone resonance assignments of free RPP30 (25 kDa) indicate that the protein is well structured in solution, with a secondary structure matching that observed in a closely related crystal structure. Chemical shift perturbations upon the addition of Pop5 (14 kDa) reveal its binding surface on RPP30. ITC experiments confirm a net 1 : 1 stoichiometry for this tight protein-protein interaction and exhibit complex isotherms, indicative of higher-order binding. Indeed, light scattering and size exclusion chromatography data reveal the complex to exist as a 78 kDa heterotetramer with two copies each of Pop5 and RPP30. These results will inform future efforts to elucidate the functional role of the Pop5-RPP30 complex in RNase P assembly and catalysis.

Calculation of chemical shift anisotropy in proteins

Individual peptide groups in proteins must exhibit some variation in the chemical shift anisotropy (CSA) of their constituent atoms, but not much is known about the extent or origins of this dispersion. Direct spectroscopic measurement of CSA remains technically challenging, and theoretical methods can help to overcome these limitations by estimating shielding tensors for arbitrary structures. Here we use an automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach to compute (15)N, (13)C' and (1)H chemical shift tensors for human ubiquitin and the GB1 and GB3 fragments of staphylococcal protein G. The average and range of variation of the anisotropies is in good agreement with experimental estimates from solid-state NMR, and the variation among residues is somewhat smaller than that estimated from solution-state measurements. Hydrogen-bond effects account for much of the variation, both between helix and sheet regions, and within elements of secondary structure, but other effects (including variations in torsion angles) may play a role as well.

An environment-dependent structural switch underlies the regulation of carnitine palmitoyltransferase 1A

The enzyme carnitine palmitoyltransferase 1 (CPT1), which is anchored in the outer mitochondrial membrane (OMM), controls the rate-limiting step in fatty acid β-oxidation in mammalian tissues. It is inhibited by malonyl-CoA, the first intermediate of fatty acid synthesis, and it responds to OMM curvature and lipid characteristics, which reflect long term nutrient/hormone availability. Here, we show that the N-terminal regulatory domain (N) of CPT1A can adopt two complex amphiphilic structural states, termed Nα and Nβ, that interchange in a switch-like manner in response to offered binding surface curvature. Structure-based site-directed mutageneses of native CPT1A suggest Nα to be inhibitory and Nβ to be noninhibitory, with the relative Nα/Nβ ratio setting the prevalent malonyl-CoA sensitivity of the enzyme. Based on the amphiphilic nature of N and molecular modeling, we propose malonyl-CoA sensitivity to be coupled to the properties of the OMM by Nα-OMM associations that alter the Nα/Nβ ratio. For enzymes residing at the membrane-water interface, this constitutes an integrative regulatory mechanism of exceptional sophistication.

The A53T mutation is key in defining the differences in the aggregation kinetics of human and mouse α-synuclein

Despite a 95% sequence similarity, the aggregation of human and mouse α-synuclein is remarkably different, as the human form is slower than the mouse form in forming fibrils but is associated with Parkinson's disease in both humans and transgenic mice. Here, the amino acid code underlying these differences is investigated by comparing the lag times, growth rates, and secondary structure propensities of a systematic series of eight human-mouse chimeras. Fluorescence analysis of these variants shows that the A53T substitution dominates the growth kinetics, while the lag phase is affected by a combination of the A53T and S87N substitutions. The secondary structure propensities derived from an NMR chemical shift analysis of the monomeric forms of the human-mouse variants enable us to establish a link between the changes in the conformational properties in the region of position 53 upon mutation and the corresponding changes in growth rates. These results suggest that the presence of an alanine residue at position 53 may be an evolutionary adaptation to minimize Parkinson's disease in humans and indicates that effective drug development efforts may be directed to target this N-terminal region of the sequence.

Comparative NMR analysis of an 80-residue G protein-coupled receptor fragment in two membrane mimetic environments

Fragments of integral membrane proteins have been used to study the physical chemical properties of regions of transporters and receptors. Ste2p(G31-T110) is an 80-residue polypeptide which contains a portion of the N-terminal domain, transmembrane domain 1 (TM1), intracellular loop 1, TM2 and part of extracellular loop 1 of the α-factor receptor (Ste2p) from Saccharomyces cerevisiae. The structure of this peptide was previously determined to form a helical hairpin in lyso-palmitoylphosphatidyl-glycerol micelles (LPPG) [1]. Herein, we perform a systematic comparison of the structure of this protein fragment in micelles and trifluoroethanol (TFE):water in order to understand whether spectra recorded in organic:aqueous medium can facilitate the structure determination in a micellar environment. Using uniformly labeled peptide and peptide selectively protonated on Ile, Val and Leu methyl groups in a perdeuterated background and a broad set of 3D NMR experiments we assigned 89% of the observable atoms. NOEs and chemical shift analysis were used to define the helical regions of the fragment. Together with constraints from paramagnetic spin labeling, NOEs were used to calculate a transiently folded helical hairpin structure for this peptide in TFE:water. Correlation of chemical shifts was insufficient to transfer assignments from TFE:water to LPPG spectra in the absence of further information.

Nonnative interactions in the FF domain folding pathway from an atomic resolution structure of a sparsely populated intermediate: an NMR relaxation dispersion study

Several all-helical single-domain proteins have been shown to fold rapidly (microsecond time scale) to a compact intermediate state and subsequently rearrange more slowly to the native conformation. An understanding of this process has been hindered by difficulties in experimental studies of intermediates in cases where they are both low-populated and only transiently formed. One such example is provided by the on-pathway folding intermediate of the small four-helix bundle FF domain from HYPA/FBP11 that is populated at several percent with a millisecond lifetime at room temperature. Here we have studied the L24A mutant that has been shown previously to form nonnative interactions in the folding transition state. A suite of Carr-Purcell-Meiboom-Gill relaxation dispersion NMR experiments have been used to measure backbone chemical shifts and amide bond vector orientations of the invisible folding intermediate that form the input restraints in calculations of atomic resolution models of its structure. Despite the fact that the intermediate structure has many features that are similar to that of the native state, a set of nonnative contacts is observed that is even more extensive than noted previously for the wild-type (WT) folding intermediate. Such nonnative interactions, which must be broken prior to adoption of the native conformation, explain why the transition from the intermediate state to the native conformer (millisecond time scale) is significantly slower than from the unfolded ensemble to the intermediate and why the L24A mutant folds more slowly than the WT.

Probing the self-association, intermolecular contacts, and folding propensity of amelogenin

Amelogenins are an intrinsically disordered protein family that plays a major role in the development of tooth enamel, one of the most highly mineralized materials in nature. Monomeric porcine amelogenin possesses random coil and residual secondary structures, but it is not known which sequence regions would be conformationally attractive to potential enamel matrix targets such as other amelogenins (self-assembly), other matrix proteins, cell surfaces, or biominerals. To address this further, we investigated recombinant porcine amelogenin (rP172) using "solvent engineering" techniques to simultaneously promote native-like structure and induce amelogenin oligomerization in a manner that allows identification of intermolecular contacts between amelogenin molecules. We discovered that in the presence of 2,2,2-trifluoroethanol (TFE) significant folding transitions and stabilization occurred primarily within the N- and C-termini, while the polyproline Type II central domain was largely resistant to conformational transitions. Seven Pro residues (P2, P127, P130, P139, P154, P157, P162) exhibited conformational response to TFE, and this indicates these Pro residues act as folding enhancers in rP172. The remaining Pro residues resisted TFE perturbations and thus act as conformational stabilizers. We also noted that TFE induced rP172 self-association via the formation of intermolecular contacts involving P4-H6, V19-P33, and E40-T58 regions of the N-terminus. Collectively, these results confirm that the N- and C-termini of amelogenin are conformationally responsive and represent potential interactive sites for amelogenin-target interactions during enamel matrix mineralization. Conversely, the Pro, Gln central domain is resistant to folding and this may have important functional significance for amelogenin.

Measurement of protein unfolding/refolding kinetics and structural characterization of hidden intermediates by NMR relaxation dispersion

Detailed understanding of protein function and malfunction hinges on the ability to characterize transiently populated states and the transitions between them. Here, we use (15)N, , and (13)CO NMR R(2) relaxation dispersion to investigate spontaneous unfolding and refolding events of native apomyoglobin. Above pH 5.0, dispersion is dominated by processes involving fluctuations of the F-helix region, which is invisible in NMR spectra. Measurements of R(2) dispersion for residues contacted by the F-helix region in the native (N) structure reveal a transient state formed by local unfolding of helix F and undocking from the protein core. A similar state was detected at pH 4.75-4.95 and determined to be an on-pathway intermediate (I1) in a linear three-state unfolding scheme (N&lrarr2;I1&lrarr2;MG) leading to a transiently populated molten globule (MG) state. The slowest steps in unfolding and refolding are N → I1 (36 s(-1)) and MG → I1 (26 s(-1)), respectively. Differences in chemical shift between N and I1 are very small, except in regions adjacent to helix F, showing that their core structures are similar. Chemical shift changes between the N and MG states, obtained from R(2) dispersion, reveal that the transient MG state is structurally similar to the equilibrium MG observed previously at high temperature and low pH. Analysis of MG state chemical shifts shows the location of residual helical structure in the transient intermediate and identifies regions that unfold or rearrange into nonnative structure during the N → MG transition. The experiments also identify regions of energetic frustration that "crack" during unfolding and impede the refolding process.

Measuring 1HN temperature coefficients in invisible protein states by relaxation dispersion NMR spectroscopy

A method based on the Carr-Purcell-Meiboom-Gill relaxation dispersion experiment is presented for measuring the temperature coefficients of amide proton chemical shifts of low populated 'invisible' protein states that exchange with a 'visible' ground state on the millisecond time-scale. The utility of the approach is demonstrated with an application to an I58D mutant of the Pfl6 Cro protein that undergoes exchange between the native, folded state and a cold denatured, unfolded conformational ensemble that is populated at a level of 6% at 2.5°C. A wide distribution of amide temperature coefficients is measured for the unfolded state. The distribution is centered about -5.6 ppb/K, consistent with an absence of intra-molecular hydrogen bonds, on average. However, the large range of values (standard deviation of 2.1 ppb/K) strongly supports the notion that the unfolded state of the protein is not a true random coil polypeptide chain.

Secondary Structure Assignment of Amyloid-β Peptide Using Chemical Shifts

The distinct conformational dependence of chemical shifts caused by α-helices and β-sheets renders NMR chemical shift analysis a powerful tool for the structural determination of proteins. However, the time scale of NMR experiments can make a secondary structure assignment of highly flexible peptides or proteins, which may be converting between conformational substates, problematic. For instance the amyloid-β monomer, according to NMR chemical shifts, adopts a predominately random coil structure in aqueous solution (with <3% α-helical content). Molecular dynamics simulations, on the other hand, suggest that α-helical content can be significant (10-25%). In this paper, we explore the possible reasons for this discrepancy and show that the different results from experiments and theory are not necessarily mutually exclusive but may reflect a general problem of secondary structure assignment of conformationally flexible biomolecules.

Backbone assignment and dynamics of human α-synuclein in viscous 2 M glucose solution

We present (1)H(N), (15)N, (13)C(α), (13)C(β) and (13)C' assignments and (15)N transverse relaxation rates (R(2)) of a Parkinson's disease-related intrinsically disordered protein, α-synuclein, in the presence of 2 M (360 g/l) glucose solution.

Advances in automated NMR protein structure determination

Around half of all protein structures solved nowadays using solution-state nuclear magnetic resonance (NMR) spectroscopy have been because of automated data analysis. The pervasiveness of computational approaches in general hides, however, a more nuanced view in which the full variety and richness of the field appears. This review is structured around a comparison of methods associated with three NMR observables: classical nuclear Overhauser effect (NOE) constraint gathering in contrast with more recent chemical shift and residual dipole coupling (RDC) based protocols. In each case, the emphasis is placed on the latest research, covering mainly the past 5 years. By describing both general concepts and representative programs, the objective is to map out a field in which--through the very profusion of approaches--it is all too easy to lose one's bearings.

Solid-state NMR characterization of gas vesicle structure

Gas vesicles are gas-filled buoyancy organelles with walls that consist almost exclusively of gas vesicle protein A (GvpA). Intact, collapsed gas vesicles from the cyanobacterium Anabaena flos-aquae were studied by solid-state NMR spectroscopy, and most of the GvpA sequence was assigned. Chemical shift analysis indicates a coil-α-β-β-α-coil peptide backbone, consistent with secondary-structure-prediction algorithms, and complementary information about mobility and solvent exposure yields a picture of the overall topology of the vesicle subunit that is consistent with its role in stabilizing an air-water interface.

Using NMR chemical shifts as structural restraints in molecular dynamics simulations of proteins

We introduce a procedure to determine the structures of proteins by incorporating NMR chemical shifts as structural restraints in molecular dynamics simulations. In this approach, the chemical shifts are expressed as differentiable functions of the atomic coordinates and used to compute forces to generate trajectories that lead to the reduction of the differences between experimental and calculated chemical shifts. We show that this strategy enables the folding of a set of proteins with representative topologies starting from partially denatured initial conformations without the use of additional experimental information. This method also enables the straightforward combination of chemical shifts with other standard NMR restraints, including those derived from NOE, J-coupling, and residual dipolar coupling measurements. We illustrate this aspect by calculating the structure of a transiently populated excited state conformation from chemical shift and residual dipolar coupling data measured by relaxation dispersion NMR experiments.

An introduction to NMR-based approaches for measuring protein dynamics

Proteins are inherently flexible at ambient temperature. At equilibrium, they are characterized by a set of conformations that undergo continuous exchange within a hierarchy of spatial and temporal scales ranging from nanometers to micrometers and femtoseconds to hours. Dynamic properties of proteins are essential for describing the structural bases of their biological functions including catalysis, binding, regulation and cellular structure. Nuclear magnetic resonance (NMR) spectroscopy represents a powerful technique for measuring these essential features of proteins. Here we provide an introduction to NMR-based approaches for studying protein dynamics, highlighting eight distinct methods with recent examples, contextualized within a common experimental and analytical framework. The selected methods are (1) Real-time NMR, (2) Exchange spectroscopy, (3) Lineshape analysis, (4) CPMG relaxation dispersion, (5) Rotating frame relaxation dispersion, (6) Nuclear spin relaxation, (7) Residual dipolar coupling, (8) Paramagnetic relaxation enhancement. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Validation of archived chemical shifts through atomic coordinates

The public archives containing protein information in the form of NMR chemical shift data at the BioMagResBank (BMRB) and of 3D structure coordinates at the Protein Data Bank are continuously expanding. The quality of the data contained in these archives, however, varies. The main issue for chemical shift values is that they are determined relative to a reference frequency. When this reference frequency is set incorrectly, all related chemical shift values are systematically offset. Such wrongly referenced chemical shift values, as well as other problems such as chemical shift values that are assigned to the wrong atom, are not easily distinguished from correct values and effectively reduce the usefulness of the archive. We describe a new method to correct and validate protein chemical shift values in relation to their 3D structure coordinates. This method classifies atoms using two parameters: the per-atom solvent accessible surface area (as calculated from the coordinates) and the secondary structure of the parent amino acid. Through the use of Gaussian statistics based on a large database of 3220 BMRB entries, we obtain per-entry chemical shift corrections as well as Z scores for the individual chemical shift values. In addition, information on the error of the correction value itself is available, and the method can retain only dependable correction values. We provide an online resource with chemical shift, atom exposure, and secondary structure information for all relevant BMRB entries (http://www.ebi.ac.uk/pdbe/nmr/vasco) and hope this data will aid the development of new chemical shift-based methods in NMR. Proteins 2010. © 2010 Wiley-Liss, Inc.

A combinatorial NMR and EPR approach for evaluating the structural ensemble of partially folded proteins

Partially folded proteins, characterized as exhibiting secondary structure elements with loose or absent tertiary contacts, represent important intermediates in both physiological protein folding and pathological protein misfolding. To aid in the characterization of the structural state(s) of such proteins, a novel structure calculation scheme is presented that combines structural restraints derived from pulsed EPR and NMR spectroscopy. The methodology is established for the protein alpha-synuclein (alphaS), which exhibits characteristics of a partially folded protein when bound to a micelle of the detergent sodium lauroyl sarcosinate (SLAS). By combining 18 EPR-derived interelectron spin label distance distributions with NMR-based secondary structure definitions and bond vector restraints, interelectron distances were correlated and a set of theoretical ensemble basis populations was calculated. A minimal set of basis structures, representing the partially folded state of SLAS-bound alphaS, was subsequently derived by back-calculating correlated distance distributions. A surprising variety of well-defined protein-micelle interactions was thus revealed in which the micelle is engulfed by two differently arranged antiparallel alphaS helices. The methodology further provided the population ratios between dominant ensemble structural states, whereas limitation in obtainable structural resolution arose from spin label flexibility and residual uncertainties in secondary structure definitions. To advance the understanding of protein-micelle interactions, the present study concludes by showing that, in marked contrast to secondary structure stability, helix dynamics of SLAS-bound alphaS correlate with the degree of protein-induced departures from free micelle dimensions.

Dynamic interactions between clathrin and locally structured elements in a disordered protein mediate clathrin lattice assembly

Assembly of clathrin lattices is mediated by assembly/adaptor proteins that contain domains that bind lipids or membrane-bound cargo proteins and clathrin binding domains (CBDs) that recruit clathrin. Here, we characterize the interaction between clathrin and a large fragment of the CBD of the clathrin assembly protein AP180. Mutational, NMR chemical shift, and analytical ultracentrifugation analyses allowed us to precisely define two clathrin binding sites within this fragment, each of which is found to bind weakly to the N-terminal domain of the clathrin heavy chain (TD). The locations of the two clathrin binding sites are consistent with predictions from sequence alignments of previously identified clathrin binding elements and, by extension, indicate that the complete AP180 CBD contains ∼12 degenerate repeats, each containing a single clathrin binding site. Sequence and circular dichroism analyses have indicated that the AP180 CBD is predominantly unstructured and our NMR analyses confirm that this is largely the case for the AP180 fragment characterized here. Unexpectedly, unlike the many proteins that undergo binding-coupled folding upon interaction with their binding partners, the AP180 fragment is similarly unstructured in its bound and free states. Instead, we find that this fragment exhibits localized β-turn-like structures at the two clathrin binding sites both when free and when bound to clathrin. These observations are incorporated into a model in which weak binding by multiple, pre-structured clathrin binding elements regularly dispersed throughout a largely unstructured CBD allows efficient recruitment of clathrin to endocytic sites and dynamic assembly of the clathrin lattice.

Solution structure and phospholipid interactions of the isolated voltage-sensor domain from KvAP

Voltage-sensor domains (VSDs) are specialized transmembrane segments that confer voltage sensitivity to many proteins such as ion channels and enzymes. The activities of these domains are highly dependent on both the chemical properties and the physical properties of the surrounding membrane environment. To learn about VSD-lipid interactions, we used nuclear magnetic resonance spectroscopy to determine the structure and phospholipid interface of the VSD from the voltage-dependent K(+) channel KvAP (prokaryotic Kv from Aeropyrum pernix). The solution structure of the KvAP VSD solubilized within phospholipid micelles is similar to a previously determined crystal structure solubilized by a nonionic detergent and complexed with an antibody fragment. The differences observed include a previously unidentified short amphipathic α-helix that precedes the first transmembrane helix and a subtle rigid-body repositioning of the S3-S4 voltage-sensor paddle. Using (15)N relaxation experiments, we show that much of the VSD, including the pronounced kink in S3 and the S3-S4 paddle, is relatively rigid on the picosecond-to-nanosecond timescale. In contrast, the kink in S3 is mobile on the microsecond-to-millisecond timescale and may act as a hinge in the movement of the paddle during channel gating. We characterized the VSD-phospholipid micelle interactions using nuclear Overhauser effect spectroscopy and showed that the micelle uniformly coats the KvAP VSD and approximates the chemical environment of a phospholipid bilayer. Using paramagnetically labeled phospholipids, we show that bilayer-forming lipids interact with the S3 and S4 helices more strongly than with S1 and S2.

A Transient and Low-Populated Protein-Folding Intermediate at Atomic Resolution

Proteins can sample conformational states that are critical for function but are seldom detected directly because of their low occupancies and short lifetimes. In this work, we used chemical shifts and bond-vector orientation constraints obtained from nuclear magnetic resonance relaxation dispersion spectroscopy, in concert with a chemical shift–based method for structure elucidation, to determine an atomic-resolution structure of an “invisible” folding intermediate of a small protein module: the FF domain. The structure reveals non-native elements preventing formation of the native conformation in the carboxyl-terminal part of the protein. This is consistent with the kinetics of folding in which a well-structured intermediate forms rapidly and then rearranges slowly to the native state. The approach introduces a general strategy for structure determination of low-populated and transiently formed protein states.

SPARTA+: a modest improvement in empirical NMR chemical shift prediction by means of an artificial neural network

NMR chemical shifts provide important local structural information for proteins and are key in recently described protein structure generation protocols. We describe a new chemical shift prediction program, SPARTA+, which is based on artificial neural networking. The neural network is trained on a large carefully pruned database, containing 580 proteins for which high-resolution X-ray structures and nearly complete backbone and (13)C(beta) chemical shifts are available. The neural network is trained to establish quantitative relations between chemical shifts and protein structures, including backbone and side-chain conformation, H-bonding, electric fields and ring-current effects. The trained neural network yields rapid chemical shift prediction for backbone and (13)C(beta) atoms, with standard deviations of 2.45, 1.09, 0.94, 1.14, 0.25 and 0.49 ppm for delta(15)N, delta(13)C', delta(13)C(alpha), delta(13)C(beta), delta(1)H(alpha) and delta(1)H(N), respectively, between the SPARTA+ predicted and experimental shifts for a set of eleven validation proteins. These results represent a modest but consistent improvement (2-10%) over the best programs available to date, and appear to be approaching the limit at which empirical approaches can predict chemical shifts.

NMR structural and dynamical investigation of the isolated voltage-sensing domain of the potassium channel KvAP: implications for voltage gating

The structure and dynamics of the isolated voltage-sensing domain (VSD) of the archaeal potassium channel KvAP was studied by high-resolution NMR. The almost complete backbone resonance assignment and partial side-chain assignment of the (2)H,(13)C,(15)N-labeled VSD were obtained for the protein domain solubilized in DPC/LDAO (2:1) mixed micelles. Secondary and tertiary structures of the VSD were characterized using secondary chemical shifts and NOE contacts. These data indicate that the spatial structure of the VSD solubilized in micelles corresponds to the structure of the domain in an open state of the channel. NOE contacts and secondary chemical shifts of amide protons indicate the presence of tightly bound water molecule as well as hydrogen bond formation involving an interhelical salt bridge (Asp62-R133) that stabilizes the overall structure of the domain. The backbone dynamics of the VSD was studied using (15)N relaxation measurements. The loop regions S1-S2 and S2-S3 were found mobile, while the S3-S4 loop (voltage-sensor paddle) was found stable at the ps-ns time scale. The moieties of S1, S2, S3, and S4 helices sharing interhelical contacts (at the level of the Asp62-R133 salt bridge) were observed in conformational exchange on the micros-ms time scale. Similar exchange-induced broadening of characteristic resonances was observed for the VSD solubilized in the membrane of lipid-protein nanodiscs composed of DMPC, DMPG, and POPC/DOPG lipids. Apparently, the observed interhelical motions represent an inherent property of the VSD of the KvAP channel and can play an important role in the voltage gating.

Enhanced conformational space sampling improves the prediction of chemical shifts in proteins

A biased-potential molecular dynamics simulation method, accelerated molecular dynamics (AMD), was combined with the chemical shift prediction algorithm SHIFTX to calculate ¹H^N, ¹⁵N, ¹³Cα, ¹³Cβ, and ¹³C′ chemical shifts of the ankyrin repeat protein IκBα (residues 67−206), the primary inhibitor of nuclear factor κ-B (NF-κB). Free-energy-weighted molecular ensembles were generated over a range of acceleration levels, affording systematic enhancement of the conformational space sampling of the protein. We have found that the predicted chemical shifts, particularly for the ¹⁵N, ¹³Cα, and ¹³Cβ nuclei, improve substantially with enhanced conformational space sampling up to an optimal acceleration level. Significant improvement in the predicted chemical shift data coincides with those regions of the protein that exhibit backbone dynamics on longer time scales. Interestingly, the optimal acceleration level for reproduction of the chemical shift data has previously been shown to best reproduce the experimental residual dipolar coupling (RDC) data for this system, as both chemical shift data and RDCs report on an ensemble and time average in the millisecond range.

The presence of a single N-terminal histidine residue enhances the fusogenic properties of a Membranotropic peptide derived from herpes simplex virus type 1 glycoprotein H

Herpes simplex virus type 1 (HSV-1)-induced membrane fusion remains one of the most elusive mechanisms to be deciphered in viral entry. The structure resolution of glycoprotein gB has revealed the presence of fusogenic domains in this protein and pointed out the key role of gB in the entry mechanism of HSV-1. A second putative fusogenic glycoprotein is represented by the heterodimer comprising the membrane-anchored glycoprotein H (gH) and the small secreted glycoprotein L, which remains on the viral envelope in virtue of its non-covalent interaction with gH. Different domains scattered on the ectodomain of HSV-1 gH have been demonstrated to display membranotropic characteristics. The segment from amino acid 626 to 644 represents the most fusogenic region identified by studies with synthetic peptides and model membranes. Herein we have identified the minimal fusogenic sequence present on gH. An enlongation at the N terminus of a single histidine (His) has proved to profoundly increase the fusogenic activity of the original sequence. Nuclear magnetic resonance (NMR) studies have shown that the addition of the N-terminal His contributes to the formation and stabilization of an alpha-helical domain with high fusion propensity.

Molecular determinants of the pKa values of Asp and Glu residues in staphylococcal nuclease

Prior computational studies of the acid-unfolding behavior of staphylococcal nuclease (SNase) suggest that the pK(a) values of its carboxylic groups are difficult to reproduce with electrostatics calculations with continuum methods. To examine the molecular determinants of the pK(a) values of carboxylic groups in SNase, the pK(a) values of all 20 Asp and Glu residues were measured with multidimensional and multinuclear NMR spectroscopy in an acid insensitive variant of SNase. The crystal structure of the protein was obtained to describe the microenvironments of the carboxylic groups. Fourteen Asp and Glu residues titrate with relatively normal pK(a) values that are depressed by less than 1.1 units relative to the normal pK(a) of Asp and Glu in water. Only six residues have pK(a) values shifted by more than 1.5 units. Asp-21 has an unusually high pK(a) of 6.5, which is probably the result of interactions with other carboxylic groups at the active site. The most perturbed pK(a) values appear to be governed by hydrogen bonding and not by Coulomb interactions. The pK(a) values calculated with standard continuum electrostatics methods applied to static structures are more depressed than the measured values because Coulomb effects are exaggerated in the calculations. The problems persist even when the protein is treated with the dielectric constant of water. This can be interpreted to imply that structural relaxation is an important determinant of the pK(a) values; however, no major pH-sensitive conformational reorganization of the backbone was detected using NMR spectroscopy.

DANGLE: A Bayesian inferential method for predicting protein backbone dihedral angles and secondary structure

This paper introduces DANGLE, a new algorithm that employs Bayesian inference to estimate the likelihood of all possible values of the backbone dihedral angles phi and psi for each residue in a query protein, based on observed chemical shifts and the conformational preferences of each amino acid type. The method provides robust estimates of phi and psi within realistic boundary ranges, an indication of the degeneracy in the relationship between shift measurements and conformation at each site, and faithful secondary structure state assignments. When a simple degeneracy-based filtering procedure is applied, DANGLE offers an ideal compromise between accuracy and coverage when compared with other shift-based dihedral angle prediction methods. In addition, per residue analysis of shift/structure degeneracy has potential to be a useful new approach for studying the properties of unfolded proteins, with sufficient sensitivity to identify regions of residual structure in the acid denatured state of apomyoglobin.

Evaluating the intrinsic cysteine redox-dependent states of the A-chain of human insulin using NMR spectroscopy, quantum chemical calculations, and mass spectrometry

Previous functional studies have proposed that solution-phase loading of human insulin A-chain peptides into cell surface Class II molecules may be limited by the redox state of intrinsic cysteine residues within the A-chain peptide. T cell functional studies of a human insulin A-chain analogue (KR A1-15) comprised of residues 1-15 of the A-chain peptide as well as an amino-terminal lysine-arginine extension have been carried out in a reducing environment. These data suggest that free thiol moieties within this peptide may participate in major histocompatibility complex (MHC) II/peptide interactions. Two-dimensional (1)H NMR spectroscopy data partnered with quantum chemical calculations identified that KR A1-15 exists in conformational flux sampling heterogeneous redox-dependent conformations including: one reduced and two oxidized states. These findings were further supported by mass spectrometry analysis of this peptide that confirmed the presence of a redox state dependent conformational equilibrium. Interestingly, the presence of a free thiol ((1)H(gamma)) resonance for cysteine 8 in the oxidized state supports the existence of the third redox-dependent conformation represented as a mixed disulfide conformation. We believe these data support the presence of a redox-dependent mechanism for regulating the activity of human insulin and provide a better understanding of redox chemistry that may be extended to other protein systems.

Automated protein NMR structure determination in solution

The main drawback of protein NMR spectroscopy today is still the extensive amount of time required for solving a single structure. The main bottleneck in this respect is the manual evaluation of the experimental spectra. A clear solution to this challenge is the development of automated methods for this purpose. At the current stage of development, this goal has been almost or in a few cases fully reached for favorable cases such as well-behaved, stably folding smaller proteins below the 25 kDa range. For larger and/or more difficult molecules, the input of a human expert is still required. However, even here, automated routines will substantially speed up the structure determination process. In this report, we will summarize recent developments in this field and especially emphasize practical aspects important for a successful automated protein structure determination in solution. An important aspect closely related to structure determination is structure validation. Therefore, we devote a section to automated approaches for this topic.

4D prediction of protein (1)H chemical shifts

A 4D approach for protein (1)H chemical shift prediction was explored. The 4th dimension is the molecular flexibility, mapped using molecular dynamics simulations. The chemical shifts were predicted with a principal component model based on atom coordinates from a database of 40 protein structures. When compared to the corresponding non-dynamic (3D) model, the 4th dimension improved prediction by 6-7%. The prediction method achieved RMS errors of 0.29 and 0.50 ppm for Halpha and HN shifts, respectively. However, for individual proteins the RMS errors were 0.17-0.34 and 0.34-0.65 ppm for the Halpha and HN shifts, respectively. X-ray structures gave better predictions than the corresponding NMR structures, indicating that chemical shifts contain invaluable information about local structures. The (1)H chemical shift prediction tool 4DSPOT is available from http://www.uku.fi/kemia/4dspot .

Isotope labeling methods for studies of excited protein states by relaxation dispersion NMR spectroscopy

The utility of nuclear magnetic resonance (NMR) spectroscopy as a tool for the study of biomolecular structure and dynamics has benefited from the development of facile labeling methods that incorporate NMR active probes at key positions in the molecule. Here we describe a protocol for the labeling of proteins that facilitates their study using a technique that is sensitive to millisecond conformational exchange processes. The samples necessary for an analysis of exchange dynamics are discussed, using the Abp1p SH3 domain from Saccharomyces cerevisiae as an example. For this system, the time frame for production of each sample, including in vitro refolding, is about 80 h. The samples so produced facilitate the measurement of accurate chemical shifts of low populated, invisible conformers that are part of the exchange pathway. The accuracy of the methodology has been established experimentally and the chemical shifts that are obtained provide important restraints in structure calculations of the excited state.

NMR-based characterization of a refolding intermediate of beta2-microglobulin labeled using a wheat germ cell-free system

In patients with dialysis-related amyloidosis, beta2-microglobulin (beta2-m) is a major structural component of amyloid fibrils. It has been suggested that the partial unfolding of beta2-m is a prerequisite to the formation of amyloid fibrils, and that the folding intermediate trapped by the non-native trans-Pro32 isomer leads to the formation of amyloid fibrils. Although clarifying the structure of this refolding intermediate by high resolution NMR spectroscopy is important, this has been made difficult by the limited lifetime of the intermediate. Here, we studied the structure of the refolding intermediate using a combination of amino acid selective labeling with wheat germ cell-free protein synthesis and NMR techniques. The HSQC spectra of beta2-ms labeled selectively at either phenylalanine, leucine, or valine enabled us to monitor the structures of the refolding intermediate. The results suggested that the refolding intermediate has an overall fold and cores similar to the native structure, but contains disordered structures around Pro32. The fluctuation of the beta-sheet regions especially the last half of the betaB strand and the first half of the betaE strand, both suggested to be important for amyloidogenicity, may transform beta2-m into an amyloidogenic structure.

Conformational averaging in structural biology: issues, challenges and computational solutions

Most experimental methods in structural biology provide time- and ensemble-averaged signals and, consequently, molecular structures based on such signals often exhibit only idealized, average features. Second, most experimental signals are only indirectly related to real, molecular geometries, and solving a structure typically involves a complicated procedure, which may not always result in a unique solution. To what extent do such conformationally-averaged, non-linear experimental signals and structural models derived from them accurately represent the underlying microscopic reality? Are there some structural motifs that are actually artificially more likely to be "seen" in an experiment simply due to the averaging artifact? Finally, what are the practical consequences of ignoring the averaging effects when it comes to functional and mechanistic implications that we try to glean from experimentally-based structural models? In this review, we critically address the work that has been aimed at studying such questions. We summarize the details of experimental methods typically used in structural biology (most notably nuclear magnetic resonance, X-ray crystallography and different types of spectroscopy), discuss their individual susceptibility to conformational (motional) averaging, and review several theoretical approaches, most importantly molecular dynamics simulations that are increasingly being used to aid experimentalists in interpreting structural biology experiments.

NMR structure of the transmembrane and cytoplasmic domains of human CD4 in micelles

The human cluster determinant 4 (CD4) is a type I transmembrane glycoprotein involved in T-cell signalling. It is expressed primarily on the surface of T helper cells but also on subsets of memory and regulatory T lymphocytes (CD4(+) cells). It serves as a coreceptor in T-cell receptor recognition of MHC II antigen complexes. Besides its cellular functions, CD4 serves as the main receptor for human immunodeficiency virus type I (HIV-1). During T-cell infection, the CD4 extracellular domain is bound by HIV-1 gp120, the viral surface glycoprotein, which triggers a number of conformational changes ultimately resulting in virion entry of the cell. Subsequently, CD4 is downregulated in infected cells by multiple strategies that involve direct interactions of the HIV-1 proteins VpU and Nef with the cytoplasmic part of CD4. In the present work, we describe the NOE-based solution structure of the transmembrane and cytoplasmic domains of the cystein-free variant of CD4 (CD4mut) in dodecylphosphocholine (DPC) micelles. Furthermore, we have characterized micelle-inserted CD4mut by paramagentic relaxation enhancement (PRE) agents and (1)H-(15)N heteronuclear NOE data. CD4mut features a stable and well-defined transmembrane helix from M372 to V395 buried in the micellar core and a cytoplasmic helix ranging from A404 to L413. Experimental data suggest the amphipathic cytoplasmic helix to be in close contact with the micellar surface. The role of the amphipathic helix and its interaction with the micellar surface is discussed with respect to the biological function of the full-length CD4 protein.

Pup, a prokaryotic ubiquitin-like protein, is an intrinsically disordered protein

Pup (prokaryotic ubiquitin-like protein) from Mycobacterium tuberculosis is the first ubiquitin-like protein identified in non-eukaryotic cells. Although different ubiquitin-like proteins from eukaryotes share low sequence similarity, their 3D (three-dimensional) structures exhibit highly conserved typical ubiquitin-like folds. Interestingly, our studies reveal that Pup not only shares low sequence similarity, but also presents a totally distinguished structure compared with other ubiquitin-like superfamily proteins. Diverse structure predictions combined with CD and NMR spectroscopic studies all demonstrate that Pup is an intrinsically disordered protein. Moreover, 1H-15N NOE (nuclear Overhauser effect) data and CSI (chemical shift index) analyses indicate that there is a residual secondary structure at the C-terminus of Pup. In M. tuberculosis, Mpa (mycobacterium proteasomal ATPase) is the regulatory cap ATPase of the proteasome that interacts with Pup and brings the substrates to the proteasome for degradation. In the present paper, SPR (surface plasmon resonance) and NMR perturbation studies imply that the C-terminus of Pup, ranging from residues 30 to 59, binds to Mpa probably through a hydrophobic interface. In addition, phylogenetic analysis clearly shows that the Pup family belongs to a unique and divergent evolutionary branch, suggesting that it is the most ancient and deeply branched family among ubiquitin-like proteins. This might explain the structural distinction between Pup and other ubiquitin-like superfamily proteins.

Quantitative determination of site-specific conformational distributions in an unfolded protein by solid-state nuclear magnetic resonance

Solid-state nuclear magnetic resonance (NMR) techniques are used to investigate the structure of the 35-residue villin headpiece subdomain (HP35) in folded, partially denatured, and fully denatured states. Experiments are carried out in frozen glycerol/water solutions, with chemical denaturation by guanidine hydrochloride (GdnHCl). Without GdnHCl, two-dimensional solid-state (13)C NMR spectra of samples prepared with uniform (13)C labeling of selected residues show relatively sharp cross-peaks at chemical shifts that are consistent with the known three-helix bundle structure of folded HP35. At high GdnHCl concentrations, most cross-peaks broaden and shift, qualitatively indicating disruption of the folded structure and development of static conformational disorder in the frozen denatured state. Conformational distributions at one residue in each helical segment are probed quantitatively with three solid-state NMR techniques that provide independent constraints on backbone varphi and psi torsion angles in samples with sequential pairs of carbonyl (13)C labels. Without GdnHCl, the combined data are well fit by alpha-helical conformations. At [GdnHCl]=4.5 M, corresponding to the approximate denaturation midpoint, the combined data are well fit by a combination of alpha-helical and partially extended conformations at each site, but with a site-dependent population ratio. At [GdnHCl]=7.0 M, corresponding to the fully denatured state, the combined data are well fit by a combination of partially extended and polyproline II conformations, again with a site-dependent population ratio. Two entirely different models for conformational distributions lead to nearly the same best-fit distributions, demonstrating the robustness of these conclusions. This work represents the first quantitative investigation of site-specific conformational distributions in partially folded and unfolded states of a protein by solid-state NMR.

Deuterium isotope effects on 15N backbone chemical shifts in proteins

Quantum mechanical calculations are presented that predict that one-bond deuterium isotope effects on the (15)N chemical shift of backbone amides of proteins, (1)Delta(15)N(D), are sensitive to backbone conformation and hydrogen bonding. A quantitative empirical model for (1)Delta(15)N(D) including the backbone dihedral angles, Phi and Psi, and the hydrogen bonding geometry is presented for glycine and amino acid residues with aliphatic side chains. The effect of hydrogen bonding is rationalized in part as an electric-field effect on the first derivative of the nuclear shielding with respect to N-H bond length. Another contributing factor is the effect of increased anharmonicity of the N-H stretching vibrational state upon hydrogen bonding, which results in an altered N-H/N-D equilibrium bond length ratio. The N-H stretching anharmonicity contribution falls off with the cosine of the N-H...O bond angle. For residues with uncharged side chains a very good prediction of isotope effects can be made. Thus, for proteins with known secondary structures, (1)Delta(15)N(D) can provide insights into hydrogen bonding geometries.