Flow-Induced β-Hairpin Folding of the Glycoprotein Ibα β-Switch

Flow-induced shear has been identified as a regulatory driving force in blood clotting. Shear induces β-hairpin folding of the glycoprotein Ibα β-switch which increases affinity for binding to the von Willebrand factor, a key step in blood clot formation and wound healing. Through 2.1-μs molecular dynamics simulations, we investigate the kinetics of flow-induced β-hairpin folding. Simulations sampling different flow velocities reveal that under flow, β-hairpin folding is initiated by hydrophobic collapse, followed by interstrand hydrogen-bond formation and turn formation. Adaptive biasing force simulations are employed to determine the free energy required for extending the unfolded β-switch from a loop to an elongated state. Lattice and freely jointed chain models illustrate how the folding rate depends on the entropic and enthalpic energy, the latter controlled by flow. The results reveal that the free energy landscape of the β-switch has two stable conformations imprinted on it, namely, loop and hairpin—with flow inducing a transition between the two.     

Macromolecular Crowding Modulates Folding Mechanism of α/β Protein Apoflavodoxin

Protein dynamics in cells may be different from those in dilute solutions in vitro, because the environment in cells is highly concentrated with other macromolecules. This volume exclusion because of macromolecular crowding is predicted to affect both equilibrium and kinetic processes involving protein conformational changes. To quantify macromolecular crowding effects on protein folding mechanisms, we investigated the folding energy landscape of an α/β protein, apoflavodoxin, in the presence of inert macromolecular crowding agents, using in silico and in vitro approaches. By means of coarse-grained molecular simulations and topology-based potential interactions, we probed the effects of increased volume fractions of crowding agents (ϕ_c) as well as of crowding agent geometry (sphere or spherocylinder) at high ϕ_c. Parallel kinetic folding experiments with purified Desulfovibro desulfuricans apoflavodoxin in vitro were performed in the presence of Ficoll (sphere) and Dextran (spherocylinder) synthetic crowding agents. In conclusion, we identified the in silico crowding conditions that best enhance protein stability, and discovered that upon manipulation of the crowding conditions, folding routes experiencing topological frustrations can be either enhanced or relieved. Our test-tube experiments confirmed that apoflavodoxin''s time-resolved folding path is modulated by crowding agent geometry. Macromolecular crowding effects may be a tool for the manipulation of protein-folding and function in living cells.     

Relationships between Structural Dynamics and Functional Kinetics in Oligomeric Membrane Receptors

Recent efforts to broaden understanding of the molecular mechanisms of membrane receptors in signal transduction make use of rate-equilibrium free-energy relationships (REFERs), previously applied to chemical reactions, enzyme kinetics, and protein folding. For oligomeric membrane receptors, we distinguish between a), the Leffler parameter α_L, to characterize the global transition state for the interconversion between conformations; and b), the Fersht parameter, ϕ_F, to assign the degree of progression of individual residue positions at the transition state. For both α_L and ϕ_F, insights are achieved by using harmonic energy profiles to reflect the dynamic nature of proteins, as illustrated with single-channel results reported for normal and mutant nicotinic receptors. We also describe new applications of α_L based on published results. For large-conductance calcium-activated potassium channels, data are satisfactorily fit with an α_L value of 0.65, in accord with REFERs. In contrast, results reported for the flip conformational state of glycine and nicotinic receptors are in disaccord with REFERs, since they yield α_L values outside the usual limits of 0–1. Concerning published ϕ_F values underlying the conformational wave hypothesis for nicotinic receptors, we note that interpretations may be complicated by variations in the width of harmonic energy profiles.     

Resolving the Negative Potential Side (n-side) Water-accessible Proton Pathway of F-type ATP Synthase by Molecular Dynamics Simulations

The rotation of F₁F_o-ATP synthase is powered by the proton motive force across the energy-transducing membrane. The protein complex functions like a turbine; the proton flow drives the rotation of the c-ring of the transmembrane F_o domain, which is coupled to the ATP-producing F₁ domain. The hairpin-structured c-protomers transport the protons by reversible protonation/deprotonation of a conserved Asp/Glu at the outer transmembrane helix (TMH). An open question is the proton transfer pathway through the membrane at atomic resolution. The protons are thought to be transferred via two half-channels to and from the conserved cAsp/Glu in the middle of the membrane. By molecular dynamics simulations of c-ring structures in a lipid bilayer, we mapped a water channel as one of the half-channels. We also analyzed the suppressor mutant cP24D/E61G in which the functional carboxylate is shifted to the inner TMH of the c-protomers. Current models concentrating on the “locked” and “open” conformations of the conserved carboxylate side chain are unable to explain the molecular function of this mutant. Our molecular dynamics simulations revealed an extended water channel with additional water molecules bridging the distance of the outer to the inner TMH. We suggest that the geometry of the water channel is an important feature for the molecular function of the membrane part of F₁F_o-ATP synthase. The inclination of the proton pathway isolates the two half-channels and may contribute to a favorable clockwise rotation in ATP synthesis mode.     

Compressing the free energy range of substructure stabilities in iso‐1‐cytochrome c

Evolutionary conservation of substructure architecture between yeast iso-1-cytochrome c and the well-characterized horse cytochrome c is studied with limited proteolysis, the alkaline conformational transition and global unfolding with guanidine-HCl. Mass spectral analysis of limited proteolysis cleavage products for iso-1-cytochrome c show that its least stable substructure is the same as horse cytochrome c. The limited proteolysis data yield a free energy of 3.8 ± 0.4 kcal mol⁻¹ to unfold the least stable substructure compared with 5.05 ± 0.30 kcal mol⁻¹ for global unfolding of iso-1-cytochrome c. Thus, substructure stabilities of iso-1-cytochrome c span only ∼1.2 kcal mol⁻¹ compared with ∼8 kcal mol⁻¹ for horse cytochrome c. Consistent with the less cooperative folding thus expected for the horse protein, the guanidine-HCl m-values are ∼3 kcal mol⁻¹M⁻¹ versus ∼4.5 kcal mol⁻¹M⁻¹ for horse versus yeast cytochrome c. The tight free energy spacing of the yeast cytochrome c substructures suggests that its folding has more branch points than for horse cytochrome c. Studies on a variant of iso-1-cytochrome c with an H26N mutation indicate that the least and most stable substructures unfold sequentially and the two least stable substructures unfold independently as for horse cytochrome c. Thus, important aspects of the substructure architecture of horse cytochrome c, albeit compressed energetically, are preserved evolutionally in yeast iso-1-cytochrome c.     

Energetic signatures of single base bulges: thermodynamic consequences and biological implications

DNA bulges are biologically consequential defects that can arise from template-primer misalignments during replication and pose challenges to the cellular DNA repair machinery. Calorimetric and spectroscopic characterizations of defect-containing duplexes reveal systematic patterns of sequence-context dependent bulge-induced destabilizations. These distinguishing energetic signatures are manifest in three coupled characteristics, namely: the magnitude of the bulge-induced duplex destabilization (ΔΔG_Bulge); the thermodynamic origins of ΔΔG_Bulge (i.e. enthalpic versus entropic); and, the cooperativity of the duplex melting transition (i.e. two-state versus non-two state). We find moderately destabilized duplexes undergo two-state dissociation and exhibit ΔΔG_Bulge values consistent with localized, nearest neighbor perturbations arising from unfavorable entropic contributions. Conversely, strongly destabilized duplexes melt in a non-two-state manner and exhibit ΔΔG_Bulge values consistent with perturbations exceeding nearest-neighbor expectations that are enthalpic in origin. Significantly, our data reveal an intriguing correlation in which the energetic impact of a single bulge base centered in one strand portends the impact of the corresponding complementary bulge base embedded in the opposite strand. We discuss potential correlations between these bulge-specific differential energetic profiles and their overall biological implications in terms of DNA recognition, repair and replication.     

Apparent two-state tendamistat folding is a sequential process along a defined route

The small all-beta-sheet protein tendamistat folds and unfolds rapidly in apparent two-state reactions. Kinetic measurements of two tendamistat variants under various solvent conditions reveal, however, that folding occurs in at least two sequential steps through a metastable obligatory intermediate. Depending on the solvent conditions either step can become rate limiting. The activation parameters indicate that the first step represents an enthalpic barrier whereas the second step is an entropic barrier at 25 degrees C. Our results suggest that initial non-specific collapse precedes formation of native secondary and tertiary structure in tendamistat folding. This points at a distinct route in tendamistat folding and indicates that partially folded metastable intermediates might play an important role in the mechanism of apparent two-state folding.     

Specificity of RSG-1.2 peptide binding to RRE-IIB RNA element of HIV-1 over Rev peptide is mainly enthalpic in origin

Rev is an essential HIV-1 regulatory protein which binds to the Rev responsive element (RRE) present within the env gene of HIV-1 RNA genome. This binding facilitates the transport of the RNA to the cytoplasm, which in turn triggers the switch between viral latency and active viral replication. Essential components of this complex have been localized to a minimal arginine rich Rev peptide and stem IIB region of RRE. A synthetic peptide known as RSG-1.2 binds with high binding affinity and specificity to the RRE-IIB than the Rev peptide, however the thermodynamic basis of this specificity has not yet been addressed. The present study aims to probe the thermodynamic origin of this specificity of RSG-1.2 over Rev Peptide for RRE-IIB. The temperature dependent melting studies show that RSG-1.2 binding stabilizes the RRE structure significantly (ΔT _m = 4.3°C), in contrast to Rev binding. Interestingly the thermodynamic signatures of the binding have also been found to be different for both the peptides. At pH 7.5, RSG-1.2 binds RRE-IIB with a K_a = 16.2±0.6×10⁷ M⁻¹ where enthalpic change ΔH = −13.9±0.1 kcal/mol is the main driving force with limited unfavorable contribution from entropic change TΔS = −2.8±0.1 kcal/mol. A large part of ΔH may be due to specific stacking between U72 and Arg15. In contrast binding of Rev (K_a = 3.1±0.4×10⁷ M⁻¹) is driven mainly by entropy (ΔH = 0 kcal/mol and TΔS = 10.2±0.2 kcal/mol) which arises from major conformational changes in the RNA upon binding.     

Assembly and activation of a kinase ribozyme

RNA activities can be regulated by modulating the relative energies of all conformations in a folding landscape; however, it is often unknown precisely how peripheral elements perturb the overall landscape in the absence of discrete alternative folds (inactive ensemble). This work explores the effects of sequence and secondary structure in governing kinase ribozyme activity. Kin.46 catalyzes thiophosphoryl transfer from ATPγS onto the 5′ hydroxyl of polynucleotide substrates, and is regulated 10,000-fold by annealing an effector oligonucleotide to form activator helix P4. Transfer kinetics for an extensive series of ribozyme variants identified several dispensable internal single-stranded segments, in addition to a potential pseudoknot at the active site between segments J1/4 and J3/2 that is partially supported by compensatory rescue. Standard allosteric mechanisms were ruled out, such as formation of discrete repressive structures or docking P4 into the rest of the ribozyme via backbone 2′ hydroxyls. Instead, P4 serves both to complete an important structural element (100-fold contribution to the reaction relative to a P4-deleted variant) and to mitigate nonspecific, inhibitory effects of the single-stranded tail (an additional 100-fold contribution to the apparent rate constant, k_obs). Thermodynamic activation parameters ΔH^‡ and ΔS^‡, calculated from the temperature dependence of k_obs, varied with tail length and sequence. Inhibitory effects of the unpaired tail are largely enthalpic for short tails and are both enthalpic and entropic for longer tails. These results refine the structural view of this kinase ribozyme and highlight the importance of nonspecific ensemble effects in conformational regulation by peripheral elements.     

Probing Position-Dependent Diffusion in Folding Reactions Using Single-Molecule Force Spectroscopy

Folding of proteins and nucleic acids involves a diffusive search over a multidimensional conformational energy landscape for the minimal-energy structure. When examining the projection of conformational motions onto a one-dimensional reaction coordinate, as done in most experiments, the diffusion coefficient D is generally position dependent. However, it has proven challenging to measure such position-dependence experimentally. We investigated the position-dependence of D in the folding of DNA hairpins as a simple model system in two ways: first, by analyzing the round-trip time to return to a given extension in constant-force extension trajectories measured by force spectroscopy, and second, by analyzing the fall time required to reach a given extension in force jump measurements. These methods yielded conflicting results: the fall time implied a fairly constant D, but the round-trip time implied variations of over an order of magnitude. Comparison of experiments with computational simulations revealed that both methods were strongly affected by experimental artifacts inherent to force spectroscopy measurements, which obscured the intrinsic position-dependence of D. Lastly, we applied Kramers’s theory to the kinetics of hairpins with energy barriers located at different positions along the hairpin stem, as a crude probe of D at different stem positions, and we found that D did not vary much as the barrier was moved along the reaction coordinate. This work underlines the difficulties faced when trying to deduce position-dependent diffusion coefficients from experimental folding trajectories.     

Conformational Transitions of Subunit ? in ATP Synthase from Thermophilic Bacillus PS3

Subunit ɛ of bacterial and chloroplast F_OF₁-ATP synthase is responsible for inhibition of ATPase activity. In Bacillus PS3 enzyme, subunit ɛ can adopt two conformations. In the “extended”, inhibitory conformation, its two C-terminal α-helices are stretched along subunit γ. In the “contracted”, noninhibitory conformation, these helices form a hairpin. The transition of subunit ɛ from an extended to a contracted state was studied in ATP synthase incorporated in Bacillus PS3 membranes at 59°C. Fluorescence energy resonance transfer between fluorophores introduced in the C-terminus of subunit ɛ and in the N-terminus of subunit γ was used to follow the conformational transition in real time. It was found that ATP induced the conformational transition from the extended to the contracted state (half-maximum transition extent at 140 μM ATP). ADP could neither prevent nor reverse the ATP-induced conformational change, but it did slow it down. Acid residues in the DELSEED region of subunit β were found to stabilize the extended conformation of ɛ. Binding of ATP directly to ɛ was not essential for the ATP-induced conformational change. The ATP concentration necessary for the half-maximal transition (140 μM) suggests that subunit ɛ probably adopts the extended state and strongly inhibits ATP hydrolysis only when the intracellular ATP level drops significantly below the normal value.     

The stability of Taq DNA polymerase results from a reduced entropic folding penalty; identification of other thermophilic proteins with similar folding thermodynamics

The thermal stability of Taq DNA polymerase is well known, and is the basis for its use in PCR. A comparative thermodynamic characterization of the large fragment domains of Taq (Klentaq) and E. coli (Klenow) DNA polymerases has been performed by obtaining full Gibbs‐Helmholtz stability curves of the free energy of folding (ΔG) versus temperature. This analysis provides the temperature dependencies of the folding enthalpy and entropy (ΔH and ΔS), and the heat capacity (ΔC_p) of folding. If increased or enhanced non‐covalent bonding in the native state is responsible for enhanced thermal stabilization of a protein, as is often proposed, then an enhanced favourable folding enthalpy should, in general, be observed for thermophilic proteins. However, for the Klenow–Klentaq homologous pair, the folding enthalpy (ΔH_fold) of Klentaq is considerably less favorable than that of Klenow at all temperatures. In contrast, it is found that Klentaq's extreme free energy of folding (ΔG_fold) originates from a significantly reduced entropic penalty of folding (ΔS_fold). Furthermore, the heat capacity changes upon folding are similar for Klenow and Klentaq. Along with this new data, comparable extended analysis of available thermodynamic data for 17 other mesophilic–thermophilic protein pairs (where enough applicable thermodynamic data exists) shows a similar pattern in seven of the 18 total systems. When analyzed with this approach, the more familiar “reduced ΔC_p mechanism” for protein thermal stabilization (observed in a different six of the 18 systems) frequently manifests as a temperature dependent shift from enthalpy driven stabilization to a reduced‐entropic‐penalty model. Proteins 2014; 82:785–793. © 2013 Wiley Periodicals, Inc.     

Sequential barriers and an obligatory metastable intermediate define the apparent two-state folding pathway of the ubiquitin-like PB1 domain of NBR1

The 90-residue N-terminal Phox and Bem1p (PB1) domain of NBR1 forms an α/β ubiquitin-like fold. Kinetic analysis using stopped-flow fluorescence reveals two-state kinetics; however, nonlinear effects in the denaturant dependence of the unfolding data demonstrate changes in the position of the rate-limiting barrier along the folding coordinate as the folding conditions change. The kinetics of wt-PB1 and several mutants show that this curvature is consistent with a single-pathway mechanism involving sequential transition states (TS1 and TS2) separated by a transiently populated high-energy intermediate, rather than movement of the transition state on a broad energy plateau. We show that the two transition states within the sequential model represent structurally and thermodynamically distinct species. TS1 is a collapsed state (α_TS1 = 0.71) with a large enthalpic barrier to formation that is rate-limiting under conditions that strongly favour folding. TS2 is highly native-like (α_TS2 = 0.93) and represents a late entropic barrier to formation of the native state. In support of the sequential transition state mechanism, we show that the G62A helix 2 substitution stabilises TS1 and the intermediate to such an extent that the latter becomes significantly populated, leading to the observation of a fast kinetic phase representing the initial U → I transition, with TS2 (α_TS2 = 0.87) becoming rate-limiting. The folding rate is not retarded by populating an intermediate, which would be expected for a misfold state, but is accelerated, suggesting that the I state is productive and on-pathway. The results show that the apparent two-state folding of the wt-PB1 domain occurs along a well-defined pathway involving structurally and thermodynamically distinct sequential transition states and an obligatory metastable intermediate that represents a productive local minimum in the energy landscape that increases the efficiency of barrier crossing through favourable effects on the entropy of activation.     

Thermodynamics of a diffusional protein folding reaction

The folding reactions of several proteins are well described as diffusional barrier crossing processes, which suggests that they should be analyzed by Kramers' rate theory rather than by transition state theory. For the cold shock protein Bc-Csp from Bacillus caldolyticus, we measured stability and folding kinetics, as well as solvent viscosity as a function of temperature and denaturant concentration. Our analysis indicates that diffusional folding reactions can be treated by transition state theory, provided that the temperature and denaturant dependence of the solvent viscosity is properly accounted for, either at the level of the measured rate constants or of the calculated activation parameters. After viscosity correction the activation barriers for folding become less enthalpic and more entropic. The transition from an enthalpic to an entropic folding barrier with increasing temperature is, however, apparent in the data before and after this correction. It is a consequence of the negative activation heat capacity of refolding, which is independent of solvent viscosity. Bc-Csp and its mesophilic homolog Bs-CspB from Bacillus subtilis differ strongly in stability but show identical enthalpic and entropic barriers to refolding. The increased stability of Bc-Csp originates from additional enthalpic interactions that are established after passage through the activated state. As a consequence, the activation enthalpy of unfolding is increased relative to Bs-CspB.     

Microsecond Barrier-Limited Chain Collapse Observed by Time-Resolved FRET and SAXS

It is generally held that random-coil polypeptide chains undergo a barrier-less continuous collapse when the solvent conditions are changed to favor the fully folded native conformation. We test this hypothesis by probing intramolecular distance distributions during folding in one of the paradigms of folding reactions, that of cytochrome c. The Trp59-to-heme distance was probed by time-resolved Förster resonance energy transfer in the microsecond time range of refolding. Contrary to expectation, a state with a Trp59–heme distance close to that of the guanidinium hydrochloride (GdnHCl) denatured state is present after ~ 27 μs of folding. A concomitant decrease in the population of this state and an increase in the population of a compact high-FRET (Förster resonance energy transfer) state (efficiency > 90%) show that the collapse is barrier limited. Small-angle X-ray scattering (SAXS) measurements over a similar time range show that the radius of gyration under native favoring conditions is comparable to that of the GdnHCl denatured unfolded state. An independent comprehensive global thermodynamic analysis reveals that marginally stable partially folded structures are also present in the nominally unfolded GdnHCl denatured state. These observations suggest that specifically collapsed intermediate structures with low stability in rapid equilibrium with the unfolded state may contribute to the apparent chain contraction observed in previous fluorescence studies using steady-state detection. In the absence of significant dynamic averaging of marginally stable partially folded states and with the use of probes sensitive to distance distributions, barrier-limited chain contraction is observed upon transfer of the GdnHCl denatured state ensemble to native-like conditions.     

Accommodating Discontinuities in Dimeric Left-Handed Coiled Coils in ATP Synthase External Stalks

ATP synthases from coupling membranes are complex rotary motors that convert the energy of proton gradients across coupling membranes into the chemical potential of the β-γ anhydride bond of ATP. Proton movement within the ring of c subunits localized in the F₀-sector drives γ and ɛ rotation within the F₁α₃β₃ catalytic core where substrates are bound and products are released. An external stalk composed of homodimeric subunits b₂ in Escherichia coli or heterodimeric bb′ in photosynthetic synthases connects F₀ subunit a with F₁ subunits δ and most likely α. The external stalk resists rotation, and is of interest both functionally and structurally. Hypotheses that the external stalk contributes to the overall efficiency of the reaction through elastic coupling of rotational substeps, and that stalks form staggered, right-handed coiled coils, are investigated here. We report on different structures that accommodate heptad discontinuities with either local or global underwinding. Analyses of the knob-and-hole packing of the E. coli b₂ and Synechocystis bb′ stalks strongly support the possibility that these proteins can adopt conventional left-handed coiled coils.     

Observation of Continuous Contraction and a Metastable Misfolded State during the Collapse and Folding of a Small Protein

To obtain proper insight into how structure develops during a protein folding reaction, it is necessary to understand the nature and mechanism of the polypeptide chain collapse reaction, which marks the initiation of folding. Here, the time-resolved fluorescence resonance energy transfer technique, in which the decay of the fluorescence light intensity with time is used to determine the time evolution of the distribution of intra-molecular distances, has been utilized to study the folding of the small protein, monellin. It is seen that when folding begins, about one-third of the protein molecules collapse into a molten globule state (I_MG), from which they relax by continuous further contraction to transit to the native state. The larger fraction gets trapped into a metastable misfolded state. Exit from this metastable state occurs via collapse to the lower free energy I_MG state. This exit is slow, on a time scale of seconds, because of activation energy barriers. The trapped misfolded molecules as well as the I_MG molecules contract continuously and slowly as structure develops. A phenomenological model of Markovian evolution of the polymer chain undergoing folding, incorporating these features, has been developed, which fits well the experimentally observed time evolution of distance distributions. The observation that the “wrong turn” to the misfolded state has not been eliminated by evolution belies the common belief that the folding of functional protein sequences is very different from that of a random heteropolymer, and the former have been selected by evolution to fold quickly.     

Kinetic and thermodynamic studies of peptidyltransferase in ribosomes from the extreme thermophile Thermus thermophilus

Throughout evolution, emerging organisms survived by adapting existing biochemical processes to new reaction conditions. Simple protein enzymes balanced changes in structural stability with changes that permitted optimal catalysis by adjustments in both entropic and enthalpic contributions to the free energy of activation for the reaction. Study of adaptive mechanisms by large multicomponent enzymes such as the ribosome has been largely unexplored. Here we have determined the kinetic and thermodynamic parameters of peptidyltransferase in ribosomes from the extreme thermophile Thermus thermophilus. Activity of thermophilic enzymes can be assayed over a wide range of temperatures, enabling one to measure accurate catalytic rates and determine enthalpic and entropic contributions to the free energy of activation of the reaction. Differences in the reaction conditions used here and in published studies on mesophilic ribosomes prevent direct comparison, but our data on Thermus ribosomes suggest that these ribosomes have adapted to changing environments using the same strategies as simple protein enzymes, balancing stability and flexibility without loss of catalytic rate. This strategy must be a very ancient process, perhaps first used by primitive ribosomes in the RNA World.     

Folding is not required for bilayer insertion: replica exchange simulations of an alpha-helical peptide with an explicit lipid bilayer

We implement the replica exchange molecular dynamics algorithm to study the interactions of a model peptide (WALP-16) with an explicitly represented DPPC membrane bilayer. We observe the spontaneous, unbiased insertion of WALP-16 into the DPPC bilayer and its folding into an alpha-helix with a transbilayer orientation. The free energy surface suggests that the insertion of the peptide into the DPPC bilayer precedes secondary structure formation. Although the peptide has some propensity to form a partially helical structure in the interfacial region of the DPPC/water system, this state is not a productive intermediate but rather an off-pathway trap for WALP-16 insertion. Equilibrium simulations show that the observed insertion/folding pathway mirrors the potential of mean force (PMF). Calculation of the enthalpic and entropic contributions to this PMF show that the surface bound conformation of WALP-16 is significantly lower in energy than other conformations, and that the insertion of WALP-16 into the bilayer without regular secondary structure is enthalpically unfavorable by 5-10 kcal/mol/residue. The observed insertion/folding pathway disagrees with the dominant conceptual model, which is that a surface-bound helix is an obligatory intermediate for the insertion of alpha-helical peptides into lipid bilayers. In our simulations, the observed insertion/folding pathway is favored because of a large (>100 kcal/mol) increase in system entropy that occurs when the unstructured WALP-16 peptide enters the lipid bilayer interior. The insertion/folding pathway that is lowest in free energy depends sensitively on the near cancellation of large enthalpic and entropic terms. This suggests the possibility that intrinsic membrane peptides may have a diversity of insertion/folding behaviors depending on the exact system of peptide and lipid under consideration.