A study of core domains, and the core domain-domain interaction of cytochrome c fragment complex

A look back at the molten globule state of proteins: thermodynamic aspects

Interest in protein folding intermediates lies in their significance to protein folding pathways. The molten globule (MG) state is one such intermediate lying on the kinetic (and sometimes thermodynamic) pathway between native and unfolded states. Development of our qualitative and quantitative understanding of the MG state can provide deeper insight into the folding pathways and hence potentially facilitate solution of the protein folding problem. An extensive look at literature suggests that most studies into protein MG states have been largely qualitative. Attempts to obtain quantitative insights into MG states have involved application of high-sensitivity calorimetry (differential scanning calorimetry and isothermal titration calorimetry). This review addresses the progress made in this direction by discussing the knowledge gained to date, along with the future promise of calorimetry, in providing quantitative information on the structural features of MG states. Particular attention is paid to the question of whether such states share common structural features or not. The difference in the nature of the transition from the MG state to the unfolded state, in terms of cooperativity, has also been addressed and discussed.

Protein reconstitution and three-dimensional domain swapping: benefits and constraints of covalency

The phenomena of protein reconstitution and three-dimensional domain swapping reveal that highly similar structures can be obtained whether a protein is comprised of one or more polypeptide chains. In this review, we use protein reconstitution as a lens through which to examine the range of protein tolerance to chain interruptions and the roles of the primary structure in related features of protein structure and folding, including circular permutation, natively unfolded proteins, allostery, and amyloid fibril formation. The results imply that noncovalent interactions in a protein are sufficient to specify its structure under the constraints imposed by the covalent backbone.

The 40s ?-loop plays a critical role in the stability and the alkaline conformational transition of cytochrome c

The structural and redox properties of a non-covalent complex reconstituted upon mixing two non-contiguous fragments of horse cytochrome c, the residues 1-38 heme-containing N-fragment with the residues 57-104 C-fragment, have been investigated. With respect to native cyt c, the complex lacks a segment of 18 residues, corresponding, in the native protein, to an omega (Omega)-loop region. The fragment complex shows compact structure, native-like alpha-helix content but a less rigid atomic packing and reduced stability with respect to the native protein. Structural heterogeneity is observed at pH 7.0, involving formation of an axially misligated low-spin species and consequent partial displacement of Met80 from the sixth coordination position of the heme-iron. Spectroscopic data suggest that a lysine (located in the Met80-containing loop, namely Lys72, Lys73, or Lys79) replaces the methionine residue. The residues 1-38/57-104 fragment complex shows an unusual biphasic alkaline titration characterized by a low (p K(a1)=6.72) and a high p K(a)-associated state transition (p K(a2)=8.56); this behavior differs from that of native cyt c, which shows a monophasic alkaline transition (p K(a)=8.9). The data indicate that the 40s Omega-loop plays an important role in the stability of cyt c and in ensuring a correct alkaline conformational transition of the protein.

Reducing the computational complexity of protein folding via fragment folding and assembly

Understanding, and ultimately predicting, how a 1-D protein chain reaches its native 3-D fold has been one of the most challenging problems during the last few decades. Data increasingly indicate that protein folding is a hierarchical process. Hence, the question arises as to whether we can use the hierarchical concept to reduce the practically intractable computational times. For such a scheme to work, the first step is to cut the protein sequence into fragments that form local minima on the polypeptide chain. The conformations of such fragments in solution are likely to be similar to those when the fragments are embedded in the native fold, although alternate conformations may be favored during the mutual stabilization in the combinatorial assembly process. Two elements are needed for such cutting: (1) a library of (clustered) fragments derived from known protein structures and (2) an assignment algorithm that selects optimal combinations to "cover" the protein sequence. The next two steps in hierarchical folding schemes, not addressed here, are the combinatorial assembly of the fragments and finally, optimization of the obtained conformations. Here, we address the first step in a hierarchical protein-folding scheme. The input is a target protein sequence and a library of fragments created by clustering building blocks that were generated by cutting all protein structures. The output is a set of cutout fragments. We briefly outline a graph theoretic algorithm that automatically assigns building blocks to the target sequence, and we describe a sample of the results we have obtained.

Comparison of protein fragments identified by limited proteolysis and by computational cutting of proteins

Here we present a comparison between protein fragments produced by limited proteolysis and those identified by computational cutting based on the building block folding model. The principles upon which the two methods are based are different. Limited proteolysis of natively folded proteins occurs at flexible sites and never at the level of chain segments of regular secondary structure such as alpha-helices. Therefore, the targets for limited proteolysis are locally unfolded regions. In contrast, the computational cutting algorithm considers the compactness of the fragments, their nonpolar buried surface area, and their isolatedness, that is, the surface area which was buried prior to the cutting and becomes exposed subsequently. Despite the different criteria, there is an overall correspondence between sites or regions of limited proteolysis with those identified by computational cutting. The computational cutting method has been applied to several model proteins for which detailed limited proteolysis data are available, namely apomyoglobin, cytochrome c, ribonuclease A, alpha-lactalbumin, and thermolysin. As expected, more cuts are obtained computationally than experimentally and the agreement is better when a number of proteolytic enzymes are used. For example, cytochrome c is cleaved by thermolysin at 56-57, 45-46, and at 80-81, and by proteinase K at 48-49 and 50-51. Incubation of the noncovalent and native-like complex of cytochrome c fragments 1-56 and 57-104 with proteinase K yielded the gapped protein species 1-48/57-104 and finally 1-40/57-104. Computational cutting of cytochrome c reproduced the major experimental observations, with cuts at 47, 64-65 or 65-66 and 80-81 and an unstable 32-47 region not assigned to any building block. The next step, not addressed in this work, is to probe the ability of the generated fragments to fold independently. Since both the computational algorithm and limited proteolysis attempt to dissect the protein folding problem, the general agreement between the two procedures is gratifying. This consistency allows us to propose the use of limited proteolysis to produce protein fragments that can adopt an independent folding and, therefore, to study folding intermediates. The results of the present study appear to validate the building block folding model and are in line with the proposal that protein folding is a hierarchical process, where parts constituting local minima of energy fold first, with their subsequent association and mutual stabilization to finally yield the global fold.

Principles of docking: An overview of search algorithms and a guide to scoring functions

The docking field has come of age. The time is ripe to present the principles of docking, reviewing the current state of the field. Two reasons are largely responsible for the maturity of the computational docking area. First, the early optimism that the very presence of the "correct" native conformation within the list of predicted docked conformations signals a near solution to the docking problem, has been replaced by the stark realization of the extreme difficulty of the next scoring/ranking step. Second, in the last couple of years more realistic approaches to handling molecular flexibility in docking schemes have emerged. As in folding, these derive from concepts abstracted from statistical mechanics, namely, populations. Docking and folding are interrelated. From the purely physical standpoint, binding and folding are analogous processes, with similar underlying principles. Computationally, the tools developed for docking will be tremendously useful for folding. For large, multidomain proteins, domain docking is probably the only rational way, mimicking the hierarchical nature of protein folding. The complexity of the problem is huge. Here we divide the computational docking problem into its two separate components. As in folding, solving the docking problem involves efficient search (and matching) algorithms, which cover the relevant conformational space, and selective scoring functions, which are both efficient and effectively discriminate between native and non-native solutions. It is universally recognized that docking of drugs is immensely important. However, protein-protein docking is equally so, relating to recognition, cellular pathways, and macromolecular assemblies. Proteins function when they are bound to other molecules. Consequently, we present the review from both the computational and the biological points of view. Although large, it covers only partially the extensive body of literature, relating to small (drug) and to large protein-protein molecule docking, to rigid and to flexible. Unfortunately, when reviewing these, a major difficulty in assessing the results is the non-uniformity in the formats in which they are presented in the literature. Consequently, we further propose a way to rectify it here.

Protein interactions leading to conformational changes monitored by limited proteolysis: apo form and fragments of horse cytochrome c

Proteolysis experiments have been used to monitor the conformational transitions from an unfolded to a folded state occurring when the apo form of horse cytochrome c (cyt c) binds the heme moiety or when two fragments of cyt c form a native-like 1:1 complex. Proteinase K was used as a proteolytic probe, in view of the fact that the broad substrate specificity of this protease allows digestion at many sites along a polypeptide chain. The rather unfolded apo form of cyt c binds heme with a concomitant conformational transition to a folded species characterized by an enhanced content of helical secondary structure. While the holoprotein is fully resistant to proteolytic digestion and the apoprotein is digested to small peptides, the noncovalent complex of the apoprotein and heme exhibits an intermediate resistance to proteolysis, in agreement with the fact that the more folded structure of the complex makes the protein substrate more resistant to proteolysis. The noncovalent native-like complex of the two fragments 1-56 and 57-104 of cyt c, covering the entire polypeptide chain of 104 residues of the protein, is rather resistant to proteolysis, while the individual fragments are easily digested. Fragment 57-104 is fast degraded to several peptides, while fragment 1-56 is slowly degraded stepwise from its C-terminal end, leading initially mostly to fragments 1-48 and 1-40 and, at later stages of proteolysis, fragments 1-38, 1-35, 1-33, and 1-31. Thus, proteolysis data indicate that the heme containing fragment 1-56 has a rather compact core and a C-terminal flexible tail. Upon prolonged incubation of the complex of fragments 1-56 and 57-104 (nicked cyt c) with proteinase K, a chain segment is removed from the nicked protein, leading to a gapped protein complex of fragments of 1-48 and 57-104 and, on further digestion, fragments 1-40 and 57-104. Of interest, the chain segment being removed by proteolysis of the complex matches the omega-loop which is evolutionarily removed in cyt c of microbial origin. Overall, rates and/or resistance to proteolysis correlates well with the extent of folding of the protein substrates, as deduced from circular dichroism measurements. Thus, our results underscore the utility of proteolytic probes for analyzing conformational and dynamic features of proteins. Finally, a specific interest of the cyt c fragment system herewith investigated resides in the fact that the fragments are exactly the exon products of the cyt c gene.

A study of the influence of the hydrophobic core residues of yeast iso-2-cytochrome c on phosphate binding: a probe of the hydrophobic core-surface charge interactions

To gain insight into the role of hydrophobic core-surface charge interactions in stabilizing cytochrome c, we investigated the influence of hydrophobic core residues on phosphate binding by mutating residues in yeast iso-2-cytochrome c to those corresponding to iso-l-cytochrome c in various combinations. Heat transition of ultraviolet CD was followed as a function of pH in the presence and absence of phosphate. Thermodynamic parameters were deduced. It was found that the I20V/V43A/M98L mutation in the hydrophobic core, whose locations are remote from the putative phosphate sites, modulates phosphate interactions. The modulation is pH dependent. The I20V/ M98L and V43A mutation effects are nonadditive. The results lead to a model analogous to that of Tsao, Evans, and Wennerstrom, where a domain associated with the ordered hydrophobic core is sensitive to the fields generated by the surface charges. Such an explanation would be in accord with the observed difference in thermal stability between iso-2 and horse cytochromes c.

In vivo assembly of aspartate transcarbamoylase from fragmented and circularly permuted catalytic polypeptide chains

Previous studies on Escherichia coli aspartate transcarbamoylase (ATCase) demonstrated that active, stable enzyme was formed in vivo from complementing polypeptides of the catalytic (c) chain encoded by gene fragments derived from the pyrBI operon. However, the enzyme lacked the allosteric properties characteristic of wild-type ATCase. In order to determine whether the loss of homotropic and heterotropic properties was attributable to the location of the interruption in the polypeptide chain rather than to the lack of continuity, we constructed a series of fragmented genes so that the breaks in the polypeptide chains would be dispersed in different domains and diverse regions of the structure. Also, analogous molecules containing circularly permuted c chains with altered termini were constructed for comparison with the ATCase molecules containing fragmented c chains. Studies were performed on four sets of ATCase molecules containing cleaved c chains at positions between residues 98 and 99, 121 and 122, 180 and 181, and 221 and 222; the corresponding circularly permuted chains had N termini at positions 99, 122, 181, and 222. All of the ATCase molecules containing fragmented or circularly permuted c chains exhibited the homotropic and heterotropic properties characteristic of the wild-type enzyme. Hill coefficients (n(H:)) and changes in them upon the addition of ATP and CTP were similar to those observed with wild-type ATCase. In addition, the conformational changes revealed by the decrease in sedimentation coefficient upon the addition of a bisubstrate analog were virtually identical to that for the wild-type enzyme. Differential scanning calorimetry showed that neither the breakage of the polypeptide chains nor the newly formed covalent bond between the termini in the wild-type enzyme had a significant impact on the thermal stability of the assembled dodecamers. The studies demonstrate that continuity of the polypeptide chain within structural domains is not essential for the assembly, activity, and allosteric properties of ATCase.

Thermodynamics of the reconstitution of tuna cytochrome c from two peptide fragments

Two peptide fragments from tuna cytochrome c (cyt c), N-fragment (residues 1-44 containing the heme) and C-fragment (residues 45-103), combine to form a 1:1 fragment complex. This was clearly proved by ion-spray mass spectrometry. It was found from CD and NMR spectra that the structure of the fragment complex formed is similar to that of an intact cyt c, although each isolated fragment itself is unstructured. Binding constants and enthalpies upon the complex formation were directly observed by isothermal titration calorimetry. Thermodynamic parameters (deltaG(o)b, deltaHb, deltaS(o)b, and deltaC(b)p)) associated with the complex formation were determined at various pHs and temperatures. DeltaHb was found to be almost independent of pH values. The change in heat capacity accompanying the complex formation (deltaC(b)p) was directly determined from the temperature dependence of deltaHb. In addition, the change in heat capacity and enthalpy upon tuna cyt c unfolding were determined by differential scanning calorimetry. Thermodynamic parameters for the unfolding/dissociation process of the fragment complex were compared with those for cyt c unfolding at pH 3.9 and 303 K. In a comparison of two unfolding processes, the heat capacity change of each was very close to the other, while both the unfolding enthalpy and entropy of the fragment complex were larger than those of tuna cyt c. These thermodynamic data suggest that the internal interactions between polar groups (hydrogen bonding) and nonpolar groups (van der Waals interactions) are preserved in the fragment complex as well as in the native state of cyt c.

Complementation of peptide fragments of the single domain protein chymotrypsin inhibitor 2

Chymotrypsin inhibitor 2 (CI2) folds kinetically as a single domain protein. It has been shown that elements of native secondary structure do not significantly form in fragments as the 64 residue protein is progressively increased in length from its N terminus, until at least 60 residues are present. Here, we analyse peptides of increasing length from the C terminus and find that native-like structure is not present even in the largest, fragment (7-64). We have examined sets of peptides of the form (1 - x) and ((x + 1)-64) to detect complementation. The only pair that readily complements and gives native-like structure is (1-40) and (41-64), where cleavage occurs in the protease-binding loop of CI2. But, all the pairs of peptides (1 - x) + (41-64) complement for x > 40, as do all pairs of (1-40) + (x-64), where x < 40. The resultant complexes appear to be equivalent to (1-40). (41-64) with the overlapping sequence being unstructured. Thus, the folding of CI2 is extremely co-operative, and interactions have to be made between subdomains (1-40) and (41-64). This is consistent with the mechanism proposed for the folding pathway of intact CI2 in which a diffuse nucleus is formed in the transition state between the alpha-helix in the N-terminal region of the protein and conserved hydrophobic contacts in the C-terminal region of the polypeptide. It is with these protein design features that CI2 can be an effective protease inhibitor.

Identification of cooperative folding units in a set of native proteins

Cooperative unfolding penalties are calculated by statistically evaluating an ensemble of denatured states derived from native structures. The ensemble of denatured states is determined by dividing the native protein into short contiguous segments and defining all possible combinations of native, i.e., interacting, and non-native, i.e., non-interacting, segments. We use a novel knowledge-based scoring function, derived from a set of non-homologous proteins in the Protein Data Bank, to describe the interactions among residues. This procedure is used for the structural identification of cooperative folding cores for four globular proteins: bovine pancreatic trypsin inhibitor, horse heart cytochrome c, French bean plastocyanin, and staphylococcal nuclease. The theoretical folding units are shown to correspond to regions that exhibit enhanced stability against denaturation as determined from experimental hydrogen exchange protection factors. Using a sequence similarity score for related sequences, we show that, in addition to residues necessary for enzymatic function, those amino acids comprising structurally important folding cores are also preferentially conserved during evolution. This implies that the identified folding cores may be part of an array of fundamental structural folding units.

Protein folding intermediates: native-state hydrogen exchange

The hydrogen exchange behavior of native cytochrome c in low concentrations of denaturant reveals a sequence of metastable, partially unfolded forms that occupy free energy levels reaching up to the fully unfolded state. The step from one form to another is accomplished by the unfolding of one or more cooperative units of structure. The cooperative units are entire omega loops or mutually stabilizing pairs of whole helices and loops. The partially unfolded forms detected by hydrogen exchange appear to represent the major intermediates in the reversible, dynamic unfolding reactions that occur even at native conditions and thus may define the major pathway for cytochrome c folding.

Automatic recognition of hydrophobic clusters and their correlation with protein folding units

A method is described to objectively identify hydrophobic clusters in proteins of known structure. Clusters are found by examining a protein for compact groupings of side chains. Compact clusters contain seven or more residues, have an average of 65% hydrophobic residues, and usually occur in protein interiors. Although smaller clusters contain only side-chain moieties, larger clusters enclose significant portions of the peptide backbone in regular secondary structure. These clusters agree well with hydrophobic regions assigned by more intuitive methods and many larger clusters correlate with protein domains. These results are in striking contrast with the clustering algorithm of J. Heringa and P. Argos (1991, J Mol Biol 220:151-171). That method finds that clusters located on a protein's surface are not especially hydrophobic and average only 3-4 residues in size. Hydrophobic clusters can be correlated with experimental evidence on early folding intermediates. This correlation is optimized when clusters with less than nine hydrophobic residues are removed from the data set. This suggests that hydrophobic clusters are important in the folding process only if they have enough hydrophobic residues.

Kinetic mechanism of cytochrome c folding: involvement of the heme and its ligands

The covalently attached heme and its axial ligands not only are essential for the structure and function of cytochrome c but they also play an important role in the folding process. Under typical denaturing conditions (concentrated guanidine hydrochloride or urea near pH 7), one of the axial ligands, His 18, remains bound to the oxidized heme iron, but the second ligand, Met 80, is replaced by a non-native histidine ligand (His 26 or His 33 in horse cytochrome c). Using quenched-flow and NMR methods, hydrogen exchange rates were measured for several individual amide protons in guanidine-denatured horse cytochrome c. The observation of a single highly protected (140-fold) backbone amide, that of His 18, suggests the presence of a persistent H-bond consistent with heme ligation of the His 18 side chain in the unfolded state. Heme absorbance changes induced by rapid acidification of oxidized cytochrome c in 4.5 M guanidine hydrochloride from pH 7.8 to 4.6 or below exhibit two kinetic phases with rates of 110 and 25 s-1, attributed to the dissociation of non-native histidine ligands from the heme in the unfolded state. The kinetics of folding from guanidine-denatured cytochrome c under a variety of initial and final conditions was investigated by stopped-flow methods, using tryptophan fluorescence as a conformational probe and Soret absorbance as a probe for the ligation state of the heme. A fast kinetic phase (80 s-1) accompanied by a major decrease in fluorescence and a minor absorbance change coincides with the formation of a partially folded intermediate with interacting chain termini detected in earlier pulsed NH exchange measurements [Roder, H., Elöve, G. A., & Englander, S. W. (1988) Nature 335, 700]. At neutral pH, an intermediate kinetic phase (1.8 s-1) accounts for 78% of the absorbance change and 47% of the fluorescence change. In contrast, the folding kinetics at pH 5 is dominated by the fast phase, and the amplitude of the intermediate phase is reduced to approximately 10%. The pH-dependent amplitude changes show titration behavior with an apparent pK of approximately 5.7, consistent with the protonation of a single histidine residue. The intermediate phase can also be suppressed by the addition of 20 mM imidazole. Since both of these conditions interfere with histidine ligation, the intermediate kinetic phase is attributed to the presence of a non-native histidine ligand (His 26 or His 33) that can become trapped in a partially folded intermediate.(ABSTRACT TRUNCATED AT 400 WORDS)

A noncovalent peptide complex as a model for an early folding intermediate of cytochrome c

Horse heart cytochrome c is one of a small number of proteins for which the folding pathway has been elucidated in structural detail by pulsed hydrogen exchange and NMR. Those studies indicated that a partially folded intermediate with interacting N- and C-terminal helices is formed at an early stage of folding when most of the chain is still disordered. This report describes a peptide model for this early intermediate, consisting of a noncovalent complex between a heme-containing N-terminal fragment (residues 1-38) and a synthetic peptide corresponding to the C-terminal helix (residues 87-104). Far-UV circular dichroism and proton NMR indicate that the isolated peptides are largely disordered, but when combined, they form a flexible, yet tightly bound complex with enhanced helical structure. These results emphasize the importance of interactions between marginally stable elements of secondary structure in forming tertiary subdomains in protein folding.

Reconstitution of active catalytic trimer of aspartate transcarbamoylase from proteolytically cleaved polypeptide chains

Treatment of the catalytic (C) trimer of Escherichia coli aspartate transcarbamoylase (ATCase) with alpha-chymotrypsin by a procedure similar to that used by Chan and Enns (1978, Can. J. Biochem. 56, 654-658) has been shown to yield an intact, active, proteolytically cleaved trimer containing polypeptide fragments of 26,000 and 8,000 MW. Vmax of the proteolytically cleaved trimer (CPC) is 75% that of the wild-type C trimer, whereas Km for aspartate and Kd for the bisubstrate analog, N-(phosphonacetyl)-L-aspartate, are increased about 7- and 15-fold, respectively. CPC trimer is very stable to heat denaturation as shown by differential scanning microcalorimetry. Amino-terminal sequence analyses as well as results from electrospray ionization mass spectrometry indicate that the limited chymotryptic digestion involves the rupture of only a single peptide bond leading to the production of two fragments corresponding to residues 1-240 and 241-310. This cleavage site involving the bond between Tyr 240 and Ala 241 is in a surface loop known to be involved in intersubunit contacts between the upper and lower C trimers in ATCase when it is in the T conformation. Reconstituted holoenzyme comprising two CPC trimers and three wild-type regulatory (R) dimers was shown by enzyme assays to be devoid of the homotropic and heterotropic allosteric properties characteristic of wild-type ATCase. Moreover, sedimentation velocity experiments demonstrate that the holoenzyme reconstituted from CPC trimers is in the R conformation. These results indicate that the intact flexible loop containing Tyr 240 is essential for stabilizing the T conformation of ATCase. Following denaturation of the CPC trimer in 4.7 M urea and dilution of the solution, the separate proteolytic fragments re-associate to form active trimers in about 60% yield. How this refolding of the fragments, docking, and association to form trimers are achieved is not known.