Stability and folding of domain proteins

Conformational stability and multistate unfolding of poly(A)-specific ribonuclease

Poly(A)-specific ribonuclease (PARN) specifically catalyzes the degradation of the poly(A) tails of single-stranded mRNAs in a highly processive mode. PARN participates in diverse and important intracellular processes by acting as a regulator of mRNA stability and translational efficiency. In this article, the equilibrium unfolding of PARN was studied using both guanidine hydrochloride and urea as chemical denaturants. The unfolding of PARN was characterized as a multistate process, but involving dissimilar equilibrium intermediates when denatured by the two denaturants. A comparison of the spectral characteristics of these intermediates indicated that the conformational changes at low concentrations of the chemical denaturants were more likely to be rearrangements of the tertiary and quaternary structures. In particular, an inactive molten globule-like intermediate was identified to exist as soluble non-native oligomers, and the formation of the oligomers was modulated by electrostatic interactions. An active dimeric intermediate unique to urea-induced unfolding was characterized to have increased regular secondary structures and modified tertiary structures, implying that additional regular structures could be induced by environmental stresses. The dissimilarity in the unfolding pathways induced by guanidine hydrochloride and urea suggest that electrostatic interactions play an important role in PARN stability and regulation. The appearance of multiple intermediates with distinct properties provides the structural basis for the multilevel regulation of PARN by conformational changes.

The partially folded homodimeric intermediate of Escherichia coli aspartate aminotransferase contains a "molten interface" structure

The role of intersubunit side chain-side chain interactions in the stability of the Escherichia coli aspartate aminotransferase (eAATase) homodimer was investigated by directed mutagenesis at 10 different interface contacts. The urea-mediated unfolding pathway of this enzyme proceeds through the formation of a dimeric intermediate, D*, that retains only 40% of the native enzyme secondary structure as judged by circular dichroism. Disruption of any single intersubunit interaction results in a >2.6 kcal mol(-1) decrease in native state stability, independent of its location or nature. However, the stability of D* with respect to U, the unfolded monomer, is the same for all mutants. The stability of the eAATase interface cannot be ascribed to the contribution of a few hot spots, or to the accumulation of a large number of weak interactions, but only to the presence of multiple important and interconnected interactions. It is proposed that a "molten interface" structure, flexible enough to accommodate point mutations, accounts for the stability of D*. Nuclei of tertiary structure, which are not involved in native intersubunit contacts, likely provide a scaffold for the unstructured interface of D*. Such a scaffold would account for the cooperative unfolding of the intermediate.

Consequences of domain insertion on the stability and folding mechanism of a protein

SlyD, the sensitive-to-lysis protein from Escherichia coli, consists of two domains. They are not arranged successively along the protein chain, but one domain, the "insert-in-flap" (IF) domain, is inserted internally as a guest into a surface loop of the host domain, which is a prolyl isomerase of the FK506 binding protein (FKBP) type. We used SlyD as a model to elucidate how such a domain insertion affects the stability and folding mechanism of the host and the guest domain. For these studies, the two-domain protein was compared with a single-domain variant SlyDDeltaIF, SlyD* without the chaperone domain (residues 1-69 and 130-165) in which the IF domain was removed and replaced by a short loop, as present in human FKBP12. Equilibrium unfolding and folding kinetics followed an apparent two-state mechanism in the absence and in the presence of the IF domain. The inserted domain decreased, however, the stability of the host domain in the transition region and decelerated its refolding reaction by about 10-fold. This originates from the interruption of the chain connectivity by the IF domain and its inherent instability. To monitor folding processes in this domain selectively, a Trp residue was introduced as fluorescent probe. Kinetic double-mixing experiments revealed that, in intact SlyD, the IF domain folds and unfolds about 1000-fold more rapidly than the FKBP domain, and that it is strongly stabilized when linked with the folded FKBP domain. The unfolding limbs of the kinetic chevrons of SlyD show a strong downward curvature. This deviation from linearity is not caused by a transition-state movement, as often assumed, but by the accumulation of a silent unfolding intermediate at high denaturant concentrations. In this kinetic intermediate, the FKBP domain is still folded, whereas the IF domain is already unfolded.

Studying the folding of multidomain proteins

There have been relatively few detailed comprehensive studies of the folding of protein domains (or modules) in the context of their natural covalently linked neighbors. This is despite the fact that a significant proportion of the proteome consists of multidomain proteins. In this review we highlight some key experimental investigations of the folding of multidomain proteins to draw attention to the difficulties that can arise in analyzing such systems. The evidence suggests that interdomain interactions can significantly affect stability, folding, and unfolding rates. However, preliminary studies suggest that folding pathways are unaffected-to this extent domains can be truly considered to be independent folding units. Nonetheless, it is clear that interactions between domains cannot be ignored, in particular when considering the effects of mutations.

Construction of recB-recD genetic fusion and functional analysis of RecBDC fusion enzyme in Escherichia coli

Background

recD, located between recB and argA, encodes the smallest polypeptide (60 kDa) of the heterotrimeric enzyme RecBCD in Escherichia coli. RecD is a 5'-3' helicase and is required for the nuclease activity of RecBCD and for tight binding to dsDNA ends. Here, we have tested the hypothesis that RecD regulates the structure and activities of RecBCD, including RecA loading.

Results

To characterize its regulatory functions, recD was genetically fused to recB through deletion and substitution mutations. The recB-recD fusion led to a decreased amount of the heterotrimer. Both fusion mutants proved to be recombination proficient, viable and resistant to DNA damaging agents, and to have DNA unwinding, ATP-dependent dsDNA exonuclease and Chi genetic activities.

Conclusion

Our findings suggest that the recB-recD fusion may form a RecBD fusion protein and therefore affect RecD assembly, but this does not change the three-dimensional structure of the heterotrimer.

The role of glycosylation and domain interactions in the thermal stability of human angiotensin-converting enzyme

The N and C domains of somatic angiotensin-converting enzyme (sACE) differ in terms of their substrate specificity, inhibitor profiling, chloride dependency and thermal stability. The C domain is thermally less stable than sACE or the N domain. Since both domains are heavily glycosylated, the effect of glycosylation on their thermal stability was investigated by assessing their catalytic and physicochemical properties. Testis ACE (tACE) expressed in mammalian cells, mammalian cells in the presence of a glucosidase inhibitor and insect cells yielded proteins with altered catalytic and physicochemical properties, indicating that the more complex glycans confer greater thermal stabilization. Furthermore, a decrease in tACE and N-domain N-glycans using site-directed mutagenesis decreased their thermal stability, suggesting that certain N-glycans have an important effect on the protein's thermodynamic properties. Evaluation of the thermal stability of sACE domain swopover and domain duplication mutants, together with sACE expressed in insect cells, showed that the C domain contained in sACE is less dependent on glycosylation for thermal stabilization than a single C domain, indicating that stabilizing interactions between the two domains contribute to the thermal stability of sACE and are decreased in a C-domain-duplicating mutant.

Subdomain competition, cooperativity, and topological frustration in the folding of CheY

The folding of multidomain proteins often proceeds in a hierarchical fashion with individual domains folding independent of one another. A large single-domain protein, however, can consist of multiple modules whose folding may be autonomous or interdependent in ways that are unclear. We used coarse-grained simulations to explore the folding landscape of the two-subdomain bacterial response regulator CheY. Thermodynamic and kinetic characterization shows the landscape to be highly analogous to the four-state landscape reported for another two-subdomain protein, T4 lysozyme. An on-pathway intermediate structured in the more stable nucleating subdomain was observed, as were transient states frustrated in off-pathway contacts prematurely structured in the weaker subdomain. Local unfolding, or backtracking, was observed in the frustrated state before the native conformation could be reached. Nonproductive frustration was attributable to competition for van der Waals contacts between the two subdomains. In an accompanying article, stopped-flow kinetic measurements support an off-pathway burst-phase intermediate, seemingly consistent with our prediction of early frustration in the folding landscape of CheY. Comparison of the folding mechanisms for CheY, T4 lysozyme, and interleukin-1 beta leads us to postulate that subdomain competition is a general feature of large single-domain proteins with multiple folding modules.

Statistical analysis of interface similarity in crystals of homologous proteins

Many proteins function as homo-oligomers and are regulated via their oligomeric state. For some proteins, the stoichiometry of homo-oligomeric states under various conditions has been studied using gel filtration or analytical ultracentrifugation experiments. The interfaces involved in these assemblies may be identified using cross-linking and mass spectrometry, solution-state NMR, and other experiments. However, for most proteins, the actual interfaces that are involved in oligomerization are inferred from X-ray crystallographic structures using assumptions about interface surface areas and physical properties. Examination of interfaces across different Protein Data Bank (PDB) entries in a protein family reveals several important features. First, similarities in space group, asymmetric unit size, and cell dimensions and angles (within 1%) do not guarantee that two crystals are actually the same crystal form, containing similar relative orientations and interactions within the crystal. Conversely, two crystals in different space groups may be quite similar in terms of all the interfaces within each crystal. Second, NMR structures and an existing benchmark of PDB crystallographic entries consisting of 126 dimers as well as larger structures and 132 monomers were used to determine whether the existence or lack of common interfaces across multiple crystal forms can be used to predict whether a protein is an oligomer or not. Monomeric proteins tend to have common interfaces across only a minority of crystal forms, whereas higher-order structures exhibit common interfaces across a majority of available crystal forms. The data can be used to estimate the probability that an interface is biological if two or more crystal forms are available. Finally, the Protein Interfaces, Surfaces, and Assemblies (PISA) database available from the European Bioinformatics Institute is more consistent in identifying interfaces observed in many crystal forms compared with the PDB and the European Bioinformatics Institute's Protein Quaternary Server (PQS). The PDB, in particular, is missing highly likely biological interfaces in its biological unit files for about 10% of PDB entries.

GANDivAWeb: a web server for detecting early folding units ("foldons") from protein 3D structures

Background

It has long been known that small regions of proteins tend to fold independently and are then stabilized by interactions between these distinct subunits or modules. Such units, also known as autonomous folding units (AFUs) or"foldons" play a key role in protein folding. A knowledge of such early folding units has diverse applications in protein engineering as well as in developing an understanding of the protein folding process. Such AFUs can also be used as model systems in order to study the structural organization of proteins.

Results

In an earlier work, we had utilized a global network partitioning algorithm to identify modules in proteins. We had shown that these modules correlate well with AFUs. In this work, we have developed a webserver, GANDivAWeb, to identify early folding units or "foldons" in networks using the algorithm described earlier. The website has three functionalities: (a) It is able to display information on the modularity of a database of 1420 proteins used in the original work, (b) It can take as input an uploaded PDB file, identify the modules using the GANDivA algorithm and email the results back to the user and (c) It can take as input an uploaded PDB file and a results file (obtained from functionality (b)) and display the results using the embedded viewer. The results include the module decomposition of the protein, plots of cartoon representations of the protein colored by module identity and connectivity as well as contour plots of the hydrophobicity and relative accessible surface area (RASA) distributions.

Conclusion

We believe that the GANDivAWeb server, will be a useful tool for scientists interested in the phenomena of protein folding as well as in protein engineering. Our tool not only provides a knowledge of the AFUs through a natural graph partitioning approach but is also able to identify residues that are critical during folding. It is our intention to use this tool to study the topological determinants of protein folding by analyzing the topological changes in proteins over the unfolding/folding pathways.

Kinetically driven refolding of the hyperstable EBNA1 origin DNA-binding dimeric beta-barrel domain into amyloid-like spherical oligomers

The Epstein-Barr nuclear antigen 1 (EBNA1) is essential for DNA replication and episome segregation of the viral genome, and participates in other gene regulatory processes of the Epstein-Barr virus in benign and malignant diseases related to this virus. Despite the participation of other regions of the protein in evading immune response, its DNA binding, dimeric beta-barrel domain (residues 452-641) is necessary and sufficient for the main functions. This domain has an unusual topology only shared by another viral origin binding protein (OBP), the E2 DNA binding domain of papillomaviruses. Both the amino acid and DNA target sequences are completely different for these two proteins, indicating a link between fold conservation and function. In this work we investigated the folding and stability of the DNA binding domain of EBNA1 OBP and found it is extremely resistant to chemical, temperature, and pH denaturation. The thiocyanate salt of guanidine is required for obtaining a complete transition to a monomeric unfolded state. The unfolding reaction is extremely slow and shows a marked uncoupling between tertiary and secondary structure, indicating the presence of intermediate species. The Gdm.SCN unfolded protein refolds to fully soluble and spherical oligomeric species of 1.2 MDa molecular weight, with identical fluorescence centre of spectral mass but different intensity and different secondary structure. The refolded spherical oligomers are substantially less stable than the native recombinant dimer. In keeping with the substantial structural rearrangement in the oligomers, the spherical oligomers do not bind DNA, indicating that the DNA binding site is either disrupted or participates in the oligomerization interface. The puzzling extreme stability of a dimeric DNA binding domain from a protein from a human infecting virus in addition to a remarkable kinetically driven folding where all molecules do not return to the most stable original species suggests a co-translational and directional folding of EBNA1 in vivo, possibly assisted by folding accessory proteins. Finally, the oligomers bind Congo red and thioflavin-T, both characteristic of repetitive beta-sheet elements of structure found in amyloids and their soluble precursors. The stable nature of the "kinetically trapped" oligomers suggest their value as models for understanding amyloid intermediates, their toxic nature, and the progress to amyloid fibers in misfolding diseases. The possible role of the EBNA1 spherical oligomers in the virus biology is discussed.

Network scaling invariants help to elucidate basic topological principles of proteins

The structural architecture of proteins continues to be an area of active research. Despite the difference in models dealing with the way proteins fold into their tertiary structures, it is recognized that small regions of proteins tend to fold independently and are then stabilized by interactions between these distinct subunits. However, there are a number of different definitions of what comprises an independent subunit. In the belief that an unequivocal definition of a domain must be based on the most fundamental property of protein 3D structure, namely, the adjacency matrix of inter-residues contact, we adopt a network representation of the protein. In this work, we used a well-established, global method for identifying modules in networks, without any specific reference to the kind of network being analyzed. The algorithm converges toward the maximization of the modularity of the given protein network and, in doing so, allows the representation of the residues of the protein in terms of their intramodule degree, z, and participation coefficient, P. We demonstrate that the labeling of residues in terms of these invariants allows for information-rich representations of the studied proteins as well as to sketch a new way to link sequence, structure, and the dynamical properties of proteins. We discovered a strong invariant character of protein molecules in terms of P/z characterization, pointing to a common topological design of all protein structures. This invariant representation, applied to different protein systems, enabled us to identify the possible functional role of high P/z residues during the folding process. Additionally, we observe a hierarchical behavior of protein structural organization that provides a sequence-secondary-tertiary structure link. The discovery of similar and repeatable scaling laws at different level of definitions going from hydrophobicity patterning along the sequence up to the size of an autonomous folding unit (AFU) and general contact distribution of the entire molecule suggest a hierarchical-like behavior of protein architecture. This implies the possibility to select different privileged scales of observation for deriving useful information on protein systems.

Exploring subdomain cooperativity in T4 lysozyme I: structural and energetic studies of a circular permutant and protein fragment

Small proteins are generally observed to fold in an apparent two-state manner. Recently, however, more sensitive techniques have demonstrated that even seemingly single-domain proteins are actually made up of smaller subdomains. T4 lysozyme is one such protein. We explored the relative autonomy of its two individual subdomains and their contribution to the overall stability of T4 lysozyme by examining a circular permutation (CP13*) that relocates the N-terminal A-helix, creating subdomains that are contiguous in sequence. By determining the high-resolution structure of CP13* and characterizing its energy landscape using native state hydrogen exchange (NSHX), we show that connectivity between the subdomains is an important determinant of the energetic cooperativity but not structural integrity of the protein. The circular permutation results in a protein more easily able to populate a partially unfolded form in which the C-terminal subdomain is folded and the N-terminal subdomain is unfolded. We also created a fragment model of this intermediate and demonstrate using X-ray crystallography that its structure is identical to the corresponding residues in the full-length protein with the exception of a small network of hydrophobic interactions. In sum, we conclude that the C-terminal subdomain dominates the energetics of T4 lysozyme folding, and the A-helix serves an important role in coupling the two subdomains.

Napin from Brassica juncea: thermodynamic and structural analysis of stability

The napin from Brassica juncea, oriental mustard, is highly thermostable, proteolysis resistant and allergenic in nature. It consists of two subunits - one small (29 amino acid residues) and one large (86 amino acids residues) - held together by disulfide bonds. The thermal unfolding of napin has been followed by differential scanning calorimetry (DSC) and circular dichroism (CD) measurements. The thermal unfolding is characterized by a three state transition with T(M1) and T(M2) at 323.5 K and 335.8 K, respectively; DeltaC(P1) and DeltaC(P2) are 2.05 kcal mol(-1) K(-1) and 1.40 kcal mol(-1) K(-1), respectively. In the temperature range 310-318 K, the molecule undergoes dimerisation. Isothermal equilibrium unfolding by guanidinium hydrochloride also follows a three state transition, N <_-_-> I <_-_-> U with DeltaG(1H2O) and DeltaG(2H2O) values of 5.2 kcal mol(-1) and 5.1 kcal mol(-1) at 300 K, respectively. Excess heat capacity values obtained, are similar to those obtained from DSC measurements. There is an increase in hydrodynamic radius from 20 A to 35.0 A due to unfolding by guanidinium hydrochloride. In silico alignment of sequences of napin has revealed that the internal repeats (40%) spanning residues 31 to 60 and 73 to 109 are conserved in all Brassica species. The internal repeats may contribute to the greater stability of napin. A thorough understanding of the structure and stability of these proteins is essential before they can be exploited for genetic improvements for nutrition.

Multiple Intermediate Conformations of Jack Bean Urease at Low pH: Anion-induced Refolding

Structural and functional characteristics of jack bean urease (JBU), a hexameric enzyme having identical subunits, were investigated under neutral as well as acidic conditions by using CD, fluorescence, ANS binding and enzyme activity measurements. At low pH and low ionic strength, JBU exists in a partially unfolded state (U(A)-state), having predominantly beta structure and no tertiary interactions along with a strong ANS binding. Addition of salts like NaCl, KCl and Na(2)SO(4) to the U(A)-state induces refolding resulting in structural propensities similar to that of native hexamer. Moreover, at low concentrations, GuHCl behaves like an anion by inducing refolding of the U(A)-state. The anion-induced refolded state (I(A)-state) is more stable than U(A)-state and the stability is nearly equal to that of the native protein against chemical-induced and thermal denaturation. Overall, these observations support a model of protein folding for a multimeric protein where certain conformations (ensembles of substates) of low energy prevail and populated under non-native conditions with different stability.

Assessing the accuracy of ancestral protein reconstruction methods

The phylogenetic inference of ancestral protein sequences is a powerful technique for the study of molecular evolution, but any conclusions drawn from such studies are only as good as the accuracy of the reconstruction method. Every inference method leads to errors in the ancestral protein sequence, resulting in potentially misleading estimates of the ancestral protein's properties. To assess the accuracy of ancestral protein reconstruction methods, we performed computational population evolution simulations featuring near-neutral evolution under purifying selection, speciation, and divergence using an off-lattice protein model where fitness depends on the ability to be stable in a specified target structure. We were thus able to compare the thermodynamic properties of the true ancestral sequences with the properties of “ancestral sequences” inferred by maximum parsimony, maximum likelihood, and Bayesian methods. Surprisingly, we found that methods such as maximum parsimony and maximum likelihood that reconstruct a “best guess” amino acid at each position overestimate thermostability, while a Bayesian method that sometimes chooses less-probable residues from the posterior probability distribution does not. Maximum likelihood and maximum parsimony apparently tend to eliminate variants at a position that are slightly detrimental to structural stability simply because such detrimental variants are less frequent. Other properties of ancestral proteins might be similarly overestimated. This suggests that ancestral reconstruction studies require greater care to come to credible conclusions regarding functional evolution. Inferred functional patterns that mimic reconstruction bias should be reevaluated.

It is now possible to apply computational methods to known current protein sequences to recreate the sequences of ancestral proteins. By synthesising these proteins and measuring their properties in the laboratory, we can gain much information about the nature of evolution, better understand how proteins change and adapt over time, and develop insights into the environments of ancient organisms. Unfortunately, the accuracy of these reconstructions is difficult to evaluate. We simulate protein evolution using a simplified computational model and apply the various reconstruction methods to the sequences that arise from our simulations. Because we have the complete record of the evolutionary history, we can evaluate the reconstruction accuracy directly. We demonstrate that the reconstruction procedures in common use may have a bias toward overestimating the properties of these ancestral proteins, opposite to what has been assumed previously. An alternative method of creating these sequences is presented, Bayesian sampling, that can eliminate this bias and provide more robust conclusions.

Structural diversity of protein segments follows a power-law distribution

The local structures of protein segments were classified and their distribution was analyzed to explore the structural diversity of proteins. Representative proteins were divided into short segments using a sliding L-residue window. Each set of local structures consisting of consecutive 1-31 amino acids was classified using a single-pass clustering method. The results demonstrate that the local structures of proteins are very unevenly distributed in the protein universe. The distribution of local structures of relatively long segments shows a power-law behavior that is formulated well by Zipf's law, implying that a protein structure possesses recursive and fractal characteristics. The degree of effective conformational freedom per residue as well as the structure entropy per residue decreases gradually with an increasing value of L and then converges to constant values. This suggests that the number of protein conformations resides within the range between 1.2L and 1.5L and that 10- to 20-residue segments are already proteinlike in terms of their structural diversity.

Dynamic protein domains: identification, interdependence, and stability

Existing methods of domain identification in proteins usually provide no information about the degree of domain independence and stability. However, this information is vital for many areas of protein research. The recently developed hierarchical clustering of correlation patterns (HCCP) technique provides machine-based domain identification in a computationally simple and physically consistent way. Here we present the modification of this technique, which not only allows determination of the most plausible number of dynamic domains but also makes it possible to estimate the degree of their independence (the extent of correlated motion) and stability (the range of environmental conditions, where domains remain intact). With this technique we provided domain assignments and calculated intra- and interdomain correlations and interdomain energies for >2500 test proteins. It is shown that mean intradomain correlation of motions can serve as a quantitative criterion of domain independence, and the HCCP stability gap is a measure of their stability. Our data show that the motions of domains with high stability are usually independent. In contrast, the domains with moderate stability usually exhibit a substantial degree of correlated motions. It is shown that in multidomain proteins the domains are most stable if they are of similar size, and this correlates with the observed abundance of such proteins.

Interdomain side-chain interactions in human gammaD crystallin influencing folding and stability

Human gammaD crystallin (HgammaD-Crys) is a two domain, beta-sheet eye lens protein that must remain soluble throughout life for lens transparency. Single amino acid substitutions of HgammaD-Crys are associated with juvenile-onset cataracts. Features of the interface between the two domains conserved among gamma-crystallins are a central six-residue hydrophobic cluster, and two pairs of interacting residues flanking the cluster. In HgammaD-Crys these pairs are Gln54/Gln143 and Arg79/Met147. We previously reported contributions of the hydrophobic cluster residues to protein stability. In this study alanine substitutions of the flanking residue pairs were constructed and analyzed. Equilibrium unfolding/refolding experiments at 37 degrees C revealed a plateau in the unfolding/refolding transitions, suggesting population of a partially folded intermediate with a folded C-terminal domain (C-td) and unfolded N-terminal domain (N-td). The N-td was destabilized by substituting residues from both domains. In contrast, the C-td was not significantly affected by substitutions of either domain. Refolding rates of the N-td were significantly decreased for mutants of either domain. In contrast, refolding rates of the C-td were similar to wild type for mutants of either domain. Therefore, domain interface residues of the folded C-td probably nucleate refolding of the N-td. We suggest that these residues stabilize the native state by shielding the central hydrophobic cluster from solvent. Glutamine and methionine side chains are among the residues covalently damaged in aged and cataractous lenses. Such damage may generate partially unfolded, aggregation- prone conformations of HgammaD-Crys that could be significant in cataract.

Stability domains, substrate-induced conformational changes, and hinge-bending motions in a psychrophilic phosphoglycerate kinase. A microcalorimetric study

The cold-active phosphoglycerate kinase from the Antarctic bacterium Pseudomonas sp. TACII18 exhibits two distinct stability domains in the free, open conformation. It is shown that these stability domains do not match the structural N- and C-domains as the heat-stable domain corresponds to about 80 residues of the C-domain, including the nucleotide binding site, whereas the remaining of the protein contributes to the main heat-labile domain. This was demonstrated by spectroscopic and microcalorimetric analyses of the native enzyme, of its mutants, and of the isolated recombinant structural domains. It is proposed that the heat-stable domain provides a compact structure improving the binding affinity of the nucleotide, therefore increasing the catalytic efficiency at low temperatures. Upon substrate binding, the enzyme adopts a uniformly more stable closed conformation. Substrate-induced stability changes suggest that the free energy of ligand binding is converted into an increased conformational stability used to drive the hinge-bending motions and domain closure.

Conformational change in the C-terminal domain is responsible for the initiation of creatine kinase thermal aggregation

Protein conformational changes may be associated with particular properties such as its function, transportation, assembly, tendency to aggregate, and potential cytotoxicity. In this research, the conformational change that is responsible for the fast destabilization and aggregation of rabbit muscle creatine kinase (EC 2.7.3.2) induced by heat was studied by intrinsic fluorescence and infrared spectroscopy. A pretransitional change of the tryptophan microenvironments was found from the intrinsic fluorescence spectra. A further analysis of the infrared spectra using quantitative second-derivative and two-dimensional correlation analysis indicated that the changes of the beta-sheet structures in the C-terminal domain and the loops occurred before the formation of intermolecular cross-beta-sheet structures and the unfolding of alpha-helices. These results suggested that the pretransitional conformational changes in the active site and the C-terminal domain might result in the modification of the domain-domain interactions and the formation of an inactive dimeric form that was prone to aggregate. Our results highlighted the fact that some minor conformational changes, which were usually negligible or undetectable by normal methods, might play a crucial role in protein stability and aggregation. Our results also suggested that the changes in domain-domain interactions, but not the dissociation of the dimer, might play a crucial role in the thermal denaturation and aggregation of this dimeric two-domain protein.

Instabilities of proteins: theoretical aspects, degradation products and methods for their detection

Recombinant DNA technology has led to a significant increase in the number of peptide and protein based pharmaceuticals, giving a new approach to combat poorly controlled diseases. This particular development has been reached in the last two decades. However, proteins are highly susceptible of physical and chemical degradation resulting in a decrease or complete loss of biological activities. Reasons for their physical and chemical instabilities and the methods for their examination, become a challenge for the pharmaceutical scientists for successful development of stabile protein - based pharmaceuticals. The stability of protein - based pharmaceuticals is significant in terms of their pharmaceutical quality and biological activity. In addition, a right choice of suitable analytical methods is needed in order to detect an early formation of degradation products or modified forms.

Lipid binding region (2303-2332) is involved in aggregation of recombinant human FVIII (rFVIII)

Factor VIII (FVIII) is a multi-domain protein that is important in the clotting cascade. Its deficiency causes Hemophilia A, a bleeding disorder. The unfolding of protein domains can lead to physical instability such as aggregation, and hinder their use in replacement therapy. It has been shown that the aggregation of rFVIIII is initiated by small fluctuations in the protein's tertiary structure (Grillo et al., 2001, Biochemistry 40:586-595). We have investigated the domain(s) involved in the initiation of aggregation using circular dichroism (CD), size exclusion chromatography (SEC), fluorescence anisotropy, domain specific antibody binding, and clotting activity studies. The studies indicated that aggregation may be initiated as a result of conformational change in the C2 domain encompassing the lipid-binding region (2303-2332). The presence of O-phospho-L-Serine (OPLS), which binds to the lipid-binding region of FVIII, prevented aggregation of the protein.

Domain coupling in a multimodular cellobiohydrolase CbhA fromClostridium thermocellum

Cellobiohydrolase A (CbhA) from Clostridium thermocellum is composed of an N-terminal carbohydrate-binding domain 4 (CBD4), an immunoglobulin-like domain (Ig), a glycoside hydrolase 9 (GH9), X1(1) and X1(2) domains, a CBD3, and a dockerin domain. All domains, except the Ig, bind Ca2+. The following constructs were made: X1(2), X1(1)X1(2), CBD3, X1(1)X1(2)-CBD3, Ig, GH9, Ig-GH9, Ig-GH9-X1(1)X1(2), and Ig-GH9-X1(1)X1(2)-CBD3. Interactions between domains in (1) buffer, (2) with Ca2+, or (3) ethylenediaminetetraacetic acid (EDTA) were studied by differential scanning calorimetry. Thermal unfoldings of all constructs were irreversible. Calcium increased T(d) and cooperativity of unfolding. Multi-domain constructs exhibited more cooperative unfolding in buffer and in the presence of EDTA than did individual domains. They denatured by mechanism simpler than expected from their modular architecture. The results indicate that domain coupling in thermophilic proteins constitutes a significant stabilizing factor.

Contributions of hydrophobic domain interface interactions to the folding and stability of human gammaD-crystallin

Human gammaD-crystallin (HgammaD-Crys) is a monomeric eye lens protein composed of two highly homologous beta-sheet domains. The domains interact through interdomain side chain contacts forming two structurally distinct regions, a central hydrophobic cluster and peripheral residues. The hydrophobic cluster contains Met43, Phe56, and Ile81 from the N-terminal domain (N-td) and Val132, Leu145, and Val170 from the C-terminal domain (C-td). Equilibrium unfolding/refolding of wild-type HgammaD-Crys in guanidine hydrochloride (GuHCl) was best fit to a three-state model with transition midpoints of 2.2 and 2.8 M GuHCl. The two transitions likely corresponded to sequential unfolding/refolding of the N-td and the C-td. Previous kinetic experiments revealed that the C-td refolds more rapidly than the N-td. We constructed alanine substitutions of the hydrophobic interface residues to analyze their roles in folding and stability. After purification from E. coli, all mutant proteins adopted a native-like structure similar to wild type. The mutants F56A, I81A, V132A, and L145A had a destabilized N-td, causing greater population of the single folded domain intermediate. Compared to wild type, these mutants also had reduced rates for productive refolding of the N-td but not the C-td. These data suggest a refolding pathway where the domain interface residues of the refolded C-td act as a nucleating center for refolding of the N-td. Specificity of domain interface interactions is likely important for preventing incorrect associations in the high protein concentrations of the lens nucleus.

Time-dependent changes in the denatured state(s) influence the folding mechanism of an all beta-sheet protein

Newt fibroblast growth factor (nFGF-1) is an approximately 15-kDa all beta-sheet protein devoid of disulfide bonds. Urea-induced equilibrium unfolding of nFGF-1, monitored by steady state fluorescence and far-UV circular dichroism spectroscopy, is cooperative with no detectable intermediate(s). Urea-induced unfolding of nFGF-1 is reversible, but the percentage of the protein recovered in the native state depends on the time of incubation of the protein in the denaturant. The yield of the protein in the native state decreases with the increase in time of incubation in the denaturant. The failure of the protein to refold to its native state is not due to trivial chemical reactions that could possibly occur upon prolonged incubation in the denaturant. (1)H-(15)N heteronuclear single quantum coherence (HSQC) spectra, limited proteolytic digestion, and fluorescence data suggest that the misfolded state(s) of nFGF-1 has structural features resembling that of the denatured state(s). GroEL, in the presence of ATP, is observed to rescue the protein from being trapped in the misfolded state(s). (1)H-(15)N HSQC data of nFGF-1, acquired in the denatured state(s) (in 8 m urea), suggest that the protein undergoes subtle time-dependent structural changes in the denaturant. To our knowledge, this report for the first time demonstrates that the commitment to adapt unproductive pathways leading to protein misfolding/aggregation occurs in the denatured state ensemble.

Pro-sequence and Ca2+-binding: Implications for Folding and Maturation of Ntn-hydrolase Penicillin Amidase from E.coli

Penicillin amidase (PA) is a bacterial periplasmic enzyme synthesized as a pre-pro-PA precursor. The pre-sequence mediates membrane translocation. The intramolecular pro-sequence is expressed along with the A and B chains but is rapidly removed in an autocatalytic manner. In extensive studies we show here that the pro-peptide is required for the correct folding of PA. Pro-PA and PA unfold via a biphasic transition that is more pronounced in the case of PA. According to size-exclusion chromatography and limited proteolysis experiments, the inflection observed in the equilibrium unfolding curves corresponds to an intermediate in which the N-terminal domain (A-chain) still possesses native-like topology, whereas the B-chain is unfolded to a large extent. In a series of in vitro experiments with a slow processing mutant pro-PA, we show that the pro-sequence in cis functions as a folding catalyst and accelerates the folding rate by seven orders of magnitude. In the absence of the pro-domain the PA refolds to a stable inactive molten globule intermediate that has native-like secondary but little tertiary structure. The pro-sequence of the homologous Alcaligenes faecalis PA can facilitate the folding of the hydrolase domain of Escherichia coli PA when added in trans (as a separate polypeptide chain). The isolated pro-sequence has a random structure in solution. However, difference circular dichroism spectra of native PA and native PA with pro-peptide added in trans suggest that the pro-sequence adopts an alpha-helical conformation in the context of the mature PA molecule. Furthermore, our results establish that Ca2+, found in the crystal structure, is not directly involved in the folding process. The cation shifts the equilibrium towards the native state and facilitates the autocatalytic processing of the pro-peptide.

Preparation and characterization of geodin. A betagamma-crystallin-type protein from a sponge

Geodin is a protein encoded by a sponge gene homologous to genes from the betagamma-crystallins superfamily. The interest for this crystallin-type protein stems from the phylogenesis of porifera, commonly called sponges, the earliest divergence event in the history of metazoans. Here we report the preparation of geodin as a recombinant protein from Escherichia coli, its characterization through physico-chemical analyses, and a model of its 3D structure based on homology modelling. Geodin is a monomeric protein of about 18 kDa, with an all-beta structure, as all other crystallins in the superfamily, but more prone to unfold in the presence of chemical denaturants, when compared with other homologues from the superfamily. Its thermal unfolding, studied by far- and near-CD, and by calorimetry, is described by a two-state model. Geodin appears to be structurally similar in many respects to the bacterial protein S crystallin, with which it also shares a significant, albeit more modest stabilizing effect exerted by calcium ions. These results suggest that the crystallin-type structural scaffold, employed in the evolution of bacteria and moulds, was successfully recruited very early in the evolution of metazoa.

Protein folding stabilizing time measurement: A direct folding process and three-dimensional random walk simulation

Protein particles undergo Brownian motion and collisions in solution. The diffusive collisions may lead to aggregation. For proteins to fold successfully the process has to occur quickly and before significant collision takes place. The speed of protein folding was deduced by studying the correlation time of a lysozyme refolding process from autocorrelation function analysis of the mean collision time and aggregation/soluble ratio of protein. It is a measure of time before which an aggregate can be formed and also is the time measure for a protein to fold into a stable state. We report on the protein folding stabilizing time of a lysozyme system to be 25.5-27.5 micros (<+/-4%) between 295 and 279K via direct folding experimental studies, supported by a three-dimensional random walk simulation of diffusion-limited aggregation model. Aggregation is suppressed when the protein is folded to a stable form. Spontaneous folding and diffusion-limited aggregation are antagonistic in nature. Meanwhile, the resultant aggresome, suggested by Raman and mass spectroscopy, may be formed by cross-linkages of disulfide bonds and hydrophobic interactions.

Probing folding and fluorescence quenching in human gammaD crystallin Greek key domains using triple tryptophan mutant proteins

Human gammaD crystallin (HgammaD-Crys), a major component of the human eye lens, is a 173-residue, primarily beta-sheet protein, associated with juvenile and mature-onset cataracts. HgammaD-Crys has four tryptophans, with two in each of the homologous Greek key domains, which are conserved throughout the gamma-crystallin family. HgammaD-Crys exhibits native-state fluorescence quenching, despite the absence of ligands or cofactors. The tryptophan absorption and fluorescence quenching may influence the lens response to ultraviolet light or the protection of the retina from ambient ultraviolet damage. To provide fluorescence reporters for each quadrant of the protein, triple mutants, each containing three tryptophan-to-phenylalanine substitutions and one native tryptophan, have been constructed and expressed. Trp 42-only and Trp 130-only exhibited fluorescence quenching between the native and denatured states typical of globular proteins, whereas Trp 68-only and Trp 156-only retained the anomalous quenching pattern of wild-type HgammaD-Crys. The three-dimensional structure of HgammaD-Crys shows Tyr/Tyr/His aromatic cages surrounding Trp 68 and Trp 156 that may be the source of the native-state quenching. During equilibrium refolding/unfolding at 37 degrees C, the tryptophan fluorescence signals indicated that domain I (W42-only and W68-only) unfolded at lower concentrations of GdnHCl than domain II (W130-only and W156-only). Kinetic analysis of both the unfolding and refolding of the triple-mutant tryptophan proteins identified an intermediate along the HgammaD-Crys folding pathway with domain I unfolded and domain II intact. This species is a candidate for the partially folded intermediate in the in vitro aggregation pathway of HgammaD-Crys.

Probing the unfolding region of ribonuclease A by site-directed mutagenesis

Ribonuclease A contains two exposed loop regions, around Ala20 and Asn34. Only the loop around Ala20 is sufficiently flexible even under native conditions to allow cleavage by nonspecific proteases. In contrast, the loop around Asn34 (together with the adjacent beta-sheet around Thr45) is the first region of the ribonuclease A molecule that becomes susceptible to thermolysin and trypsin under unfolding conditions. This second region therefore has been suggested to be involved in early steps of unfolding and was designated as the unfolding region of the ribonuclease A molecule. Consequently, modifications in this region should have a great impact on the unfolding and, thus, on the thermodynamic stability. Also, if the Ala20 loop contributes to the stability of the ribonuclease A molecule, rigidification of this flexible region should stabilize the entire protein molecule. We substituted several residues in both regions without any dramatic effects on the native conformation and catalytic activity. As a result of their remarkably differing stability, the variants fell into two groups carrying the mutations: (a) A20P, S21P, A20P/S21P, S21L, or N34D; (b) L35S, L35A, F46Y, K31A/R33S, L35S/F46Y, L35A/F46Y, or K31A/R33S/F46Y. The first group showed a thermodynamic and kinetic stability similar to wild-type ribonuclease A, whereas both stabilities of the variants in the second group were greatly decreased, suggesting that the decrease in DeltaG can be mainly attributed to an increased unfolding rate. Although rigidification of the Ala20 loop by introduction of proline did not result in stabilization, disturbance of the network of hydrogen bonds and hydrophobic interactions that interlock the proposed unfolding region dramatically destabilized the ribonuclease A molecule.

Enhanced Ligand Affinity for Receptors in which Components of the Binding Site Are Independently Mobile

Using calmodulin antagonism as a model, it is demonstrated that, under circumstances in which binding sites are motionally independent, it is possible to create bifunctional ligands that bind with significant affinity enhancement over their monofunctional counterparts. Suitable head groups were identified by using a semiquantitative screen of monofunctional tryptophan analogs. Two bifunctional ligands, which contained two copies of the highest-affinity head group tethered by rigid linkers, were synthesized. The bifunctional ligands bound to calmodulin with a stoichiometry of 1:1 and with an affinity enhancement over their monofunctional counterparts; the latter bound with a stoichiometry of 2:1 ligand:protein. A lower limit to the effective concentrations of the domains of calmodulin relative to each other (0.2-2 mM) was determined. A comparable effective concentration was achieved for bifunctional ligands based on higher-affinity naphthalene sulphonamide derivatives.

Ageing and vision: structure, stability and function of lens crystallins

The alpha-, beta- and gamma-crystallins are the major protein components of the vertebrate eye lens, alpha-crystallin as a molecular chaperone as well as a structural protein, beta- and gamma-crystallins as structural proteins. For the lens to be able to retain life-long transparency in the absence of protein turnover, the crystallins must meet not only the requirement of solubility associated with high cellular concentration but that of longevity as well. For proteins, longevity is commonly assumed to be correlated with long-term retention of native structure, which in turn can be due to inherent thermodynamic stability, efficient capture and refolding of non-native protein by chaperones, or a combination of both. Understanding how the specific interactions that confer intrinsic stability of the protein fold are combined with the stabilizing effect of protein assembly, and how the non-specific interactions and associations of the assemblies enable the generation of highly concentrated solutions, is thus of importance to understand the loss of transparency of the lens with age. Post-translational modification can have a major effect on protein stability but an emerging theme of the few studies of the effect of post-translational modification of the crystallins is one of solubility and assembly. Here we review the structure, assembly, interactions, stability and post-translational modifications of the crystallins, not only in isolation but also as part of a multi-component system. The available data are discussed in the context of the establishment, the maintenance and finally, with age, the loss of transparency of the lens. Understanding the structural basis of protein stability and interactions in the healthy eye lens is the route to solve the enormous medical and economical problem of cataract.

Effect of dimer dissociation on activity and thermostability of the alpha-glucuronidase from Geobacillus stearothermophilus: dissecting the different oligomeric forms of family 67 glycoside hydrolases

The oligomeric organization of enzymes plays an important role in many biological processes, such as allosteric regulation, conformational stability and thermal stability. alpha-Glucuronidases are family 67 glycosidases that cleave the alpha-1,2-glycosidic bond between 4-O-methyl-D-glucuronic acid and xylose units as part of an array of hemicellulose-hydrolyzing enzymes. Currently, two crystal structures of alpha-glucuronidases are available, those from Geobacillus stearothermophilus (AguA) and from Cellvibrio japonicus (GlcA67A). Both enzymes are homodimeric, but surprisingly their dimeric organization is different, raising questions regarding the significance of dimerization for the enzymes' activity and stability. Structural comparison of the two enzymes suggests several elements that are responsible for the different dimerization organization. Phylogenetic analysis shows that the alpha-glucuronidases AguA and GlcA67A can be classified into two distinct subfamilies of bacterial alpha-glucuronidases, where the dimer-forming residues of each enzyme are conserved only within its own subfamily. It seems that the different dimeric forms of AguA and GlcA67A represent the two alternative dimeric organizations of these subfamilies. To study the biological significance of the dimerization in alpha-glucuronidases, we have constructed a monomeric form of AguA by mutating three of its interface residues (W328E, R329T, and R665N). The activity of the monomer was significantly lower than the activity of the wild-type dimeric AguA, and the optimal temperature for activity of the monomer was around 35 degrees C, compared to 65 degrees C of the wild-type enzyme. Nevertheless, the melting temperature of the monomeric protein, 72.9 degrees C, was almost identical to that of the wild-type, 73.4 degrees C. It appears that the dimerization of AguA is essential for efficient catalysis and that the dissociation into monomers results in subtle conformational changes in the structure which indirectly influence the active site region and reduce the activity. Structural and mechanistic explanations for these effects are discussed.

Do Proteins Always Benefit from a Stability Increase? Relevant and Residual Stabilisation in a Three-state Protein by Charge Optimisation

The vast majority of our knowledge on protein stability arises from the study of simple two-state models. However, proteins displaying equilibrium intermediates under certain conditions abound and it is unclear whether the energetics of native/intermediate equilibria is well represented in current knowledge. We consider here that the overall conformational stability of three-state proteins is made of a "relevant" term and a "residual" one, corresponding to the free energy differences of the native to intermediate (N-to-I) and intermediate to denatured (I-to-D) equilibria, respectively. The N-to-I free energy difference is considered to be the relevant stability because protein-unfolding intermediates are likely devoid of biological activity. We use surface charge optimisation to first increase the overall (N-to-D) stability of a model three-state protein (apoflavodoxin) and then investigate whether the stabilisation obtained is realised into relevant or into residual stability. Most of the mutations designed from electrostatic calculations or from simple sequence conservation analysis produce large increases in the overall stability of the protein. However, in most cases, this simply leads to similarly large increases of the residual stability. Two mutations, nevertheless, show a different trend and increase the relevant stability of the protein substantially. When all the mutations are mapped onto the structure of the apoflavodoxin thermal-unfolding intermediate (obtained independently by equilibrium phi-analysis and NMR) they cluster perfectly so that the mutations increasing the relevant stability appear in the small unstructured region of the intermediate and the others in the native-like region. This illustrates the need for specific investigation of N-to-I equilibria and the structure of protein intermediates, and indicates that it is possible to rationally stabilise a protein against partial unfolding once the structure of the intermediate conformation is known, even if at low resolution.

The structural dynamics of myoglobin

Conformational fluctuations in proteins were initially invoked to explain the observation that diffusion of small ligands through the matrix is a global phenomenon. Small globular proteins contain internal cavities that play a role not only in matrix dynamics but also in controlling function, tracing a pathway for the diffusion of the ligand to and from the active site. This is the main point addressed in this Review, which presents pertinent information obtained on myoglobin (Mb). Mb, a simple globular heme protein which binds reversibly oxygen and other ligands. The bond between the heme Fe(II) and gaseous ligands can be photodissociated by a laser pulse, generating a non-equilibrium population of protein structures that relaxes on a picosecond to millisecond time range. This process is associated with migration of the ligand to internal cavities of the protein, which are known to bind xenon. Some of the results obtained by laser photolysis, molecular dynamics simulations, and X-ray diffraction of intermediate states of wild-type and mutant myoglobins are summarized. The extended relaxation of the globin moiety directly observed by Laue crystallography reflects re-equilibration among conformational substates known to play an essential role in controlling protein function.

10 Residue Folded Peptide Designed by Segment Statistics

We have designed a peptide termed chignolin, consisting of only 10 amino acid residues (GYDPETGTWG), on the basis of statistics derived from more than 10,000 protein segments. The peptide folds into a unique structure in water and shows a cooperative thermal transition, both of which may be hallmarks of a protein. Also, the experimentally determined beta-hairpin structure was very close to what we had targeted. The performance of the short peptide not only implies that the methodology employed here can contribute toward development of novel techniques for protein design, but it also yields insights into the raison d'etre of an autonomous element involved in a natural protein. This is of interest for the pursuit of folding mechanisms and evolutionary processes of proteins.

Kinetics of heat- and acidification-induced aggregation of firefly luciferase

The general approach to analysis of the kinetics of protein aggregation registered by the turbidimetric method has been elaborated. The terminal part of the kinetic curves is analyzed using a theoretical equation connecting the derivative of the apparent absorbance (A) with respect to time (dA/dt) and A (t is time). This analysis allows the limiting value of A at t--> infinity (A(lim)) and the order of aggregation with respect to protein (n) to be calculated. Approach proposed was applied to analysis of thermal and acidification-induced aggregation of firefly luciferase. In both cases the A(lim) value is a linear function of the protein concentration. The terminal part of the kinetic curves of thermal aggregation follows the first-order kinetics (n=1), whereas the kinetics of acidification-induced aggregation are characterized by the value of n higher than unity (n=1.29). The mechanism of nucleation-dependent aggregation has been discussed.

Fast folding of the two-domain semliki forest virus capsid protein explains co-translational proteolytic activity

The capsid protein of Semliki Forest virus constitutes the N-terminal part of a large viral polyprotein. It consists of an unstructured basic segment (residues 1-118) and a 149 residue serine protease module (SFVP, residues 119-267) comprised of two beta-barrel domains. Previous in vivo and in vitro translation experiments have demonstrated that SFVP folds co-translationally during synthesis of the viral polyprotein and rapidly cleaves itself off the nascent chain. To test whether fast co-translation folding of SFVP is an intrinsic property of the polypeptide chain or whether folding is accelerated by cellular components, we investigated spontaneous folding of recombinant SFVP in vitro. The results show that the majority of unfolded SFVP molecules fold faster than any previously studied two-domain protein (tau=50 ms), and that folding of the N-terminal domain precedes structure formation of the C-terminal domain. This shows that co-translational folding of SFVP does not require additional cellular components and suggests that rapid folding is the result of molecular evolution towards efficient virus biogenesis.

Folding and association of an extremely stable dimeric protein from Sulfolobus islandicus

ORF56 is a plasmid-encoded protein from Sulfolobus islandicus, which probably controls the copy number of the pRN1 plasmid by binding to its own promotor. The protein showed an extremely high stability in denaturant, heat, and pH-induced unfolding transitions, which can be well described by a two-state reaction between native dimers and unfolded monomers. The homodimeric character of native ORF56 was confirmed by analytical ultracentrifugation. Far-UV circular dichroism and fluorescence spectroscopy gave superimposable denaturant-induced unfolding transitions and the midpoints of both heat as well as denaturant-induced unfolding depend on the protein concentration supporting the two-state model. This model was confirmed by GdmSCN-induced unfolding monitored by heteronuclear 2D NMR spectroscopy. Chemical denaturation was accomplished by GdmCl and GdmSCN, revealing a Gibbs free energy of stabilization of -85.1 kJ/mol at 25 degrees C. Thermal unfolding was possible only above 1 M GdmCl, which shifted the melting temperature (t(m)) below the boiling point of water. Linear extrapolation of t(m) to 0 M GdmCl yielded a t(m) of 107.5 degrees C (5 microM monomer concentration). Additionally, ORF56 remains natively structured over a remarkable pH range from pH 2 to pH 12. Folding kinetics were followed by far-UV CD and fluorescence after either stopped-flow or manual mixing. All kinetic traces showed only a single phase and the two probes revealed coincident folding rates (k(f), k(u)), indicating the absence of intermediates. Apparent first-order refolding rates depend linearly on the protein concentration, whereas the unfolding rates do not. Both lnk(f) and lnk(u) depend linearly on the GdmCl concentration. Together, folding and association of homodimeric ORF56 are concurrent events. In the absence of denaturant ORF56 refolds fast (7.0 x 10(7)M(-1)s(-1)) and unfolds extremely slowly (5.7 year(-1)). Therefore, high stability is coupled to a slow unfolding rate, which is often observed for proteins of extremophilic organisms.

A proline switch controls folding and domain interactions in the gene-3-protein of the filamentous phage fd

The amino-terminal domains N1 and N2 of the gene-3-protein of phage fd form a bilobal structural and functional entity that protrudes from the phage tip. Domain N2 initiates the infection of Escherichia coli by binding to the F pilus. This binding results in the dissociation of the two domains and allows N1 to interact with the TolA receptor at the cell surface. The refolding of the N1-N2 fragment begins with the folding of domain N1, which takes a few milliseconds, followed by the folding of domain N2, which is complete within five minutes. The subsequent domain assembly is unusually slow and shows a time-constant of 6200 s at 25 degrees C. We found that the rate of this reaction is controlled by the trans to cis isomerization of the Gln212-Pro213 bond in the hinge subdomain of N2, a region that provides many interactions between N1 and N2 in the gene-3-protein. The substitution of Pro213 by Gly accelerated domain association 30-fold and revealed that the folding of the two individual domains and their assembly are indeed sequential steps in the refolding of the gene-3-protein. In the course of infection, the domains must separate to expose the binding site for TolA on domain N1. The kinetic block of domain reassembly caused by Pro213 isomerization could ensure that after the initial binding of N2 to the F pilus the open state persists until N1 and TolA are close enough for their mutual interaction. Pro213 isomerization might thus serve as a slow conformational switch in the function of the gene-3-protein.

A simple procedure to weight empirical potentials in a fitness function so as to optimize its performance in ab initio protein-folding problem

In an approach to the protein folding problem by a Genetic Algorithm, the fitness function plays a critical role. Empirical potentials are generally used to build the fitness function, and they must be weighted to obtain a valuable one. The weights are generally found by the comparison with a set of misfolded structures (decoys), but a dependence of the obtained fitness generally arises on the used decoys. Here we describe a general procedure to find out, from a given set of potentials, their better linear combination that could either identify the wild structure or prove their powerlessness. We use topological considerations over the hyperspace of the potentials, and a multiple linear inequalities solver. The iterated method flows through the following steps: it determines a direction in the hyperspace of the potentials, which identifies the native structure as a vertex among a set of misfolded decoys. A multiple linear inequalities solver obtains the direction. A Genetic Algorithm, tailored to the specific problem, uses the fitness function defined by this direction and generally reaches a new structure better than the experimental one, which is added to the ensemble. The decoys so generated are not dependent on a deterministic criterion. This iterative procedure can be stopped either by identifying an effective fitness function or by proving the impossibility of its achievement. In order to test the method under the hardest conditions, we choose numerous and heterogeneous quantities as components of the fitness function. This method could be a useful tool for the scientific community in order to test any fitness proposed and to recognize the most important components on which it is built.

The folding mechanism of a two-domain protein: folding kinetics and domain docking of the gene-3 protein of phage fd

The gene-3 protein (G3P) of filamentous phages is essential for the infection of Escherichia coli. The carboxy-terminal domain anchors this protein in the phage coat, whereas the two amino-terminal domains N1 and N2 protrude from the phage surface. We analyzed the folding mechanism of the two-domain fragment N1-N2 of G3P (G3P(*)) and the interplay between folding and domain assembly. For this analysis, a variant of G3P(*) was used that contained four stabilizing mutations (IIHY-G3P(*)). The observed refolding kinetics extend from 10 ms to several hours. Domain N1 refolds very rapidly (with a time constant of 9.4 ms at 0.5 M guanidinium chloride, 25 degrees C) both as a part of IIHY-G3P(*) and as an isolated protein fragment. The refolding of domain N2 is slower and involves two reactions with time constants of seven seconds and 42 seconds. These folding reactions of the individual domains are followed by a very slow, spectroscopically silent docking process, which shows a time constant of 6200 seconds. This reaction was detected by a kinetic unfolding assay for native molecules. Before docking, N1 and N2 unfold fast and independently, after docking they unfold slowly in a correlated fashion. A high energy barrier is thus created by domain docking, which protects G3P kinetically against unfolding. The slow domain docking is possibly important for the infection of E.coli by the phage. Upon binding to the F pilus, the N2 domain separates from N1 and the binding site for TolA on domain N1 is exposed. Since domain reassembly is so slow, this binding site remains accessible until pilus retraction has brought N1 close to TolA on the bacterial surface.

Calcium and domain interactions contribute to the thermostability of domains of the multimodular cellobiohydrolase, CbhA, a subunit of the Clostridium thermocellum cellulosome

Each of three internal domains of multi-modular cellobiohydrolase CbhA from Clostridium thermocellum, X1(1), X1(2) (previously designated as fibronectin type 3-like modules, Fn3(1) and Fn3(2)) and family 3 carbohydrate-binding module (CBM3) binds 1 mol of Ca(2+). Structures and thermal stabilities of X1(1), X1(2), CBM3, X1(1)X1(2), and X1(1)X1(2)-CBM3 containing Ca(2+) (holo-proteins) and without Ca(2+) (apo-proteins) have been studied using CD spectroscopy. All domains are beta-proteins with irregular far-UV CD spectra due to the aromatic side chain contributions. The positive signal at 294 nm in the near-UV CD spectrum of X1(1) lacking a tryptophan residue might be attributed to the presence of aromatic clusters. Thermal denaturation of all proteins is reversible and results in the total loss of tertiary structure and preservation of significant amount of ordered secondary structure. Removal of Ca(2+) destabilizes polypeptides in a different way and to a different extent. It decreases the melting temperature ( T (m)) (by 20 degrees C) and co-operativity of thermal transition of X1(1), increases the number of transitions and lowers the co-operativity of unfolding of CBM3, and slightly decreases T (m)s (2.4-4.2 degrees C) of X1(2), X1(1)X1(2), and X1(1)X1(2)-CBM3. Transitions of X1(1)X1(2) and X1(1)X1(2)-CBM3 follow a two-state model regardless of the presence of Ca(2+). X1(1) is strongly stabilized in the apo-X1(1)X1(2) and apo-X1(1)X1(2)-CBM3 as they display T (m)s similar to those of individual and combined holo-modules. Observed CD spectra of X1(1)X1(2) and X1(1)X1(2)-CBM3 differ from those calculated as the simple weighted sum of individual modules. These differences are more prominent in spectra of apo-proteins. The results indicate the presence of inter-domain interactions in CbhA. Holo-modules, i.e. containing Ca(2+), behave essentially independently, but in the absence of Ca(2+) domain interactions are more important for the conformation of the polypeptides.

Evolutionary stabilization of the gene-3-protein of phage fd reveals the principles that govern the thermodynamic stability of two-domain proteins

The gene-3-protein (G3P) of filamentous phage is essential for their propagation. It consists of three domains. The CT domain anchors G3P in the phage coat, the N2 domain binds to the F pilus of Escherichia coli and thus initiates infection, and the N1 domain continues by interacting with the TolA receptor. Phage are thus only infective when the three domains of G3P are tightly linked, and this requirement is exploited by Proside, an in vitro selection method for proteins with increased stability. In Proside, a repertoire of variants of the protein to be stabilized is inserted between the N2 and the CT domains of G3P. Stabilized variants can be selected because they resist cleavage by a protease and thus maintain the essential linkage between the domains. The method is limited by the proteolytic stability of G3P itself. We improved the stability of G3P by subjecting the phage without a guest protein to rounds of random in vivo mutagenesis and proteolytic Proside selections. Variants of G3P with one to four mutations were selected, and the temperature at which the corresponding phage became accessible for a protease increased in a stepwise manner from 40 degrees C to almost 60 degrees C. The N1-N2 fragments of wild-type gene-3-protein and of the four selected variants were purified and their stabilities towards thermal and denaturant-induced unfolding were determined. In the biphasic transitions of these proteins domain dissociation and unfolding of N2 occur in a concerted reaction in the first step, followed by the independent unfolding of domain N1 in the second step. N2 is thus less stable than N1, and it unfolds when the interactions with N1 are broken. The strongest stabilizations were caused by mutations in domain N2, in particular in its hinge subdomain, which provides many stabilizing interactions between the N1 and N2 domains. These results reveal how the individual domains and their assembly contribute to the overall stability of two-domain proteins and how mutations are optimally placed to improve the stability of such proteins.

In vitro unfolding, refolding, and polymerization of human gammaD crystallin, a protein involved in cataract formation

Human gammaD crystallin (HgammaD-Crys), a major protein of the human eye lens, is a primary component of cataracts. This 174-residue primarily beta-sheet protein is made up of four Greek keys separated into two domains. Mutations in the human gene sequence encoding HgammaD-Crys are implicated in early-onset cataracts in children, and the mutant protein expressed in Escherichia coli exhibits properties that reflect the in vivo pathology. We have characterized the unfolding, refolding, and competing aggregation of human wild-type HgammaD-Crys as a function of guanidinium hydrochloride (GuHCl) concentration at neutral pH and 37 degrees C, using intrinsic tryptophan fluorescence to monitor in vitro folding. Wild-type HgammaD-Crys exhibited reversible refolding above 1.0 M GuHCl. The GuHCl unfolded protein was more fluorescent than its native counterpart despite the absence of metal or ion-tryptophan interactions. Aggregation of refolding intermediates of HgammaD-Crys was observed in both equilibrium and kinetic refolding processes. The aggregation pathway competed with productive refolding at denaturant concentrations below 1.0 M GuHCl, beyond the major conformational transition region. Atomic force microscopy of samples under aggregating conditions revealed the sequential appearance of small nuclei, thin protofibrils, and fiber bundles. The HgammaD-Crys fibrous aggregate species bound bisANS appreciably, indicating the presence of exposed hydrophobic pockets. The mechanism of HgammaD-Crys aggregation may provide clues to understanding age-onset cataract formation in vivo.