The folding pathway of T4 lysozyme: an on-pathway hidden folding intermediate

T4 lysozyme has two easily distinguishable but energetically coupled domains: the N and C-terminal domains. In earlier studies, an amide hydrogen/deuterium exchange pulse-labeling experiment detected a stable submillisecond intermediate that accumulates before the rate-limiting transition state. It involves the formation of structures in both the N and C-terminal regions. However, a native-state hydrogen exchange experiment subsequently detected an equilibrium intermediate that only involves the formation of the C-terminal domain. Here, using stopped-flow circular dichroism and fluorescence, amide hydrogen exchange-folding competition, and protein engineering methods, we re-examined the folding pathway of T4-lysozyme. We found no evidence for the existence of a stable folding intermediate before the rate-limiting transition state at neutral pH. In addition, using native-state hydrogen exchange-directed protein engineering, we created a mimic of the equilibrium intermediate. We found that the intermediate mimic folds with the same rate as the wild-type protein, suggesting that the equilibrium intermediate is an on-pathway intermediate that exists after the rate-limiting transition state.     

Identification of a collapsed intermediate with non-native long-range interactions on the folding pathway of a pair of Fyn SH3 domain mutants by NMR relaxation dispersion spectroscopy

Recent 15N and 13C spin-relaxation dispersion studies of fast-folding mutants of the Fyn SH3 domain have established that folding proceeds through a low-populated on-pathway intermediate (I) where the central beta-sheet is at least partially formed, but without interactions between the NH2- and COOH-terminal beta-strands that exist in the folded state (F). Initial studies focused on mutants where Gly48 is replaced; in an effort to establish whether this intermediate is a general feature of Fyn SH3 folding a series of 15N relaxation experiments monitoring the folding of Fyn SH3 mutants N53P/V55L and A39V/N53P/V55L are reported here. For these mutants as well, folding proceeds through an on-pathway intermediate with similar features to those observed for G48M and G48V Fyn SH3 domains. However, the 15N chemical shifts extracted for the intermediate indicate pronounced non-native contacts between the NH2 and COOH-terminal regions not observed previously. The kinetic parameters extracted for the folding of A39V/N53P/V55L Fyn SH3 from the three-state folding model F<-->I<-->U are in good agreement with folding and unfolding rates extrapolated to zero denaturant obtained from stopped-flow experiments analyzed in terms of a simplified two-state folding reaction. The folding of the triple mutant was studied over a wide range of temperatures, establishing that there is no difference in heat capacities between F and I states. This confirms a compact folding intermediate structure, which is supported by the 15N chemical shifts of the I state extracted from the dispersion data. The temperature-dependent relaxation data simplifies data analysis because at low temperatures (< 25 degrees C) the unfolded state (U) is negligibly populated relative to I and F. A comparison between parameters extracted at low temperatures where the F<-->I exchange model is appropriate with those from the more complex, three-state model at higher temperatures has been used to validate the protocol for analysis of three-site exchange relaxation data.     

(1)H/(15)N heteronuclear NMR spectroscopy shows four dynamic domains for phospholamban reconstituted in dodecylphosphocholine micelles

We report the backbone dynamics of monomeric phospholamban in dodecylphosphocholine micelles using (1)H/(15)N heteronuclear NMR spectroscopy. Phospholamban is a 52-amino acid membrane protein that regulates Ca-ATPase in cardiac muscle. Phospholamban comprises three structural domains: a transmembrane domain from residues 22 to 52, a connecting loop from 17 to 21, and a cytoplasmic domain from 1 to 16 that is organized in an "L"-shaped structure where the transmembrane and the cytoplasmic domain form an angle of approximately 80 degrees (Zamoon et al., 2003; Mascioni et al., 2002). T(1), T(2), and (1)H/(15)N nuclear Overhauser effect values measured for the amide backbone resonances were interpreted using the model-free approach of Lipari and Szabo. The results point to the existence of four dynamic domains, revealing the overall plasticity of the cytoplasmic helix, the flexible loop, and part of the transmembrane domain (residues 22-30). In addition, using Carr-Purcell-Meiboom-Gill-based experiments, we have characterized phospholamban dynamics in the micros-ms timescale. We found that the majority of the residues in the cytoplasmic domain, the flexible loop, and the first ten residues of the transmembrane domain undergo dynamics in the micros-ms range, whereas minimal dynamics were detected for the transmembrane domain. Hydrogen/deuterium exchange factors measured at different temperatures support the existence of slow motion in both the loop and the cytoplasmic helix. We propose that these dynamic properties are critical factors in the biomolecular recognition of phospholamban by Ca-ATPase and other interacting proteins such as protein kinase A and protein phosphatase 1.     

Probing the flexibility of the DsbA oxidoreductase from Vibrio cholerae--a 15N - 1H heteronuclear NMR relaxation analysis of oxidized and reduced forms of DsbA

We have determined the structure of the reduced form of the DsbA oxidoreductase from Vibrio cholerae. The reduced structure shows a high level of similarity to the crystal structure of the oxidized form and is typical of this class of enzyme containing a thioredoxin domain with an inserted alpha-helical domain. Proteolytic and thermal stability measurements show that the reduced form of DsbA is considerably more stable than the oxidized form. NMR relaxation data have been collected and analyzed using a model-free approach to probe the dynamics of the reduced and oxidized states of DsbA. Akaike's information criteria have been applied both in the selection of the model-free models and the diffusion tensors that describe the global motions of each redox form. Analysis of the dynamics reveals that the oxidized protein shows increased disorder on the pico- to nanosecond and micro- to millisecond timescale. Many significant changes in dynamics are located either close to the active site or at the insertion points between the domains. In addition, analysis of the diffusion data shows there is a clear difference in the degree of interdomain movement between oxidized and reduced DsbA with the oxidized form being the more rigid. Principal components analysis has been employed to indicate possible concerted movements in the DsbA structure, which suggests that the modeled interdomain motions affect the catalytic cleft of the enzyme. Taken together, these data provide compelling evidence of a role for dynamics in the catalytic cycle of DsbA.     

Analysis of internal motions of interleukin-13 variant associated with severe bronchial asthma using (15)N NMR relaxation measurements

The single nucleotide polymorphism interleukin-13 (IL-13) R110Q is associated with severe bronchial asthma because its lower affinity leads to the augmentation of local IL-13 concentration, resulting in an increase in the signal transduction via IL-13R. Since the mutation site does not directly bind to IL-13Ralpha2, we carried out NMR relaxation analyses of the wild-type IL-13 and IL-13-R110Q in order to examine whether the R110Q mutation affects the internal motions in IL-13 molecules. The results showed that the internal motion in the micro- to millisecond time scale on helix D, which is suggested to be important for the interaction between IL-13 and IL-13Ralpha2, is increased in IL-13-R110Q compared with that in the wild-type IL-13. It therefore appears that the difference in the internal motions on helix D between the wild-type IL-13 and IL-13-R110Q may be involved in their affinity differences with IL-13Ralpha2.     

1H(C) and 1H(N) total NOE correlations in a single 3D NMR experiment. 15N and 13C time-sharing in t1 and t2 dimensions for simultaneous data acquisition

Simultaneous data acquisition in time-sharing (TS) multi-dimensional NMR experiments has been shown an effective means to reduce experimental time, and thus to accelerate structure determination of proteins. This has been accomplished by spin evolution time-sharing of the X and Y heteronuclei, such as ¹⁵N and ¹³C, in one of the time dimensions. In this work, we report a new 3D TS experiment, which allows simultaneous ¹³C and ¹⁵N spin labeling coherence in both t ₁ and t ₂ dimensions to give four NOESY spectra in a single 3D experiment. These spectra represent total NOE correlations between ¹H_N and ¹H_C resonances. This strategy of double time-sharing (2TS) results in an overall four-fold reduction in experimental time compared with its conventional counterpart. This 3D 2TS CN-CN-H HSQC-NOESY-HSQC pulse sequence also demonstrates improvements in water suppression, ¹⁵N spectral resolution and sensitivity, which were developed based on 2D TS CN-H HSQC and 3D TS H-CN-H NOESY-HSQC experiments. Combining the 3D TS and the 3D 2TS NOESY experiments, NOE assignment ambiguities and errors are considerably reduced. These results will be useful for rapid protein structure determination to complement the effort of discerning the functions of diverse genomic proteins.     

Improved sensitivity in indirect monitoring of chemical shifts of proton-heteronuclear spin pairs (1H-13C and 1H-15N) in 3D and 4D NMR spectroscopy

In three-dimensional and four-dimensional experiments on doubly labelled proteins not only heteronuclear (¹³C or ¹⁵N) but also proton (¹H) frequencies are often indirectly monitored, rather than being directly observed. In this communication we show how in these experiments by overlaying ¹H and heteronuclear evolutions one can obtain decreased apparent relaxation rates of ¹H signals, yielding improved sensitivity. The new method applies to spin pairs like ¹H-¹⁵N, as in amide groups, or ¹H-¹³C, as in methine groups of alpha or aromatic systems.     

Architecture of the hepatitis C virus E1 glycoprotein transmembrane domain studied by NMR

Oligomerization of hepatitis C viral envelope proteins E1 and E2 is essential to virus fusion and assembly. Although interactions within the transmembrane (TM) domains of these glycoproteins have proven contributions to the E1/E2 heterodimerization process and consequent infectivity, there is little structural information on this entry mechanism. Here, as a first step towards our long-term goal of understanding the interaction between E1 and E2 TM-domains, we have expressed, purified and characterized E1-TM using structural biomolecular NMR methods. An MBP-fusion expression system yielded sufficient quantities of pure E1-TM, which was solubilized in two membrane-mimicking environments, SDS- and LPPG-micelles, affording samples amenable to NMR studies. Triple resonance assignment experiments and relaxation measurements provided information on the secondary structure and global fold of E1-TM in these environments. In SDS micelles E1-TM adopts a helical conformation, with helical stretches at residues 354–363 and 371–379 separated by a more flexible segment of residues 364–370. In LPPG micelles a helical conformation was observed for residues 354–377 with greater flexibility in the 366–367 dyad, suggesting LPPG provides a more native environment for the peptide. Replacement of key positively charged residue K370 with an alanine did not affect the secondary structure of E1-TM but did change the relative positioning within the micelle of the two helices. These results lay the foundation for structure determination of E1-TM and a molecular understanding of how E1-TM flexibility enhances its interaction with E2-TM during heterodimerization and membrane fusion.     

Solution structure of the MID1 B-box2 CHC(D/C)C(2)H(2) zinc-binding domain: insights into an evolutionarily conserved RING fold

The B-box type 2 domain is a prominent feature of a large and growing family of RING, B-box, coiled-coil (RBCC) domain-containing proteins and is also present in more than 1500 additional proteins. Most proteins usually contain a single B-box2 domain, although some proteins contain tandem domains consisting of both type 1 and type 2 B-boxes, which actually share little sequence similarity. Recently, we determined the solution structure of B-box1 from MID1, a putative E3 ubiquitin ligase that is mutated in X-linked Opitz G/BBB syndrome, and showed that it adopted a betabetaalpha RING-like fold. Here, we report the tertiary structure of the B-box2 (CHC(D/C)C(2)H(2)) domain from MID1 using multidimensional NMR spectroscopy. This MID1 B-box2 domain consists of a short alpha-helix and a structured loop with two short anti-parallel beta-strands and adopts a tertiary structure similar to the B-box1 and RING structures, even though there is minimal primary sequence similarity between these domains. By mutagenesis, ESI-FTICR and ICP mass spectrometry, we show that the B-box2 domain coordinates two zinc atoms with a 'cross-brace' pattern: one by Cys175, His178, Cys195 and Cys198 and the other by Cys187, Asp190, His204, and His207. Interestingly, this is the first case that an aspartic acid is involved in zinc atom coordination in a zinc-finger domain, although aspartic acid has been shown to coordinate non-catalytic zinc in matrix metalloproteinases. In addition, the finding of a Cys195Phe substitution identified in a patient with X-linked Opitz GBBB syndrome supports the importance of proper zinc coordination for the function of the MID1 B-box2 domain. Notably, however, our structure differs from the only other published B-box2 structure, that from XNF7, which was shown to coordinate one zinc atom. Finally, the similarity in tertiary structures of the B-box2, B-box1 and RING domains suggests these domains have evolved from a common ancestor.     

Biophysical characterization of structural properties and folding of interleukin-1 receptor antagonist

Structural properties and folding of interleukin-1 receptor antagonist (IL-1ra), a therapeutically important cytokine with a symmetric beta-trefoil topology, are characterized using optical spectroscopy, high-resolution NMR, and size-exclusion chromatography. Spectral contributions of two tryptophan residues, Trp17 and Trp120, present in the wild-type protein, have been determined from mutational analysis. Trp17 dominates the emission spectrum of IL-1ra, while Trp120 is quenched presumably by the nearby cysteine residues in both folded and unfolded states. The same Trp17 gives rise to two characteristic negative peaks in the aromatic CD. Urea denaturation of the wild-type protein is probed by measuring intrinsic and extrinsic (binding of 1-anilinonaphthalene-8-sulfonic acid) fluorescence, near- and far-UV CD, and 1D and 2D ((1)H-(15)N heteronuclear single quantum coherence (HSQC)) NMR. Overall, the data suggest an essentially two-state equilibrium denaturation mechanism with small, but detectable structural changes within the pretransition region. The majority of the (1)H-(15)N HSQC cross-peaks of the folded state show only a limited chemical shift change as a function of the denaturant concentration. However, the amide cross-peak of Leu31 demonstrates a significant urea dependence that can be fitted to a two-state binding model with a dissociation constant of 0.95+/-0.04 M. This interaction has at least a five times higher affinity than reported values for nonspecific urea binding to denatured proteins and peptides, suggesting that the structural context around Leu31 stabilizes the protein-urea interaction. A possible role of denaturant binding in inducing the pretransition changes in IL-1ra is discussed. Urea unfolding of wild-type IL-1ra is sufficiently slow to enable HPLC separation of folded and unfolded states. Quantitative size-exclusion chromatography has provided a hydrodynamic view of the kinetic denaturation process. Thermodynamic stability and unfolding kinetics of IL-1ra resemble those of structurally and evolutionary close IL-1beta, suggesting similarity of their free energy landscapes.     

Global fold and backbone dynamics of the hepatitis C virus E2 glycoprotein transmembrane domain determined by NMR

E1 and E2 are two hepatitis C viral envelope glycoproteins that assemble into a heterodimer that is essential for membrane fusion and penetration into the target cell. Both extracellular and transmembrane (TM) glycoprotein domains contribute to this interaction, but study of TM–TM interactions has been limited because synthesis and structural characterization of these highly hydrophobic segments present significant challenges. In this NMR study, by successful expression and purification of the E2 transmembrane domain as a fusion construct we have determined the global fold and characterized backbone motions for this peptide incorporated in phospholipid micelles. Backbone resonance frequencies, relaxation rates and solvent exposure measurements concur in showing this domain to adopt a helical conformation, with two helical segments spanning residues 717–726 and 732–746 connected by an unstructured linker containing the charged residues D728 and R730 involved in E1 binding. Although this linker exhibits increased local motions on the ps timescale, the dominating contribution to its relaxation is the global tumbling motion with an estimated correlation time of 12.3 ns. The positioning of the helix–linker–helix architecture within the mixed micelle was established by paramagnetic NMR spectroscopy and phospholipid-peptide cross relaxation measurements. These indicate that while the helices traverse the hydrophobic interior of the micelle, the linker lies closer to the micelle perimeter to accommodate its charged residues. These results lay the groundwork for structure determination of the E1/E2 complex and a molecular understanding of glycoprotein heterodimerization.