Conformational transition states of a beta-hairpin peptide between the ordered and disordered conformations in explicit water

The conformational transition states of a beta-hairpin peptide in explicit water were identified from the free energy landscapes obtained from the multicanonical ensemble, using an enhanced conformational sampling calculation. The beta-hairpin conformations were significant at 300 K in the landscape, and the typical nuclear Overhauser effect signals were reproduced, consistent with the previously reported experiment. In contrast, the disordered conformations were predominant at higher temperatures. Among the stable conformations at 300 K, there were several free energy barriers, which were not visible in the landscapes formed with the conventional parameters. We identified the transition states around the saddle points along the putative folding and unfolding paths between the beta-hairpin and the disordered conformations in the landscape. The characteristic features of these transition states are the predominant hydrophobic contacts and the several hydrogen bonds among the side-chains, as well as some of the backbone hydrogen bonds. The unfolding simulations at high temperatures, 400 K and 500 K, and their principal component analyses also provided estimates for the transition state conformations, which agreed well with those at 400 K and 500 K deduced from the current free energy landscapes at 400 K and 500 K, respectively. However, the transition states at high temperatures were much more widely distributed on the landscape than those at 300 K, and their conformations were different.     

Design and NMR conformational study of a beta-sheet peptide based on Betanova and WW domains

A good approach to test our current knowledge on formation of protein beta-sheets is de novo protein design. To obtain a three-stranded beta-sheet mini-protein, we have built a series of chimeric peptides by taking as a template a previously designed beta-sheet peptide, Betanova-LLM, and incorporating N- and/or C-terminal extensions taken from WW domains, the smallest natural beta-sheet domain that is stable in absence of disulfide bridges. Some Betanova-LLM strand residues were also substituted by those of a prototype WW domain. The designed peptides were cloned and expressed in Escherichia coli. The ability of the purified peptides to adopt beta-sheet structures was examined by circular dichroism (CD). Then, the peptide showing the highest beta-sheet population according to the CD spectra, named 3SBWW-2, was further investigated by 1H and 13C NMR. Based on NOE and chemical shift data, peptide 3SBWW-2 adopts a well defined three-stranded antiparallel beta-sheet structure with a disordered C-terminal tail. To discern between the contributions to beta-sheet stability of strand residues and the C-terminal extension, the structural behavior of a control peptide with the same strand residues as 3SBWW-2 but lacking the C-terminal extension, named Betanova-LYYL, was also investigated. beta-Sheet stability in these two peptides, in the parent Betanova-LLM and in WW-P, a prototype WW domain, decreased in the order WW-P > 3SBWW-2 > Betanova-LYYL > Betanova-LLM. Conclusions about the contributions to beta-sheet stability were drawn by comparing structural properties of these four peptides.     

Cation-pi interactions studied in a model coiled-coil peptide

Cation-pi interactions between aromatic amino acids and the positively charged residues lysine and arginine have been proposed to play an important role in stabilizing protein structure. We have used a peptide that adopts a coiled coil structure as a model system to evaluate the energetic contribution of cation-pi interactions to protein folding. Peptides were designed in which phenylalanine, tyrosine, and tryptophan were placed at a solvent-exposed position of the helix, one turn removed from an arginine residue that could provide a favorable cation-pi interaction. Only the arginine-phenylalanine pairing provided significant stabilization of the peptide structure and it appears that hydrophobic packing, rather than the cation-pi effect, is more likely to be responsible for the stability of this peptide. We conclude that any stabilizing effect of cation-pi interactions in these peptides is much smaller than that predicted from computational studies.     

Investigation of the nature of the methionine-pi interaction in beta-hairpin peptide model systems

There are frequent contacts between aromatic rings and sulfur atoms in proteins. However, it is unclear to what degree this putative interaction is stabilizing and what the nature of the interaction is. We have investigated the aryl-sulfur interaction by placing a methionine residue diagonal to an aromatic ring on the same face of a beta-hairpin, which places the methionine side chain in close proximity to the aryl side chain. The methionine (Met)-aryl interaction was compared with an equivalent hydrophobic and cation-pi interaction in the context of the beta-hairpin. The interaction between phenylalanine (Phe), tryptophan (Trp), or cyclohexylalanine (Cha) and Met stabilized the beta-hairpin by -0.3 to -0.5 kcal mole(-1), as determined by double-mutant cycles. The peptides were subjected to thermal denaturations that suggest a hydrophobic driving force for the interactions between Met and Trp or Cha. The observed interaction of Met or norleucine (Nle) with Trp or Cha are quite similar, implying a hydrophobic driving force for the Met-pi interaction. However, the thermodynamic data suggest that there may be some differences between the interaction of Met with Trp and Phe and that there may be a small thermodynamic component to the Met...Phe interaction.     

Thermodynamics and stability of a beta-sheet complex: molecular dynamics simulations on simplified off-lattice protein models

We have performed discontinuous molecular dynamics simulations of the thermodynamics and stability of a tetrameric beta-sheet complex that contains four identical four-stranded antiparallel beta-sheet peptides. The potential used in the simulation is a hybrid Go-type potential characterized by the bias gap parameter g, an artificial measure of the preference of a model protein for its native state, and the intermolecular contact parameter eta, which measures the ratio of intermolecular to intramolecular native attractions. Despite the simplicity of the model, a complex set of thermodynamic transitions for the beta-sheet complex is revealed that shows there are three distinct oligomer (partially ordered, ordered, and highly ordered beta-sheet complex) states and four noninteracting monomers phases. The thermodynamic properties of the three oligomer states strongly depend on both the size of the intermolecular contact parameter eta and the temperature. The partially ordered beta-sheet complex is made up of four ordered globules and is observed at intermediate to large eta at high temperatures. The ordered beta-sheet complex contains four native beta-sheets and is located at small to intermediate eta at low temperatures in the phase diagram. The highly ordered beta-sheet complex has fully-stiff beta-sheet strands, the same as the global energy minimum structure, and is observed for all eta at low temperatures.     

EstB from Burkholderia gladioli: a novel esterase with a beta-lactamase fold reveals steric factors to discriminate between esterolytic and beta-lactam cleaving activity

Esterases form a diverse class of enzymes of largely unknown physiological role. Because many drugs and pesticides carry ester functions, the hydrolysis of such compounds forms at least one potential biological function. Carboxylesterases catalyze the hydrolysis of short chain aliphatic and aromatic carboxylic ester compounds. Esterases, D-alanyl-D-alanine-peptidases (DD-peptidases) and beta-lactamases can be grouped into two distinct classes of hydrolases with different folds and topologically unrelated catalytic residues, the one class comprising of esterases, the other one of beta-lactamases and DD-peptidases. The chemical reactivities of esters and beta-lactams towards hydrolysis are quite similar, which raises the question of which factors prevent esterases from displaying beta-lactamase activity and vice versa. Here we describe the crystal structure of EstB, an esterase isolated from Burkholderia gladioli. It shows the protein to belong to a novel class of esterases with homology to Penicillin binding proteins, notably DD-peptidase and class C beta-lactamases. Site-directed mutagenesis and the crystal structure of the complex with diisopropyl-fluorophosphate suggest Ser75 within the "beta-lactamase" Ser-x-x-Lys motif to act as catalytic nucleophile. Despite its structural homology to beta-lactamases, EstB shows no beta-lactamase activity. Although the nature and arrangement of active-site residues is very similar between EstB and homologous beta-lactamases, there are considerable differences in the shape of the active site tunnel. Modeling studies suggest steric factors to account for the enzyme's selectivity for ester hydrolysis versus beta-lactam cleavage.     

Solution structure of the PWWP domain of the hepatoma-derived growth factor family

Among the many PWWP-containing proteins, the largest group of homologous proteins is related to hepatoma-derived growth factor (HDGF). Within a well-conserved region at the extreme N-terminus, HDGF and five HDGF-related proteins (HRPs) always have a PWWP domain, which is a module found in many chromatin-associated proteins. In this study, we determined the solution structure of the PWWP domain of HDGF-related protein-3 (HRP-3) by NMR spectroscopy. The structure consists of a five-stranded beta-barrel with a PWWP-specific long loop connecting beta2 and beta3 (PR-loop), followed by a helical region including two alpha-helices. Its structure was found to have a characteristic solvent-exposed hydrophobic cavity, which is composed of an abundance of aromatic residues in the beta1/beta2 loop (beta-beta arch) and the beta3/beta4 loop. A similar ligand binding cavity occurs at the corresponding position in the Tudor, chromo, and MBT domains, which have structural and probable evolutionary relationships with PWWP domains. These findings suggest that the PWWP domains of the HDGF family bind to some component of chromatin via the cavity.     

Solution structure of termicin,an antimicrobial peptide from the termite Pseudacanthotermes spiniger

The solution structure of termicin from hemocytes of the termite Pseudacanthotermes spiniger was determined by proton two-dimensional nuclear magnetic resonance spectroscopy and molecular modeling techniques. Termicin is a cysteine-rich antifungal peptide also exhibiting a weak antibacterial activity. The global fold of termicin consists of an alpha-helical segment (Phe4-Gln14) and a two-stranded (Phe19-Asp25 and Gln28-Phe33) antiparallel beta-sheet forming a "cysteine stabilized alphabeta motif" (CSalphabeta) also found in antibacterial and antifungal defensins from insects and from plants. Interestingly, termicin shares more structural similarities with the antibacterial insect defensins and with MGD-1, a mussel defensin, than with the insect antifungal defensins such as drosomycin and heliomicin. These structural comparisons suggest that global fold alone does not explain the difference between antifungals and antibacterials. The antifungal properties of termicin may be related to its marked hydrophobicity and its amphipatic structure as compared to the antibacterial defensins. [SWISS-PROT accession number: Termicin (P82321); PDB accession number: 1MM0.]     

Ultrastructural localization of integrin subunits beta4 and alpha3 within the migrating epithelial tongue of in vivo human wounds

Subsequent to wounding, keratinocytes must quickly restore barrier function. In vitro wound models have served to elucidate mechanisms of epithelial closure and key roles for integrins alpha6beta4 and alpha3beta1. To extrapolate in vitro data to in vivo human tissues, we used ultrathin cryomicrotomy to simultaneously observe tissue ultrastructure and immunogold localization in unwounded skin and acute human cutaneous wounds. Localization of the beta4 integrin subunit in unwounded skin shows dominant hemidesmosomal association and minor basal keratinocyte lateral filopodic cell-cell expression. After wounding, beta4 dominantly localized to cytokeratin-rich regions (trailing edge hemidesmosomes) and minor association with lamellipodia (leading edge). beta4 colocalizes with alpha3 within filopodia juxtaposed to wound matrix, and increased concentrations of beta4 were found in cytoplasmic vesicles within basal keratinocytes of the migrating tongue. alpha3 integrin subunit dominantly localized to filopodia within basal keratinocyte lateral cell-cell interfaces in unwounded skin and both cell-cell and cell-matrix filopodic interactions in wounded skin. This study indicates that beta4 interacts with the extracellular environment through both stable and transient interactions and may be managed through a different endosomal trafficking pathway than alpha3. alpha3 integrin, despite its ability to respond to alternate ligands after wounding, does so through a single structure, the filopodia.     

Enzyme specificity under dynamic control II: Principal component analysis of α-lytic protease using global and local solvent boundary conditions

The contributions of conformational dynamics to substrate specificity have been examined by the application of principal component analysis to molecular dynamics trajectories of alpha-lytic protease. The wild-type alpha-lytic protease is highly specific for substrates with small hydrophobic side chains at the specificity pocket, while the Met190-->Ala binding pocket mutant has a much broader specificity, actively hydrolyzing substrates ranging from Ala to Phe. Based on a combination of multiconformation analysis of cryo-X-ray crystallographic data, solution nuclear magnetic resonance (NMR), and normal mode calculations, we had hypothesized that the large alteration in specificity of the mutant enzyme is mainly attributable to changes in the dynamic movement of the two walls of the specificity pocket. To test this hypothesis, we performed a principal component analysis using 1-nanosecond molecular dynamics simulations using either a global or local solvent boundary condition. The results of this analysis strongly support our hypothesis and verify the results previously obtained by in vacuo normal mode analysis. We found that the walls of the wild-type substrate binding pocket move in tandem with one another, causing the pocket size to remain fixed so that only small substrates are recognized. In contrast, the M190A mutant shows uncoupled movement of the binding pocket walls, allowing the pocket to sample both smaller and larger sizes, which appears to be the cause of the observed broad specificity. The results suggest that the protein dynamics of alpha-lytic protease may play a significant role in defining the patterns of substrate specificity. As shown here, concerted local movements within proteins can be efficiently analyzed through a combination of principal component analysis and molecular dynamics trajectories using a local solvent boundary condition to reduce computational time and matrix size.     

Long-range side-chain-main-chain interactions play crucial roles in stabilizing the (betaalpha)8 barrel motif of the alpha subunit of tryptophan synthase

The role of hither-to-fore unrecognized long-range hydrogen bonds between main-chain amide hydrogens and polar side chains on the stability of a well-studied (betaalpha)8, TIM barrel protein, the alpha subunit of tryptophan synthase (alphaTS), was probed by mutational analysis. The F19-D46 and I97-D124 hydrogen bonds link the N terminus of a beta-strand with the C terminus of the succeeding antiparallel alpha-helix, and the A103-D130 hydrogen bond links the N terminus of an alpha-helix with the C terminus of the succeeding antiparallel beta-strand, forming clamps for the respective betaalpha or alphabeta hairpins. The individual replacement of these aspartic acid side chains with alanine leads to what appear to be closely related partially folded structures with significantly reduced far-UV CD ellipticity and thermodynamic stability. Comparisons with the effects of eliminating another main-chain-side-chain hydrogen bond, G26-S33, and two electrostatic side-chain-side-chain hydrogen bonds, D38-H92 and D112-H146, all in the same N-terminal folding unit of alphaTS, demonstrated a unique role for the clamp interactions in stabilizing the native barrel conformation. Because neither the asparagine nor glutamic acid variant at position 46 can completely reproduce the spectroscopic, thermodynamic, or kinetic folding properties of aspartic acid, both size and charge are crucial to its unique role in the clamp hydrogen bond. Kinetic studies suggest that the three clamp hydrogen bonds act in concert to stabilize the transition state leading to the fully folded TIM barrel motif.     

The role of side chain conformational flexibility in surface recognition by Tenebrio molitor antifreeze protein

Two-dimensional nuclear magnetic resonance spectroscopy was used to investigate the flexibility of the threonine side chains in the beta-helical Tenebrio molitor antifreeze protein (TmAFP) at low temperatures. From measurement of the (3)J(alphabeta) (1)H-(1)H scalar coupling constants, the chi(1) angles and preferred rotamer populations can be calculated. It was determined that the threonines on the ice-binding face of the protein adopt a preferred rotameric conformation at near freezing temperatures, whereas the threonines not on the ice-binding face sample many rotameric states. This suggests that TmAFP maintains a preformed ice-binding conformation in solution, wherein the rigid array of threonines that form the AFP-ice interface matches the ice crystal lattice. A key factor in binding to the ice surface and inhibition of ice crystal growth appears to be the close surface-to-surface complementarity between the AFP and crystalline ice, and the lack of an entropic penalty associated with freezing out motions in a flexible ligand.     

Energy landscape of a peptide consisting of alpha-helix, 3(10)-helix, beta-turn, beta-hairpin, and other disordered conformations

The energy landscape of a peptide [Ace-Lys-Gln-Cys-Arg-Glu-Arg-Ala-Nme] in explicit water was studied with a multicanonical molecular dynamics simulation, and the AMBER parm96 force field was used for the energy calculation. The peptide was taken from the recognition helix of the DNA-binding protein, c-MYB: A rugged energy landscape was obtained, in which the random-coil conformations were dominant at room temperature. The CD spectra of the synthesized peptide revealed that it is in the random state at room temperature. However, the 300 K canonical ensemble, Q(300K), contained alpha-helix, 3(10)-helix, beta-turn, and beta-hairpin structures with small but notable probabilities of existence. The complete alpha-helix, imperfect alpha-helix, and random-coil conformations were separated from one another in the conformational space. This means that the peptide must overcome energy barriers to form the alpha-helix. The overcoming process may correspond to the hydrogen-bond rearrangements from peptide-water to peptide-peptide interactions. The beta-turn, imperfect 3(10)-helix, and beta-hairpin structures, among which there are no energy barriers at 300 K, were embedded in the ensemble of the random-coil conformations. Two types of beta-hairpin with different beta-turn regions were observed in Q(300K). The two beta-hairpin structures may have different mechanisms for the beta-hairpin formation. The current study proposes a scheme that the random state of this peptide consists of both ordered and disordered conformations. In contrast, the energy landscape obtained from the parm94 force field was funnel like, in which the peptide formed the helical conformation at room temperature and random coil at high temperature.     

Minimum model for the alpha-helix-beta-hairpin transition in proteins

We present a minimal model for proteins, which is able to capture the structural conversion between the alpha-helix and beta-hairpin. In most regimes of the parameter space, the model produces a stable structure at a low temperature; in a few limited regimes of the parameter space, the model displays an beta-hairpin transition as the physical conditions vary. These variations include a perturbation on hydrogen bonding propensity at the middle of the modeled chain, or the change of the hydrophobicity of a designated pair along the chain. Using Monte Carlo simulations, we demonstrate the structural conversion by means of state diagrams, heat capacity maps, and free energy maps.     

Elongation of the BH8 β-hairpin peptide: Electrostatic interactions in β-hairpin formation and stability

An elongated version of the de novo designed beta-hairpin peptide, BH8, has allowed us to gain insight into the role of electrostatic interactions in beta-hairpin stability. A Lys-Glu electrostatic pair has been introduced by adding a residue at the beginning and at the end of the N-terminal and C-terminal strands, respectively, of the beta-hairpin structure, in both orientations. The two resulting peptides and controls having Ala residues at these positions and different combinations of Ala with Lys, or Glu residues, have been analyzed by nuclear magnetic resonance (NMR), under different pH and ionic strength conditions. All of the NMR parameters, in particular the conformational shift analysis of Calpha protons and the coupling constants, (3)J(HNalpha), correlate well and the population estimates are in reasonable agreement among the different methods used. In the most structured peptides, we find an extension of the beta-hairpin structure comprising the two extra residues. Analysis of the pH and salt dependence shows that ionic pairs contribute to beta-hairpin stability. The interaction is electrostatic in nature and can be screened by salt. There is also an important salt-independent contribution of negatively charged groups to the stability of this family of beta-hairpin peptides.     

Energy landscape and dynamics of the beta-hairpin G peptide and its isomers: Topology and sequences

We have investigated free energy landscape [MM/PBSA + normal modes entropy] of permutations in the G peptide (41-56) from the protein G B1 domain by studying six isomers corresponding to moving the hydrophobic cluster along the beta-strands (toward the turn: T1, AGEWTYDDKTFTVTET; T2, GEDTWDYATFTVTKTE; T3, GEDDWTYATFTVTKTE; toward the end: E1, WTYDDAGETKTFTVT; E2, WEYTGDDATKTETFTV; E3, WTYEGDDATKTETFTV). The free energy terms include molecular mechanics energy, Poisson-Boltzmann electrostatic solvation energy, surface area solvation energy, and conformational entropy estimated by using normal mode analysis. From the wild type to T1, then T3, and finally T2, we see a progressively changing energy landscape, toward a less stable beta-hairpin structure. Moving the hydrophobic cluster outside toward the end region causes a greater change in the energy landscape. alpha-Helical instead of a beta-hairpin structure was the most stable form for the E2 isomer. However, no matter how much the sequence changes, for all variants studied, ideal "native" beta-hairpin topologies remain as minima (regardless of whether global or local) in the energy landscape. In general, we find that the energy landscape is dependent on the hydrophobic cluster topology and on the sequence. Our present study indicates that the key is the relative conformational energies of the different conformations. Changes in the sequence strongly modulate the relative stabilities of topologically similar regions in the energy landscape, rather than redefine the topology space. This finding is consistent with a population redistribution in the process of protein folding. The limited variation of topological space, compared with the number of possible sequence changes, may relate to the observation that the number of known protein folds are far less than the sequential allowance.     

Structure-function studies on the active site of the coelenterazine-dependent luciferase from Renilla

Renilla luciferase (RLUC) is a versatile tool for gene expression assays and in vivo biosensor applications, but its catalytic mechanism remains to be elucidated. RLUC is evolutionarily related to the alpha/beta hydrolase family. Its closest known homologs are bacterial dehalogenases, raising the question of how a protein with a hydrolase fold can function as a decarboxylating oxygenase. Molecular docking simulations with the coelenterazine substrate against an RLUC homology model as well as a recently determined RLUC crystal structure were used to build hypotheses to identify functionally important residues, which were subsequently tested by site-directed mutagenesis, heterologous expression, and bioluminescence emission spectroscopy. The data highlighted two triads of residues that are critical for catalysis. The putative catalytic triad residues D120, E144, and H285 bear only limited resemblance to those found in the active site of aequorin, a coelenterazine-utilizing photoprotein, suggesting that the reaction scheme employed by RLUC differs substantially from the one established for aequorin. The role of H285 in catalysis was further supported by inhibition using diethylpyrocarbonate. Multiple substitutions of N53, W121, and P220--three other residues implicated in product binding in the homologous dehalogenase Sphingomonas LinB--also supported their involvement in catalysis. Together with luminescence spectra, our data lead us to propose that the conserved catalytic triad of RLUC is directly involved in the decarboxylation reaction of coelenterazine to produce bioluminescence, while the other active-site residues are used for binding of the substrate.     

Unfolding crystallins: the destabilizing role of a beta-hairpin cysteine in betaB2-crystallin by simulation and experiment

The thermodynamic and kinetic stabilities of the eye lens family of betagamma-crystallins are important factors in the etiology of senile cataract. They control the chance of proteins unfolding, which can lead to aggregation and loss of transparency. betaB2-Crystallin orthologs are of low stability and comprise two typical betagamma-crystallin domains, although, uniquely, the N-terminal domain has a cysteine in one of the conserved folded beta-hairpins. Using high-temperature (500 K) molecular dynamics simulations with explicit solvent on the N-terminal domain of rodent betaB2-crystallin, we have identified in silico local flexibility in this folded beta-hairpin. We have shown in vitro using two-domain human betaB2-crystallin that replacement of this cysteine with a more usual aromatic residue (phenylalanine) results in a gain in conformational stability and a reduction in the rate of unfolding. We have used principal components analysis to visualize and cluster the coordinates from eight separate simulated unfolding trajectories of both the wild-type and the C50F mutant N-terminal domains. These data, representing fluctuations around the native well, show that although the mutant and wild-type appear to behave similarly over the early time period, the wild type appears to explore a different region of conformational space. It is proposed that the advantage of having this low-stability cysteine may be correlated with a subunit-exchange mechanism that allows betaB2-crystallin to interact with a range of other beta-crystallin subunits.     

Characterization of two potentially universal turn motifs that shape the repeated five-residues fold--crystal structure of a lumenal pentapeptide repeat protein from Cyanothece 51142

The genome of the diurnal cyanobacterium Cyanothece sp. PCC 51142 has recently been sequenced and observed to contain 35 pentapeptide repeat proteins (PRPs). These proteins, while present throughout the prokaryotic and eukaryotic kingdoms, are most abundant in cyanobacteria. The sheer number of PRPs in cyanobacteria coupled with their predicted location in every cellular compartment argues for important, yet unknown, physiological and biochemical functions. To gain biochemical insights, the crystal structure for Rfr32, a 167-residue PRP with an N-terminal 29-residue signal peptide, was determined at 2.1 A resolution. The structure is dominated by 21 tandem pentapeptide repeats that fold into a right-handed quadrilateral beta-helix, or Rfr-fold, as observed for the tandem pentapeptide repeats in the only other PRP structure, the mycobacterial fluoroquinoline resistance protein MfpA from Mycobacterium tuberculosis. Sitting on top of the Rfr-fold are two short, antiparallel alpha-helices, bridged with a disulfide bond, that perhaps prevent edge-to-edge aggregation at the C terminus. Analysis of the main-chain (Phi,Psi) dihedral orientations for the pentapeptide repeats in Rfr32 and MfpA makes it possible to recognize the structural details for the two distinct types of four-residue turns adopted by the pentapeptide repeats in the Rfr-fold. These turns, labeled type II and type IV beta-turns, may be universal motifs that shape the Rfr-fold in all PRPs.     

Cofactor effects on the protein folding reaction: acceleration of alpha-lactalbumin refolding by metal ions

About 30% of proteins require cofactors for their proper folding. The effects of cofactors on the folding reaction have been investigated with alpha-lactalbumin as a model protein and metal ions as cofactors. Metal ions accelerate the refolding of alpha-lactalbumin by lessening the energy barrier between the molten globule state and the transition state, mainly by decreasing the difference of entropy between the two states. These effects are linked to metal ion binding to the protein in the native state. Hence, relationships between the metal affinities for the intermediate states and those for the native state are observed. Some residual specificity for the calcium ion is still observed in the molten globule state, this specificity getting closer in the transition state to that of the native state. The comparison between kinetic and steady-state data in association with the Phi value method indicates the binding of the metal ions on the unfolded state of alpha-lactalbumin. Altogether, these results provide insight into cofactor effects on protein folding. They also suggest new possibilities to investigate the presence of residual native structures in the unfolded state of protein and the effects of such structures on the protein folding reaction and on protein stability.