A polar, solvent-exposed residue can be essential for native protein structure

BACKGROUND: A large energy gap between the native state and the non-native folded states is required for folding into a unique three-dimensional structure. The features that define this energy gap are not well understood, but can be addressed using de novo protein design. Previously, alpha(2)D, a dimeric four-helix bundle, was designed and shown to adopt a native-like conformation. The high-resolution solution structure revealed that this protein adopted a bisecting U motif. Glu7, a solvent-exposed residue that adopts many conformations in solution, might be involved in defining the unique three-dimensional structure of alpha(2)D. RESULTS: A variety of hydrophobic and polar residues were substituted for Glu7 and the dynamic and thermodynamic properties of the resulting proteins were characterized by analytical ultracentrifugation, circular dichroism spectroscopy, and nuclear magnetic resonance spectroscopy. The majority of substitutions at this solvent-exposed position had little affect on the ability to fold into a dimeric four-helix bundle. The ability to adopt a unique conformation, however, was profoundly modulated by the residue at this position despite the similar free energies of folding of each variant. CONCLUSIONS: Although Glu7 is not involved directly in stabilizing the native state of alpha(2)D, it is involved indirectly in specifying the observed fold by modulating the energy gap between the native state and the non-native folded states. These results provide experimental support for hypothetical models arising from lattice simulations of protein folding, and underscore the importance of polar interfacial residues in defining the native conformations of proteins.     

Protein refolding versus aggregation: computer simulations on an intermediate-resolution protein model.

Computer simulations are performed on a system of eight model peptide chains to study how the competition between protein refolding and aggregation affects the optimal conditions for refolding of four-helix bundles. The discontinuous molecular dynamics algorithm is utilized along with an intermediate-resolution protein model that we developed for this work. Physically, the model is much more detailed than any model used to date for simulations of protein aggregation. Each model residue consists of a detailed, three-bead backbone and a simplified, single-bead side-chain. Excluded volume, hydrogen bond, and hydrophobic interactions are modeled with discontinuous (i.e. hard-sphere and square-well) potentials. Simulations efficiently sample conformational space, and complete folding trajectories from random initial configurations to two four-helix bundles are possible within two days on a single processor workstation. Folding of the bundles follows two main pathways, one through a trimeric intermediate and the other through an intermediate with two dimers. The proportion of trajectories that follow each route is significantly different for the eight-peptide system in this work than in a previously studied four-peptide system, which yields one four-helix bundle, suggesting, as our previous simulations have, that protein folding properties are strongly influenced by the presence of other proteins. Folding of the bundles is optimal within a fixed temperature range, with the high-temperature boundary a function of the complexity of the protein (or oligomer) to be folded and the low-temperature boundary a function of the complexity of the protein's environment. Above the optimal temperature range for folding, the model chains tend to unfold; below the optimal range, the model chains tend to aggregate. As has been seen previously, aggregates have substantial levels of native secondary structure, suggesting that aggregates are composed largely of partially folded intermediates, not denatured chains.     

Tuning the strength of a bacterial N-end rule degradation signal

The N-end rule is a degradation pathway conserved from bacteria to mammals that links a protein's stability in vivo to the identity of its N-terminal residue. In Escherichia coli, the components of this pathway directly responsible for protein degradation are the ClpAP protease and its adaptor ClpS. We recently demonstrated that ClpAP is able to recognize N-end motifs in the absence of ClpS although with significantly reduced substrate affinity. In this study, a systematic sequence analysis reveals new features of N-end rule degradation signals. To achieve specificity, recognition of an N-end motif by the protease-adaptor complex uses both the identity of the N-terminal residue and a free alpha-amino group. Acidic residues near the first residue decrease substrate affinity, demonstrating that the identity of adjacent residues can affect recognition although significant flexibility is tolerated. However, shortening the distance between the N-end residue and the stably folded portion of a protein prevents degradation entirely, indicating that an N-end signal alone is not always sufficient for degradation. Together, these data define in vitro the sequence and structural requirements for the function of bacterial N-end signals.     

Populating partially unfolded forms by hydrogen exchange-directed protein engineering

The native-state hydrogen exchange of a redesigned apocytochrome b(562) suggests that at least two partially unfolded forms (PUFs) exist for this four-helix bundle protein under native conditions. The more stable PUF has the N-terminal helix unfolded. To verify the conclusion further and obtain more detailed structural information about this PUF, five hydrophobic core residues in the N-terminal helix were mutated to Gly and Asp to destabilize the native state selectively and populate the PUF for structural studies. The secondary structure and the backbone dynamics of this mutant were characterized using multidimensional NMR. Consistent with the prediction, the N-terminal region of the mutant was found to be unfolded while other parts of the proteins remained folded. These results suggest that native-state hydrogen exchange-directed protein engineering can be a useful approach to populating partially unfolded forms for detailed structural studies.     

Evaluation of the protein interfaces that form an intermolecular four-helix bundle as studied by computer simulation

Rational design of protein surface is important for creating higher order protein structures, but it is still challenging. In this study, we designed in silico the several binding interfaces on protein surfaces that allow a de novo protein–protein interaction to be formed. We used a computer simulation technique to find appropriate amino acid arrangements for the binding interface. The protein–protein interaction can be made by forming an intermolecular four-helix bundle structure, which is often found in naturally occurring protein subunit interfaces. As a model protein, we used a helical protein, YciF. Molecular dynamics simulation showed that a new protein–protein interaction is formed depending on the number of hydrophobic and charged amino acid residues present in the binding surfaces. However, too many hydrophobic amino acid residues present in the interface negatively affected on the binding. Finally, we found an appropriate arrangement of hydrophobic and charged amino acid residues that induces a protein–protein interaction through an intermolecular four-helix bundle formation.     

A designed four-alpha-helix bundle that binds the volatile general anesthetic halothane with high affinity

The structural features of volatile anesthetic binding sites on proteins are being examined with the use of a defined model system consisting of a four-alpha-helix bundle scaffold with a hydrophobic core. Previous work has suggested that introducing a cavity into the hydrophobic core improves anesthetic binding affinity. The more polarizable methionine side chain was substituted for a leucine, in an attempt to enhance the dispersion forces between the ligand and the protein. The resulting bundle variant has an improved affinity (K(d) = 0.20 +/- 0.01 mM) for halothane binding, compared with the leucine-containing bundle (K(d) = 0.69 +/- 0.06 mM). Photoaffinity labeling with (14)C-halothane reveals preferential labeling of the W15 residue in both peptides, supporting the view that fluorescence quenching by bound anesthetic reports both the binding energetics and the location of the ligand in the hydrophobic core. The rates of amide hydrogen exchange were similar for the two bundles, suggesting that differences in binding affinity were not due to changes in protein stability. Binding of halothane to both four-alpha-helix bundle proteins stabilized the native folded conformations. Molecular dynamics simulations of the bundles illustrate the existence of the hydrophobic core, containing both W15 residues. These results suggest that in addition to packing defects, enhanced dispersion forces may be important in providing higher affinity anesthetic binding sites. Alternatively, the effect of the methionine substitution on halothane binding energetics may reflect either improved access to the binding site or allosteric optimization of the dimensions of the binding pocket. Finally, preferential stabilization of folded protein conformations may represent a fundamental mechanism of inhaled anesthetic action.     

Repacking of hydrophobic residues in a stable mutant of apocytochrome b562 selected by phage-display and proteolysis

To test a hydrophobic core-directed protein design approach, we previously have used phage-display and proteolysis to select stably folded proteins from a library of mutants of apocytochrome b562. The consensus sequence of the selected mutants has hydrophilic residues at two of the three positions that are designed to form a hydrophobic core. To understand this unexpected result, we determined the high-resolution structure of one of the selected mutants using multi-dimensional nuclear magnetic resonance (NMR). The structure shows that the two hydrophilic residues in the consensus sequence were on the surface of the structure. Instead, two of their neighboring hydrophobic residues reorganized their side-chain conformations and formed the hydrophobic core. This result suggests that the hydrophobic core-directed protein design by phage-display and proteolysis is a valid method in general but alternative hydrophobic packing needs to be considered in the initial design. The unexpected repacking of the hydrophobic residues also highlights the plastic nature of protein structures.     

Solution structure and dynamics of LuxU from Vibrio harveyi, a phosphotransferase protein involved in bacterial quorum sensing

The marine bacterium Vibrio harveyi controls its bioluminescence by a process known as quorum sensing. In this process, autoinducer molecules are detected by membrane-bound sensor kinase/response regulator proteins (LuxN and LuxQ) that relay a signal via a series of protein phosphorylation reactions to another response regulator protein, LuxO. Phosphorylated LuxO indirectly represses the expression of the proteins responsible for bioluminescence. Integral to this quorum sensing process is the function of the phosphotransferase protein, LuxU. LuxU acts to shuttle the phosphate from the membrane-bound proteins, LuxN and LuxQ, to LuxO. LuxU is a 114 amino acid residue monomeric protein. Solution NMR was used to determine the three-dimensional structure of LuxU. LuxU contains a four-helix bundle topology with the active-site histidine residue (His58) located on alpha-helix C and exposed to solution. The active site represents a cluster of positively charged residues located on an otherwise hydrophobic protein face. NMR spin-relaxation experiments identify a collection of flexible residues localized on the same region of LuxU as His58. The studies described here represent the first structural characterization of an isolated, monomeric bacterial phosphotransferase protein.     

Functional specificity conferred by the unique plasticity of fully alpha-helical Ras and Rho GAPs

Structural comparisons of the two GTPase activating proteins (GAPs) p120 and p50 in complex with Ras and Rho, respectively, allowed us to decipher the functional role of specific structural features, such as helix alpha8c of p120 and helix A1 of p50, necessary for small GTPase recognition. We identified important residues that may be critical for stabilization of the GAP/GTPase binary complexes. Detection of topohydrophobic positions (positions which are most often occupied by hydrophobic amino acids within a family of protein domains) conserved between the two GAP families led to the characterization of a common flexible four-helix bundle. Altogether, these data are consistent with a rearrangement of several helices around a common core, which strongly supports the assumption that p50 and p120 GAPs derive from a unique fold. Considered as a whole, the remarkable plasticity of GAPs appears to be a means used by nature to accurately confer functional specificity.     

Fluorous effect in proteins: de novo design and characterization of a four-alpha-helix bundle protein containing hexafluoroleucine

Several studies have demonstrated that proteins incorporating fluorinated analogues of hydrophobic amino acids such as leucine and valine into their hydrophobic cores exhibit increased stability toward thermal denaturation and unfolding by guanidinium chloride. However, estimates for the increase in the thermodynamic stability of a protein (DeltaDeltaG(unfold)) afforded by the substitution of a hydrophobic amino acid with its fluorinated analogue vary quite significantly. To address this, we have designed a peptide that adopts an antiparallel four-helix bundle structure in which the hydrophobic core is packed with leucine, and investigated the effects of substituting the central two layers of the core with L-5,5,5,5',5',5'-hexafluoroleucine (hFLeu). We find that DeltaDeltaG(unfold) is increased by 0.3 kcal/mol per hFLeu residue. This is in good agreement with the predicted increase in DeltaDeltaG(unfold) of 0.4 kcal/mol per residue arising from the increased hydrophobicity of the hFLeu side chain, which we determined experimentally from partitioning measurements on hFLeu and leucine. The increased stability of this fluorinated protein may therefore be ascribed to simple hydrophobic effects, rather than specific "fluorous" interactions between the hFLeu residues.     

Conformation of the coat protein of filamentous bacteriophage Pf1 determined by neutron diffraction from magnetically oriented gels of specifically deuterated virions

The structure of filamentous bacteriophage Pf1 has been studied using neutron diffraction from magnetically oriented gels of native and valine-deuterated phage. Neutron diffraction intensities were measured to approximately 8 A resolution along the equator and first six layer-lines, and differences due to the deuterated valine residues were apparent. Analysis of equatorial data indicate that one valine residue is located at a radius of about 13 A, three are in the hydrophobic center of the protein coat at an average of about 22 A radius, and one is near the outer surface of the virion at about 28 A radius. Analysis of the three-dimensional data was initiated using the rod model for the alpha-helices of the coat protein derived from earlier X-ray diffraction studies. This model was refined against the neutron diffraction intensities from native phage to obtain a phase set that was used to calculate a difference map between the valine-deuterated and native phage. The difference map exhibits peaks that correspond to the positions of the five valine residues in the coat protein. From the amino acid sequence and the alpha-helical conformation of the coat protein, the five valine residues can be unambiguously assigned to the difference peaks. This assignment indicates that the two alpha-helices of the coat protein are parallel to one another, connected by a short stretch of non-helical peptide. The valine positions also indicate that the helical surface lattice of the phage particle is right-handed.     

Binding of small molecules to cavity forming mutants of a de novo designed protein

A central goal of protein design is to devise novel proteins for applications in biotechnology and medicine. Many applications, including those focused on sensing and catalysis will require proteins that recognize and bind to small molecules. Here, we show that stably folded α-helical proteins isolated from a binary patterned library of designed sequences can be mutated to produce binding sites capable of binding a range of small aromatic compounds. Specifically, we mutated two phenylalanine side chains to alanine in the known structure of de novo protein S-824 to create buried cavities in the core of this four-helix bundle. The parental protein and the Phe→Ala variants were exposed to mixtures of compounds, and selective binding was assessed by saturation transfer difference NMR. The affinities of benzene and a number of its derivatives were determined by pulse field gradient spin echo NMR, and several of the compounds were shown to bind the mutated protein with micromolar dissociation constants. These studies suggest that stably folded de novo proteins from binary patterned libraries are well-suited as scaffolds for the design of binding sites.     

Design and synthesis of de novo cytochromes c

Natural c-type cytochromes are characterized by the consensus Cys-X-X-Cys-His heme-binding motif (where X is any amino acid) by which the heme is covalently attached to protein by the addition of the sulfhydryl groups of two cysteine residues to the vinyl groups of the heme. In this work, the consensus sequence was used for the heme-binding site of a designed four-helix bundle, and the apoproteins with either a histidine residue or a methionine residue positioned at the sixth coordination site were synthesized and reacted with iron protoporphyrin IX (protoheme) under mild reducing conditions in vitro. These polypeptides bound one heme per helix-loop-helix monomer via a single thioether bond and formed four-helix bundle dimers in the holo forms as designed. They exhibited visible absorption spectra characteristic of c-type cytochromes, in which the absorption bands shifted to lower wavelengths in comparison with the b-type heme binding intermediates of the same proteins. Unexpectedly, the designed cytochromes c with bis-His-coordinated heme iron exhibited oxidation-reduction potentials similar to those of their b-type intermediates, which have no thioether bond. Furthermore, the cytochrome c with His and Met residues as the axial ligands exhibited redox potentials increased by only 15-30 mV in comparison with the cytochrome with the bis-His coordination. These results indicate that highly positive redox potentials of natural cytochromes c are not only due to the heme covalent structure, including the Met ligation, but also due to noncovalent and hydrophobic environments surrounding the heme. The covalent attachment of heme to the polypeptide in natural cytochromes c may contribute to their higher redox potentials by reducing the thermodynamic stability of the oxidized forms relatively against that of the reduced forms without the loss of heme.     

Solvent-exposed residues in the Tet repressor (TetR) four-helix bundle contribute to subunit recognition and dimer stability

Dimerization specificity of Tet repressor (TetR) can be altered by changes in the core of the four-helix bundle that mediates protein-protein recognition. We demonstrate here that the affinity of subunit interaction depends also on the solvent-exposed residues at positions 128 and 179'-184', which interact across the dimerization surface. TetR(B) and (D), two naturally occurring sequence variants, differ at position 128 with respect to the monomer-monomer distances in the crystal structures and the charge of the amino acids, being glutamate in TetR(B) and arginine in TetR(D). In vivo analysis of chimeric TetR(B/D) variants revealed that the single E128R exchange does not alter the dimerization specificity of TetR(B) to the one of TetR(D). When combined with specificity mutations in alpha10, it is, however, able to increase dimerization efficiency of the TetR(B/D) chimera with TetR(D). A loss of contact analysis revealed a positive interaction between Arg-128 and residues located at positions 179'-184' of the second monomer. We constructed a hyperstable TetR(B) variant by replacing residues 128 and 179-184 by the respective TetR(D) sequence. These results establish that in addition to a region in the hydrophobic core residues at the solvent-exposed periphery of the dimerization surface participate in protein-protein recognition in the TetR four-helix bundle.     

A method for predicting protein structure from sequence

BACKGROUND: The ability to predict the native conformation of a globular protein from its amino-acid sequence is an important unsolved problem of molecular biology. We have previously reported a method in which reduced representations of proteins are folded on a lattice by Monte Carlo simulation, using statistically-derived potentials. When applied to sequences designed to fold into four-helix bundles, this method generated predicted conformations closely resembling the real ones. RESULTS: We now report a hierarchical approach to protein-structure prediction, in which two cycles of the above-mentioned lattice method (the second on a finer lattice) are followed by a full-atom molecular dynamics simulation. The end product of the simulations is thus a full-atom representation of the predicted structure. The application of this procedure to the 60 residue, B domain of staphylococcal protein A predicts a three-helix bundle with a backbone root mean square (rms) deviation of 2.25-3 A from the experimentally determined structure. Further application to a designed, 120 residue monomeric protein, mROP, based on the dimeric ROP protein of Escherichia coli, predicts a left turning, four-helix bundle native state. Although the ultimate assessment of the quality of this prediction awaits the experimental determination of the mROP structure, a comparison of this structure with the set of equivalent residues in the ROP dime- crystal structure indicates that they have a rms deviation of approximately 3.6-4.2 A. CONCLUSION: Thus, for a set of helical proteins that have simple native topologies, the native folds of the proteins can be predicted with reasonable accuracy from their sequences alone. Our approach suggest a direction for future work addressing the protein-folding problem.     

pH-induced conformational change in an alpha-helical coiled-coil is controlled by His residues in the hydrophobic core

An alpha-helical coiled-coil structure is one of the basic structural units in proteins. Hydrophilic residues at the hydrophobic positions in the coiled-coil structure play important roles in structures and functions of natural proteins. We reported here a peptide that formed a triple stranded alpha-helical coiled-coil showing the pH-dependent structural change. The peptide was designed to have two His residues at the hydrophobic positions of the center of the coiled-coil structure. The peptide folded into a triple stranded coiled-coil at neutral pH, while it unfolded at acidic pH. This construct is useful to create a protein that the structure or function is controlled by pH.     

Engineering, biophysical characterisation and binding properties of a soluble mutant form of annexin A2 domain IV that adopts a partially folded conformation

The four approximately 75-residue domains (repeats) that constitute the annexin core structure all possess an identical five-alpha-helix bundle topology, but the physico-chemical properties of the isolated domains are different. Domain IV of the annexins has previously been expressed only as inclusion bodies, resistant to solubilisation. Analysis of the conserved, exposed hydrophobic residues of the four annexin domains reveals that domain IV contains the largest number of hydrophobic residues involved in interfacial contacts with the other domains. We designed five constructs of domain IV of annexin A2 in which several interfacial hydrophobic residues were substituted by hydrophilic residues. The mutant domain, in which all fully exposed hydrophobic interfacial residues were substituted, was isolated as a soluble protein. Circular dichroism measurements indicate that it harbours a high content of alpha-helical secondary structure and some tertiary structure. The CD-monitored (lambda=222 nm) thermal melting profile suggests a weak cooperative transition. Nuclear magnetic resonance (1H-15N) correlation spectroscopy reveals heterogeneous line broadening and an intermediate spectral dispersion. These properties are indicative of a partially folded protein in which some residues are in a fairly structured conformation, whereas others are in an unfolded state. This conclusion is corroborated by 1-anilinonaphthalene-8-sulfonate fluorescence (ANS) analyses. Surface plasmon resonance measurements also indicate that this domain binds heparin, a known ligand of domain IV in the full-length annexin A2, although with lower affinity.