Solution NMR structures of IgG binding domains with artificially evolved high levels of sequence identity but different folds

We describe here the solution NMR structures of two IgG binding domains with highly homologous sequences but different three-dimensional structures. The proteins, G311 and A219, are derived from the IgG binding domains of their wild-type counterparts, protein G and protein A, respectively. Through a series of site-directed mutations and phage display selections, the sequences of G311 and A219 were designed to converge to a point of high-level sequence identity while keeping their respective wild-type tertiary folds. Structures of both artificially evolved sequences were determined by NMR spectroscopy. The main chain fold of G311 can be superimposed on the wild-type alpha/beta protein G structure with a backbone rmsd of 1.4 A, and the A219 structure can be overlaid on the wild-type three-alpha-helix protein A fold also with a backbone rmsd of 1.4 A. The structure of G311, in particular, accommodates a large number of mutational changes without undergoing a change in the overall fold of the main chain. The structural differences are maintained despite a high level (59%) of sequence identity. These proteins serve as starting points for further experiments that will probe basic concepts of protein folding and conformational switching.     

The different folding behavior of insulin and insulin-like growth factor 1 is mainly controlled by their B-chain/domain

Although insulin and insulin-like growth factor 1 (IGF-1) share homologous sequence, similar tertiary structure, weakly overlapped biological activity, and a common ancestor, the two highly homologous sequences encode different folding behavior: insulin folds into one unique stable tertiary structure while IGF-1 folds into two disulfide isomers with similar thermodynamic stability. To further elucidate the molecular mechanism of their different folding behavior, we prepared two single-chain hybrids of insulin and IGF-1, Ins(A)/IGF-1(B) and Ins(B)/IGF-1(A), as well as a mini-IGF-1 by means of protein engineering and studied their structure as well as folding behavior. Both mini-IGF-1 and Ins(A)/IGF-1(B) fold into two thermodynamically stable disulfide isomers in vivo and in vitro just like that of IGF-1, while Ins(B)/IGF-1(A) folds into one unique thermodynamically stable tertiary structure in vivo and in vitro just like that of insulin. So we deduce that the different folding behavior of insulin and IGF-1 is mainly controlled by their B-chain/domain. By V8 endoproteinase digestion and circular dichroism analysis, as well as insulin receptor binding assay, we deduce that Ins(B)/IGF-1(A), isomer 2 of mini-IGF-1, and isomer 2 of Ins(A)/IGF-1(B) adopt native IGF-1/insulin-like three-dimensional structure with native disulfides, while isomer 1 of mini-IGF-1 and isomer 1 of Ins(A)/IGF-1(B) adopt the swap IGF-1-like three-dimensional structure with swap disulfides.     

Folding mechanisms of proteins with high sequence identity but different folds

The problem of how a protein folds from a linear chain of amino acids to the three-dimensional structure necessary for function is often investigated using proteins with a low degree of sequence identity that adopt different folds. The design of pairs of proteins with a high degree of sequence identity but different folds offers the opportunity for a complementary study; in two highly similar sequences, which residues are the most important in directing folding to a particular structure? Here we use molecular dynamics simulations to characterize the folding-unfolding pathways of a pair of proteins designed by Bryan and co-workers [Alexander, P. A., et al. (2005) Biochemistry 44, 14045-14054; He, Y. N., et al. (2005) Biochemistry 44, 14055-14061]. Despite being 59% identical, the two protein sequences fold to two different structures. The first sequence folds to the alpha+beta protein G structure and the second to the all-alpha-helical protein A structure. We show that the final protein structure is determined early along the folding pathway. In folding to the protein G structure, the single alpha-helix (alpha1) and the beta3-beta4 turn fold early. Formation of the hairpin turn essentially prevents folding to helical structure in this region of the protein. This early structure is then consolidated by formation of long-range hydrophobic interactions between alpha1 and the beta3-beta4 turn. The protein A sequence differs both in the residues that form the beta3-beta4 turn and also in many of the residues that form the early hydrophobic interactions in the protein G structure. Instead, in the protein A sequence, a more hierarchical mechanism is observed, with helices folding before many of the tertiary interactions are formed. We find that small, but critical, sequence differences determine the topology of the protein early along the folding pathway, which help to explain the process by which one fold can evolve into another.     

The sequence determinant causing different folding behaviors of insulin and insulin-like growth factor-1

Although insulin and insulin-like growth factor-1 (IGF-1) belong to the insulin superfamily and share highly homologous sequences, similar tertiary structure, and a common ancestor molecule, amphioxus insulin-like peptide, they have different folding behaviors: IGF-1 folds into two thermodynamically stable tertiary structures (native and swap forms), while insulin folds into one unique stable structure. To further understand which part of the sequence determines their different folding behavior, based on previous reports from the laboratory, two peptide models, [B9A][1-4]porcine insulin precursor (PIP) and [B10E][1-4]PIP, were constructed. The plasmids encoding the peptides were transformed into yeast cells for expression of the peptides; the results showed that the former peptide was expressed as single component, while the latter was expressed as a mixture of two components (isomer 1 and isomer 2). The expression results together with studies of circular dichoism, disulfide rearrangement, and refolding lead us to deduce that isomer 1 corresponds to the swap form and the isomer 2 corresponds to the native form. We further demonstrate that the sequence 1-4 plus B9 of IGF-1 B-domain can make PIP fold into two structures, while sequence 1-5 of insulin B-chain can make IGF-1 fold into one unique structure. In other words, it is the IGF-1 B-domain sequence that 1-4 allows IGF-1 folding into two thermodynamically stable tertiary structures; this sequence plus its residue B9E can change PIP folding behavior from folding into one unique structure to two thermodynamically stable structures, like that of IGF-1.     

The two cysteine-rich head domains of minicollagen from Hydra nematocysts differ in their cystine framework and overall fold despite an identical cysteine sequence pattern

Synthetic replicates of naturally occurring cysteine-rich peptides such as hormones, neurotransmitters, growth factors, enzyme inhibitors, defensins and toxins often can be oxidatively folded in high yields to their native structure in simple redox buffers. Thereby, identical cysteine patterns in the sequence were found to generate identical disulfide connectivities and homologous spatial structures despite significant variability in the non-cysteine positions. Minicollagen-1 from the nematocysts of Hydra is a trimeric protein that contains cysteine-rich domains at the N and C termini, which are involved in the assembly of an intermolecular disulfide network. Determination of the three-dimensional structures of peptides corresponding to the N-terminal and C-terminal domains by NMR spectroscopy revealed a remarkable exception from the general rule. Despite an identical cysteine pattern, the two domains of minicollagen-1 form different disulfide bridges and exhibit distinctly different folds, both of which are not found in the current structural databases. To our knowledge, this is the first case where two relatively short peptides with the abundant cysteine residues in identical sequence positions fold uniquely and with high yields into defined, but differing, structures. Therefore, the cysteine-rich domains of minicollagen constitute ideal model systems for studies of the interplay between folding and oxidation in proteins.     

The thermodynamic stability of insulin disulfides is not affected by the C-domain of insulin-like growth factor 1

Both Insulin and insulin-like growth factor 1 are members of insulin superfamily. They share homologous primary and tertiary structure as well as weakly overlapping biological activity. However, their folding behavior is different: insulin and its recombinant precursor (PIP) fold into one unique tertiary structure, while IGF-1 folds into two disulfides isomers with similar thermodynamic stability. To elucidate the molecular mechanism of their different folding behavior, we prepared a singlechain hybrid of insulin and IGF-1, [B10Glu]Ins/IGF-1(C), and studied its folding behavior compared with that of PIP and IGF-1. We also separated a major non-native disulfides isomer of the hybrid and studied its refolding. The data showed that the C-domain of IGF-1 did not affect the folding thermodynamics of insulin, that is, the primary structure of the hybrid encoded only one thermodynamically stable disulfides linkage. However, the folding kinetics of insulin was affected by the C-domain of IGF-1.     

Ion concentration and temperature dependence of DNA binding: comparison of PurR and LacI repressor proteins.

Purine repressor (PurR) binding to specific DNA is enhanced by complexing with purines, whereas lactose repressor (LacI) binding is diminished by interaction with inducer sugars despite 30% identity in their protein sequences and highly homologous tertiary structures. Nonetheless, in switching from low- to high-affinity DNA binding, these proteins undergo a similar structural change in which the hinge region connecting the DNA and effector binding domains folds into an alpha-helix and contacts the DNA minor groove. The differences in response to effector for these proteins should be manifest in the polyelectrolyte effect which arises from cations displaced from DNA by interaction with positively charged side chains on a protein and is quantitated by measurement of DNA binding affinity as a function of ion concentration. Consistent with structural data for these proteins, high-affinity operator DNA binding by the PurR-purine complex involved approximately 15 ion pairs, a value significantly greater than that for the corresponding state of LacI (approximately 6 ion pairs). For both proteins, however, conversion to the low-affinity state results in a decrease of approximately 2-fold in the number of cations released per dimeric DNA binding site. Heat capacity changes (DeltaC(p)) that accompany DNA binding, derived from buried apolar surface area, coupled folding, and restriction of motional freedom of polar groups in the interface, also reflect the differences between these homologous repressor proteins. DNA binding of the PurR-guanine complex is accompanied by a DeltaC(p) (-2.8 kcal mol(-1) K(-1)) more negative than that observed previously for LacI (-0.9 to -1.5 kcal mol(-1) K(-1)), suggesting that more extensive protein folding and/or enhanced structural rigidity may occur upon DNA binding for PurR compared to DNA binding for LacI. The differences between these proteins illustrate plasticity of function despite high-level sequence and structural homology and undermine efforts to predict protein behavior on the basis of such similarities.     

The folding nucleus of a fibronectin type III domain is composed of core residues of the immunoglobulin-like fold

To identify the contacts that stabilise the rate-limiting transition state for folding of FNfn10 (the tenth fnIII domain of human fibronectin), 42 mutants have been analysed at 29 positions across this domain. An anomalous response to mutation means that structure formation in the A, B and G strands cannot be evaluated by this method. In all the residues analysed, phi-values are fractional and no completely structured region is observed. The analysis reveals that hydrophobic residues from the central strands of the beta-sandwich form a large core of interactions in the transition state. Br?nsted analysis shows that the stabilisation energy from the amino acid side-chains in the transition state is approximately 40 % of that in the native state. The protein folds by a nucleation-condensation mechanism, and tertiary interactions within the core make up the folding nucleus. Local interactions, in turns and loops, are apparently much less significant. Comparison with an homologous domain from human tenascin (TNfn3), shows that FNfn10 has a more extended, structured transition state spanning three different "layers" of the beta-sandwich. The results support the hypothesis that interactions in the common structural core guide the folding of these domains.     

Abp1p and Fyn SH3 domains fold through similar low-populated intermediate states

Src homology 3 (SH3) domains are small modules that are thought to fold via a two-state mechanism, without the accumulation of significant populations of intermediate states. Relaxation dispersion NMR studies of the folding of G48V and G48M mutants of the Fyn SH3 domain have established that, at least for these modules, folding proceeds through the formation of a transient on-pathway intermediate with an equilibrium population of 1-2% that can be readily detected [Korzhnev, D. M., et al. (2004) Nature 430, 586-590]. To investigate the generality of this result, we present an (15)N relaxation dispersion NMR study of a pair of additional SH3 domains, including a G48V mutant of a stabilized Abp1p SH3 domain that shares 36% sequence identity with the Fyn SH3 module, and a A39V/N53P/V55L mutant Fyn SH3 domain. A transient folding intermediate is detected for both of the proteins studied here, and the dispersion data are well fit to a folding model of the form F <--> I <--> U, where F, I, and U correspond to folded, intermediate, and unfolded states, respectively. The temperature dependencies of the folding/unfolding rate constants were obtained so that the thermodynamic properties of each of F, I, and U could be established. The detection of I states in folding pathways of all SH3 domains examined to date via relaxation dispersion NMR spectroscopy indicates that such intermediates may well be a conserved feature in the folding of such domains in general but that their transient nature along with their low population makes detection difficult using more well-established approaches to the study of folding.     

From Discrete to Continuum Models of Three-Dimensional Deformations in Epithelial Sheets

Metamorphic proteins, including proteins with high levels of sequence identity but different folds, are exceptions to the long-standing rule-of-thumb that proteins with as little as 30% sequence identity adopt the same fold. Which topologies can be bridged by these highly identical sequences remains an open question. Here we bridge two 3-α-helix bundle proteins with two radically different folds. Using a straightforward approach, we engineered the sequences of one subdomain within maltose binding protein (MBP, α/β/α-sandwich) and another within outer surface protein A (OspA, β-sheet) to have high sequence identity (80 and 77%, respectively) with engineered variants of protein G (G_A, 3-α-helix bundle). Circular dichroism and nuclear magnetic resonance spectra of all engineered variants demonstrate that they maintain their native conformations despite substantial sequence modification. Furthermore, the MBP variant (80% identical to G_A) remained active. Thermodynamic analysis of numerous G_A and MBP variants suggests that the key to our approach involved stabilizing the modified MBP and OspA subdomains via external interactions with neighboring substructures, indicating that subdomain interactions can stabilize alternative folds over a broad range of sequence variation. These findings suggest that it is possible to bridge one fold with many other topologies, which has implications for protein folding, evolution, and misfolding diseases.     

人IgG四亚类对噬菌体展示Ig结合蛋白单结构域随机组合文库的体外进化筛选

金黄色葡萄球菌蛋白A(Staphylococcal protein A,SpA)和链球菌蛋白G(Streptococcal protein G,SpG)是细菌产生的特异结合宿主抗体的细菌免疫球蛋白结合蛋白(Immunoglobulin(Ig)-binding proteins,IBPs)的代表分子。SpA和SpG均包含由多个序列高度同源的结合结构域重复组成的抗体结合区,各单结构域都具有完全的结合IgG的功能。为研究这些单结构域随机组合能否产生具有新结合特性的组合分子,将SpA的A、B、C、D、E以及SpG的B2、B3共7个单结合结构域随机组合构建成噬菌体展示文库后,应用人IgG1、2、3、4为诱饵分子对该文库进行4轮筛选,获得了SpA天然分子中不存在的单结构域排列组合分子D-C。在筛选过程中,阴性对照噬菌体的逐渐减少、展示两个结构域以上的噬菌体比例不断增多,尤其是D-C组合的选择性富集和其随机连接肽的严格筛选都显示了筛选的有效性和D-C组合的重要性。噬菌体ELISA进一步证实D-C与人IgG四亚类的结合能力远强于天然SpA分子。该研究应用分子进化技术首次获得了一种与人IgG四亚类具有结合优势的新型组合分子D-C,不仅可为IgG纯化、制备、检测等方面的应用提供新的候选分子,还为细菌IBP结构功能的进一步研究提供新的手段。     

Structure-based classification of 45 FK506-binding proteins

The FK506-binding proteins (FKBPs) are a unique group of chaperones found in a wide variety of organisms. They perform a number of cellular functions including protein folding, regulation of cytokines, transport of steroid receptor complexes, nucleic acid binding, histone assembly, and modulation of apoptosis. These functions are mediated by specific domains that adopt distinct tertiary conformations. Using the Threading/ASSEmbly/Refinement (TASSER) approach, tertiary structures were predicted for a total of 45 FKBPs in 23 species. These models were compared with previously characterized FKBP solution structures and the predicted structures were employed to identify groups of homologous proteins. The resulting classification may be utilized to infer functional roles of newly discovered FKBPs. The three-dimensional conformations revealed that this family may have undergone several modifications throughout evolution, including loss of N- and C-terminal regions, duplication of FKBP domains as well as insertions of entire functional motifs. Docking simulations suggest that additional sequence segments outside FKBP domains may modulate the binding affinity of FKBPs to immunosuppressive drugs. The docking models also indicate the presence of a helix-loop-helix (HLH) region within a subset of FKBPs, which may be responsible for the interaction between this group of proteins and nucleic acids.     

Molecular recognition by the Candida albicans group I intron: tertiary interactions with an imino G.A pair facilitate binding of the 5' exon and lower the KM for guanosine

A G.A pair at position -5 in the P1 helix of the Candida albicans ribozyme contributes to tertiary binding of the 5' exon substrate [Disney, M. D., Haidaris, C. G., and Turner, D. H. (2001) Biochemistry 40, 6507-6519]. Here, the G in the G.A pair is replaced with inosine (I) in both semisynthetic ribozymes and oligonucleotide mimics of the internal guide sequence. Comparisons of oligonucleotide binding affinity for these and other sequences indicate that the G.A pair is in an imino conformation where the exocyclic amine of G contributes approximately 1.4 kcal/mol to tertiary interactions that help dock the ribozyme's P1 helix. Furthermore, replacement of the G.A pair with a G-C pair produces less favorable interactions with the 2'-hydroxyl group at the -3 position and a less favorable K(M) for pG in a ribozyme-catalyzed transesterification reaction. These results are also consistent with the G.A pair promoting docking of the P1 helix into the catalytic core. Evidently, tertiary interactions with the exocyclic amino group of a G in a single G.A pair can increase the equilibrium constant for tertiary folding of RNA by roughly 10-fold at 37 degrees C. Results with a G.U or G.G pair replacing the G.A pair at the -5 position suggest similar tertiary interactions with these pairs.     

Conserved residues and their role in the structure, function, and stability of acyl-coenzyme A binding protein

In the family of acyl-coenzyme A binding proteins, a subset of 26 sequence sites are identical in all eukaryotes and conserved throughout evolution of the eukaryotic kingdoms. In the context of the bovine protein, the importance of these 26 sequence positions for structure, function, stability, and folding has been analyzed using single-site mutations. A total of 28 mutant proteins were analyzed which covered 17 conserved sequence positions and three nonconserved positions. As a first step, the influence of the mutations on the protein folding reaction has been probed, revealing a folding nucleus of eight hydrophobic residues formed between the N- and C-terminal helices [Kragelund, B. B., et al. (1999) Nat. Struct. Biol. (In press)]. To fully analyze the role of the conserved residues, the function and the stability have been measured for the same set of mutant proteins. Effects on function were measured by the extent of binding of the ligand dodecanoyl-CoA using isothermal titration calorimetry, and effects on protein stability were measured with chemical denaturation followed by intrinsic tryptophan and tyrosine fluorescence. The sequence sites that have been conserved for direct functional purposes have been identified. These are Phe5, Tyr28, Tyr31, Lys32, Lys54, and Tyr73. Binding site residues are mainly polar or charged residues, and together, four of these contribute approximately 8 kcal mol-1 of the total free energy of binding of 11 kcal mol-1. The sequence sites conserved for stability of the structure have likewise been identified and are Phe5, Ala9, Val12, Leu15, Leu25, Tyr28, Lys32, Gln33, Tyr73, Val77, and Leu80. Essentially, all of the conserved residues that maintain the stability are hydrophobic residues at the interface of the helices. Only one conserved polar residue, Gln33, is involved in stability. The results indicate that conservation of residues in homologous proteins may result from a summed optimization of an effective folding reaction, a stable native protein, and a fully active binding site. This is important in protein design strategies, where optimization of one of these parameters, typically function or stability, may influence any of the others markedly.     

Solvent protection of the hammerhead ribozyme in the ground state: evidence for a cation-assisted conformational change leading to catalysis

Tertiary folding of the hammerhead ribozyme has been analyzed by hydroxyl radical footprinting. Three hammerhead constructs with distinct noncore sequences, connectivities, and catalytic properties show identical protection patterns, in which conserved core residues (G5, A6, U7, G8, and A9) and the cleavage site (C17, G1.1, and U1.2) become reproducibly protected from nucleolytic attack by radicals. Metal ion titrations show that all protections appear together, suggesting a single folding event to a common tertiary structure, rather than an ensemble of different folds. The apparent binding constants for folding and catalysis by Mg(2+) are lower than those for Li(+) by 3 orders of magnitude, but in each case the protected sites are identical. For both Mg(2+) and Li(+), the ribozyme folds into the protected tertiary structure at significantly lower cation concentrations than those required for cleavage. The sites of protection include all of the sites of reduced solvent accessibility calculated from two different crystal structures, including both core and noncore nucleotides. In addition, experimentally observed protected sites include additional sequences adjacent to those predicted by the crystal structures, suggesting that the solution structure may be folded into a more compact shape. A 2'-deoxy substitution at G5 abolishes all protection, indicating that the 2'-OH is essential for folding. Together, these results support a model in which low concentrations of metal ions fold the ribozyme into a stable ground state tertiary structure that is similar to the crystallographic structures, and higher concentrations of metal ions support a transient conformational change into the transition state for catalysis. These data do not themselves address the issue as to whether a large- or small-scale conformational change is required for catalysis.     

Folding of a synthetic parallel beta-roll protein

Recently, the design of beta-sheet proteins and concomitant folding studies have attracted increasing attention. A unique natural all-beta domain occurs in a family of cytolytic bacterial toxins, the so-called RTX toxins. This domain consists of a variable number (about 6-45) of tandem repeats of a glycine-rich nine-residue motif with the consensus sequence GGXGXDX(L/I/F)X. The analysis of the three-dimensional structure of alkaline protease from Pseudomonas aeruginosa which possesses six of these repeats revealed that they fold into a novel 'parallel beta-roll' where calcium is bound within the turns connecting the beta-strands. A 75-mer peptide of the sequence NH(2)-WLS-[GGSGNDNLS](8)-COOH was chemically synthesised. Circular dichroism spectroscopy showed that this polypeptide folds in the presence of Ca(2+) and polyethylene glycol into a beta-structure which is presumably identical with the parallel beta-roll. This synthetic beta-roll behaves similarly to the isolated beta-roll domains from Escherichia coli haemolysin or Bordetella pertussis cyclolysin in terms of calcium binding and polymerisation behaviour.     

A comparison of the GroE chaperonin requirements for sequentially and structurally homologous malate dehydrogenases: the importance of folding kinetics and solution environment

Escherichia coli malate dehydrogenase (EcMDH) and its eukaryotic counterpart, porcine mitochondrial malate dehydrogenase (PmMDH), are highly homologous proteins with significant sequence identity (60%) and virtually identical native structural folds. Despite this homology, EcMDH folds rapidly and efficiently in vitro and does not seem to interact with GroE chaperonins at physiological temperatures (37 degrees C), whereas PmMDH folds much slower than EcMDH and requires these chaperonins to fold to the native state at 37 degrees C. Double jump experiments indicate that the slow folding behavior of PmMDH is not limited by proline isomerization. Although the folding enhancer glycerol (<5 m) does not alter the renaturation kinetics of EcMDH, it dramatically accelerates the spontaneous renaturation of PmMDH at all temperatures tested. Kinetic analysis of PmMDH renaturation with increasing glycerol concentrations suggests that this osmolyte increases the on-pathway kinetics of the monomer folding to assembly-competent forms. Other osmolytes such as trimethylamine N-oxide, sucrose, and betaine also reactivate PmMDH at nonpermissive temperatures (37 degrees C). Glycerol jump experiments with preformed GroEL.PmMDH complexes indicate that the shift between stringent (requires ATP and GroES) and relaxed (only requires ATP) complex conformations is rapid (<3-5 s). The similarity in irreversible misfolding kinetics of PmMDH measured with glycerol or the activated chaperonin complex (GroEL.GroES.ATP) suggests that these folding aids may influence the same step in the PmMDH folding reaction. Moreover, the interactions between glycerol-induced PmMDH folding intermediates and GroEL.GroES.ATP are diminished. Our results support the notion that the protein folding kinetics of sequentially and structurally homologous proteins, rather than the structural fold, dictates the GroE chaperonin requirement.     

Assembly of transmembrane helices of simple polytopic membrane proteins from sequence conservation patterns

The transmembrane (TM) domains of most membrane proteins consist of helix bundles. The seemingly simple task of TM helix bundle assembly has turned out to be extremely difficult. This is true even for simple TM helix bundle proteins, i.e., those that have the simple form of compact TM helix bundles. Herein, we present a computational method that is capable of generating native-like structural models for simple TM helix bundle proteins having modest numbers of TM helices based on sequence conservation patterns. Thus, the only requirement for our method is the presence of more than 30 homologous sequences for an accurate extraction of sequence conservation patterns. The prediction method first computes a number of representative well-packed conformations for each pair of contacting TM helices, and then a library of tertiary folds is generated by overlaying overlapping TM helices of the representative conformations. This library is scored using sequence conservation patterns, and a subsequent clustering analysis yields five final models. Assuming that neighboring TM helices in the sequence contact each other (but not that TM helices A and G contact each other), the method produced structural models of Calpha atom root-mean-square deviation (CA RMSD) of 3-5 A from corresponding crystal structures for bacteriorhodopsin, halorhodopsin, sensory rhodopsin II, and rhodopsin. In blind predictions, this type of contact knowledge is not available. Mimicking this, predictions were made for the rotor of the V-type Na(+)-adenosine triphosphatase without such knowledge. The CA RMSD between the best model and its crystal structure is only 3.4 A, and its contact accuracy reaches 55%. Furthermore, the model correctly identifies the binding pocket for sodium ion. These results demonstrate that the method can be readily applied to ab initio structure prediction of simple TM helix bundle proteins having modest numbers of TM helices.     

Thermodynamics of a beta-hairpin structure: evidence for cooperative formation of folding nucleus

To elucidate early nucleation stages in protein folding, multi-probed thermodynamic characterization was applied to the beta-hairpin structural formation of G-peptide, which is a C-terminal fragment of the B1 domain of streptococcal protein G. The segment corresponding to the sequence of G-peptide is believed to act as a nucleus during the folding process of the B1 domain. In spite of the broad thermal transition of G-peptide, nuclear magnetic resonance (NMR) melting measurements combined with our original analytical theory enabled us to obtain the thermodynamic properties of the beta-hairpin formation with considerable accuracy. Additionally, all the thermodynamic properties determined by every NMR probe on both the main-chain and the side-chains were quite similar, and also comparable to the values that were independently determined by calorimetric analysis of G-peptide. These results demonstrate that G-peptide folds cooperatively throughout the molecule. In other words, the formation of the beta-hairpin is interpreted as the fashion of a first-order phase transition between two states without any distinguishable intermediates. This cooperative formation of the short linear peptide consisting of only 16 residues provides insight into not only the first folding events of the B1 domain, but also the general principles of proteins in terms of structural hierarchy, stability and folding mechanism.