Energetic coupling between native-state prolyl isomerization and conformational protein folding

In folded proteins, prolyl peptide bonds are usually thought to be either trans or cis because only one of the isomers can be accommodated in the native folded protein. For the N-terminal domain of the gene-3 protein of the filamentous phage fd (N2 domain), Pro161 resides at the tip of a beta hairpin and was found to be cis in the crystal structure of this protein. Here we show that Pro161 exists in both the cis and the trans conformations in the folded form of the N2 domain. We investigated how conformational folding and prolyl isomerization are coupled in the unfolding and refolding of N2 domain. A combination of single-mixing and double-mixing unfolding and refolding experiments showed that, in unfolded N2 domain, 7% of the molecules contain a cis-Pro161 and 93% of the molecules contain a trans-Pro161. During refolding, the fraction of molecules with a cis-Pro161 increases to 85%. This implies that 10.3 kJ mol(-1) of the folding free energy was used to drive this 75-fold change in the Pro161 cis/trans equilibrium constant during folding. The stabilities of the forms with the cis and the trans isomers of Pro161 and their folding kinetics could be determined separately because their conformational folding is much faster than the prolyl isomerization reactions in the native and the unfolded proteins. The energetic coupling between conformational folding and Pro161 isomerization is already fully established in the transition state of folding, and the two isomeric forms are thus truly native forms. The folding kinetics are well described by a four-species box model, in which the N2 molecules with either isomer of Pro161 can fold to the native state and in which cis/trans isomerization occurs in both the unfolded and the folded proteins.     

Folding and association of the antibody domain CH3: prolyl isomerization preceeds dimerization.

The simplest naturally occurring model system for studying immunoglobulin folding and assembly is the non-covalent homodimer formed by the C-terminal domains (CH3) of the heavy chains of IgG. Here, we describe the structure of recombinant CH3 dimer as determined by X-ray crystallography and an analysis of the folding pathway of this protein.Under conditions where prolyl isomerization does not contribute to the folding kinetics, formation of the beta-sandwich structure is the rate-limiting step. beta-Sheet formation of CH3 is a slow process, even compared to other antibody domains, while the subsequent association of the folded monomers is fast. After long-time denaturation, the majority of the unfolded CH3 molecules reaches the native state in two serial reactions, involving the re-isomerization of the Pro35-peptide bond to the cis configuration. The species with the wrong isomer accumulate as a monomeric intermediate. Importantly, the isomerization to the correct cis configuration is the prerequisite for dimerization of the CH3 domain. In contrast, in the Fab fragment of the same antibody, prolyl isomerization occurs after dimerization demonstrating that within one protein, comprised of highly homologous domains, both the kinetics of beta-sandwich formation and the stage at which prolyl isomerization occurs during the folding process can be completely different.     

Kinetic mechanism and catalysis of a native-state prolyl isomerization reaction

Folding of tendamistat is a rapid two-state process for the majority of the unfolded molecules. In fluorescence-monitored refolding kinetics about 8% of the unfolded molecules fold slowly (lambda=0.083s(-1)), limited by peptidyl-prolyl cis-trans isomerization. This is significantly less than expected from the presence of three trans prolyl-peptide bonds in the native state. In interrupted refolding experiments we detected an additional very slow folding reaction (lambda=0.008s(-1) at pH 2) with an amplitude of about 12%. This reaction is caused by the interconversion of a highly structured intermediate to native tendamistat. The intermediate has essentially native spectroscopic properties and about 2% of it remain populated in equilibrium after folding is complete. Catalysis by human cyclophilin 18 identifies this very slow reaction as a prolyl isomerization reaction. This shows that prolyl-isomerases are able to efficiently catalyze native state isomerization reactions, which allows them to influence biologically important regulatory conformational transitions. Folding kinetics of the proline variants P7A, P9A, P50A and P7A/P9A show that the very slow reaction is due to isomerization of the Glu6-Pro7 and Ala8-Pro9 peptide bonds, which are located in a region that makes strong backbone and side-chain interactions to both beta-sheets. In the P50A variant the very slow isomerization reaction is still present but native state heterogeneity is not observed any more, indicating a long-range destabilizing effect on the alternative native state relative to N. These results enable us to include all prolyl and non-prolyl peptide bond isomerization reactions in the folding mechanism of tendamistat and to characterize the kinetic mechanism and the energetics of a native-state prolyl isomerization reaction.     

The Folding Pathway of the Antibody VL Domain

Antibodies are modular proteins consisting of domains that exhibit a β-sandwich structure, the so-called immunoglobulin fold. Despite structural similarity, differences in folding and stability exist between different domains. In particular, the variable domain of the light chain V_L is unusual as it is associated with misfolding diseases, including the pathologic assembly of the protein into fibrillar structures. Here, we have analysed the folding pathway of a V_L domain with a view to determine features that may influence the relationship between productive folding and fibril formation. The V_L domain from MAK33 (murine monoclonal antibody of the subtype κ/IgG1) has not previously been associated with fibrillisation but is shown here to be capable of forming fibrils. The folding pathway of this V_L domain is complex, involving two intermediates in different pathways. An obligatory early molten globule-like intermediate with secondary structure but only loose tertiary interactions is inferred. The native state can then be formed directly from this intermediate in a phase that can be accelerated by the addition of prolyl isomerases. However, an alternative pathway involving a second, more native-like intermediate is also significantly populated. Thus, the protein can reach the native state via two distinct folding pathways. Comparisons to the folding pathways of other antibody domains reveal similarities in the folding pathways; however, in detail, the folding of the V_L domain is striking, with two intermediates populated on different branches of the folding pathway, one of which could provide an entry point for molecules diverted into the amyloid pathway.     

Combination of the human prolyl isomerase FKBP12 with unrelated chaperone domains leads to chimeric folding enzymes with high activity

Folding enzymes often use distinct domains for the binding of substrate proteins ("chaperone domains") and for the catalysis of slow folding reactions such as disulfide formation or prolyl isomerization. The human prolyl isomerase FKBP12 is a small single-domain protein without a chaperone domain. Its very low folding activity could previously be increased by inserting the chaperone domain from the homolog SlyD (sensitive-to-lysis protein D) of Escherichia coli. We now inserted three unrelated chaperone domains into human FKBP12: the apical domain of the chaperonin GroEL from E. coli, the chaperone domain of protein disulfide isomerase from yeast, or the chaperone domain of SurA from the periplasm of E. coli. All three conveyed FKBP12 with a high affinity for unfolded proteins and increased its folding activity. Substrate binding and release of the chimeric folding enzymes were found to be very fast. This allows rapid substrate transfer from the chaperone domain to the catalytic domain and ensures efficient rebinding of protein chains that were unable to complete folding. The advantage of having separate sites, first for generic protein binding and then for specific catalysis, explains why our construction of the artificial folding enzymes with foreign chaperone domains was successful.     

Effect of multiple prolyl isomerization reactions on the stability and folding kinetics of the notch ankyrin domain: experiment and theory

Studies on the folding kinetics of the Notch ankyrin domain have demonstrated that the major refolding phase is slow, the minor refolding phase is limited by the isomerization of prolyl peptide bonds, and that unfolding is multiexponential. Here, we explore the relationship between prolyl isomerization and folding heterogeneity using a combination of experiment and simulation. Proline residues were replaced with alanine, both singly and in various combinations. These destabilizing substitutions combine to eliminate the minor refolding phase, although unfolding heterogeneity persists even when all seven proline residues are replaced. To test whether prolyl isomerization influences the major refolding phase, we modeled folding and prolyl isomerization as a system of sequential reactions. Simulations that use rate constants of the major folding phase of the Notch ankyrin domain to represent intrinsic folding indicate that even with seven prolyl isomerization reactions, only two significant phases should be observed, and that the fast observed phase provides a good approximation of the intrinsic folding in the absence of prolyl isomerization. These results indicate that the major refolding phase of the Notch ankyrin domain reflects an intrinsically slow folding transition, rather than coupling of fast folding events with slow prolyl isomerization steps. This is consistent with the observation that the single observed refolding phase of a construct in which all proline residues are replaced remains slow. Finally, the simulation fails to produce a second unfolding phase at high urea concentrations, indicating that prolyl isomerization does not play a role in the three-state mechanism that leads to this heterogeneity.     

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

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

Multiple parallel-pathway folding of proline-free Staphylococcal nuclease

When a protein exhibits complex kinetics of refolding, we often ascribe the complexity to slow isomerization events in the denatured protein, such as cis/trans isomerization of peptidyl prolyl bonds. Does the complex folding kinetics arise only from this well-known reason? Here, we have investigated the refolding of a proline-free variant of staphylococcal nuclease by stopped-flow, double-jump techniques, to examine the folding reactions without the slow prolyl isomerizations. As a result, the protein folds into the native state along at least two accessible parallel pathways, starting from a macroscopically single denatured-state ensemble. The presence of intermediates on the individual folding pathways has revealed the existence of multiple parallel pathways, and is characterized by multi-exponential folding kinetics with a lag phase. Therefore, a "single" amino acid sequence can fold along the multiple parallel pathways. This observation in staphylococcal nuclease suggests that the multiple folding may be more general than we have expected, because the multiple parallel-pathway folding cannot be excluded from proteins that show simpler kinetics.     

Nonprolyl cis peptide bonds in unfolded proteins cause complex folding kinetics

Folding of tendamistat, an inhibitor of alpha-amylase, is a fast two-state process accompanied by two minor slow reactions, which were assigned to prolyl isomerization. In a proline-free variant, 5% of the molecules still fold slowly with a rate constant of 2.5 s(-1). This reaction is caused by a slow equilibrium between two populations of unfolded molecules. The time constant for this equilibration process, its sensitivity to LiCl and its temperature dependence identify it as a cis-trans isomerization of nonprolyl peptide bonds. Although nonprolyl peptide bonds have the cis conformation populating only approximately 0.15% in unfolded proteins, their large number generates a significant fraction of slow-folding molecules. This emphasizes that heterogeneous populations in an unfolded protein can induce complex folding kinetics on various time scales.     

Generation of a Highly Active Folding Enzyme by Combining a Parvulin-Type Prolyl Isomerase from SurA with an Unrelated Chaperone Domain

Parvulins are small prolyl isomerases and serve as catalytic domains of folding enzymes. SurA (survival protein A) from the periplasm of Escherichia coli consists of an inactive (Par1) and an active (Par2) parvulin domain as well as a chaperone domain. In the absence of the chaperone domain, the folding activity of Par2 is virtually abolished. We created a chimeric protein by inserting the chaperone domain of SlyD, an unrelated folding enzyme from the FKBP family, into a loop of the isolated Par2 domain of SurA. This increased its folding activity 450-fold to a value higher than the activity of SurA, in which Par2 is linked with its natural chaperone domain. In the presence of both the natural and the foreign chaperone domain, the folding activity of Par2 was 1500-fold increased. Related and unrelated chaperone domains thus are similarly efficient in enhancing the folding activity of the prolyl isomerase Par2. A sequence analysis of various chaperone domains suggests that clusters of exposed methionine residues in mobile chain regions might be important for a generic interaction with unfolded protein chains. This binding is highly dynamic to allow frequent transfer of folding protein chains between chaperone and catalytic domains.     

Folding and oxidation of the antibody domain C(H)3

The non-covalent homodimer formed by the C-terminal domains of the IgG1 heavy chains (C(H)3) is the simplest naturally occurring model system for studying immunoglobulin folding and assembly. In the native state, the intrachain disulfide bridge, which connects a three-stranded and a four-stranded beta-sheet is buried in the hydrophobic core of the protein. Here, we show that the disulfide bridge is not required for folding and association, since the reduced C(H)3 domain folds to a dimer with defined secondary and tertiary structure. However, the thermodynamic stability of the reduced C(H)3 dimer is much lower than that of the oxidized state. This allows the formation of disulfide bonds either concomitant with folding (starting from the reduced, denatured state) or after folding (starting from the reduced dimer). The analysis of the two processes revealed that, under all conditions investigated, one of the cysteine residues, Cys 86, reacts preferentially with oxidized glutathione to a mixed disulfide that subsequently interacts with the less-reactive second thiol group of the intra-molecular disulfide bond. For folded C(H)3, the second step in the oxidation process is slow. In contrast, starting from the unfolded and reduced protein, the oxidation reaction is faster. However, the overall folding reaction of C(H)3 during oxidative folding is a slow process. Especially, dimerization is slow, compared to the association starting from the denatured oxidized state. This deceleration may be due to misfolded conformations trapped by the disulfide bridge.     

Molecular Determinants of a Native-State Prolyl Isomerization

Prolyl cis/trans isomerizations determine the rates of many protein-folding reactions, and they can serve as molecular switches and timers. The energy required to shift the prolyl cis/trans equilibrium during these processes originates from conformational reactions that are linked structurally and energetically with prolyl isomerization. We used the N2 domain of the gene-3-protein of phage fd to elucidate how such an energetic linkage develops in the course of folding. The Asp160-Pro161 bond at the tip of a β hairpin of N2 is cis in the crystal structure, but in fact, it exists as a mixture of conformers in folded N2. During refolding, about 10 kJ mol^− 1 of conformational energy becomes available for a 75-fold shift of the cis/trans equilibrium constant at Pro161, from 7/93 in the unfolded to 90/10 in the folded form. We combined single- and double-mixing kinetic experiments with a mutational analysis to identify the structural origin of this proline shift energy and to elucidate the molecular path for the transfer of this energy to Pro161. It originates largely, if not entirely, from the two-stranded β sheet at the base of the Pro161 hairpin. The two strands improve their stabilizing interactions when Pro161 is cis, and this stabilization is propagated to Pro161, because the connector peptides between the β strands and Pro161 are native-like folded when Pro161 is cis. In the presence of a trans-Pro161, the connector peptides are locally unfolded, and thus, Pro161 is structurally and energetically uncoupled from the β sheet. Such interrelations between local folding and prolyl isomerization and the potential modulation by prolyl isomerases might also be used to break and reestablish slow communication pathways in proteins.     

Reversible unfolding and refolding behavior of a monomeric aldolase from Staphylococcus aureus.

Thermal and GdmCl-induced unfolding transitions of aldolase from Staphylococcus aureus are reversible under a variety of solvent conditions. Analysis of the transitions reveals that no partially folded intermediates can be detected under equilibrium conditions. The stability of the enzyme is very low with a delta G0 value of -9 +/- 2 kJ/mol at 20 degrees C. The kinetics of unfolding and refolding of aldolase are complex and comprise at least one fast and two slow reactions. This complexity arises from prolyl isomerization reactions in the unfolded chain, which are kinetically coupled to the actual folding reaction. Comparison with model calculations shows that at least two prolyl peptide bonds give rise to the observed slow folding reactions of aldolase and that all of the involved bonds are presumably in the trans conformation in the native state. The rate constant of the actual folding reaction is fast with a relaxation time of about 15 s at the midpoint of the folding transition at 15 degrees C. The data presented on the folding and stability of aldolase are comparable to the properties of much smaller proteins. This might be connected with the simple and highly repetitive tertiary structure pattern of the enzyme, which belongs to the group of alpha/beta barrel proteins.     

Proline isomerization and the slow folding reactions of the alpha subunit of tryptophan synthase from Escherichia coli

Previous studies on the refolding of the alpha subunit of tryptophan synthase from Escherichia coli assigned two slow refolding phases to rate-limiting isomerizations of two 'essential' proline residues, one in each of the two domains of the protein (Matthews, C.R., Crisanti, M.M., Manz, J.T. and Gepner, G.L. (1983) Biochemistry 22, 1445-1452). The double-jump experiment (Brandts, J.F., Halvorson, H.R. and Brennan, M. (1975) Biochemistry 14, 4953-4963) was used to further investigate this phenomenon. The reaction assigned to the carboxyl domain is consistent with the proline isomerization hypothesis. The amino domain process is more rapid than expected for proline isomerization and may reflect another type of slow folding reaction. The results permit a further refinement of the folding model for the alpha subunit and demonstrate the existence of a third unfolded species whose folding is not limited by either of these two reactions.     

Osmolyte effects on kinetics of FKBP12 C22A folding coupled with prolyl isomerization

Unfolding and refolding kinetics of human FKBP12 C22A were monitored by fluorescence emission over a wide range of urea concentration in the presence and absence of protecting osmolytes glycerol, proline, sarcosine and trimethylamine-N-oxide (TMAO). Unfolding is well described by a mono-exponential process, while refolding required a minimum of two exponentials for an adequate fit throughout the urea concentration range considered. The bi-exponential behavior resulted from complex coupling between protein folding, and prolyl isomerization in the denatured state in which the urea-dependent rate constant for folding was greater than, equal to, and less than the rate constants for prolyl isomerization within the urea concentration range of zero to five molar. Amplitudes and the observed folding and unfolding rate constants were fitted to a reversible three-state model composed of two sequential steps involving the native state and a folding-competent denatured species thermodynamically linked to a folding-incompetent denatured species. Excellent agreement between thermodynamic parameters for FKBP12 C22A folding calculated from the kinetic parameters and those obtained directly from equilibrium denaturation assays provides strong support for the applicability of the mechanism, and provides evidence that FKBP12 C22A folding/unfolding is two-state, with prolyl isomer heterogeneity in the denatured ensemble. Despite the chemical diversity of the protecting osmolytes, they all exhibit the same kinetic behavior of increasing the rate constant of folding and decreasing the rate constant for unfolding. Osmolyte effects on folding/unfolding kinetics are readily explained in terms of principles established in understanding osmolyte effects on protein stability. These principles involve the osmophobic effect, which raises the Gibbs energy of the denatured state due to exposure of peptide backbone, thereby increasing the folding rate. This effect also plays a key role in decreasing the unfolding rate when, as is often the case, the activated complex exposes more backbone than is exposed in the native state.     

Interaction of trigger factor with the ribosome

The trigger factor of Escherichia coli is a prolyl isomerase and a chaperone. It interacts with the ribosome and affects the folding of newly formed protein chains. Therefore, the dynamics of the interactions of trigger factor with the ribosome and with unfolded protein chains should be tailored for this function. Previously, we had found that binding of unfolded proteins to trigger factor is fast and that the lifetime of the complex between these two components is only about 100 ms. Here, we have labeled the trigger factor in its amino-terminal, ribosome-binding domain with a fluorescent dye and investigated how it interacts with the ribosome. We found that this association, as well as the dissociation of the complex, are rather slow processes. The average lifetime of the complex is about 30 seconds (at 20 degrees C). The strong differences in the dynamics of the interactions of trigger factor with the ribosome and with protein substrates might ensure that, on the one hand, trigger factor remains bound to the ribosome while a protein chain is being synthesized, and, on the other hand, allows it to scan the newly formed protein for prolyl bonds that need catalysis of isomerization.     

Cooperation of the prolyl isomerase and chaperone activities of the protein folding catalyst SlyD

The SlyD (sensitive to lysis D) protein of Escherichia coli is a folding enzyme with a chaperone domain and a prolyl isomerase domain of the FK506 binding protein type. Here we investigated how the two domains and their interplay are optimized for function in protein folding. Unfolded protein molecules initially form a highly dynamic complex with the chaperone domain of SlyD, and they are then transferred to the prolyl isomerase domain. The turnover number of the prolyl isomerase site is very high and guarantees that, after transfer, prolyl peptide bonds in substrate proteins are isomerized very rapidly. The Michaelis constant of catalyzed folding reflects the substrate affinity of the chaperone domain, and the turnover number is presumably determined by the rate of productive substrate transfer from the chaperone to the prolyl isomerase site and by the intrinsic propensity of the refolding protein chain to leave the active site with the native prolyl isomer. The efficiency of substrate transfer is high because dissociation from the chaperone site is very fast and because the two sites are close to each other. Protein molecules that left the prolyl isomerase site with an incorrect prolyl isomer can rapidly be re-bound by the chaperone domain because the association rate is very high as well.     

Evaluation of similarities in the cis/trans isomerase function of trigger factor and DnaK

Two functionally redundant enzymes, trigger factor and the hsp70 chaperone DnaK, have been found to assist de novo protein folding in E coli. Trigger factor is a peripheral peptidyl prolyl cis/trans isomerase (PPIase) of the large subunit of the ribosome. In contrast, DnaK displays two catalytic features: the secondary amide peptide bond cis/trans isomerase (APIase) function supplemented by the ATPase site. APIases accelerate the cis/trans isomerization of nonprolyl peptide bonds. Both enzymes have affinity for an unfolded polypeptide chain. The diminished low temperature cell viability in the presence of trigger factor variants with impaired PPlase activity indicates that the enhancement of folding rates plays a crucial role in protein folding in vivo. For the DnaK-mediated increase in the folding yield in vitro, the minimal model for APlase catalysis involves the catalyzed partitioning of a rapidly formed folding intermediate as could be inferred from the DnaK/DnaJ/GrpE/ATP-assisted refolding of GdmCl-denatured luciferase. Using three different peptide bond cis/trans isomerization assays in vitro, we could show that there is no overlapping substrate specificity of trigger factor and DnaK. We propose that only if trigger factor recruits supplementing molecules is it capable of exhibiting functional complementarity with DnaK in protein folding.     

Coupled kinetic traps in cytochrome c folding: His-heme misligation and proline isomerization

The effect of His-heme misligation on folding has been investigated for a triple mutant of yeast iso-2 cytochrome c (N26H,H33N,H39K iso-2). The variant contains a single misligating His residue at position 26, a location at which His residues are found in several cytochrome c homologues, including horse, tuna, and yeast iso-1. The amplitude for fast phase folding exhibits a strong initial pH dependence. For GdnHCl unfolded protein at an initial pH<5, the observed refolding at final pH 6 is dominated by a fast phase (tau(2f)=20 ms, alpha(2f)=90 %) that represents folding in the absence of misligation. For unfolded protein at initial pH 6, folding at final pH 6 occurs in a fast phase of reduced amplitude (alpha(2f) approximately 20 %) but the same rate (tau(2f)=20 ms), and in two slower phases (tau(m)=6-8 seconds, alpha(m) approximately 45 %; and tau(1b)=16-20 seconds, alpha(1b) approximately 35 %). Double jump experiments show that the initial pH dependence of the folding amplitudes results from a slow pH-dependent equilibrium between fast and slow folding species present in the unfolded protein. The slow equilibrium arises from coupling of the His protonation equilibrium to His-heme misligation and proline isomerization. Specifically, Pro25 is predominantly in trans in the unligated low-pH unfolded protein, but is constrained in a non-native cis isomerization state by His26-heme misligation near neutral pH. Refolding from the misligated unfolded form proceeds slowly due to the large energetic barrier required for proline isomerization and displacement of the misligated His26-heme ligand.     

The major outer membrane protein of Fusobacterium nucleatum (FomA) folds and inserts into lipid bilayers via parallel folding pathways

Membrane protein insertion and folding was studied for the major outer membrane protein of Fusobacterium nucleatum (FomA), which is a voltage-dependent general diffusion porin. The transmembrane domain of FomA forms a beta-barrel that is predicted to consist of 14 beta-strands. Here, unfolded FomA is shown to insert and fold spontaneously and quantitatively into phospholipid bilayers upon dilution of the denaturant urea, which was shown previously only for outer membrane protein A (OmpA) of Escherichia coli. Folding of FomA is demonstrated by circular dichroism and fluorescence spectroscopy, by SDS-polyacrylamide gel electrophoresis, and by single-channel recordings. Refolded FomA had a single-channel conductance of 1.1 nS at 1 M KCl, in agreement with the conductance of FomA isolated from membranes in native form. In contrast to OmpA, which forms a smaller eight-stranded beta-barrel domain, folding kinetics of the larger FomA were slower and provided evidence for parallel folding pathways of FomA into lipid bilayers. Two pathways were observed independent of membrane thickness with two different lipid bilayers, which were either composed of dicapryl phosphatidylcholine or dioleoyl phosphatidylcholine. This is the first observation of parallel membrane insertion and folding pathways of a beta-barrel membrane protein from an unfolded state in urea into lipid bilayers. The kinetics of both folding pathways depended on the chain length of the lipid and on temperature with estimated activation energies of 19 kJ/mol (dicapryl phosphatidylcholine) and 70 kJ/mol (dioleoyl phosphatidylcholine) for the faster pathways.