Exploring the folding free energy surface of a three-helix bundle protein

The multidimensional free energy surface for a small fast folding helical protein is explored based on first-principle calculations. The model represents the 46-residue segment from fragment B of staphylococcal protein A. The relationship between collapse and tertiary structure formation, and the order of collapse and secondary structure formation, are investigated. We find that the initial collapse process gives rise to a transition state with about 30% of the native tertiary structure and 50–70% of the native helix content. We also observe two distinct distributions of native helix in this collapsed state (Rg ≈ 12 Å), one with about 20% of the native helical hydrogen bonds, the other with near 70%. The former corresponds to a local minimum. The barrier from this metastable state to the native state is about 2 kBT. In the latter case, folding is essentially a downhill process involving topological assembly. In addition, the order of formation of secondary structure among the three helices is examined. We observe cooperative formation of the secondary structure in helix I and helix II. Secondary structure in helix III starts to form following the formation of certain secondary structure in both helix I and helix II. Comparisons of our results with those from theory and experiment are made.

Folding thermodynamics of a model three-helix-bundle protein

The calculated folding thermodynamics of a simple off-lattice three-helix-bundle protein model under equilibrium conditions shows the experimentally observed protein transitions: a collapse transition, a disordered-to-ordered globule transition, a globule to native-state transition, and the transition from the active native state to a frozen inactive state. The cooperativity and physical origin of the various transitions are explored with a single “optimization” parameter and characterized with the Lindemann criterion for liquid versus solid-state dynamics. Below the folding temperature, the model has a simple free energy surface with a single basin near the native state; the surface is similar to that calculated from a simulation of the same three-helix-bundle protein with an all-atom representation [Boczko, E. M. & Brooks III, C. L. (1995) Science 269, 393–396].

Cytochrome c′ folding triggered by electron transfer: Fast and slow formation of four-helix bundles

Reduced (FeII) Rhodopseudomonas palustris cytochrome c′ (Cyt c′) is more stable toward unfolding ([GuHCl]1/2 = 2.9(1) M) than the oxidized (FeIII) protein ([GuHCl]1/2 = 1.9(1) M). The difference in folding free energies (ΔΔGf° = 70 meV) is less than half of the difference in reduction potentials of the folded protein (100 mV vs. NHE) and a free heme in aqueous solution (≈−150 mV). The spectroscopic features of unfolded FeII–Cyt c′ indicate a low-spin heme that is axially coordinated to methionine sulfur (Met-15 or Met-25). Time-resolved absorption measurements after CO photodissociation from unfolded FeII(CO)–Cyt c′ confirm that methionine can bind to the ferroheme on the microsecond time scale [kobs = 5(2) × 104 s−1]. Protein folding was initiated by photoreduction (two-photon laser excitation of NADH) of unfolded FeIII–Cyt c′ ([GuHCl] = 2.02–2.54 M). Folding kinetics monitored by heme absorption span a wide time range and are highly heterogeneous; there are fast-folding (≈103 s−1), intermediate-folding (102–101 s−1), and slow-folding (10−1 s−1) populations, with the last two likely containing methionine-ligated (Met-15 or Met-25) ferrohemes. Kinetics after photoreduction of unfolded FeIII–Cyt c′ in the presence of CO are attributable to CO binding [1.4(6) × 103 s−1] and FeII(CO)–Cyt c′ folding [2.8(9) s−1] processes; stopped-flow triggered folding of FeIII–Cyt c′ (which does not contain a protein-derived sixth ligand) is adequately described by a single kinetics phase with an estimated folding time constant of ≈4 ms [ΔGf° = −33(3) kJ mol−1] at zero denaturant.

Brownian dynamics simulations of protein folding: Access to milliseconds time scale and beyond

Protein folding occurs on a time scale ranging from milliseconds to minutes for a majority of proteins. Computer simulation of protein folding, from a random configuration to the native structure, is nontrivial owing to the large disparity between the simulation and folding time scales. As an effort to overcome this limitation, simple models with idealized protein subdomains, e.g., the diffusion–collision model of Karplus and Weaver, have gained some popularity. We present here new results for the folding of a four-helix bundle within the framework of the diffusion–collision model. Even with such simplifying assumptions, a direct application of standard Brownian dynamics methods would consume 10,000 processor-years on current supercomputers. We circumvent this difficulty by invoking a special Brownian dynamics simulation. The method features the calculation of the mean passage time of an event from the flux overpopulation method and the sampling of events that lead to productive collisions even if their probability is extremely small (because of large free-energy barriers that separate them from the higher probability events). Using these developments, we demonstrate that a coarse-grained model of the four-helix bundle can be simulated in several days on current supercomputers. Furthermore, such simulations yield folding times that are in the range of time scales observed in experiments.

Channel formation by antiapoptotic protein Bcl-2

Bcl-2 is the prototypical member of a large family of apoptosis-regulating proteins, consisting of blockers and promoters of cell death. The three-dimensional structure of a Bcl-2 homologue, Bcl-XL, suggests striking similarity to the pore-forming domains of diphtheria toxin and the bacterial colicins, prompting exploration of whether Bcl-2 is capable of forming pores in lipid membranes. Using chloride efflux from KCl-loaded unilamellar lipid vesicles as an assay, purified recombinant Bcl-2 protein exhibited pore-forming activity with properties similar to those of the bacterial toxins, diphtheria toxin, and colicins, i.e., dependence on low pH and acidic lipid membranes. In contrast, a mutant of Bcl-2 lacking the two core hydrophobic α-helices (helices 5 and 6), predicted to be required for membrane insertion and channel formation, produced only nonspecific effects. In planar lipid bilayers, where detection of single channels is possible, Bcl-2 formed discrete ion-conducting, cation-selective channels, whereas the Bcl-2 (Δh5, 6) mutant did not. The most frequent conductance observed (18 ± 2 pS in 0.5 M KCl at pH 7.4) is consistent with a four-helix bundle structure arising from Bcl-2 dimers. However, larger channel conductances (41 ± 2 pS and 90 ± 10 pS) also were detected with progressively lower occurrence, implying the step-wise formation of larger oligomers of Bcl-2 in membranes. These findings thus provide biophysical evidence that Bcl-2 forms channels in lipid membranes, suggesting a novel function for this antiapoptotic protein.

Crystal structure of truncated human apolipoprotein A-I suggests a lipid-bound conformation

The structure of truncated human apolipoprotein A-I (apo A-I), the major protein component of high density lipoprotein, has been determined at 4-Å resolution. The crystals comprise residues 44–243 (exon 4) of apo A-I, a fragment that binds to lipid similarly to intact apo A-I and that retains the lipid-bound conformation even in the absence of lipid. The molecule consists almost entirely of a pseudo-continuous, amphipathic α-helix that is punctuated by kinks at regularly spaced proline residues; it adopts a shape similar to a horseshoe of dimensions 125 × 80 × 40 Å. Four molecules in the asymmetric unit associate via their hydrophobic faces to form an antiparallel four-helix bundle with an elliptical ring shape. Based on this structure, we propose a model for the structure of apo A-I bound to high density lipoprotein.

Conserved structural features of the synaptic fusion complex: SNARE proteins reclassified as Q- and R-SNAREs

SNARE [soluble NSF (N-ethylmaleimide-sensitive fusion protein) attachment protein receptor] proteins are essential for membrane fusion and are conserved from yeast to humans. Sequence alignments of the most conserved regions were mapped onto the recently solved crystal structure of the heterotrimeric synaptic fusion complex. The association of the four α-helices in the synaptic fusion complex structure produces highly conserved layers of interacting amino acid side chains in the center of the four-helix bundle. Mutations in these layers reduce complex stability and cause defects in membrane traffic even in distantly related SNAREs. When syntaxin-4 is modeled into the synaptic fusion complex as a replacement of syntaxin-1A, no major steric clashes arise and the most variable amino acids localize to the outer surface of the complex. We conclude that the main structural features of the neuronal complex are highly conserved during evolution. On the basis of these features we have reclassified SNARE proteins into Q-SNAREs and R-SNAREs, and we propose that fusion-competent SNARE complexes generally consist of four-helix bundles composed of three Q-SNAREs and one R-SNARE.

Comparative assessment of virulence of recombinant vaccinia viruses expressing IL-2 and IL-15 in immunodeficient mice

IL-2 and -15 belong to the four α-helix bundle family of cytokines and display a spectrum of overlapping immune functions because of shared signal transducing receptor components of the IL-2 receptor complex. However, recent evidence suggests a nonredundant unique role for IL-15 in the establishment and perhaps maintenance of peripheral natural killer (NK) cell populations in vivo. To explore the contribution of locally released IL-15 on peripheral NK-cell-mediated innate immune responses, we generated a recombinant vaccinia virus that expresses IL-15 and evaluated the course of vaccinial disease in athymic nude mice. Coexpression of IL-15 resulted in the attenuation of virulence of vaccinia virus, and mice inoculated with 105 plaque-forming units or less resolved the infection successfully. In contrast, mice inoculated with a similar dose of the control vaccinia virus failed to eliminate the virus and died of generalized vaccinial disease. Enhanced expression of IL-12 and IFN-γ as well as induction of chemokines were evident in the mice inoculated with IL-15-expressing vaccinia virus in addition to an increase in NK cells in the spleen. However, in this model system, the degree of attenuation in viral virulence attained with coexpression of IL-15 was much less than that achieved with coexpression of IL-2, suggesting that the peripheral NK-cell-mediated events are more responsive to IL-2 than to IL-15.

The crystal structure of the immunity protein of colicin E7 suggests a possible colicin-interacting surface.

The immunity protein of colicin E7 (ImmE7) can bind specifically to the DNase-type colicin E7 and inhibit its bactericidal activity. Here we report the 1.8-angstrom crystal structure of the ImmE7 protein. This is the first x-ray structure determined in the superfamily of colicin immunity proteins. The ImmE7 protein consists of four antiparallel alpha-helices, folded in a topology similar to the architecture of a four-helix bundle structure. A region rich in acidic residues is identified. This negatively charged area has the greatest variability within the family of DNase-type immunity proteins; thus, it seems likely that this area is involved in specific binding to colicin. Based on structural, genetic, and kinetic data, we suggest that all the DNase-type immunity proteins, as well as colicins, share a "homologous-structural framework" and that specific interaction between a colicin and its cognate immunity protein relies upon how well these two proteins' charged residues match on the interaction surface, thus leading to specific immunity of the colicin.

Alternate protein frameworks for molecular recognition.

In an effort to determine whether proteins with structures other than the immunoglobulin fold can be used to mimic the ligand binding properties of antibodies, we generated a library from the four-helix bundle protein cytochrome b562 in which the two loops were randomized. Panning of this library against the bovine serum albumin (BSA) conjugate of N-methyl-p-nitrobenzylamine derivative 1 by phage display methods yielded cytochromes in which residues Trp-20, Arg-21, and Ser-22 in loop A and Arg-83 and Trp-84 in loop B were conserved. The individual mutants, which fold into native-like structure, bind selectively to the BSA-1 conjugate with micromolar dissociation constants (Kd), in comparison to a monoclonal antibody that binds selectively to 1 with a Kd of 290 nM. These and other antibody-like receptors may prove useful as therapeutic agents or as reagents for both intra- and extracellular studies.

Structure of the Ets-1 pointed domain and mitogen-activated protein kinase phosphorylation site

The Pointed (PNT) domain and an adjacent mitogen-activated protein (MAP) kinase phosphorylation site are defined by sequence conservation among a subset of ets transcription factors and are implicated in two regulatory strategies, protein interactions and posttranslational modifications, respectively. By using NMR, we have determined the structure of a 110-residue fragment of murine Ets-1 that includes the PNT domain and MAP kinase site. The Ets-1 PNT domain forms a monomeric five-helix bundle. The architecture is distinct from that of any known DNA- or protein-binding module, including the helix-loop-helix fold proposed for the PNT domain of the ets protein TEL. The MAP kinase site is in a highly flexible region of both the unphosphorylated and phosphorylated forms of the Ets-1 fragment. Phosphorylation alters neither the structure nor monomeric state of the PNT domain. These results suggest that the Ets-1 PNT domain functions in heterotypic protein interactions and support the possibility that target recognition is coupled to structuring of the MAP kinase site.

Structure and function in rhodopsin: Rhodopsin mutants with a neutral amino acid at E134 have a partially activated conformation in the dark state*

The Glu-134–Arg-135 residues in rhodopsin, located near the cytoplasmic end of the C helix, are involved in G protein binding, or activation, or both. Furthermore, the charge-neutralizing mutation Glu-134 to Gln-134 produces hyperactivity in the activated state and produces constitutive activity in opsin. The Glu/Asp-Arg charge pair is highly conserved in equivalent positions in other G protein-coupled receptors. To investigate the structural consequences of charge-neutralizing mutations at Glu-134 and Arg-135 in rhodopsin, single spin-labeled side chains were introduced at sites in the cytoplasmic domains of helices C (140), E (227), F (250), or G (316) to serve as “molecular sensors” of the local helix bundle conformation. In each of the spin-labeled rhodopsins, a Gln substitution was introduced at either Glu-134 or Arg-135, and the electron paramagnetic resonance spectrum of the spin label was used to monitor the structural response of the helix bundle. The results indicate that a Gln substitution at Glu-134 induces a photoactivated conformation around helices C and G even in the dark state, an observation of potential relevance to the hyperactivity and constitutive activity of the mutant. In contrast, little change is induced in helix F, which has been shown to undergo a dominant motion upon photoactivation. This result implies that the multiple helix motions accompanying photoactivation are not strongly coupled and can be induced to take place independently. Gln substitution at Arg-135 produces only minor structural changes in the dark- or light-activated conformation, suggesting that this residue is not a determinant of structure in the regions investigated, although it may be functionally important.

Ferredoxin and ferredoxin–heme maquettes

A 16-amino acid residue peptide derived from a consensus motif of natural ferredoxins incorporates a tetranuclear iron sulfur cluster under physiological conditions. Successful assembly of the [4Fe–4S]2+/1+ cluster within a monomeric peptide was demonstrated using size exclusion chromatography, UV-visible, visible CD, and cryogenic EPR spectroscopies. The robustness of [4Fe–4S]2+/1+ formation was tested using peptides with either the ligating cysteine exchanged for alanine or with the intervening amino acids replaced by glycine. The small size of the peptide allows for modular incorporation into more complex protein structures. In one larger structure, we describe a tetra-α-helix bundle that self-assembles both iron–sulfur clusters and hemes, thereby demonstrating feasibility for the general synthesis of maquettes containing multiple, juxtaposed redox cofactors. This is a motif common to the catalytic sites of native oxidoreductases.

Staphylococcal protein A: Unfolding pathways, unfolded states, and differences between the B and E domains

We present multiple native and denaturation simulations of the B and E domains of the three-helix bundle protein A, totaling 60 ns. The C-terminal helix (H3) consistently denatures later than either of the other two helices and contains residual helical structure in the denatured state. These results are consistent with experiments suggesting that the isolated H3 fragment is more stable than H1 and H2 and that H3 forms early in folding. Interestingly, the denatured state of the B domain is much more compact than that of the E domain. This sequence-dependent effect on the dimensions of the denatured state and the lack of correlation with structure suggest that the radius of gyration can be a misleading reaction coordinate for unfolding/folding. Various unfolding and refolding events are observed in the denaturation simulations. In some cases, the transitions are facilitated through interactions with other portions of the protein—contact-assisted helix formation. In the native simulations, the E domain is very stable: after 6 ns, the Cα root-mean-square deviation from the starting structure is less than 1.4 Å. In contrast, the native state of the B domain deviates more and its inter-helical angles fluctuate. In apparent contrast, we note that the B domain is thermodynamically more stable than the E domain. The simulations suggest that the increased stability of the B domain may be due to heightened mobility, and therefore entropy, in the native state and decreased mobility/entropy in the more compact denatured state.

Detecting the conformational change of transmembrane signaling in a bacterial chemoreceptor by measuring effects on disulfide cross-linking in vivo.

Transmembrane signaling by bacterial chemoreceptors is thought to involve relative movement among the four transmembrane helices of the homodimer. We assayed that movement by measuring effects of ligand occupancy on rates of oxidative cross-linking between cysteines introduced into neighboring helices of the transmembrane domain of chemoreceptor Trg from Escherichia coli. Measurements were done on chemoreceptors in their native environment, intact cells that were motile and chemotactically responsive. Receptor occupancy did not appear to cause drastic rearrangement of the four-helix structure since, among 67 cysteine pairs tested, the same 19 exhibited oxidative cross-linking in the presence or absence of saturating chemoattractant. However, occupancy did cause subtle changes that were detected as effects on rates of cross-linking. Among the seven disulfides appropriate for measurements of initial rates of formation, ligand occupancy had significant and different effects on all three cross-links that connected the two helices within a subunit but had minimal effects on the four that spanned the packing interface between subunits. This constitutes direct evidence that the conformational change of transmembrane signaling involves significant movement within a subunit and minimal movement between subunits, a pattern deduced from several previous studies and now documented directly. Among possible modes of movement between the two helices of a subunit, axial sliding of one helix relative to the other was the conformational change that best accounted for the observed effects on cross-linking.