De novo design of helical bundles as models for understanding protein folding and function

Hot Spot Occlusion from Bulk Water: a Comprehensive Study of the Complex between the Lysozyme HEL and the Antibody FVD1.3

Alanine scanning of protein-protein interfaces has shown that there are some residues in the protein-protein interfaces, responsible for most of the binding free energy, which are called hot spots. Hot spots tend to exist in densely packed central clusters, and a hypothesis has been proposed that considers that inaccessibility to the solvent must be a necessary condition to define a residue as a binding hot spot. This O-ring hypothesis is mainly based on the analysis of the accessible surface area (ASA) of 23 static, crystallographic structures of protein complexes. It is known, however, that protein flexibility allows for temporary exposures of buried interfacial groups, and even though the ASA provides a general trend of the propensity for hydration, protein/solvent-specific interactions or hydrogen bonding cannot be considered here. Therefore, a microscopic level, atomistic picture of hot spot solvation is needed to support the O-ring hypothesis. In this study, we began by applying a computational alanine-scanning mutagenesis technique, which reproduces the experimental results and allows for decomposing the binding free energy difference in its different energetic factors. Subsequently, we calculated the radial distribution function and residence times of the water molecules near the hot/warm spots to study the importance of the water environment around those energetically important amino acid residues. This study shows that within a flexible, dynamic protein framework, the warm/hot spot residues are, indeed, kept sheltered from the bulk solvent during the whole simulation, which allows a better interacting microenvironment.

Analysis of peptide design in four-, five-, and six-helix bundle template assembled synthetic protein molecules

Four-, five-, and six-helix bundle template assembled synthetic proteins (TASPs) have been synthesized using disulfide bonds between cavitand templates and peptides, and characterized in terms of stability and structural specificity. The peptide sequence (CGGGEELLKKLEE LLKKG) used was originally designed for a four-helix bundle. The TASPs were analyzed using CD spectroscopy, chemical denaturation studies, NMR spectroscopy, sedimentation equilibria studies, and hydrophobic dye binding studies to determine the effect of a single peptide sequence when incorporated into bundles with different numbers of helices. If the design was indeed idealized for a four-helix bundle, then the five- and six-helix bundles should be less stable and manifest lower conformational specificity. The TASPs all demonstrated high stability and cooperative unfolding. However, the four-helix bundle was found to be significantly more stable and nativelike compared to the five- and six-helix bundles. This suggests that the peptide sequence is specific to the four-helix bundle, as designed. This result demonstrates the ability to design de novo proteins with specified structure, not just generic stability.

Selective intercalation of six ligand molecules in a self-assembled triple helix

The addition of a ligand molecule to an artificial self-assembled triple helix leads to the selective intercalation of two hydrogen-bonded trimers in specific binding pockets. Furthermore, the triple helix suffers large conformational rearrangements in order to accommodate the ligand molecules in a highly organized manner.

Synthesis of functionalized de novo designed 8–16 kDa model proteins towards metal ion-binding and esterase activity

De novo design and total chemical synthesis of proteins provides a powerful approach for biological and biophysical studies with the ability to prepare artificial proteins with tailored properties, potentially of importance for biophysical studies, material science, nanobioscience, and as molecular probes. In this paper, the previously developed concept of carbohydrates as templates is employed in the de novo design of model proteins (artificial helix bundles) termed 'carboproteins'. The 4-alpha-helix bundle is a macromolecular structure, where four amphiphilic alpha-helical peptide strands form a hydrophobic core. Here this structure is modified towards achieving metal ion-binding and catalytic activity. We report: (i) test of directional effects from different tetravalent carbohydrate templates, (ii) synthesis and evaluation of carboproteins functionalized with phenol, pyridyl or imidazolyl moieties as potential ligands for metal ion-binding as well as for catalysis. Our results include: (i) support of our previous 'controversial' finding that for some carboproteins the degree of alpha-helicity depends on the template, i.e., that there is, to some extent, a controlling effect from the template, (ii) demonstration of binding of Cu(ii) to tetra-functional carboproteins by electrospray ionization-time of flight-mass spectrometry (ESI-TOF-MS), UV-VIS absorption spectroscopy and size exclusion chromatography-inductively coupled plasma-mass spectrometry (SEC-ICP-MS); (iii) a kinetic investigation of the esterase activity.

Characterization of de novo four-helix bundles by molecular dynamics simulations

We have investigated the structure and dynamics of three cavitand-based four-helix bundles (caviteins) by computer simulation. In these systems, designed de novo, each of the four helices contain the identical basis sequence EELLKKLEELLKKG (N1). Each cavitein consists of a rigid macrocycle (cavitand) with four aryl linkages, to each of which is connected an N1 peptide by means of a linker peptide. The three caviteins studied here differ only in the linker peptide, which consist of one, two, or three glycine residues. Previous experimental work has shown that these systems exhibit very different behavior in terms of stability and oligomerization states despite the small differences in the linker peptide. Given that to date no three-dimensional structure is available for these caviteins, we have undertaken a series of molecular dynamics (MD) simulations in explicit water to try to rationalize the large differences in the experimentally observed behavior of these systems. Our results provide insight, for the first time, into why and how the cavitein with a single glycine linker forms dimers. In addition, our results indicate why although the two- and three-glycine-linked caviteins have similar stabilities, they have different native-like characteristics: the cavitein with three glycines can form a supercoiled helix, whereas the one with two glycines cannot. These findings may provide a useful guide in the rational de novo design of novel proteins with finely tunable structures and functions in the future.

The side chain interaction index as a tool for predicting fast-folding elements and the structure and stability of engineered peptides

The side chain interaction index (SCII) is a method of calculating the propensity for short-range interactions among side chains within a peptide sequence. Here, it is shown that the SCII values of secondary structure elements that have been shown to fold early and independently cluster separately from those of structures that fold later and/or are dependent on long-range interactions. In addition, the SCII values of engineered peptides that spontaneously adopt a particular desired fold in solution are significantly different from those of engineered peptides that fail to exhibit a stable conformation. Thus, the SCII, as a measure of local structural stability, constitutes a useful tool in folding prediction and in protein/peptide engineering. A program that allows rapid calculation of SCII values is presented.

Theodor Bücher Lecture. Metabolomics, modelling and machine learning in systems biology - towards an understanding of the languages of cells. Delivered on 3 July 2005 at the 30th FEBS Congress and the 9th IUBMB conference in Budapest

The newly emerging field of systems biology involves a judicious interplay between high-throughput 'wet' experimentation, computational modelling and technology development, coupled to the world of ideas and theory. This interplay involves iterative cycles, such that systems biology is not at all confined to hypothesis-dependent studies, with intelligent, principled, hypothesis-generating studies being of high importance and consequently very far from aimless fishing expeditions. I seek to illustrate each of these facets. Novel technology development in metabolomics can increase substantially the dynamic range and number of metabolites that one can detect, and these can be exploited as disease markers and in the consequent and principled generation of hypotheses that are consistent with the data and achieve this in a value-free manner. Much of classical biochemistry and signalling pathway analysis has concentrated on the analyses of changes in the concentrations of intermediates, with 'local' equations - such as that of Michaelis and Menten v=(Vmax x S)/(S+K m) - that describe individual steps being based solely on the instantaneous values of these concentrations. Recent work using single cells (that are not subject to the intellectually unsupportable averaging of the variable displayed by heterogeneous cells possessing nonlinear kinetics) has led to the recognition that some protein signalling pathways may encode their signals not (just) as concentrations (AM or amplitude-modulated in a radio analogy) but via changes in the dynamics of those concentrations (the signals are FM or frequency-modulated). This contributes in principle to a straightforward solution of the crosstalk problem, leads to a profound reassessment of how to understand the downstream effects of dynamic changes in the concentrations of elements in these pathways, and stresses the role of signal processing (and not merely the intermediates) in biological signalling. It is this signal processing that lies at the heart of understanding the languages of cells. The resolution of many of the modern and postgenomic problems of biochemistry requires the development of a myriad of new technologies (and maybe a new culture), and thus regular input from the physical sciences, engineering, mathematics and computer science. One solution, that we are adopting in the Manchester Interdisciplinary Biocentre (http://www.mib.ac.uk/) and the Manchester Centre for Integrative Systems Biology (http://www.mcisb.org/), is thus to colocate individuals with the necessary combinations of skills. Novel disciplines that require such an integrative approach continue to emerge. These include fields such as chemical genomics, synthetic biology, distributed computational environments for biological data and modelling, single cell diagnostics/bionanotechnology, and computational linguistics/text mining.

A designed well-folded monomeric four-helix bundle protein prepared by Fmoc solid-phase peptide synthesis and native chemical ligation

The design and total chemical synthesis of a monomeric native-like four-helix bundle protein is presented. The designed protein, GTD-Lig, consists of 90 amino acids and is based on the dimeric structure of the de novo designed helix-loop-helix GTD-43. GTD-Lig was prepared by the native chemical ligation strategy and the fragments (45 residues long) were synthesized by applying standard fluorenylmethoxycarbonyl (Fmoc) chemistry. The required peptide-thioester fragment was prepared by anchoring the free gamma-carboxy group of Fmoc-Glu-allyl to the solid phase. After chain elongation the allyl moiety was orthogonally removed and the resulting carboxy group was functionalized with a glycine-thioester followed by standard trifluoroacetic acid (TFA) cleavage to produce the unprotected peptide-thioester. The structure of the synthetic protein was examined by far- and near-UV circular dichroism (CD), sedimentation equilibrium ultracentrifugation, and NMR and fluorescence spectroscopy. The spectroscopic methods show a highly helical and native-like monomeric protein consistent with the design. Heat-induced unfolding was studied by tryptophan absorbance and far-UV CD. The thermal unfolding of GTD-Lig occurs in two steps; a cooperative transition from the native state to an intermediate state and thereafter by noncooperative melting to the unfolded state. The intermediate exhibits the properties of a molten globule such as a retained native secondary structure and a compact hydrophobic core. The thermodynamics of GuHCl-induced unfolding were evaluated by far-UV CD monitoring and the unfolding exhibited a cooperative transition that is well-fitted by a two-state mechanism from the native to the unfolded state. GTD-Lig clearly shows the characteristics of a native protein with a well-defined structure and typical unfolding transitions. The design and synthesis presented herein is of general applicability for the construction of large monomeric proteins.

Intelligent design: the de novo engineering of proteins with specified functions

One of the principal successes of de novo protein design has been the creation of small, robust protein-cofactor complexes which can serve as simplified models, or maquettes, of more complicated multicofactor protein complexes commonly found in nature. Different maquettes, generated by us and others, recreate a variety of aspects, or functional elements, recognized as parts of natural enzyme function. The current challenge is to both expand the palette of functional elements and combine and/or integrate them in recreating familiar enzyme activities or generating novel catalysis in the simplest protein scaffolds.

Protein oligomerization: how and why

A large fraction of cellular proteins are oligomeric. Protein oligomerization may often be an advantageous feature from the perspective of protein evolution and has probably evolved by a variety of mechanisms. The study of protein oligomerization may provide insights into the early protein environment and the evolution of modern proteins. Oligomeric mini-proteins, short peptides with discrete protein-like structures, may serve as valuable models for understanding features of protein oligomerization.

Design of amphiphilic protein maquettes: controlling assembly, membrane insertion, and cofactor interactions

We have designed polypeptides combining selected lipophilic (LP) and hydrophilic (HP) sequences that assemble into amphiphilic (AP) alpha-helical bundles to reproduce key structure characteristics and functional elements of natural membrane proteins. The principal AP maquette (AP1) developed here joins 14 residues of a heme binding sequence from a structured diheme-four-alpha-helical bundle (HP1), with 24 residues of a membrane-spanning LP domain from the natural four-alpha-helical M2 channel of the influenza virus, through a flexible linking sequence (GGNG) to make a 42 amino acid peptide. The individual AP1 helices (without connecting loops) assemble in detergent into four-alpha-helical bundles as observed by analytical ultracentrifugation. The helices are oriented parallel as indicated by interactions typical of adjacent hemes. AP1 orients vectorially at nonpolar-polar interfaces and readily incorporates into phospholipid vesicles with >97% efficiency, although most probably without vectorial bias. Mono- and diheme-AP1 in membranes enhance functional elements well established in related HP analogues. These include strong redox charge coupling of heme with interior glutamates and internal electric field effects eliciting a remarkable 160 mV splitting of the redox potentials of adjacent hemes that leads to differential heme binding affinities. The AP maquette variants, AP2 and AP3, removed heme-ligating histidines from the HP domain and included heme-ligating histidines in LP domains by selecting the b(H) heme binding sequence from the membrane-spanning d-helix of respiratory cytochrome bc(1). These represent the first examples of AP maquettes with heme and bacteriochlorophyll binding sites located within the LP domains.

Design of a heterospecific, tetrameric, 21-residue miniprotein with mixed alpha/beta structure

The study of short, autonomously folding peptides, or "miniproteins," is important for advancing our understanding of protein stability and folding specificity. Although many examples of synthetic alpha-helical structures are known, relatively few mixed alpha/beta structures have been successfully designed. Only one mixed-secondary structure oligomer, an alpha/beta homotetramer, has been reported thus far. In this report, we use structural analysis and computational design to convert this homotetramer into the smallest known alpha/beta-heterotetramer. Computational screening of many possible sequence/structure combinations led efficiently to the design of short, 21-residue peptides that fold cooperatively and autonomously into a specific complex in solution. A 1.95 A crystal structure reveals how steric complementarity and charge patterning encode heterospecificity. The first- and second-generation heterotetrameric miniproteins described here will be useful as simple models for the analysis of protein-protein interaction specificity and as structural platforms for the further elaboration of folding and function.

Kinase activation through dimerization by human SH2-B

The isoforms of SH2-B, APS, and Lnk form a family of signaling proteins that have been described as activators, mediators, or inhibitors of cytokine and growth factor signaling. We now show that the three alternatively spliced isoforms of human SH2-B readily homodimerize in yeast two-hybrid and cellular transfections assays, and this is mediated specifically by a unique domain in its amino terminus. Consistent with previous reports, we further show that the SH2 domains of SH2-B and APS bind JAK2 at Tyr813. These findings suggested a model in which two molecules of SH2-B or APS homodimerize with their SH2 domains bound to two JAK2 molecules, creating heterotetrameric JAK2-(SH2-B)2-JAK2 or JAK2-(APS)2-JAK2 complexes. We further show that APS and SH2-B isoforms heterodimerize. At lower levels of SH2-B or APS expression, dimerization approximates two JAK2 molecules to induce transactivation. At higher relative concentrations of SH2-B or APS, kinase activation is blocked. SH2-B or APS homodimerization and SH2-B/APS heterodimerization thus provide direct mechanisms for activating and inhibiting JAK2 and other kinases from the inside of the cell and for potentiating or attenuating cytokine and growth factor receptor signaling when ligands are present.

Symmetry and frustration in protein energy landscapes: a near degeneracy resolves the Rop dimer-folding mystery

Protein folding has become one of the best understood biochemical reactions from a kinetic viewpoint. The funneled energy landscape, a consequence of the minimal frustration achieved by evolution in sequences, explains how most proteins fold efficiently and robustly to their functional structure and allows robust prediction of folding kinetics. The folding of Rop (repressor of primer) dimer is exceptional because some of its mutants with a redesigned hydrophobic core both fold and unfold much faster than the WT protein, which seems to conflict with a simple funneled energy landscape for which topology mainly determines the kinetics. We propose that the mystery of Rop folding can be unraveled by assuming a double-funneled energy landscape on which there are two basins that correspond to distinct but related topological structures. Because of the near symmetry of the molecule, mutations can cause a conformational switch to a nearly degenerate yet distinct topology or lead to a mixture of both topologies. The topology predicted to have the lower free-energy barrier height for folding was further found by all-atom modeling to give a better structural fit for those mutants with the extreme folding and unfolding rates. Thus, the non-Hammond effects can be understood within energy-landscape theory if there are in fact two different but nearly degenerate structures for Rop. Mutations in symmetric and regular structures may give rise to frustration and thus result in degeneracy.

Artificial metalloenzymes: proteins as hosts for enantioselective catalysis

Enantioselective catalysis is one of the most efficient ways to synthesize high-added-value enantiomerically pure organic compounds. As the subtle details which govern enantioselection cannot be reliably predicted or computed, catalysis relies more and more on a combinatorial approach. Biocatalysis offers an attractive, and often complementary, alternative for the synthesis of enantiopure products. From a combinatorial perspective, the potential of directed evolution techniques in optimizing an enzyme's selectivity is unrivaled. In this review, attention is focused on the construction of artificial metalloenzymes for enantioselective catalytic applications. Such systems are shown to combine properties of both homogeneous and enzymatic kingdoms. This review also includes our recent research results and implications in the development of new semisynthetic metalloproteins for the enantioselective hydrogenation of N-protected dehydro-amino acids.

The design of coiled-coil structures and assemblies

Protein design allows sequence-to-structure relationships in proteins to be examined and, potentially, new protein structures and functions to be made to order. To succeed, however, the protein-design process requires reliable rules that link protein sequence to structure?function. Although our present understanding of coiled-coil folding and assembly is not complete, through numerous bioinformatics and experimental studies there are now sufficient rules to allow confident design attempts of naturally observed and even novel coiled-coil motifs. This review summarizes the current design rules for coiled coils, and describes some of the key successful coiled-coil designs that have been created to date. The designs range from those for relatively straightforward, naturally observed structures-including parallel and antiparallel dimers, trimers and tetramers, all of which have been made as homomers and heteromers-to more exotic structures that expand the repertoire of Nature's coiled-coil structures. Examples in the second bracket include a probe that binds a cancer-associated coiled-coil protein; a tetramer with a right-handed supercoil; sticky-ended coiled coils that self-assemble to form fibers; coiled coils that switch conformational state; a three-component two-stranded coiled coil; and an antiparallel dimer that directs fragment complementation of larger proteins. Some of the more recent examples show an important development in the field; namely, new designs are being created with function as well as structure in mind. This will remain one of the key challenges in coiled-coil design in the next few years. Other challenges that lie ahead include the need to discover more rules for coiled-coil prediction and design, and to implement these in prediction and design algorithms. The considerable success of coiled-coil design so far bodes well for this, however. It is likely that these challenges will be met and surpassed.

Helix packing motif common to the crystal structures of two undecapeptides containing dehydrophenylalanine residues: Implications for the de novo design of helical bundle super secondary structural modules

De novo designed peptide based super secondary structures are expected to provide scaffolds for the incorporation of functional sites as in proteins. Self-association of peptide helices of similar screw sense, mediated by weak interactions, has been probed by the crystal structure determination of two closely related peptides: Ac-Gly1-Ala2-Delta Phe3-Leu4-Val5-DeltaPhe6-Leu7-Val8-DeltaPhe9-Ala10-Gly11-NH2 (I) and Ac-Gly1-Ala2-DeltaPhe3-Leu4-Ala5-DeltaPhe6-Leu7-Ala8-DeltaPhe9-Ala10-Gly11-NH2 (II). The crystal structures determined to atomic resolution and refined to R factors 8.12 and 4.01%, respectively, reveal right-handed 3(10)-helical conformations for both peptides. CD has also revealed the preferential formation of right-handed 3(10)-helical conformations for both molecules. Our aim was to critically analyze the packing of the helices in the solid state with a view to elicit clues for the design of super secondary structural motifs such as two, three, and four helical bundles based on helix-helix interactions. An important finding is that a packing motif could be identified common to both the structures, in which a given peptide helix is surrounded by six other helices reminiscent of transmembrane seven helical bundles. The outer helices are oriented either parallel or antiparallel to the central helix. The helices interact laterally through a combination of N--H...O, C--H...O, and C--H...pi hydrogen bonds. Layers of interacting leucine residues are seen in both peptide crystal structures. The packing of the peptide helices in the solid state appears to provide valuable leads for the design of super secondary structural modules such as two, three, or four helix bundles by connecting adjacent antiparallel helices through suitable linkers such as tetraglycine segments.

Ab initio prediction of the three-dimensional structure of a de novo designed protein: A double-blind case study

Ab initio structure prediction and de novo protein design are two problems at the forefront of research in the fields of structural biology and chemistry. The goal of ab initio structure prediction of proteins is to correctly characterize the 3D structure of a protein using only the amino acid sequence as input. De novo protein design involves the production of novel protein sequences that adopt a desired fold. In this work, the results of a double-blind study are presented in which a new ab initio method was successfully used to predict the 3D structure of a protein designed through an experimental approach using binary patterned combinatorial libraries of de novo sequences. The predicted structure, which was produced before the experimental structure was known and without consideration of the design goals, and the final NMR analysis both characterize this protein as a 4-helix bundle. The similarity of these structures is evidenced by both small RMSD values between the coordinates of the two structures and a detailed analysis of the helical packing.

Improved side-chain prediction accuracy using an ab initio potential energy function and a very large rotamer library

Accurate prediction of the placement and comformations of protein side chains given only the backbone trace has a wide range of uses in protein design, structure prediction, and functional analysis. Prediction has most often relied on discrete rotamer libraries so that rapid fitness of side-chain rotamers can be assessed against some scoring function. Scoring functions are generally based on experimental parameters from small-molecule studies or empirical parameters based on determined protein structures. Here, we describe the NCN algorithm for predicting the placement of side chains. A predominantly first-principles approach was taken to develop the potential energy function incorporating van der Waals and electrostatics based on the OPLS parameters, and a hydrogen bonding term. The only empirical knowledge used is the frequency of rotameric states from the PDB. The rotamer library includes nearly 50,000 rotamers, and is the most extensive discrete library used to date. Although the computational time tends to be longer than most other algorithms, the overall accuracy exceeds all algorithms in the literature when placing rotamers on an accurate backbone trace. Considering only the most buried residues, 80% of the total residues tested, the placement accuracy reaches 92% for chi(1), and 83% for chi(1 + 2), and an overall RMS deviation of 1 A. Additionally, we show that if information is available to restrict chi(1) to one rotamer well, then this algorithm can generate structures with an average RMS deviation of 1.0 A for all heavy side-chains atoms and a corresponding overall chi(1 + 2) accuracy of 85.0%.

Effect of chain length on coiled-coil stability: Decreasing stability with increasing chain length

The de novo design and biophysical characterization of three series of two-stranded alpha-helical coiled coils with different chain lengths are described. Our goal was to examine how increasing chain length would affect protein folding and stability when one or more heptad repeat(s) of K-A-E-A-L-E-G (gabcdef) was inserted into the central region of different coiled-coil host proteins. This heptad was designed to maintain the continuous 3-4 hydrophobic repeat of the coiled-coil host and introduce an Ala and Leu residue in the hydrophobic core at the a and d position, respectively, and a pair of stabilizing interchain ionic i to i' + 5 (g to e') interactions per heptad inserted. The secondary structures of the three series of disulfide-bridged polypeptides were studied by CD spectroscopy and their stabilities determined by chemical and thermal denaturation. The results showed that successive insertions of this heptad systematically decreased the stability of all the coiled coils studied regardless of the overall initial stability of the host coiled coil. These observations are in contrast to the generally accepted implication that the folding and stability of coiled coils are enhanced with increasing chain length. Our results imply that, in these examples where an Ala and Leu hydrophobic residue were introduced into the coiled-coil core per inserted heptad, there was still insufficient stability to overcome unfavorable entropy associated with chain length extension, even though the inserted heptad contained the most stabilizing hydrophobic residue (Leu) at position d and stabilizing ionic attractions.

A phenylalanine zipper mediates APS dimerization

The APS, SH2-B and LNK proteins are adapters that activate and modulate receptor tyrosine kinase and JAK/STAT signaling. We now show that a conserved N-terminal domain mediates APS homodimerization. We determined the crystal structure of the dimerization domain at a resolution of 1.7 A using bromide ion MAD phasing. Each molecule contributes two helices to a compact four-helix bundle having a bisecting-U topology. Its most conspicuous feature is a stack of interdigitated phenylalanine side chains at the domain core. These residues create a new motif we refer to as a 'phenylalanine zipper,' which is critical to dimerization. A newly developed bridging yeast tri-hybrid assay showed that APS dimerizes JAK2, insulin receptor and IGF1 receptor kinases using its SH2 and dimerization domains. Dimerization via the phenylalanine zipper domain provides a mechanism for activating and modulating tyrosine kinase activity even in the absence of extracellular ligands.

X-ray structure analysis of a designed oligomeric miniprotein reveals a discrete quaternary architecture

The x-ray crystal structure of an oligomeric miniprotein has been determined to a 1.2-A resolution by means of multiwavelength anomalous diffraction phasing with selenomethionine analogs that retain the biophysical characteristics of the native peptide. Peptide 1, comprising alpha and beta secondary structure elements with only 21 aa per monomer, associates as a discrete tetramer. The peptide adopts a previously uncharacterized quaternary structure in which alpha and beta components interact to form a tightly packed and well defined hydrophobic core. The structure provides insight into the origins of the unusual thermal stability of the oligomer. The miniprotein shares many characteristics of larger proteins, including cooperative folding, lack of 1-anilino-8-naphthalene sulfonate binding, and limited deuterium exchange, and possesses a buried surface area typical of native proteins.

De novo design of catalytic proteins

The de novo design of catalytic proteins provides a stringent test of our understanding of enzyme function, while simultaneously laying the groundwork for the design of novel catalysts. Here we describe the design of an O(2)-dependent phenol oxidase whose structure, sequence, and activity are designed from first principles. The protein catalyzes the two-electron oxidation of 4-aminophenol (k(cat)/K(M) = 1,500 M(-1).min(-1)) to the corresponding quinone monoimine by using a diiron cofactor. The catalytic efficiency is sensitive to changes of the size of a methyl group in the protein, illustrating the specificity of the design.

Proteins with simplified hydrophobic cores compared to other packing mutants

Efforts to design proteins with greatly reduced sequence diversity have often resulted in proteins with so-called molten globule properties. Substitutions were made at six neighboring sites in the major hydrophobic core of staphylococcal nuclease to create variants with all leucine, all isoleucine or all valine at these sites. The mutant proteins with simplified cores constructed here are quite unstable and have poorly packed cores, attested to by interaction energies. Eight related mutants with greater sequence diversity were also constructed. Comparison to these mutants and 159 other permutations of these 3 aliphatic side chains at these same 6 sites previously constructed shows that the simplified cores are not unusual in their stabilities or interaction energies. Further, crystal structures of the two mutants with the worst packing, as measured by interaction energies, showed no unusual disorder in the core. Therefore, reduction of sequence diversity is not necessarily incompatible with a single stable native structure. Other factors must also contribute to previous protein design failures.

Kinetic traps in the folding/unfolding of procaspase-1 CARD domain

We have examined the folding and unfolding of the caspase recruitment domain of procaspase-1 (CP1-CARD), a member of the alpha-helical Greek key protein family. The equilibrium folding/unfolding of CP1-CARD is described by a two-state mechanism, and the results show CP1-CARD is marginally stable with a DeltaG(H2O) of 1.1 +/- 0.2 kcal/mole and an m-value of 0.65 +/- 0.06 kcal/mole/M (10 mM Tris-HCl at pH 8.0, 1 mM DTT, 25 degrees C). Consistent with the equilibrium folding data, CP1-CARD is a monomer in solution when examined by size exclusion chromatography. Single-mixing stopped-flow refolding and unfolding studies show that CP1-CARD folds and unfolds rapidly, with no detectable slow phases, and the reactions appear to reach equilibrium within 10 msec. However, double jump kinetic experiments demonstrate the presence of an unfolded-like intermediate during unfolding. The intermediate converts to the fully unfolded conformation with a half-time of 10 sec. Interrupted refolding studies demonstrate the presence of one or more nativelike intermediates during refolding, which convert to the native conformation with a half-time of about 60 sec. Overall, the data show that both unfolding and refolding processes are slow, and the pathways contain kinetically trapped species.

A computer simulation study of water drying at the interface of protein chains

This study investigated the water drying (cavitation) in the interfacial region of two chains of a dimeric protein by nanosecond molecular dynamics simulations using explicit water representation. Separation-induced cavity of water was directly observed in the region. We evaluated the separation length scale of two chains on which the drying transition occurs, and the average number of water molecules that are expelled from the interfacial region during the transition. The obtained values can be rationalized by Kelvin equation for finite lateral size of confinement [K. Lum and A. Luzar, Phys. Rev. E 56, R6283 (1997)]. Also, we found that the drying transition is accompanied by an exponential reduction in the average hydrogen-bond number per interfacial water molecule. The results of this study may deepen the understanding of how hydrophobic interaction drives the assembly of protein chains.

Complementary nucleobase interaction enhances peptide-peptide recognition and self-replicating catalysis

The availability of the complementary interaction of nucleobases for influencing the formation of peptide architectures was explored. Nucleobases were incorporated as additional recognition elements in coiled-coil peptides by employing nucleobase amino acids (NBAs), which are artificial L-alpha-amino gamma-nucleobase-butyric acids. The effect of the base-pair interaction on intermolecular recognition between peptides was evaluated through a self-replication reaction. The self-replication reactions of the peptides with complementary base pairs such as thymine-adenine or guanine-cytosine at the g-g' heptad positions were accelerated in comparison with those of the peptides with mismatched base pairs or without nucleobases. However, thymine-adenine pairs at the e-e' positions did not enhance the self-replication. In the presence of a denaturant, the enhancement effects of complementary base pairs on the reaction disappeared. Thermal denaturation studies showed that the thymine-adenine pairs contributed to stabilization of the coiled-coil structure and that the pairs at the g-g' positions were more effective than those at the e-e' positions. The peptide-peptide interaction was reinforced by complementary nucleobase interactions appropriately arranged in the peptide structure; these led to acceleration of the self-replication reactions.

Designing proteins from the inside out

Globular proteins are characterized by the specific and tight packing of hydrophobic side-chains in the so-called "hydrophobic core." Formation of the core is key in folding, stabilization, and conformational specificity. The critical role of hydrophobic cores in maintaining the highly ordered structures present in natural proteins justifies the tremendous efforts devoted to their redesign. Both experimental and computational combinatorial-based approaches have been reported in the last years as powerful protein design tools. These manage to explore large regions of the sequence/conformational space, allowing the search for alternative protein core arrangements displaying native-like properties. The overall results obtained from core design projects have contributed significantly to our present knowledge of protein folding and function. In addition, core design has worked as a benchmark for the development of ambitious protein design projects that nowadays are allowing the de novo design of novel protein structures and functions.

Stabilizing and Destabilizing Clusters in the Hydrophobic Core of Long Two-stranded α-Helical Coiled-coils

Detailed sequence analyses of the hydrophobic core residues of two long two-stranded alpha-helical coiled-coils that differ dramatically in sequence, function, and length were performed (tropomyosin of 284 residues and the coiled-coil domain of the myosin rod of 1086 residues). Three types of regions were present in the hydrophobic core of both proteins: stabilizing clusters and destabilizing clusters, defined as three or more consecutive core residues of either stabilizing (Leu, Ile, Val, Met, Phe, and Tyr) or destabilizing (Gly, Ala, Cys, Ser, Thr, Asn, Gln, Asp, Glu, His, Arg, Lys, and Trp) residues, and intervening regions that consist of both stabilizing and destabilizing residues in the hydrophobic core but no clusters. Subsequently, we designed a series of two-stranded coiled-coils to determine what defines a destabilizing cluster and varied the length of the destabilizing cluster from 3 to 7 residues to determine the length effect of the destabilizing cluster on protein stability. The results showed a dramatic destabilization, caused by a single Leu to Ala substitution, on formation of a 3-residue destabilizing cluster (DeltaT(m) of 17-21 degrees C) regardless of the stability of the coiled-coil. Any further substitution of Leu to Ala that increased the size of the destabilizing cluster to 5 or 7 hydrophobic core residues in length had little effect on stability (DeltaT(m) of 1.4-2.8 degrees C). These results suggested that the contribution of Leu to protein stability is context-dependent on whether the hydrophobe is in a stabilizing cluster or its proximity to neighboring destabilizing and stabilizing clusters.

Combinatorial approaches to protein stability and structure

Why do proteins adopt the conformations that they do, and what determines their stabilities? While we have come to some understanding of the forces that underlie protein architecture, a precise, predictive, physicochemical explanation is still elusive. Two obstacles to addressing these questions are the unfathomable vastness of protein sequence space, and the difficulty in making direct physical measurements on large numbers of protein variants. Here, we review combinatorial methods that have been applied to problems in protein biophysics over the last 15 years. The effects of hydrophobic core composition, the most important determinant of structure and stability, are still poorly understood. Particular attention is given to core composition as addressed by library methods. Increasingly useful screens and selections, in combination with modern high-throughput approaches borrowed from genomics and proteomics efforts, are making the empirical, statistical correlation between sequence and structure a tractable problem for the coming years.

Exploring folding free energy landscapes using computational protein design

Recent advances in computational protein design have allowed exciting new insights into the sequence dependence of protein folding free energy landscapes. Whereas most previous studies have examined the sequence dependence of protein stability and folding kinetics by characterizing naturally occurring proteins and variants of these proteins that contain a small number of mutations, it is now possible to generate and characterize computationally designed proteins that differ significantly from naturally occurring proteins in sequence and/or structure. These computer-generated proteins provide insights into the determinants of protein structure, stability and folding, and make it possible to disentangle the properties of proteins that are the consequence of natural selection from those that reflect the fundamental physical chemistry of polypeptide chains.

Computational design of protein–protein interactions

Computational protein design strategies have been developed to reengineer protein-protein interfaces in an automated, generalizable fashion. In the past two years, these methods have been successfully applied to generate chimeric proteins and protein pairs with specificities different from naturally occurring protein-protein interactions. Although there are shortcomings in current approaches, both in the way conformational space is sampled and in the energy functions used to evaluate designed conformations, the successes suggest we are now entering an era in which computational methods can be used to modulate, reengineer and design protein-protein interaction networks in living cells.

De novo designed cyclic-peptide heme complexes

The structural characterization of de novo designed metalloproteins together with determination of chemical reactivity can provide a detailed understanding of the relationship between protein structure and functional properties. Toward this goal, we have prepared a series of cyclic peptides that bind to water-soluble metalloporphyrins (FeIII and CoIII). Neutral and positively charged histidine-containing peptides bind with a high affinity, whereas anionic peptides bind only weakly to the negatively charged metalloporphyrin. Additionally, it was found that the peptide becomes helical only in the presence of the metalloporphyrin. CD experiments confirm that the metalloporphyrin binds specific cyclic peptides with high affinity and with isodichroic behavior. Thermal unfolding experiments show that the complex has "native-like" properties. Finally, NMR spectroscopy produced well dispersed spectra and experimental restraints that provide a high-resolution solution structure of the complexed peptide.