Fast and one-step folding of closely and distantly related homologous proteins of a four-helix bundle family

An overview on the conversion of glycerol to value-added industrial products via chemical and biochemical routes

Glycerol is a common by-product of industrial biodiesel syntheses. Due to its properties, availability, and versatility, residual glycerol can be used as a raw material in the production of high value-added industrial inputs and outputs. In particular, products like hydrogen, propylene glycol, acrolein, epichlorohydrin, dioxalane and dioxane, glycerol carbonate, n-butanol, citric acid, ethanol, butanol, propionic acid, (mono-, di-, and triacylglycerols), cynamoil esters, glycerol acetate, benzoic acid, and other applications. In this context, the present study presents a critical evaluation of the innovative technologies based on the use of residual glycerol in different industries, including the pharmaceutical, textile, food, cosmetic, and energy sectors. Chemical and biochemical catalysts in the transformation of residual glycerol are explored, along with the factors to be considered regarding the choice of catalyst route used in the conversion process, aiming at improving the production of these industrial products.

Effects of Ligand Binding on the Energy Landscape of Acyl-CoA-Binding Protein

Binding of ligands is often crucial for function yet the effects of ligand binding on the mechanical stability and energy landscape of proteins are incompletely understood. Here, we use a combination of single-molecule optical tweezers and MD simulations to investigate the effect of ligand binding on the energy landscape of acyl-coenzyme A (CoA)-binding protein (ACBP). ACBP is a topologically simple and highly conserved four-α-helix bundle protein that acts as an intracellular transporter and buffer for fatty-acyl-CoA and is active in membrane assembly. We have previously described the behavior of ACBP under tension, revealing a highly extended transition state (TS) located almost halfway between the unfolded and native states. Here, we performed force-ramp and force-jump experiments, in combination with advanced statistical analysis, to show that octanoyl-CoA binding increases the activation free energy for the unfolding reaction of ACBP without affecting the position of the transition state along the reaction coordinate. It follows that ligand binding enhances the mechanical resistance and thermodynamic stability of the protein, without changing its mechanical compliance. Steered molecular dynamics simulations allowed us to rationalize the results in terms of key interactions that octanoyl-CoA establishes with the four α-helices of ACBP and showed that the unfolding pathway is marginally affected by the ligand. The results show that ligand-induced mechanical stabilization effects can be complex and may prove useful for the rational design of stabilizing ligands.

Global analysis of protein stability by temperature and chemical denaturation

The stability of a protein is a fundamental property that determines under which conditions, the protein is functional. Equilibrium unfolding with denaturants requires preparation of several samples and only provides the free energy of folding when performed at a single temperature. The typical sample requirement is around 0.5-1 mg of protein. If the stability of many proteins or protein variants needs to be determined, substantial protein production may be needed. Here we have determined the stability of acyl-coenzyme A binding protein at pH 5.3 and chymotrypsin inhibitor 2 at pH 3 and pH 6.25 by combined temperature and denaturant unfolding. We used a setup where tryptophan fluorescence is measured in quartz capillaries where only 10 μl is needed. Temperature unfolding of a series of 15 samples at increasing denaturant concentrations provided accurate and precise thermodynamic parameters. We find that the number of samples may be further reduced and less than 10 μg of protein in total are needed for reliable stability measurements. For assessment of stability of protein purified in small scale e.g. in micro plate format, our method will be highly applicable. The routine for fitting the experimental data is made available as a python notebook.

Structure and conformational stability of the triosephosphate isomerase from Zea mays. Comparison with the chemical unfolding pathways of other eukaryotic TIMs

We studied the structure, function and thermodynamic properties for the unfolding of the Triosephosphate isomerase (TIM) from Zea mays (ZmTIM). ZmTIM shows a catalytic efficiency close to the diffusion limit. Native ZmTIM is a dimer that dissociates upon dilution into inactive and unfolded monomers. Its thermal unfolding is irreversible with a T_m of 61.6 ± 1.4 °C and an activation energy of 383.4 ± 11.5 kJ mol^-1. The urea-induced unfolding of ZmTIM is reversible. Transitions followed by catalytic activity and spectroscopic properties are monophasic and superimposable, indicating that ZmTIM unfolds/refolds in a two-state behavior with an unfolding ΔG°_(H20) = 99.8 ± 5.3 kJ mol^-1. This contrasts with most other studied TIMs, where folding intermediates are common. The three-dimensional structure of ZmTIM was solved at 1.8 Å. A structural comparison with other eukaryotic TIMs shows a similar number of intramolecular and intermolecular interactions. Interestingly the number of interfacial water molecules found in ZmTIM is lower than those observed in most TIMs that show folding intermediates. Although with the available data, there is no clear correlation between structural properties and the number of equilibrium intermediates in the unfolding of TIM, the identification of such structural properties should increase our understanding of folding mechanisms.

Protein folding, misfolding and aggregation: The importance of two-electron stabilizing interactions

Proteins associated with neurodegenerative diseases are highly pleiomorphic and may adopt an all-α-helical fold in one environment, assemble into all-β-sheet or collapse into a coil in another, and rapidly polymerize in yet another one via divergent aggregation pathways that yield broad diversity of aggregates’ morphology. A thorough understanding of this behaviour may be necessary to develop a treatment for Alzheimer’s and related disorders. Unfortunately, our present comprehension of folding and misfolding is limited for want of a physicochemical theory of protein secondary and tertiary structure. Here we demonstrate that electronic configuration and hyperconjugation of the peptide amide bonds ought to be taken into account to advance such a theory. To capture the effect of polarization of peptide linkages on conformational and H-bonding propensity of the polypeptide backbone, we introduce a function of shielding tensors of the C^α atoms. Carrying no information about side chain-side chain interactions, this function nonetheless identifies basic features of the secondary and tertiary structure, establishes sequence correlates of the metamorphic and pH-driven equilibria, relates binding affinities and folding rate constants to secondary structure preferences, and manifests common patterns of backbone density distribution in amyloidogenic regions of Alzheimer’s amyloid β and tau, Parkinson’s α-synuclein and prions. Based on those findings, a split-intein like mechanism of molecular recognition is proposed to underlie dimerization of Aβ, tau, αS and PrP^C, and divergent pathways for subsequent association of dimers are outlined; a related mechanism is proposed to underlie formation of PrP^Sc fibrils. The model does account for: (i) structural features of paranuclei, off-pathway oligomers, non-fibrillar aggregates and fibrils; (ii) effects of incubation conditions, point mutations, isoform lengths, small-molecule assembly modulators and chirality of solid-liquid interface on the rate and morphology of aggregation; (iii) fibril-surface catalysis of secondary nucleation; and (iv) self-propagation of infectious strains of mammalian prions.

Evidence for the principle of minimal frustration in the evolution of protein folding landscapes

Theoretical and experimental studies have firmly established that protein folding can be described by a funneled energy landscape. This funneled energy landscape is the result of foldable protein sequences evolving following the principle of minimal frustration, which allows proteins to rapidly fold to their native biologically functional conformations. For a protein family with a given functional fold, the principle of minimal frustration suggests that, independent of sequence, all proteins within this family should fold with similar rates. However, depending on the optimal living temperature of the organism, proteins also need to modulate their thermodynamic stability. Consequently, the difference in thermodynamic stability should be primarily caused by differences in the unfolding rates. To test this hypothesis experimentally, we performed comprehensive thermodynamic and kinetic analyses of 15 different proteins from the thioredoxin family. Eight of these thioredoxins were extant proteins from psychrophilic, mesophilic, or thermophilic organisms. The other seven protein sequences were obtained using ancestral sequence reconstruction and can be dated back over 4 billion years. We found that all studied proteins fold with very similar rates but unfold with rates that differ up to three orders of magnitude. The unfolding rates correlate well with the thermodynamic stability of the proteins. Moreover, proteins that unfold slower are more resistant to proteolysis. These results provide direct experimental support to the principle of minimal frustration hypothesis.

The octadecaneuropeptide stimulates somatolactin release from cultured goldfish pituitary cells

The present study aimed to investigate the distribution of the octadecaneuropeptide (ODN) in the goldfish brain and to look for a possible effect of ODN on somatolactin (SL) release from pituitary cells. A discrete population of ODN-immunoreactive neurones was localised in the lateral part of the nucleus lateralis tuberis. These neurones sent projections through the neurohypophyseal tract towards the neurohypophysis, and nerve fibres were seen in the close vicinity of SL-producing cells in the pars intermedia. Incubation of cultured goldfish pituitary cells with graded concentrations of ODN (10(-9) -10(-5 ) m) induced a dose-dependent stimulation of SL-β, but not SL-α, release. ODN-evoked SL release was blocked by the metabotrophic endozepine receptor antagonist cyclo(1-8) [DLeu(5) ]OP but was not affected by the central-type benzodiazepine receptor antagonist flumazenil. ODN-induced SL release was suppressed by treatment with the phospholipase C (PLC) inhibitor U-73122 but not with the protein kinase A (PKA) inhibitor H-89. These results indicate that, in fish, ODN produced by hypothalamic neurones acts as a hypophysiotrophic neuropeptide stimulating SL release. The effect of ODN is mediated through a metabotrophic endozepine receptor positively coupled to the PLC/inositol 1,4,5-trisphosphate/protein kinase C-signalling pathway.

Analysis of the distribution of functionally relevant rare codons

Background

The substitution of rare codons with more frequent codons is a commonly applied method in heterologous gene expression to increase protein yields. However, in some cases these substitutions lead to a decrease of protein solubility or activity. To predict these functionally relevant rare codons, a method was developed which is based on an analysis of multisequence alignments of homologous protein families.

Results

The method successfully predicts functionally relevant codons in fatty acid binding protein and chloramphenicol acetyltransferase which had been experimentally determined. However, the analysis of 16 homologous protein families belonging to the α/β hydrolase fold showed that functionally rare codons share no common location in respect to the tertiary and secondary structure.

Conclusion

A systematic analysis of multisequence alignments of homologous protein families can be used to predict rare codons with a potential impact on protein expression. Our analysis showed that most genes contain at least one putative rare codon rich region. Rare codons located near to those regions should be excluded in an approach of improving protein expression by an exchange of rare codons by more frequent codons.

Acyl-CoA binding proteins; structural and functional conservation over 2000 MYA

Besides serving as essential substrates for beta-oxidation and synthesis of triacylglycerols and more complex lipids like sphingolipids and sterol esters, long-chain fatty acyl-CoA esters are increasingly being recognized as important regulators of enzyme activities and gene transcription. Acyl-CoA binding protein, ACBP, has been proposed to play a pivotal role in the intracellular trafficking and utilization of long-chain fatty acyl-CoA esters. Depletion of acyl-CoA binding protein in yeast results in aberrant organelle morphology incl. fragmented vacuoles, multi-layered plasma membranes and accumulation of vesicles of variable sizes. In contrast to synthesis and turn-over of glycerolipids, the levels of very-long-chain fatty acids, long-chain bases and ceramide are severely affected by Acb1p depletion, suggesting that Acb1p, rather than playing a general role, serves specific roles in cellular lipid metabolism.

Kinetics analysis methods for approximate folding landscapes

MOTIVATION

Protein motions play an essential role in many biochemical processes. Lab studies often quantify these motions in terms of their kinetics such as the speed at which a protein folds or the population of certain interesting states like the native state. Kinetic metrics give quantifiable measurements of the folding process that can be compared across a group of proteins such as a wild-type protein and its mutants.

RESULTS

We present two new techniques, map-based master equation solution and map-based Monte Carlo simulation, to study protein kinetics through folding rates and population kinetics from approximate folding landscapes, models called maps. From these two new techniques, interesting metrics that describe the folding process, such as reaction coordinates, can also be studied. In this article we focus on two metrics, formation of helices and structure formation around tryptophan residues. These two metrics are often studied in the lab through circular dichroism (CD) spectra analysis and tryptophan fluorescence experiments, respectively. The approximated landscape models we use here are the maps of protein conformations and their associated transitions that we have presented and validated previously. In contrast to other methods such as the traditional master equation and Monte Carlo simulation, our techniques are both fast and can easily be computed for full-length detailed protein models. We validate our map-based kinetics techniques by comparing folding rates to known experimental results. We also look in depth at the population kinetics, helix formation and structure near tryptophan residues for a variety of proteins.

AVAILABILITY

We invite the community to help us enrich our publicly available database of motions and kinetics analysis by submitting to our server: http://parasol.tamu.edu/foldingserver/.

Rapid folding and unfolding of Apaf-1 CARD

Caspase recruitment domains (CARDs) are members of the death domain superfamily and contain six antiparallel helices in an alpha-helical Greek key topology. We have examined the equilibrium and kinetic folding of the CARD of Apaf-1 (apoptotic protease activating factor 1), which consists of 97 amino acid residues, at pH 6 and pH 8. The results showed that an apparent two state equilibrium mechanism is not adequate to describe the folding of Apaf-1 CARD at either pH, suggesting the presence of intermediates in equilibrium unfolding. Interestingly, the results showed that the secondary structure is less stable than the tertiary structure, based on the transition mid-points for unfolding. Single mixing and sequential mixing stopped-flow studies showed that Apaf-1 CARD folds and unfolds rapidly and suggest a folding mechanism that contains parallel channels with two unfolded conformations folding to the native conformation. Kinetic simulations show that a slow folding phase is described by a third conformation in the unfolded ensemble that interconverts with one or both unfolded species. Overall, the native ensemble is formed rapidly upon refolding. This is in contrast to other CARDs in which folding appears to be dominated by formation of kinetic traps.

Structural and mechanistic exploration of acid resistance: kinetic stability facilitates evolution of extremophilic behavior

Kinetically stable proteins are unique in that their stability is determined solely by kinetic barriers rather than by thermodynamic equilibria. To better understand how kinetic stability promotes protein survival under extreme environmental conditions, we analyzed the unfolding behavior and determined the structure of Nocardiopsis alba Protease A (NAPase), an acid-resistant, kinetically stable protease, and compared these results with a neutrophilic homolog, alpha-lytic protease (alphaLP). Although NAPase and alphaLP have the same number of acid-titratable residues, kinetic studies revealed that the height of the unfolding free energy barrier for NAPase is less sensitive to acid than that of alphaLP, thereby accounting for NAPase's improved tolerance of low pH. A comparison of the alphaLP and NAPase structures identified multiple salt-bridges in the domain interface of alphaLP that were relocated to outer regions of NAPase, suggesting a novel mechanism of acid stability in which acid-sensitive electrostatic interactions are rearranged to similarly affect the energetics of both the native state and the unfolding transition state. An acid-stable variant of alphaLP in which a single interdomain salt-bridge is replaced with a corresponding intradomain NAPase salt-bridge shows a dramatic >15-fold increase in acid resistance, providing further evidence for this hypothesis. These observations also led to a general model of the unfolding transition state structure for alphaLP protease family members in which the two domains separate from each other while remaining relatively intact themselves. These results illustrate the remarkable utility of kinetic stability as an evolutionary tool for developing longevity over a broad range of harsh conditions.

Massive sequence perturbation of the Raf ras binding domain reveals relationships between sequence conservation, secondary structure propensity, hydrophobic core organization and stability

The contributions of specific residues to the delicate balance between function, stability and folding rates could be determined, in part by [corrected] comparing the sequences of structures having identical folds, but insignificant sequence homology. Recently, we have devised an experimental strategy to thoroughly explore residue substitutions consistent with a specific class of structure. Using this approach, the amino acids tolerated at virtually all residues of the c-Raf/Raf1 ras binding domain (Raf RBD), an exemplar of the common beta-grasp ubiquitin-like topology, were obtained and used to define the sequence determinants of this fold. Herein, we present analyses suggesting that more subtle sequence selection pressure, including propensity for secondary structure, the hydrophobic core organization and charge distribution are imposed on the Raf RBD sequence. Secondly, using the Gibbs free energies (DeltaG(F-U)) obtained for 51 mutants of Raf RBD, we demonstrate a strong correlation between amino acid conservation and the destabilization induced by truncating mutants. In addition, four mutants are shown to significantly stabilize Raf RBD native structure. Two of these mutations, including the well-studied R89L, are known to severely compromise binding affinity for ras. Another stabilized mutant consisted of a deletion of amino acid residues E104-K106. This deletion naturally occurs in the homologues a-Raf and b-Raf and could indicate functional divergence. Finally, the combination of mutations affecting five of 78 residues of Raf RBD results in stabilization of the structure by approximately 12 kJ mol(-1) (DeltaG(F-U) is -22 and -34 kJ mol(-1) for wt and mutant, respectively). The sequence perturbation approach combined with sequence/structure analysis of the ubiquitin-like fold provide a basis for the identification of sequence-specific requirements for function, stability and folding rate of the Raf RBD and structural analogues, highlighting the utility of conservation profiles as predictive tools of structural organization.

Catabolism of the octadecaneuropeptide ODN by prolyl endopeptidase: identification of an unusual cleavage site

The octadecaneuropeptide ODN (QATVGDVNTDRPGLLDLK), a biologically active fragment of diazepam-binding inhibitor, exerts a number of behavioral and neurophysiological activities. The presence of a proline residue in the sequence of ODN led us to investigate the role of proline endopeptidase (PEP) in the catabolism of this neuropeptide. The effect of PEP on the breakdown of ODN and related analogs was studied by combining RP-HPLC analysis and MALDI-TOF MS characterization. Incubation of ODN with PEP generated two products, i.e. ODN3-18 and ODN5-18 which resulted from cleavage of the Ala-Thr and Val-Gly peptide bonds. S 17092, a specific PEP inhibitor, significantly reduced the PEP-induced cleavages of ODN. Similarly, [Ala2]OP showed S 17092-sensitive post-alanine cleavage, while [pGlu1]ODN and OP (ODN11-18) were not catabolized by the enzyme. For all these peptides, cleavage of the Pro-Gly peptide bond by PEP was never observed, even after prolonged incubation times. In contrast, PEP hydrolyzed human urotensin II at the canonical post-proline site. Collectively, these data suggest that the Ala2 residue is the preferential cleavage site of ODN and that the Pro-Gly bond of ODN is not hydrolyzed by PEP. In addition, this study reveals for the first time that the endoproteolytic activity of PEP can specifically take place after a valine moiety.

Do proteins with similar folds have similar transition state structures? A diffuse transition state of the 169 residue apoflavodoxin

Apoflavodoxin from Anabaena PCC 7119 is a 169 residue globular protein of known structure and energetics. Here, we present a comprehensive Phi-value analysis to characterize the structure of its transition state. A total of 34 non-disruptive mutations are made throughout the structure and a range of Phi-values from zero to one are observed. In addition, a small set of eight aliphatic small-to-large mutations have been introduced in the hydrophobic core of the protein and they have been analyzed to investigate the feasibility of stabilizing the unfolding transition state by creating new non-native interactions. We find that the transition state of apoflavodoxin (so far the largest protein subjected to Phi-analysis) is diffuse and that it can be stabilized by unspecific hydrophobic interactions that can speed up the folding reaction. The data gathered on the apoflavodoxin transition state are compared with results from experimental studies in other proteins to revisit the relationship between the native state topology and transition state structure.

Biophysical characterization of the enzyme I of the Streptomyces coelicolor phosphoenolpyruvate:sugar phosphotransferase system

The first protein in the bacterial phosphoenolpyruvate (PEP):sugar phosphotransferase system is the homodimeric 60-kDa enzyme I (EI), which autophosphorylates in the presence of PEP and Mg2+. The conformational stability and structure of the EI from Streptomyces coelicolor, EI(sc), were explored in the absence and in the presence of its effectors by using several biophysical probes (namely, fluorescence, far-ultraviolet circular dichroism, Fourier transform infrared spectroscopy (FTIR), and differential scanning calorimetry) and computational approaches. The structure of EI(sc) was obtained by homology modeling of the isolated N- and C-terminal domains of other EI proteins. The experimental results indicate that at physiological pH, the dimeric EI(sc) had a well-folded structure; however, at low pH, EI(sc) showed a partially unfolded state with the features of a molten globule, as suggested by fluorescence, far-ultraviolet circular dichroism, FTIR, and 8-anilino-1-naphthalene-sulfonic acid binding. The thermal stability of EI(sc), in the absence of PEP and Mg2+, was maximal at pH 7. The presence of PEP and Mg2+ did not change substantially the secondary structure of the protein, as indicated by FTIR measurements. However, quenching experiments and proteolysis patterns suggest conformational changes in the presence of PEP; furthermore, the thermal stability of EI(sc) was modified depending on the effector added. Our approach suggests that thermodynamical analysis might reveal subtle conformational changes.

Evolution of the acyl-CoA binding protein (ACBP)

Acyl-CoA-binding protein (ACBP) is a 10 kDa protein that binds C12-C22 acyl-CoA esters with high affinity. In vitro and in vivo experiments suggest that it is involved in multiple cellular tasks including modulation of fatty acid biosynthesis, enzyme regulation, regulation of the intracellular acyl-CoA pool size, donation of acyl-CoA esters for beta-oxidation, vesicular trafficking, complex lipid synthesis and gene regulation. In the present study, we delineate the evolutionary history of ACBP to get a complete picture of its evolution and distribution among species. ACBP homologues were identified in all four eukaryotic kingdoms, Animalia, Plantae, Fungi and Protista, and eleven eubacterial species. ACBP homologues were not detected in any other known bacterial species, or in archaea. Nearly all of the ACBP-containing bacteria are pathogenic to plants or animals, suggesting that an ACBP gene could have been acquired from a eukaryotic host by horizontal gene transfer. Many bacterial, fungal and higher eukaryotic species only harbour a single ACBP homologue. However, a number of species, ranging from protozoa to vertebrates, have evolved two to six lineage-specific paralogues through gene duplication and/or retrotransposition events. The ACBP protein is highly conserved across phylums, and the majority of ACBP genes are subjected to strong purifying selection. Experimental evidence indicates that the function of ACBP has been conserved from yeast to humans and that the multiple lineage-specific paralogues have evolved altered functions. The appearance of ACBP very early on in evolution points towards a fundamental role of ACBP in acyl-CoA metabolism, including ceramide synthesis and in signalling.

Ultra-fast Barrier-limited Folding in the Peripheral Subunit-binding Domain Family

We have determined the solution structures, equilibrium properties and ultra-fast folding kinetics for three bacterial homologues of the peripheral subunit-binding domain (PSBD) family. The mesophilic homologue, BBL, was less stable than the thermophilic and hyper-thermophilic variants (E3BD and POB, respectively). The broad unfolding transitions of each PSBD, when probed by different techniques, were essentially superimposable, consistent with co-operative denaturation. Temperature-jump and continuous-flow fluorescence methods were used to measure the folding kinetics for E3BD, POB and BBL. E3BD folded fairly rapidly at 298K (folding half-time approximately 25 micros) and BBL and POB folded even faster (folding half-times approximately 3-5 micros). The variations in equilibrium and kinetic behaviour observed for the PSBD family resembles that of the homeodomain family, where the folding pattern changes from apparent two-state transitions to multi-state kinetics as the denatured state becomes more structured. The faster folding of POB may be a consequence of its higher propensity to form helical structure in the region corresponding to the folding nucleus of E3BD. The ultra-fast folding of BBL appears to be a consequence of residual structure in the denatured ensemble, as with engrailed homeodomain. We discuss issues concerning "one-state", downhill folding, and find no evidence for, and strong evidence against, it occurring in these PSBDs. The shorter construct used previously for BBL was destabilized significantly and the stability further perturbed by the introduction of fluorescent probes. Thermal titrations for 11 side-chains scattered around the protein, when probed by (13)C-NMR experiments, could be fit globally to a common co-operative transition.

NMR Spectroscopic Characterization of Millisecond Protein Folding by Transverse Relaxation Dispersion Measurements

The cold shock protein CspB adopts its native and functional tertiary structure on the millisecond time scale. We employed transverse relaxation NMR methods, which allow a quantitative measurement of the cooperativity of this fast folding reaction on a residue basis. Thereby, chemical exchange contributions to the transverse relaxation rate (R(2)) were observed for every residue of CspB verifying the potential of this method to identify not only local dynamics but also to characterize global events. Toward this end, the homogeneity of the transition state of folding was probed by comparing Chevron plots (i.e., dependence of the apparent folding rate on the denaturant concentration) determined by stopped-flow fluorescence with Chevron plots of six residues acquired by R(2) dispersion experiments. The coinciding results obtained for probes at different locations in the three-dimensional structure of CspB indicate the ability and significance of transverse relaxation NMR to determine Chevron plots on a residue-by-residue basis providing detailed insights on the nature of the transition state of folding.

Reversible Equilibrium Unfolding of Triosephosphate Isomerase from Trypanosoma cruzi in Guanidinium Hydrochloride Involves Stable Dimeric and Monomeric Intermediates

The reversible guanidinium hydrochloride-induced unfolding of Trypanosoma cruzi triosephosphate isomerase (TcTIM) was characterized under equilibrium conditions. The catalytic activity was followed as a native homodimeric functional probe. Circular dichroism, intrinsic fluorescence, and size-exclusion chromatography were used as secondary, tertiary, and quaternary structural probes, respectively. The change in ANS fluorescence intensity with increasing denaturant concentrations was also determined. The results show that two stable intermediates exist in the transition from the homodimeric native enzyme to the unfolded monomers: one (N(2*)) is a slightly more expanded, non-native, and active dimer, and the other is a partially expanded monomer (M) that binds ANS. Spectroscopic and activity data were used to reach a thermodynamic characterization. The results indicate that the Gibbs free energies for the partial reactions are 4.5 (N(2) <==> N(2*)), 65.8 (N(2*) <==> 2M), and 17.8 kJ/mol (M <==> U). It appears that TcTIM monomers are more stable than those found for other TIM species (except yeast TIM), where monomer stability is only marginal. These results are compared with those for the guanidinium hydrochloride-induced denaturation of TIM from different species, where despite the functional and three-dimensional similarities, a remarkable heterogeneity exists in the unfolding pathways.

Formation of Native and Non-native Interactions in Ensembles of Denatured ACBP Molecules from Paramagnetic Relaxation Enhancement Studies

Paramagnetic relaxation enhancement measurements in the denatured state of ACBP have provided distance restraints that have been used in computer simulations to determine the conformational ensembles representing the denatured states of ACBP under a variety of conditions. A detailed comparison of the residual structure in the denatured state of ACBP under these different conditions has enabled us to infer that regions in the N and C-terminal parts of the protein sequence have a high tendency to interact in the unfolded state under physiological conditions. By comparing the structural features in the denatured states with those in the transition state for folding we also provided new insights into the mechanism of formation of the native state of this protein.

Different secondary structure elements as scaffolds for protein folding transition states of two homologous four-helix bundles

Comparison of the folding processes for homologue proteins can provide valuable information about details in the interactions leading to the formation of the folding transition state. Here the folding kinetics of 18 variants of yACBP and 3 variants of bACBP have been studied by Phi-value analysis. In combination with Phi-values from previous work, detailed insight into the transition states for folding of both yACBP and bACBP has been obtained. Of the 16 sequence positions that have been studied in both yACBP and bACBP, 5 (V12, I/L27, Y73, V77, and L80) have high Phi-values and appear to be important for the transition state formation in both homologues. Y31, A34, and A69 have high Phi-values only in yACBP, while F5, A9, and I74 have high Phi-values only in bACBP. Thus, additional interactions between helices A2 and A4 appear to be important for the transition state of yACBP, whereas additional interactions between helices A1 and A4 appear to be important for the transition state of bACBP. To examine whether these differences could be assigned to different packing of the residues in the native state, a solution structure of yACBP was determined by NMR. Small changes in the packing of the hydrophobic side-chains, which strengthen the interactions between helices A2 and A4, are observed in yACBP relative to bACBP. It is suggested that different structure elements serve as scaffolds for the folding of the 2 ACBP homologues.

A Statistical Model for Predicting Protein Folding Rates from Amino Acid Sequence with Structural Class Information

Prediction of protein folding rates from amino acid sequences is one of the most important challenges in molecular biology. In this work, I have related the protein folding rates with physical-chemical, energetic and conformational properties of amino acid residues. I found that the classification of proteins into different structural classes shows an excellent correlation between amino acid properties and folding rates of two- and three-state proteins, indicating the importance of native state topology in determining the protein folding rates. I have formulated a simple linear regression model for predicting the protein folding rates from amino acid sequences along with structural class information and obtained an excellent agreement between predicted and experimentally observed folding rates of proteins; the correlation coefficients are 0.99, 0.96 and 0.95, respectively, for all-alpha, all-beta and mixed class proteins. This is the first available method, which is capable of predicting the protein folding rates just from the amino acid sequence with the aid of generic amino acid properties and structural class information.

The folding of spectrin domains II: phi-value analysis of R16

Studies on the folding of helical proteins have shown a wide range of different mechanisms and highlighted the importance of helical propensity as a factor in determining folding mechanism. Here, we contribute to this interesting field with the protein engineering phi-value analysis of the 16th domain of chicken brain alpha-spectrin, R16. The fortuitous curvature seen in the unfolding arm of the chevron plot allows us to investigate both early and late events in folding. R16 is the first two-state helical protein for which this has been possible.

Topological determinants of protein unfolding rates

For proteins that fold by two-state kinetics, the folding and unfolding processes are believed to be closely related to their native structures. In particular, folding and unfolding rates are influenced by the native structures of proteins. Thus, we focus on finding important topological quantities from a protein structure that determine its unfolding rate. After constructing graphs from protein native structures, we investigate the relationships between unfolding rates and various topological quantities of the graphs. First, we find that the correlation between the unfolding rate and the contact order is not as prominent as in the case of the folding rate and the contact order. Next, we investigate the correlation between the unfolding rate and the clustering coefficient of the graph of a protein native structure, and observe no correlation between them. Finally, we find that a newly introduced quantity, the impact of edge removal per residue, has a good overall correlation with protein unfolding rates. The impact of edge removal is defined as the ratio of the change of the average path length to the edge removal probability. From these facts, we conclude that the protein unfolding process is closely related to the protein native structure.

Analysis of the differences in the folding kinetics of structurally homologous proteins based on predictions of the gross features of residue contacts

It is a general notion that proteins with very similar three-dimensional structures would show very similar folding kinetics. However, recent studies reveal that the folding kinetic properties of some proteins contradict this thought (i.e., the members in a same protein family fold through different pathways). For example, it has been reported that some beta-proteins in the intracellular lipid-binding protein family fold through quite different pathways (Burns et al., Proteins 1998;33:107-118). Similar differences in folding kinetics are also observed in the members of the globin family (Nishimura et al., Nat Struct Biol 2000;7:679-686). In our study, we examine the possibility of predicting qualitative differences in folding kinetics of the intracellular lipid-binding proteins and two globin proteins (i.e., myoglobin and leghemoglobin). The problem is tackled by means of a contact map based on the average distance statistics between residues, the Average Distance Map (ADM), as constructed from sequence. The ADMs for the three proteins show overall similarity, but some local differences among maps are also observed. Our results demonstrate that some properties of the protein folding kinetics are consistent with local differences in the ADMs. We also discuss the general possibility of predicting folding kinetics from sequence information.

Equilibrium and kinetic folding of an alpha-helical Greek key protein domain: caspase recruitment domain (CARD) of RICK

We have characterized the equilibrium and kinetic folding of a unique protein domain, caspase recruitment domain (CARD), of the RIP-like interacting CLARP kinase (RICK) (RICK-CARD), which adopts a alpha-helical Greek key fold. At equilibrium, the folding of RICK-CARD is well described by a two-state mechanism representing the native and unfolded ensembles. The protein is marginally stable, with a DeltaG(H)()2(O) of 3.0 +/- 0.15 kcal/mol and an m-value of 1.27 +/- 0.06 kcal mol(-1) M(-1) (30 mM Tris-HCl, pH 8, 1 mM DTT, 25 degrees C). While the m-value is constant, the protein stability decreases in the presence of moderate salt concentrations (below 200 mM) and then increases at higher salt concentrations. The results suggest that electrostatic interactions are stabilizing in the native protein, and the favorable Coulombic interactions are reduced at low ionic strength. Above 200 mM salt, the results are consistent with Hofmeister effects. The unfolding pathway of RICK-CARD is complex and contains at least three non-native conformations. The refolding pathway of RICK-CARD also is complex, and the data suggest that the unfolded protein folds via two intermediate conformations prior to reaching the native state. Overall, the data suggest the presence of kinetically trapped, or misfolded, species that are on-pathway both in refolding and in unfolding.

Prediction of folding rates and transition-state placement from native-state geometry

A variety of experimental and theoretical studies have established that the folding process of monomeric proteins is strongly influenced by the topology of the native state. In particular, folding times have been shown to correlate well with the contact order, a measure of contact locality. Our investigation focuses on identifying additional topologic properties that correlate with experimentally measurable quantities, such as folding rates and transition-state placement, for both two- and three-state folders. The validation against data from 40 experiments shows that a particular topological property that measures the interdependence of contacts, termed cliquishness or clustering coefficient, can account with statistically significant accuracy both for the transition state placement and especially for folding rates. The observed correlations can be further improved by optimally combining the distinct topological information captured by cliquishness and contact order.

Optical characterization of armadillo acyl-CoA binding protein

Acyl-CoA binding protein (ACBP) and fatty acid binding protein (FABP) are intracellular transporters of activated and free fatty acids, respectively. Unlike other tissues with active lipid metabolism, armadillo Harderian gland contains much more ACBP than FABP. To characterize armadillo ACBP structure and binding properties, we produced it in Escherichia coli and carried out detailed fluorescence and circular dichroism spectroscopy studies. The K(D) for palmitoyl-CoA, measured directly by fluorescence and rotatory power, was 34+/-12 and 75+/-39 nM, respectively. The structure of armadillo ACBP appears to be very similar to that of bovine and rat liver ACBPs.

Structural features of cytochrome c' folding intermediates revealed by fluorescence energy-transfer kinetics

We employed fluorescence energy-transfer probes to investigate the polypeptide dynamics accompanying cytochrome c' folding. Analysis of fluorescence energy-transfer kinetics from wild-type Trp-72 or Trp-32 in a crystallographically characterized (1.78 A) Q1A/F32W/W72F mutant shows that there is structural heterogeneity in denatured cytochrome c'. Even at guanidine hydrochloride concentrations well beyond the unfolding transition, a substantial fraction of the polypeptides ( approximately 50%) adopts compact conformations (tryptophan-to-heme distance, approximately 25 A) in both pseudo-wild-type (Q1A) and mutant proteins. A burst phase (< or =5 ms) is revealed when stopped flow-triggered refolding is probed by tryptophan intensity: measurements on the Q1A protein show that approximately 75% of the Trp-72 fluorescence (83% for Trp-32) is quenched within the mixing deadtime, suggesting that most of the polypeptides have collapsed.

Fast folding of a four-helical bundle protein

The FK506-FKBP12 binding-domain of the kinase FRAP (FRB) forms a classic up-down four-helical bundle. The folding pathway of this protein has been investigated using a combination of equilibrium and kinetic studies. The native state of the protein is stable with respect to the unfolded state by some 7 kcal mol(-1) at pH 6.0, 10 degrees C. A kinetic analysis of unfolding and refolding rate constants as a function of chemical denaturant concentration suggests that an intermediate state may be populated during folding at low concentrations of denaturant. The presence of this intermediate state is confirmed by refolding experiments performed in the presence of the hydrophobic dye 8-anilinonaphthalene-1 sulfonate (ANS). ANS binds to the partially folded intermediate state populated during the folding of FRB and undergoes a large change in fluorescence that can be detected using stopped-flow techniques. Analysis of the kinetic data suggests that the intermediate state is compact and it may even be a misfolded species that has to partially unfold before it can reach the transition state. Folding and unfolding rate constants in water are approximately 150-200 s(-1) and 0.005-0.06 s(-1), respectively, at neutral pH and 10 degrees C. The folding of FRB is somewhat slower than for other all-helical proteins, probably as a consequence of the formation of a metastable intermediate state. The folding rate constant in the absence of any populated intermediate can be estimated to be 8800 s(-1). Despite the presence of an intermediate state, which effectively slows folding, the protein still folds rapidly with a half-life of 5 ms at 10 degrees C. The dependence of the rate constants on denaturant concentration indicates that the transition state for folding is compact with some 80% of the surface area exposed in the unfolded state buried in the transition state. Data presented for FRB is compared with kinetic data obtained for other all-helical proteins.

Folding mechanism of indole-3-glycerol phosphate synthase from Sulfolobus solfataricus: a test of the conservation of folding mechanisms hypothesis in (beta(alpha))(8) barrels

As a test of the hypothesis that folding mechanisms are better conserved than sequences in TIM barrels, the equilibrium and kinetic folding mechanisms of indole-3-glycerol phosphate synthase (sIGPS) from the thermoacidophilic archaebacterium Sulfolobus solfataricus were compared to the well-characterized models of the alpha subunit of tryptophan synthase (alphaTS) from Escherichia coli. A multifaceted approach combining urea denaturation and far-UV circular dichroism, tyrosine fluorescence total intensity, and tyrosine fluorescence anisotropy was employed. Despite a sequence identity of only 13%, a stable intermediate (I) in sIGPS was found to be similar to a stable intermediate in alphaTS in terms of its thermodynamic properties and secondary structure. Kinetic experiments revealed that the fastest detectable folding event for sIGPS involves a burst-phase (<5ms) reaction that leads directly to the stable intermediate. The slower of two subsequent phases reflects the formation/disruption of an off-pathway dimeric form of I. The faster phase reflects the conversion of I to the native state and is limited by folding under marginally stable conditions and by isomerization or rearrangement under strongly folding conditions. By contrast, alphaTS is thought to fold via an off-pathway burst-phase intermediate whose unfolding controls access to a set of four on-pathway intermediates that comprise the stable equilibrium intermediate. At least three proline isomerization reactions are known to limit their interconversions and lead to a parallel channel mechanism. The simple sequential mechanism deduced for sIGPS reflects the dominance of the on-pathway burst-phase intermediate and the absence of prolyl residues that partition the stable intermediate into kinetically distinguishable species. Comparison of the results for sIGPS and alphaTS demonstrates that the thermodynamic properties and the final steps of the folding reaction are better conserved than the early events. The initial events in folding appear to be more sensitive to the sequence differences between the two TIM barrel proteins.

Early kinetic intermediate in the folding of acyl-CoA binding protein detected by fluorescence labeling and ultrarapid mixing

Early conformational events during folding of acyl-CoA binding protein (ACBP), an 86-residue alpha-helical protein, were explored by using a continuous-flow mixing apparatus with a dead time of 70 micros to measure changes in intrinsic tryptophan fluorescence and tryptophan-dansyl fluorescence energy transfer. Although the folding of ACBP was initially described as a concerted two-state process, the tryptophan fluorescence measurements revealed a previously unresolved phase with a time constant tau = 80 micros, indicating formation of an intermediate with only slightly enhanced fluorescence of Trp-55 and Trp-58 relative to the unfolded state. To amplify this phase, a dansyl fluorophore was introduced at the C terminus by labeling an I86C mutant of ACBP with 5-IAEDANS [5-((((2-iodoacetyl)amino)ethyl)amino)naphthalene-1-sulfonic acid]. Continuous-flow refolding of guanidine HCl-denatured ACBP showed a major increase in tryptophan-dansyl fluorescence energy transfer, indicating formation of a partially collapsed ensemble of states on the 100-micros time scale. A subsequent decrease in dansyl fluorescence is attributed to intramolecular quenching of donor fluorescence on formation of the native state. The kinetic data are fully accounted for by three-state mechanisms with either on- or off-pathway intermediates. The intermediate accumulates to a maximum population of 40%, and its stability depends only weakly on denaturant concentration, which is consistent with a marginally stable ensemble of partially collapsed states with approximately 1/3 of the solvent-accessible surface buried. The findings indicate that ultrafast mixing methods combined with sensitive conformational probes can reveal transient accumulation of intermediate states in proteins with apparent two-state folding mechanisms.

Transient intermediary states with high and low folding probabilities in the apparent two-state folding equilibrium of ACBP at low pH

Measurements of the stability as a function of pH for the acyl-coenzyme A binding protein (ACBP) has shown a significant difference in the pH transition midpoint measured by NMR spectroscopy at pH 3.12 and the transition midpoint measured at pH 2.92 and 2.97 by circular dichroism and by fluorescence spectroscopy, respectively. A similar behavior has not been observed in other proteins. It is suggested that these differences arise because the population of the unfolded molecules still contains significant amounts of native like secondary and tertiary structure. NMR spectroscopy measures the concentration of the two components of the folding unfolding equilibrium individually, whereas circular dichroism and fluorescence measure the concentration of the conformations of the light-absorbing chromophores present in both the folded and the unfolded molecules. In the narrow pH range, nascent structure can be detected as the average amount of secondary structure per unfolded molecule and hydrophobic interactions in the population of unfolded molecules. These structures are not observable immediately by NMR spectroscopy; however, a chemical shift analysis of the peptide backbone (13)C chemical shift indicates strongly the existence of short-lived and transient helical structures at pH 2.3. Magnetization transfer studies have been applied to study the equilibrium between folded and unfolded ACBP near the pH transition point measured by NMR. This study has shown that there are two categories of subpopulations in the population of unfolded ACBP. One for which magnetization can be transferred to the folded form during the folding process, and one for which transfer is not observed. The molecules of the latter population of unfolded protein apparently, do not fold within the time-frame of the magnetization transfer experiment. This result suggests the existence of a subpopulation of the acid-unfolded protein molecules with a high propensity for folding. It is suggested that in this subpopulation, a particular set of native like interactions in the peptide backbone and between side-chains in the peptide chain have to be formed.

Folding rate prediction using total contact distance

Linear regression analysis found that either contact order (CO) or long-range order (LRO) parameter has a significant correlation with the logarithms of folding rates. This suggests that sequence separation per contact and total number of contacts are both important in determining the rate of folding. Here, the two factors are incorporated into a new parameter, total contact distance (TCD). Using a database of 28 two-state or weakly three-state folding proteins, TCD is found to be the most accurate among the three parameters (CO, LRO, and TCD) in terms of correlation and prediction. It provides even more accurate prediction than the best neural network results with two descriptors (contact order and stability per residue). The improvement is achieved in all three-structural classes (all alpha, beta, and mixed). The accuracy of total contact distance in predicting folding rates is essentially unchanged if "short"-ranged contacts (absolute value of i - j < or = 14) are not included in calculation. Thus, only long-range contacts with a sequence separation of more than 14 residues are important in determining the rate of folding. This is consistent with the results from the long-range order parameter. One of the significant outliers in prediction is found to be associated with the only protein in the database that involves nonlocal disulfide bonds. Removing the protein leads to a correlation coefficient of 0.89 between experimental observed and predicted folding rates in jackknife cross validation. The corresponding values for CO and LRO are 0.71 and 0.80, respectively.

Diffusion-collision model study of misfolding in a four-helix bundle protein

Proteins with complex folding kinetics will be susceptible to misfolding at some stage in the folding process. We simulate this problem by using the diffusion-collision model to study non-native kinetic intermediate misfolding in a four-helix bundle protein. We find a limit on the size of the pairwise hydrophobic area loss in non-native intermediates, such that burying above this limit creates long-lasting non-native kinetic intermediates that would disrupt folding and prevent formation of the native state. Our study of misfolding suggests a method for limiting the production of misfolded kinetic intermediates for helical proteins and could, perhaps, lead to more efficient production of proteins in bulk.

Acidic conditions stabilise intermediates populated during the folding of Im7 and Im9

The helical bacterial immunity proteins Im7 and Im9 have been shown to fold via kinetic mechanisms of differing complexity, despite having 60 % sequence identity. At pH 7.0 and 10 degrees C, Im7 folds in a three-state mechanism involving an on-pathway intermediate, while Im9 folds in an apparent two-state transition. In order to examine the folding mechanisms of these proteins in more detail, the folding kinetics of both Im7 and Im9 (at 10 degrees C in 0.4 M sodium sulphate) have been examined as a function of pH. Kinetic modelling of the folding and unfolding data for Im7 between pH 5.0 and 8.0 shows that the on-pathway intermediate is stabilised by more acidic conditions, whilst the native state is destabilised. The opposing effect of pH on the stability of these states results in a significant population of the intermediate at equilibrium at pH 6.0 and below. At pH 7.0, the folding and unfolding kinetics for Im9 can be fitted adequately by a two-state model, in accord with previous results. However, under acidic conditions there is a clear change of slope in the plot of the logarithm of the folding rate constant versus denaturant concentration, consistent with the population of one or more intermediate(s) early during folding. The kinetic data for Im9 at these pH values can be fitted to a three-state model, where the intermediate ensemble is stabilised and the native state destabilised as the pH is reduced, rationalising previous results that showed that an intermediate is not observed experimentally at pH 7.0. The data suggest that intermediate formation is a general step in immunity protein folding and demonstrate that it is necessary to explore a wide range of refolding conditions in order to show that intermediates do not form in the folding of other small, single-domain proteins.

On the relationship between protein stability and folding kinetics: a comparative study of the N-terminal domains of RNase HI, E. coli and Bacillus stearothermophilus L9

There is currently a great deal of interest in proteins that fold in a single highly cooperative step. Particular attention has been focused on elucidating the factors that govern folding rates of simple proteins. Recently, the topology of the native state has been proposed to be the most important determinant of their folding rates. Here we report a comparative study of the folding of three topologically equivalent proteins that adapt a particularly simple alpha/beta fold. The folding kinetics of the N-terminal domain of RNase HI and the N-terminal domain of the ribosomal protein L9 from Escherichia coli (eNTL9) were compared to the previously characterized N-terminal domain of L9 from Bacillus stearothermophilus (bNTL9). This 6.2 kDa protein, which is one of simplest examples of the ABCalphaD motif, folds via a two-state mechanism on the millisecond to submillisecond time scale. The RNase HI domain and bNTL9 have very similar tertiary structures but there is little similarity in primary sequence. bNTL9 and eNTL9 share the same biological function and a similar primary sequence but differ significantly in stability. Fluorescence-detected stopped-flow experiments showed that the three proteins fold in a two-state fashion. The folding rates in the absence of denaturant were found to be very different, ranging form 21 s(-1) to 790 s(-1) at 10 degrees C. The diverse folding rates appear to reflect large differences in the stability of the proteins. When compared at an isostability point, the folding rates converged to a similar value and there is a strong linear correlation between the log of the folding rate and stability for this set of proteins. These observations are consistent with the idea that stability can play an important role in dictating relative folding rates among topologically equivalent proteins.

Theoretical study of a landscape of protein folding-unfolding pathways. Folding rates at midtransition

This paper presents a new method for calculating the folding-unfolding rates of globular proteins. The method is based on solution of kinetic equations for a network of folding-unfolding pathways of the proteins. The rates are calculated in the point of thermodynamic equilibrium between the native and completely unfolded states. The method has been applied to all the proteins listed by Jackson [Jackson, S. E. (1998) Folding Des. 3, R81-R91] and some peptides. Although the studied protein chains differ by more than 1 order of magnitude in size and exhibit two- as well as three-state kinetics in water, and their folding rates cover more than 11 orders of magnitude, the theoretical estimates are reasonable close to the experimentally measured folding rates in midtransition (the correlation coefficient being as high as 0.78). This means that the presented theory (having no adjustable parameters at all) is consistent with the experimental observations.

Cytochrome c' folding triggered by electron transfer: fast and slow formation of four-helix bundles

Reduced (Fe(II)) Rhodopseudomonas palustris cytochrome c' (Cyt c') is more stable toward unfolding ([GuHCl](1/2) = 2.9(1) M) than the oxidized (Fe(III)) protein ([GuHCl](1/2) = 1.9(1) M). The difference in folding free energies (Delta Delta G(f) degrees = 70 meV) is less than half of the difference in reduction potentials of the folded protein (100 mV vs. NHE) and a free heme in aqueous solution ( approximately -150 mV). The spectroscopic features of unfolded Fe(II)-Cyt c' indicate a low-spin heme that is axially coordinated to methionine sulfur (Met-15 or Met-25). Time-resolved absorption measurements after CO photodissociation from unfolded Fe(II)(CO)-Cyt c' confirm that methionine can bind to the ferroheme on the microsecond time scale [k(obs) = 5(2) x 10(4) s(-1)]. Protein folding was initiated by photoreduction (two-photon laser excitation of NADH) of unfolded Fe(III)-Cyt c' ([GuHCl] = 2.02--2.54 M). Folding kinetics monitored by heme absorption span a wide time range and are highly heterogeneous; there are fast-folding ( approximately 10(3) s(-1)), intermediate-folding (10(2)-10(1) s(-1)), and slow-folding (10(-1) s(-1)) populations, with the last two likely containing methionine-ligated (Met-15 or Met-25) ferrohemes. Kinetics after photoreduction of unfolded Fe(III)-Cyt c' in the presence of CO are attributable to CO binding [1.4(6) x 10(3) s(-1)] and Fe(II)(CO)-Cyt c' folding [2.8(9) s(-1)] processes; stopped-flow triggered folding of Fe(III)-Cyt c' (which does not contain a protein-derived sixth ligand) is adequately described by a single kinetics phase with an estimated folding time constant of approximately 4 ms [Delta G(f) degrees = -33(3) kJ mol(-1)] at zero denaturant.

Anabaena apoflavodoxin hydrogen exchange: on the stable exchange core of the alpha/beta(21345) flavodoxin-like family

An important issue in modern protein biophysics is whether structurally homologous proteins share common stability and/or folding features. Flavodoxin is an archetypal alpha/beta protein organized in three layers: a central beta-sheet (strand order 21345) flanked by helices 1 and 5 on one side and helices 2, 3, and 4 on the opposite side. The backbone internal dynamics of the apoflavodoxin from Anabaena is analyzed here by the hydrogen exchange method. The hydrogen exchange rates indicate that 46 amide protons, distributed throughout the structure of apoflavodoxin, exchange relatively slowly at pH 7.0 (k(ex) < 10(-1) min(-1)). According to their distribution in the structure, protein stability is highest on the beta-sheet, helix 4, and on the layer formed by helices 1 and 5. The exchange kinetics of Anabaena apoflavodoxin was compared with those of the apoflavodoxin from Azotobacter, with which it shares a 48% sequence identity, and with Che Y and cutinase, two other alpha/beta (21345) proteins with no significant sequence homology with flavodoxins. Both similarities and differences are observed in the cores of these proteins. It is of interest that a cluster of a few structurally equivalent residues in the central beta-strands and in helix 5 is common to the cores.

Folding of intracellular retinol and retinoic acid binding proteins

The folding mechanisms of cellular retinol binding protein II (CRBP II), cellular retinoic acid binding protein I (CRABP I), and cellular retinoic acid binding protein II (CRABP II) were examined. These beta-sheet proteins have very similar structures and higher sequence homologies than most proteins in this diverse family. They have similar stabilities and show completely reversible folding at equilibrium with urea as a denaturant. The unfolding kinetics of these proteins were monitored during folding and unfolding by circular dichroism (CD) and fluorescence. During unfolding, CRABP II showed no intermediates, CRABP I had an intermediate with nativelike secondary structure, and CRBP II had an intermediate that lacked secondary structure. The refolding kinetics of these proteins were more similar. Each protein showed a burst-phase change in intensity by both CD and fluorescence, followed by a single observed phase by both CD and fluorescence and one or two additional refolding phases by fluorescence. The fluorescence spectral properties of the intermediate states were similar and suggested a gradual increase in the amount of native tertiary structure present for each step in a sequential path. However, the rates of folding differed by as much as 3 orders of magnitude and were slower than those expected from the contact order and topology of these proteins. As such, proteins with the same final structure may not follow the same route to the native state.

The fast folding pathway in human lysozyme and its blockage by appropriate mutagenesis: a sequential stopped-flow fluorescence study 1 1Edited by C. R. Matthews

In this work we were able to show that human lysozyme refolds along two parallel pathways: a fast path followed by 13% of the molecules that leads directly from a collapsed state to the native protein and a slow one for the remaining molecules that involves a partially unfolded intermediate state. However, in the refolding process of LYLA1, a chimera of human lysozyme which possesses the Ca2+-binding loop and helix C of bovine alpha-lactalbumin, the direct pathway is no longer accessible. This indicates that these structural elements, which are located in the interface region between the alpha- and beta-domain of the protein, and their interaction with the environment play an important role in the fast folding of the molecules. These results also shed some light on the conservation of folding patterns amongst structurally homologous proteins. In recent years it was often stated that structurally homologous proteins with high sequence identity follow the same folding pattern. Human lysozyme and LYLA1 have a sequence identity of 87%. However, we have shown that their folding patterns are different. Therefore, a high degree of sequence identity for two proteins belonging to the same family is not a guarantee for an identical folding pattern.