Molecular dynamics simulations of synthetic peptide folding

Because the time scale of protein folding is much greater than that of the widely used simulations of native structures, a detailed report of molecular dynamics simulations of folding has not been available. In this study, we Included the average solvent effect in the potential functions to simplify the calculation of the solvent effect and carried out long molecular dynamics simulations of the alanine-based synthetic peptides at 274 K. From either an extended or a randomly generated conformation, the simulations approached a helix-coil equilibrium in about 3 ns. The multiple minima problem did not prevent helix folding. The calculated helical ratio of Ac-AAQAAAAQAAAAQAAY-NH₂ was 47%, in good agreement with the circular dichroism measurement (about 50%). A helical segment with frayed ends was the most stable conformation, but the hydrophobic interaction favored the compact, distorted helix-turn-helix conformations. The transition between the two types of conformations occurred in a much larger time scale than helix propagation. The transient hydrogen bonds between the glutamine side chain and the backbone carbonyl group could reduce the free energy barrier of helix folding and unfolding. The substitution of a single alanine residue in the middle of the peptide with valine or glycine decreased the average helical ratio significantly, in agreement with experimental observations. © 1996 Wiley-Liss, Inc.     

Conformational and dynamic characterization of the molten globule state of an apomyoglobin mutant with an altered folding pathway.

Kinetic and equilibrium studies of apomyoglobin folding pathways and intermediates have provided important insights into the mechanism of protein folding. To investigate the role of intrinsic helical propensities in the apomyoglobin folding process, a mutant has been prepared in which Asn132 and Glu136 have been substituted with glycine to destabilize the H helix. The structure and dynamics of the equilibrium molten globule state formed at pH 4.1 have been examined using NMR spectroscopy. Deviations of backbone (13)C(alpha) and (13)CO chemical shifts from random coil values reveal high populations of helical structure in the A and G helix regions and in part of the B helix. However, the H helix is significantly destabilized compared to the wild-type molten globule. Heteronuclear [(1)H]-(15)N NOEs show that, although the polypeptide backbone in the H helix region is more flexible than in the wild-type protein, its motions are restricted by transient hydrophobic interactions with the molten globule core. Quench flow hydrogen exchange measurements reveal stable helical structure in the A and G helices and part of the B helix in the burst phase kinetic intermediate and confirm that the H helix is largely unstructured. Stabilization of structure in the H helix occurs during the slow folding phases, in synchrony with the C and E helices and the CD region. The kinetic and equilibrium molten globule intermediates formed by N132G/E136G are similar in structure. Although both the wild-type apomyoglobin and the mutant fold via compact helical intermediates, the structures of the intermediates and consequently the detailed folding pathways differ. Apomyoglobin is therefore capable of compensating for mutations by using alternative folding pathways within a common basic framework. Tertiary hydrophobic interactions appear to play an important role in the formation and stabilization of secondary structure in the H helix of the N132G/E136G mutant. These studies provide important insights into the interplay between secondary and tertiary structure formation in protein folding.     

Phosphorylation stabilizes the N‐termini of α‐helices

The role of phosphorylation in stabilizing the N‐termini of α‐helices is examined using computer simulations of model peptides. The models comprise either a phosphorylated or unphosphorylated serine at the helix N‐terminus, followed by nine alanines. Monte Carlo/stochastic Dynamics simulations were performed on the model helices. The simulations revealed a distinct stabilization of the helical conformation at the N‐terminus after phosphorylation. The stabilization was attributable to favorable electrostatic interactions between the phosphate and the helix backbone. However, direct helix capping by the phosphorylated sidechain was not observed. The results of the calculations are consistent with experimental evidence on the stabilization of helices by phosphates and other anions. © 1999 John Wiley & Sons, Inc. Biopoly 49: 225–233, 1999     

Conformational preferences of a short Aib/Ala‐based water‐soluble peptide as a function of temperature

The amino acid Aib predisposes a peptide to be helical with context‐dependent preference for either 3₁₀‐ or α‐ or a mixed helical conformation. Short peptides also show an inherent tendency to be unfolded. To characterize helical and unfolded states adopted by water‐soluble Aib‐containing peptides, the conformational preference of Ac‐Ala‐Aib‐Ala‐Lys‐Ala‐Aib‐Lys‐Ala‐Lys‐Ala‐Aib‐Tyr‐NH₂ was determined by CD, NMR and MD simulations as a function of temperature. Temperature‐dependent CD data indicated the contribution of two major components, each an admixture of helical and extended/polyproline II structures. Both right‐ and left‐handed helical conformations were detected from deconvolution of CD data and ¹³C NMR experiments. The presence of a helical backbone, more pronounced at the N‐terminal, and a temperature‐induced shift in α‐helix/3₁₀‐helix equilibrium, more pronounced at the C‐terminal, emerged from NMR data. Starting from polyproline II, the N‐terminal of the peptide folded into a helical backbone in MD simulations within 5 ns at 60°C. Longer simulations showed a mixed‐helical backbone to be stable over the entire peptide at 5°C while at 60°C the mixed‐helix was either stable at the N‐terminus or occurred in short stretches through out the peptide, along with a significant population of polyproline II. Our results point towards conformational heterogeneity of water‐soluble Aib‐based peptide helices and the associated subtleties. The problem of analyzing CD and NMR data of both left‐ and right‐handed helices are discussed, especially the validity of the ellipticity ratio [θ]₂₂₂/[θ]₂₀₇, as a reporter of α‐/3₁₀‐ population ratio, in right‐ and left‐handed helical mixtures. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Sequence and conformational preferences at termini of α‐helices in membrane proteins: Role of the helix environment

α‐helices are amongst the most common secondary structural elements seen in membrane proteins and are packed in the form of helix bundles. These α‐helices encounter varying external environments (hydrophobic, hydrophilic) that may influence the sequence preferences at their N and C‐termini. The role of the external environment in stabilization of the helix termini in membrane proteins is still unknown. Here we analyze α‐helices in a high‐resolution dataset of integral α‐helical membrane proteins and establish that their sequence and conformational preferences differ from those in globular proteins. We specifically examine these preferences at the N and C‐termini in helices initiating/terminating inside the membrane core as well as in linkers connecting these transmembrane helices. We find that the sequence preferences and structural motifs at capping (Ncap and Ccap) and near‐helical (N' and C') positions are influenced by a combination of features including the membrane environment and the innate helix initiation and termination property of residues forming structural motifs. We also find that a large number of helix termini which do not form any particular capping motif are stabilized by formation of hydrogen bonds and hydrophobic interactions contributed from the neighboring helices in the membrane protein. We further validate the sequence preferences obtained from our analysis with data from an ultradeep sequencing study that identifies evolutionarily conserved amino acids in the rat neurotensin receptor. The results from our analysis provide insights for the secondary structure prediction, modeling and design of membrane proteins. Proteins 2014; 82:3420–3436. © 2014 Wiley Periodicals, Inc.     

Fundamental processes of protein folding: measuring the energetic balance between helix formation and hydrophobic interactions

Theories of protein folding often consider contributions from three fundamental elements: loops, hydrophobic interactions, and secondary structures. The pathway of protein folding, the rate of folding, and the final folded structure should be predictable if the energetic contributions to folding of these fundamental factors were properly understood. alphatalpha is a helix-turn-helix peptide that was developed by de novo design to provide a model system for the study of these important elements of protein folding. Hydrogen exchange experiments were performed on selectively 15N-labeled alphatalpha and used to calculate the stability of hydrogen bonds within the peptide. The resulting pattern of hydrogen bond stability was analyzed using a version of Lifson-Roig model that was extended to include a statistical parameter for tertiary interactions. This parameter, x, represents the additional statistical weight conferred upon a helical state by a tertiary contact. The hydrogen exchange data is most closely fit by the XHC model with an x parameter of 9.25. Thus the statistical weight of a hydrophobic tertiary contact is approximately 5.8x the statistical weight for helix formation by alanine. The value for the x parameter derived from this study should provide a basis for the understanding of the relationship between hydrophobic cluster formation and secondary structure formation during the early stages of protein folding.     

Conformational features of a hexapeptide model Ac-TGAAKA-NH2 corresponding to a hydrated alpha helical segment from glyceraldehyde 3-phosphate dehydrogenase: implications for the role of turns in helix folding

Recent analysis of alpha helices in protein crystal structures, available in literature, revealed hydrated alpha helical segments in which, water molecule breaks open helix 5-->1 hydrogen bond by inserting itself, hydrogen bonds to both C=O and NH groups of helix hydrogen bond without disrupting the helix hydrogen bond, and hydrogen bonds to either C=O or NH of helix hydrogen bond. These hydrated segments display a variety of turn conformations and are thought to be 'folding intermediates' trapped during folding-unfolding of alpha helices. A role for reverse turns is implicated in the folding of alpha helices. We considered a hexapeptide model Ac-1TGAAKA6-NH2 from glyceraldehyde 3-phosphate dehydrogenase, corresponding to a hydrated helical segment to assess its role in helix folding. The sequence is a site for two 'folding intermediates'. The conformational features of the model peptide have been investigated by 1H 2D NMR techniques and quantum mechanical perturbative configuration interaction over localized orbitals (PCILO) method. Theoretical modeling largely correlates with experimental observations. Based upon the amide proton temperature coefficients, the observed d alpha n(i, i + 1), d alpha n(i, i + 2), dnn(i, i + 1), d beta n(i, i + 1) NOEs and the results from theoretical modeling, we conclude that the residues of the peptide sample alpha helical and neck regions of the Ramachandran phi, psi map with reduced conformational entropy and there is a potential for turn conformations at N and C terminal ends of the peptide. The role of reduced conformational entropy and turn potential in helix formation have been discussed. We conclude that the peptide sequence can serve as a 'folding intermediate' in the helix folding of glyceraldehyde 3-phosphate dehydrogenase.     

Non‐native α‐helix formation is not necessary for folding of lipocalin: Comparison of burst‐phase folding between tear lipocalin and β‐lactoglobulin

Tear lipocalin and β‐lactoglobulin are members of the lipocalin superfamily. They have similar tertiary structures but unusually low overall sequence similarity. Non‐native helical structures are formed during the early stage of β‐lactoglobulin folding. To address whether the non‐native helix formation is found in the folding of other lipocalin superfamily proteins, the folding kinetics of a tear lipocalin variant were investigated by stopped‐flow methods measuring the time‐dependent changes in circular dichroism (CD) spectrum and small‐angle X‐ray scattering (SAXS). CD spectrum showed that extensive secondary structures are not formed during a burst‐phase (within a measurement dead time). The SAXS data showed that the radius of gyration becomes much smaller than in the unfolded state during the burst‐phase, indicating that the molecule is collapsed during an early stage of folding. Therefore, non‐native helix formation is not general for folding of all lipocalin family members. The non‐native helix content in the burst‐phase folding appears to depend on helical propensities of the amino acid sequence. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Helix folding simulations with various initial conformations.

Using a solvent-referenced energy calculation, a 16-residue peptide with alanine side chains folded into predominantly alpha-helical conformations during constant temperature (274 K) simulations. From different initial conformations, helical conformations were reached and the multiple energy minima did not become a serious problem. Under the same conditions, the simulation did not indiscriminately fold a sequence such as polyglycine into stable helices. Interesting observations from the simulations relate to the folding mechanism. The electrostatic interactions between the successive amides favored extended conformations (or beta strands) and caused energy barriers to helix folding. beta-bends were observed as intermediates during helix nucleation. The helix propagation toward the C-terminus seemed faster than that toward the N-terminus. In helical conformations, hydrogen bond oscillation between the i,i+ 4 and the i,i+3 patterns was observed. The i,i+3 hydrogen bonds occurred more frequently during helix propagation and deformation near both ends of the helical segment.     

Helix stability of oligoglycine,oligoalanine, and oligo‐β‐alanine dodecamers reflected by hydrogen‐bond persistence

Helices are important structural/recognition elements in proteins and peptides. Stability and conformational differences between helices composed of α‐ and β‐amino acids as scaffolds for mimicry of helix recognition has become a theme in medicinal chemistry. Furthermore, helices formed by β‐amino acids are experimentally more stable than those formed by α‐amino acids. This is paradoxical because the larger sizes of the hydrogen‐bonding rings required by the extra methylene groups should lead to entropic destabilization. In this study, molecular dynamics simulations using the second‐generation force field, AMOEBA (Ponder, J.W., et al., Current status of the AMOEBA polarizable force field. J Phys Chem B, 2010. 114 (8): p. 2549–64.) explored the stability and hydrogen‐bonding patterns of capped oligo‐β‐alanine, oligoalanine, and oligoglycine dodecamers in water. The MD simulations showed that oligo‐β‐alanine has strong acceptor+2 hydrogen bonds, but surprisingly did not contain a large content of 3₁₂‐helical structures, possibly due to the sparse distribution of the 3₁₂‐helical structure and other structures with acceptor+2 hydrogen bonds. On the other hand, despite its backbone flexibility, the β‐alanine dodecamer had more stable and persistent <3.0 Å hydrogen bonds. Its structure was dominated more by multicentered hydrogen bonds than either oligoglycine or oligoalanine helices. The 3₁ (PII) helical structure, prevalent in oligoglycine and oligoalanine, does not appear to be stable in oligo‐β‐alanine indicating its competition with other structures (stacking structure as indicated by MD analyses). These differences are among the factors that shape helical structural preferences and the relative stabilities of these three oligopeptides. Proteins 2014; 82:3043–3061. © 2014 Wiley Periodicals, Inc.     

Modular design of synthetic protein mimics. Characterization of the helical conformation of a 13-residue peptide in crystals

The incorporation of alpha-aminoisobutyryl (Aib) residues into peptide sequences facilitates helical folding. Aib-containing sequences have been chosen for the design of rigid helical segments in a modular approach to the construction of a synthetic protein mimic. The helical conformation of the synthetic peptide Boc-Aib-(Val-Ala-Leu-Aib)3-OMe in crystals is established by X-ray diffraction. The 13-residue apolar peptide adopts a helical form in the crystal with seven alpha-type hydrogen bonds in the middle and 3(10)-type hydrogen bonds at either end. The helices stack in columns, zigzag rather than linear, by means of direct NH...OC head to tail hydrogen bonds. Leucyl side chains are extended on one side of the helix and valyl side chains on the other side. Water molecules form hydrogen bonds with several backbone carbonyl oxygens that also participate in alpha-helix hydrogen bonds. There is no apparent distortion of the helix caused by hydration. The space group is P2(1)2(1)2(1), with a = 9.964 (3) A, b = 20.117 (3) A, c = 39.311 (6) A, Z = 4, and dx = 1.127 g/cm3 for C64H106N13O16.1.33H2O. The final agreement factor R was 0.089 for 3667 data observed greater than 3 sigma(F) with a resolution of 0.9 A.     

Conformational studies of retro–inverso peptides: The crystal and molecular structure of the hydantoin from H-Ala-G-Ala-mGly-OBzl

The crystal structure of the hydantoin 1-[(S)-1′–aminoethylmalonyl benzyl ester]-(S)-4-methylimidazolidin-2.5-dione (1) derived from the peptide H-Ala-gAla-mGly-OBzl, Having the retro-inverso modification of the Ala-Gly bond, has been determined by x-ray diffraction analysis. The crystals are orthorhombic, space group P2₁2₁2₁ with a = 6.539. b = 14.721, c = 17.101 Å, z = 4. The structure was solved by direct methods and refined with anisotropic thermal factors to a final R value of 0.067 for the 947 observed reflections. Reversal of the Ala-Gly amide bond perturbs the folding tendency of the backbone shown by the parent peptide t-BuCO-Ala-Gly-NHiPr. The gem-diamino residue, gAla, and the malonyl moieties are found in the helical and the extended conformations, respectively. Intramolecular hydrogen bonding is not observed. The molecules in the crystal are held together by the formation of two intermolecular hydrogen bonds of the N? H?O?C type with N?O distances of 2.86 and 3.17 Å respectively. 1995 John Wiley&Sons. Inc. © 1995 John Wiley & Sons, Inc.     

Polynucleotide folding in yeast tRNAPhe: elucidation of short-, medium-, and long-range interactions of sugar-phosphate-sugar backbone and base using a "blocked" nucleotide probe

It has recently been proposed that the repeating backbone nucleotide may be regarded as consisting of two blocks of equal magnitude representable by two virtual bonds. Implicit consideration of the nucleotide (ψ,ψ) and internucleotide (ω′,ω) geometry that generate variety in polynucleotide conformations, and of the constancy of the repeating structural moieties (P-C4′ and C4′-P) independent of the above rotations, has enabled us to utilize this scheme in the study of ordered structures such as di-, oligonucleotides and, most significantly, tRNA. The polynucleotide folding dictated by short-, intermediate-, and long-range interactions in the monoclinic and orthorhombic forms is described and compared through circular plots depicting the virtual bond torsions and distance plots constructed independently for backbone as well as bases. The torsions and the bond angles associated with the virtual bonds afford a clear distinction between ordered helical segments from loops and bends of tRNA. Lower virtual bond torsions (?60° to 60°) concomitant with higher values of virtual bond angles characterize various bend regions, while torsions around 160°–210° typify ordered helical strands. The distance plot elucidates the type of interaction associated with various sub-structures (helix–helix, helix–loop, and loop–loop) that form the constituents of different structural domains. Several other features such as the manifestation of the P₁₀ loop and the approximate twofold symmetry in the tRNA molecule are conspicuous on the distance plot.     

Structure, stability and folding of the alpha-helix

Pauling first described the alpha-helix nearly 50 years ago, yet new features of its structure continue to be discovered, using peptide model systems, site-directed mutagenesis, advances in theory, the expansion of the Protein Data Bank and new experimental techniques. Helical peptides in solution form a vast number of structures, including fully helical, fully coiled and partly helical. To interpret peptide results quantitatively it is essential to use a helix/coil model that includes the stabilities of all these conformations. Our models now include terms for helix interiors, capping, side-chain interactions, N-termini and 3(10)-helices. The first three amino acids in a helix (N1, N2 and N3) and the preceding N-cap are unique, as their amide NH groups do not participate in backbone hydrogen bonding. We surveyed their structures in proteins and measured their amino acid preferences. The results are predominantly rationalized by hydrogen bonding to the free NH groups. Stabilizing side-chain-side-chain energies, including hydrophobic interactions, hydrogen bonding and polar/non-polar interactions, were measured accurately in helical peptides. Helices in proteins show a preference for having approximately an integral number of turns so that their N- and C-caps lie on the same side. There are also strong periodic trends in the likelihood of terminating a helix with a Schellman or alpha L C-cap motif. The kinetics of alpha-helix folding have been studied with stopped-flow deep ultraviolet circular dichroism using synchrotron radiation as the light source; this gives a far superior signal-to-noise ratio than a conventional instrument. We find that poly(Glu), poly(Lys) and alanine-based peptides fold in milliseconds, with longer peptides showing a transient overshoot in helix content.     

Peptide design: Crystal structure of a helical peptide module attached to a potentially nonhelical amino terminal segment

The peptide Boc-Gly-Dpg-Gly-Val-Ala-Leu-Aib-Val-Ala-Leu-OMe has been designed to examine the structural consequences of placing a short segment with a low helix propensity at the amino terminus of a helical heptapeptide module. The Gly-Dpg-Gly segment is a potential connecting element in the synthetic construction of a helix-linker-helix motif. Crystal parameters for the peptide are P2₁, a = 8.651(3) Å, b = 46.826(13) Å, c = 16.245 Å, β = 90.13(3)*, Z = 4; 2 independent molecules/asymmetric unit. The structure reveals almost identical conformations for the two independent molecules. The backbone is completely helical for residues 2–9, with one 4 → 1 hydrogen bond and six 5 → 1 hydrogen bonds. The α,α-di-n-propylglycine residue adopts a helical conformation. Gly(1) adopts an extended conformation resulting in a nonhelical N-terminus, with the Boc group swinging away from the helix. The lateral association of helices in the b axis direction is unusual in that the helix axes are directed up or down (parallel or antiparallel) by pairs: ↓↓↑↑↓↓, etc. © 1996 John Wiley & Sons, Inc.     

Global optimisation by replica exchange with scaled hybrid Hamiltonians

A global optimisation scheme based on replica-exchange molecular dynamics simulation with scaled hybrid Hamiltonians is presented and applied to fold trpzip2 peptide from extended structures with explicit water model using only eight replicas. The algorithm is shown to be capable of reproducibly optimising the peptide to structures of root mean squared deviations less than 2.5 Å with respect to all heavy atoms of NMR models. Moreover, the large amount of structural data sampled in the optimisation process enables us to provide a possible folding mechanism. The transition state ensemble is characterised by a largely formed turn and a compact packing of tryptophan 2, 4 and 9. The tryptophan 2/11 pair is found to form at a very late stage of the folding process. The first (closest to the turn) and the fourth native backbone hydrogen bonds form earlier and a picture of strict zipping up of hydrogen bonds is not observed. It is demonstrated in the present study that this global optimisation method which integrates structure prediction with approximated conformational samplings may be of help to understand the folding puzzle.     

Dynamics of hydrogen bond desolvation in protein folding

As proteins fold, a progressive structuring, immobilization and eventual exclusion of water surrounding backbone hydrogen bonds takes place. This process turns hydrogen bonds into major determinants of the folding pathway and compensates for the penalty of desolvation of the backbone polar groups. Taken as an average over all hydrogen bonds in a native fold, this extent of protection is found to be nearly ubiquitous. It is dynamically crucial, determining a constraint in the long-time limit behavior of coarse-grained ab initio simulations. Furthermore, an examination of one of the longest available (1micros) all-atom simulations with explicit solvent reveals that this average extent of protection is a constant of motion for the folding trajectory. We propose how such a stabilization is best achieved by clustering five hydrophobes around the backbone hydrogen bonds, an arrangement that yields the optimal stabilization. Our results support and clarify the view that hydrophobic surface burial should be commensurate with hydrogen-bond formation and enable us to define a basic desolvation motif inherent to structure and folding dynamics.     

Structure of the hypothetical protein Ton 1535 from Thermococcus onnurineus NA1 reveals unique structural properties by a left‐handed helical turn in normal α‐solenoid protein

The crystal structure of Ton1535, a hypothetical protein from Thermococcus onnurineus NA1, was determined at 2.3 Å resolution. With two antiparallel α‐helices in a helix‐turn‐helix motif as a repeating unit, Ton1535 consists of right‐handed coiled N‐ and C‐terminal regions that are stacked together using helix bundles containing a left‐handed helical turn. One left‐handed helical turn in the right‐handed coiled structure produces two unique structural properties. One is the presence of separated concave grooves rather than one continuous concave groove, and the other is the contribution of α‐helices on the convex surfaces of the N‐terminal region to the extended surface of the concave groove of the C‐terminal region and vice versa. Proteins 2014; 82:1072–1078. © 2013 Wiley Periodicals, Inc.     

Chain conformation in polyretropeptides III: Design of a 310 helix using α,α-dialkylated amino acids and retropeptide bonds

Computer simulations have been used to design a polypeptide with a 3₁₀ helix conformation. The study has been been performed taking advantage of the intrinsic helix forming tendency of α-Aminoisobutyric acid. In order to avoid the formation of the α helix, which is the other common helical conformation adopted by α-Aminoisobutyric acid-based peptides, retropeptide bonds have been included in the sequence. Thus, retropeptides are not able to form the intramolecular hydrogen bonding interactions characteristic of the α helix. The influences of both the peptide length and the solvent have been examined and compared with those of the polypeptide without retropeptide bonds. Proteins 29:575–582,1997. © 1997 Wiley-Liss, Inc.     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin‐binding protein G from Streptococcus. IV. Implication for the mechanism of folding of the parent protein

A 34‐residue α/β peptide [IG(28–61)], derived from the C‐terminal part of the B3 domain of the immunoglobulin binding protein G from Streptoccocus, was studied using CD and NMR spectroscopy at various temperatures and by differential scanning calorimetry. It was found that the C‐terminal part (a 16‐residue‐long fragment) of this peptide, which corresponds to the sequence of the β‐hairpin in the native structure, forms structure similar to the β‐hairpin only at T = 313 K, and the structure is stabilized by non‐native long‐range hydrophobic interactions (Val47–Val59). On the other hand, the N‐terminal part of IG(28–61), which corresponds to the middle α‐helix in the native structure, is unstructured at low temperature (283 K) and forms an α‐helix‐like structure at 305 K, and only one helical turn is observed at 313 K. At all temperatures at which NMR experiments were performed (283, 305, and 313 K), we do not observe any long‐range connectivities which would have supported packing between the C‐terminal (β‐hairpin) and the N‐terminal (α‐helix) parts of the sequence. Such interactions are absent, in contrast to the folding pathway of the B domain of protein G, proposed recently by Kmiecik and Kolinski (Biophys J 2008, 94, 726–736), based on Monte‐Carlo dynamics studies. Alternative folding mechanisms are proposed and discussed. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 469–480, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com