Global Fold Research Articles

Decoupling immunomodulatory properties from lipid binding in the α-pore-forming toxin Sticholysin II

Sticholysin II (StII), a pore-forming toxin from the marine anemone Stichodactyla helianthus, enhances an antigen-specific cytotoxic T lymphocyte (CTL) response when co-encapsulated in liposomes with a model antigen. This capacity does not depend exclusively on its pore-forming activity and is partially supported by its ability to activate Toll-like receptor 4 (TLR4) in dendritic cells, presumably by interacting with this receptor or by triggering signaling cascades upon binding to lipid membrane. In order to investigate whether the lipid binding capacity of StII is required for immunomodulation, we designed a mutant in which the aromatic amino acids from the interfacial binding site Trp110, Tyr111 and Trp114 were substituted by Ala. In the present work, we demonstrated that StII3A keeps the secondary structure composition and global folding of StII, while it loses its lipid binding and permeabilization abilities. Despite this, StII3A upregulates dendritic cells maturation markers, enhances an antigen-specific effector CD8+ T cells response and confers antitumor protection in a preventive scenario in C57BL/6 mice. Our results indicate that a mechanism independent of its lipid binding ability is involved in the immunomodulatory capacity of StII, pointing to StII3A as a promising candidate to improve the reliability of the Sts-based vaccine platform.

Exploring the Kinetics and Thermodynamics of a Novel Histidine Ammonia-Lyase from Geobacillus kaustophilus.

Histidine ammonia-lyase (HAL) plays a pivotal role in the non-oxidative deamination of L-histidine to produce trans-urocanic, a crucial process in amino acid metabolism. This study examines the cloning, purification, and biochemical characterization of a novel HAL from Geobacillus kaustophilus (GkHAL) and eight active site mutants to assess their effects on substrate binding, catalysis, thermostability, and secondary structure. The GkHAL enzyme was successfully overexpressed and purified to homogeneity. Its primary sequence displayed 40.7% to 43.7% similarity with other known HALs and shared the same oligomeric structure in solution. Kinetic assays showed that GkHAL has optimal activity at 85 °C and pH 8.5, with high thermal stability even after preincubation at high temperatures. Mutations at Y52, H82, N194, and E411 resulted in a complete loss of catalytic activity, underscoring their essential role in enzyme function, while mutations at residues Q274, R280, and F325 did not abolish activity but did reduce catalytic efficiency. Notably, mutants R280K and F325Y displayed novel activity with L-histidinamide, expanding the substrate specificity of HAL enzymes. Circular dichroism (CD) analysis showed minor secondary structure changes in the mutants but no significant effect on global GkHAL folding. These findings suggest that GkHAL could be a promising candidate for potential biotechnological applications.

Increasing the Measured Effective Quantum Volume with Zero Noise Extrapolation

Quantum volume is a full-stack benchmark for near-term quantum computers. It quantifies the largest size of a square circuit which can be executed on the target device with reasonable fidelity. Error mitigation is a set of techniques intended to remove the effects of noise present in the computation of noisy quantum computers when computing an expectation value of interest. Effective quantum volume is a proposed metric that applies error mitigation to the quantum volume protocol to evaluate the effectiveness not only of the target device but also of the error mitigation algorithm. Digital zero-noise extrapolation is an error mitigation technique that estimates the noiseless expectation value using circuit folding to amplify errors by known scale factors and then extrapolating computed expectation values to the zero-noise limit. Here we demonstrate that zero-noise extrapolation, with global and local unitary folding with fractional scale factors, in conjunction with dynamical decoupling, can increase the effective quantum volume over the vendor-measured quantum volume. Specifically, we measure the effective quantum volume of four IBM Quantum superconducting processor units, obtaining values that are larger than the vendor-measured quantum volume on each device. This is the first such increase reported.

The importance of the location of the N-terminus in successful protein folding in vivo and in vitro

Protein folding in the cell often begins during translation. Many proteins fold more efficiently cotranslationally than when refolding from a denatured state. Changing the vectorial synthesis of the polypeptide chain through circular permutation could impact functional, soluble protein expression and interactions with cellular proteostasis factors. Here, we measure the solubility and function of every possible circular permutant (CP) of HaloTag in Escherichia coli cell lysate using a gel-based assay, and in living E. coli cells via FACS-seq. We find that 78% of HaloTag CPs retain protein function, though a subset of these proteins are also highly aggregation-prone. We examine the function of each CP in E. coli cells lacking the cotranslational chaperone trigger factor and the intracellular protease Lon and find no significant changes in function as a result of modifying the cellular proteostasis network. Finally, we biophysically characterize two topologically interesting CPs in vitro via circular dichroism and hydrogen-deuterium exchange coupled with mass spectrometry to reveal changes in global stability and folding kinetics with circular permutation. For CP33, we identify a change in the refolding intermediate as compared to wild-type (WT) HaloTag. Finally, we show that the strongest predictor of aggregation-prone expression in cells is the introduction of termini within the refolding intermediate. These results, in addition to our finding that termini insertion within the conformationally restrained core is most disruptive to protein function, indicate that successful folding of circular permutants may depend more on changes in folding pathway and termini insertion in flexible regions than on the availability of proteostasis factors.

Structural Analysis of the Drosophila melanogaster GSTome.

Glutathione transferase (GST) is a superfamily of ubiquitous enzymes, multigenic in numerous organisms and which generally present homodimeric structures. GSTs are involved in numerous biological functions such as chemical detoxification as well as chemoperception in mammals and insects. GSTs catalyze the conjugation of their cofactor, reduced glutathione (GSH), to xenobiotic electrophilic centers. To achieve this catalytic function, GSTs are comprised of a ligand binding site and a GSH binding site per subunit, which is very specific and highly conserved; the hydrophobic substrate binding site enables the binding of diverse substrates. In this work, we focus our interest in a model organism, the fruit fly Drosophila melanogaster (D. mel), which comprises 42 GST sequences distributed in six classes and composing its GSTome. The goal of this study is to describe the complete structural GSTome of D. mel to determine how changes in the amino acid sequence modify the structural characteristics of GST, particularly in the GSH binding sites and in the dimerization interface. First, we predicted the 3D atomic structures of each GST using the AlphaFold (AF) program and compared them with X-ray crystallography structures, when they exist. We also characterized and compared their global and local folds. Second, we used multiple sequence alignment coupled with AF-predicted structures to characterize the relationship between the conservation of amino acids in the sequence and their structural features. Finally, we applied normal mode analysis to estimate thermal B-factors of all GST structures of D. mel. Particularly, we extracted flexibility profiles of GST and identify key residues and motifs that are systematically involved in the ligand binding/dimerization processes and thus playing a crucial role in the catalytic function. This methodology will be extended to guide the in silico design of synthetic GST with new/optimal catalytic properties for detoxification applications.

Comparative analysis of RNA 3D structure prediction methods: towards enhanced modeling of RNA-ligand interactions.

Accurate RNA structure models are crucial for designing small molecule ligands that modulate their functions. This study assesses six standalone RNA 3D structure prediction methods-DeepFoldRNA, RhoFold, BRiQ, FARFAR2, SimRNAand Vfold2, excluding web-based tools due to intellectual property concerns. We focus on reproducing the RNA structure existing in RNA-small molecule complexes, particularly on the ability to model ligand binding sites. Using a comprehensive set of RNA structures from the PDB, which includes diverse structural elements, we found that machine learning (ML)-based methods effectively predict global RNA folds but are less accurate with local interactions. Conversely, non-ML-based methods demonstrate higher precision in modeling intramolecular interactions, particularly with secondary structure restraints. Importantly, ligand-binding site accuracy can remain sufficiently high for practical use, even if the overall model quality is not optimal. With the recent release of AlphaFold 3, we included this advanced method in our tests. Benchmark subsets containing new structures, not used in the training of the tested ML methods, show that AlphaFold 3's performance was comparable to other ML-based methods, albeit with some challenges in accurately modeling ligand binding sites. This study underscores the importance of enhancing binding site prediction accuracy and the challenges in modeling RNA-ligand interactions accurately.

FRET-guided modeling of nucleic acids.

The functional diversity of RNAs is encoded in their innate conformational heterogeneity. The combination of single-molecule spectroscopy and computational modeling offers new attractive opportunities to map structural transitions within nucleic acid ensembles. Here, we describe a framework to harmonize single-molecule Förster resonance energy transfer (FRET) measurements with molecular dynamics simulations and de novo structure prediction. Using either all-atom or implicit fluorophore modeling, we recreate FRET experiments in silico, visualize the underlying structural dynamics and quantify the reaction coordinates. Using multiple accessible-contact volumes as a post hoc scoring method for fragment assembly in Rosetta, we demonstrate that FRET can be used to filter a de novo RNA structure prediction ensemble by refuting models that are not compatible with in vitro FRET measurement. We benchmark our FRET-assisted modeling approach on double-labeled DNA strands and validate it against an intrinsically dynamic manganese(II)-binding riboswitch. We show that a FRET coordinate describing the assembly of a four-way junction allows our pipeline to recapitulate the global fold of the riboswitch displayed by the crystal structure. We conclude that computational fluorescence spectroscopy facilitates the interpretability of dynamic structural ensembles and improves the mechanistic understanding of nucleic acid interactions.

Molecular dynamics of structural effects of reactive carbonyl species derivate of lipid peroxidation on bovine serum albumin

BackgroundSerum albumin is the most abundant protein in the Mammalia blood plasma at where plays a decisive role in the transport wide variety of hydrophobic ligands. BSA undergoes oxidative modifications like the carbonylation by the reactive carbonyl species (RCSs) 4-hydroxy-2-nonenal (HNE), 4 hydroxy-2-hexenal (HHE), malondialdehyde (MDA) and 4-oxo-2-nonenal (ONE), among others. The structural and functional changes induced by protein carbonylation have been associated with the advancement of neurodegenerative, cardiovascular, metabolic and cancer diseases. MethodsTo elucidate structural effects of protein carbonylation with RCSs on BSA, parameters for six new non-standard amino acids were designated and molecular dynamics simulations of its mono‑carbonylated-BSA systems were conducted in the AMBER force field. Trajectories were evaluated by RMSD, RMSF, PCA, RoG and SASA analysis. ResultsAn increase in the conformational instability for all proteins modified with local changes were observed, without significant changes on the BSA global three-dimensional folding. A more relaxed compaction level and major solvent accessible surface area for modified systems was found. Four regions of high molecular fluctuation were identified in all modified systems, being the subdomains IA and IIIB those with the most remarkable local conformational changes. Regarding essential modes of domain movements, it was evidenced that the most representatives were those related to IA subdomain, while IIIB subdomain presented discrete changes. ConclusionsRCSs induces local structural changes on mono‑carbonylated BSA. Also, this study extends our knowledge on how carbonylation by RCSs induce structural effects on proteins.

Swirling electrolyte. Part 1. Hierarchy of catastrophic transitions in steady annular flows

This study extends our previous work (McCloughan & Suslov, J. Fluid Mech., vol. 887, 2020, A23), where the existence of a saddle-node bifurcation of steady axisymmetric electrolyte flows driven by the Lorentz force in a shallow annular domain was first reported. Here we perform further weakly nonlinear analysis over a wider range of the governing parameters to demonstrate that the previously reported saddle-node bifurcation is a local feature of a global fold catastrophe, which, in turn, is a part of cusp catastrophe occurring as the thickness of the fluid layer increases. The amplitude equation characterising multiple flow solutions in the finite vicinity of catastrophe points is derived. The sensitivity of its coefficients and solutions to the distance from the catastrophe points is assessed demonstrating the robustness of the used analytical procedure. The asymptotic flow solution past the catastrophe point is subsequently obtained and its topology is explored confirming the existence of the secondary circulation in the bulk of flow (two-tori background flow structure). The latter is argued to lead to the appearance of experimentally observable vortices on the fluid surface. The rigorous justification of this conjecture is to be given in Part 2 of the study.

Large deviations and conditioning for chaotic non-invertible deterministic maps: analysis via the forward deterministic dynamics and the backward stochastic dynamics

The large deviation properties of trajectory observables for chaotic non-invertible deterministic maps as studied recently by Smith (2022 Phys. Rev. E 106 L042202) and by Gutierrez et al (2023 arXiv:2304.13754) are revisited in order to analyze in detail the similarities and the differences with the case of stochastic Markov chains. More concretely, we focus on the simplest example displaying the two essential properties of local stretching and global folding, namely the doubling map xt+1=2xt[mod1] on the real-space interval x∈[0,1[ that can also be analyzed via the decomposition x=∑l=1+∞σl2l into binary coefficients σl=0,1 . The large deviation properties of trajectory observables can be studied either via deformations of the forward deterministic dynamics or via deformations of the backward stochastic dynamics. Our main conclusions concerning the construction of the corresponding Doob canonical conditioned processes are: (i) non-trivial conditioned dynamics can be constructed only in the backward stochastic perspective where the reweighting of existing transitions is possible, and not in the forward deterministic perspective; (ii) the corresponding conditioned steady state is not smooth on the real-space interval x∈[0,1[ and can be better characterized in the binary space σl=1,2,..,+∞ . As a consequence, the backward stochastic dynamics in the binary space are also the most appropriate framework to analyze higher levels of large deviations, and we obtain the explicit large deviations at level 2 for the probability of the empirical density of long backward trajectories.

Structural and dynamic mechanisms for coupled folding and tRNA recognition of a translational T-box riboswitch

T-box riboswitches are unique riboregulators where gene regulation is mediated through interactions between two highly structured RNAs. Despite extensive structural insights, how RNA-RNA interactions drive the folding and structural transitions of T-box to achieve functional conformations remains unclear. Here, by combining SAXS, single-molecule FRET and computational modeling, we elaborate the folding energy landscape of a translational T-box aptamer consisting of stems I, II and IIA/B, which Mg2+-induced global folding and tRNA binding are cooperatively coupled. smFRET measurements reveal that high Mg2+ stabilizes IIA/B and its stacking on II, which drives the pre-docking of I and II into a competent conformation, subsequent tRNA binding promotes docking of I and II to form a high-affinity tRNA binding groove, of which the essentiality of IIA/B and S-turn in II is substantiated with mutational analysis. We highlight a delicate balance among Mg2+, the intra- and intermolecular RNA-RNA interactions in modulating RNA folding and function.

Assessment of three-dimensional RNA structure prediction in CASP15.

The prediction of RNA three-dimensional structures remains an unsolved problem. Here, we report assessments of RNA structure predictions in CASP15, the first CASP exercise that involved RNA structure modeling. Forty-two predictor groups submitted models for at least one of twelve RNA-containing targets. These models were evaluated by the RNA-Puzzles organizers and, separately, by a CASP-recruited team using metrics (GDT, lDDT) and approaches (Z-score rankings) initially developed for assessment of proteins and generalized here for RNA assessment. The two assessments independently ranked the same predictor groups as first (AIchemy_RNA2), second (Chen), and third (RNAPolis and GeneSilico, tied); predictions from deep learning approaches were significantly worse than these top ranked groups, which did not use deep learning. Further analyses based on direct comparison of predicted models to cryogenic electron microscopy (cryo-EM) maps and x-ray diffraction data support these rankings. With the exception of two RNA-protein complexes, models submitted by CASP15 groups correctly predicted the global fold of the RNA targets. Comparisons of CASP15 submissions to designed RNA nanostructures as well as molecular replacement trials highlight the potential utility of current RNA modeling approaches for RNA nanotechnology and structural biology, respectively. Nevertheless, challenges remain in modeling fine details such as noncanonical pairs, in ranking among submitted models, and in prediction of multiple structures resolved by cryo-EM or crystallography.

Structural features of thyroglobulin linked to protein trafficking.

Thyroglobulin must pass endoplasmic reticulum (ER) quality control to become secreted for thyroid hormone synthesis. Defective thyroglobulin, blocked in trafficking, can cause hypothyroidism. Thyroglobulin is a large protein (~2750 residues) spanning regions I-II-III plus a C-terminal cholinesterase-like domain. The cholinesterase-like domain functions as an intramolecular chaperone for regions I-II-III, but the folding pathway leading to successful thyroglobulin trafficking remains largely unknown. Here, informed by the recent three-dimensional structure of thyroglobulin as determined by cryo-electron microscopy, we have bioengineered three novel classes of mutants yielding three entirely distinct quality control phenotypes. Specifically, upon expressing recombinant thyroglobulin, we find that first, mutations eliminating a disulfide bond enclosing a 200-amino acid loop in region I have surprisingly little impact on the ability of thyroglobulin to fold to a secretion-competent state. Next, we have identified a mutation on the surface of the cholinesterase-like domain that has no discernible effect on regional folding yet affects contact between distinct regions and thereby triggers impairment in the trafficking of full-length thyroglobulin. Finally, we have probed a conserved disulfide in the cholinesterase-like domain that interferes dramatically with local folding, and this defect then impacts on global folding, blocking the entire thyroglobulin in the ER. These data highlight variants with distinct effects on ER quality control, inhibiting domain-specific folding; folding via regional contact; neither; or both.

The NCIWEB Server: A Novel Implementation of the Noncovalent Interactions Index for Biomolecular Systems.

It is well-known that the activity and function of proteins is strictly correlated with their secondary, tertiary, and quaternary structures. Their biological role is regulated by their conformational flexibility and global fold, which, in turn, is largely governed by complex noncovalent interaction networks. Because of the large size of proteins, the analysis of their noncovalent interaction networks is challenging, but can provide insights into the energetics of conformational changes or protein-protein and protein-ligand interactions. The noncovalent interaction (NCI) index, based on the reduced density gradient, is a well-established tool for the detection of weak contacts in biological systems. In this work, we present a web-based application to expand the use of this index to proteins, which only requires a molecular structure as input and provides a mapping of the number, type, and strength of noncovalent interactions. Structure preparation is automated and allows direct importing from the PDB database, making this server (https://nciweb.dsi.upmc.fr) accessible to scientists with limited experience in bioinformatics. A quick overview of this tool and concise instructions are presented, together with an illustrative application.

Two-Dimensional NMR Spectroscopy of the G Protein-Coupled Receptor A2AAR in Lipid Nanodiscs.

Eight hundred and twenty-six human G protein-coupled receptors (GPCRs) mediate the actions of two-thirds of the human hormones and neurotransmitters and over one-third of clinically used drugs. Studying the structure and dynamics of human GPCRs in lipid bilayer environments resembling the native cell membrane milieu is of great interest as a basis for understanding structure-function relationships and thus benefits continued drug development. Here, we incorporate the human A2A adenosine receptor (A2AAR) into lipid nanodiscs, which represent a detergent-free environment for structural studies using nuclear magnetic resonance (NMR) in solution. The [15N,1H]-TROSY correlation spectra confirmed that the complex of [u-15N, ~70% 2H]-A2AAR with an inverse agonist adopts its global fold in lipid nanodiscs in solution at physiological temperature. The global assessment led to two observations of practical interest. First, A2AAR in nanodiscs can be stored for at least one month at 4 °C in an aqueous solvent. Second, LMNG/CHS micelles are a very close mimic of the environment of A2AAR in nanodiscs. The NMR signal of five individually assigned tryptophan indole 15N-1H moieties located in different regions of the receptor structure further enabled a detailed assessment of the impact of nanodiscs and LMNG/CHS micelles on the local structure and dynamics of A2AAR. As expected, the largest effects were observed near the lipid-water interface along the intra- and extracellular surfaces, indicating possible roles of tryptophan side chains in stabilizing GPCRs in lipid bilayer membranes.

Biophysical characterization of folded state type II luciferase‐like monooxygenase

We noticed that the Priestia megaterium genome contains five Luciferase‐like monooxygenase (LLM) encoding genes, however, their functions are unknown. The objective of this work was to characterize the biophysical properties of the recombinant LLM2 from Priestia megaterium PSA10 through in vitro and in silico approaches. We successfully cloned into the pET vector system and expressed the recombinant LLM2 in Escherichia coli BL21(DE3). The recombinant LLM2 was overproduced and purified in the form of an inclusion body with a molecular weight of ±39.5 kDa when it was analyzed in 15% SDS‐PAGE. The inclusion body of recombinant LLM2 was then refolded and characterized for its biophysical properties by measuring the UV spectrum of 200 to 250 nm wavelength and determining the change of enthalpy (ΔH) and entropy (ΔS) at the melting temperature. The refolded recombinant LLM2 exhibited a strong spectrum at 205 nm, while the unfolded recombinant LLM2 did not. The Tm, ΔHTm, and ΔSTm values of the refolded recombinant LLM2 were determined to be 318.31±4.4 K, 11.76±1.3 kJ.mol‐1, and (3.74±0.48)x10‐2 kJ.mol‐1.K‐1, respectively. The predicted 3D structure of LLM2 showed that the protein contains the TIM‐barrel, resembling the common global fold of bacterial luciferases. Determination of the cofactor preference suggested that the LLM2 preferred FAD for its cofactor.

Cryo-EM reveals the structure and dynamics of a 723-residue malate synthase G

Determination of sub-100 kDa (kDa) structures by cryo-electron microscopy (EM) is a longstanding but not straightforward goal. Here, we present a 2.9-Å cryo-EM structure of a 723-amino acid apo-form malate synthase G (MSG) from Escherichia coli. The cryo-EM structure of the 82-kDa MSG exhibits the same global folding as structures resolved by crystallography and nuclear magnetic resonance (NMR) spectroscopy, and the crystal and cryo-EM structures are indistinguishable. Analyses of MSG dynamics reveal consistent conformational flexibilities among the three experimental approaches, most notably that the α/β domain exhibits structural heterogeneity. We observed that sidechains of F453, L454, M629, and E630 residues involved in hosting the cofactor acetyl-CoA and substrate rotate differently between the cryo-EM apo-form and complex crystal structures. Our work demonstrates that the cryo-EM technique can be used to determine structures and conformational heterogeneity of sub-100 kDa biomolecules to a quality as high as that obtained from X-ray crystallography and NMR spectroscopy.

De novo protein fold design through sequence-independent fragment assembly simulations

De novo protein design generally consists of two steps, including structure and sequence design. Many protein design studies have focused on sequence design with scaffolds adapted from native structures in the PDB, which renders novel areas of protein structure and function space unexplored. We developed FoldDesign to create novel protein folds from specific secondary structure (SS) assignments through sequence-independent replica-exchange Monte Carlo (REMC) simulations. The method was tested on 354 non-redundant topologies, where FoldDesign consistently created stable structural folds, while recapitulating on average 87.7% of the SS elements. Meanwhile, the FoldDesign scaffolds had well-formed structures with buried residues and solvent-exposed areas closely matching their native counterparts. Despite the high fidelity to the input SS restraints and local structural characteristics of native proteins, a large portion of the designed scaffolds possessed global folds completely different from natural proteins in the PDB, highlighting the ability of FoldDesign to explore novel areas of protein fold space. Detailed data analyses revealed that the major contributions to the successful structure design lay in the optimal energy force field, which contains a balanced set of SS packing terms, and REMC simulations, which were coupled with multiple auxiliary movements to efficiently search the conformational space. Additionally, the ability to recognize and assemble uncommon super-SS geometries, rather than the unique arrangement of common SS motifs, was the key to generating novel folds. These results demonstrate a strong potential to explore both structural and functional spaces through computational design simulations that natural proteins have not reached through evolution.

Global and local cortical folding alterations are associated with neurodevelopmental subtype in bipolar disorders: a sulcal pits analysis

BackgroundAnalyzing cortical folding may provide insight into the biological underpinnings of neurodevelopmental diseases. A neurodevelopmental subtype of bipolar disorders (BD-ND) has been characterized by the combination of early age of onset and psychotic features. We investigate potential cortical morphology differences associated with this subtype. We analyze, for the first time in bipolar disorders, the sulcal pits, the deepest points in each fold of the cerebral cortex. MethodsWe extracted the sulcal pits from anatomical MRI among 512 participants gathered from 7 scanning sites. We compared the number of sulcal pits in each hemisphere as well as their regional occurrence and depth between the BD-ND subgroup (N = 184), a subgroup without neurodevelopmental features (BD, N = 77) and a group of healthy controls (HC, N = 251). ResultsIn whole brain analysis, BD-ND group have a higher number of sulcal pits in comparison to the BD group. The local analysis revealed, after correction for multiple testing, a higher occurrence of sulcal pits in the left premotor cortex among the BD-ND subgroup compared to the BD and the HC groups. ConclusionOur findings confirm that BD-ND is associated with a specific brain morphology revealed by the analysis of sulcal pits. These markers may help to better understand neurodevelopment in mood disorder and stratify patients according to a pathophysiological hypothesis.

GB1 hairpin kinetics: capturing the folding pathway with molecular dynamics, replica exchange and optimal dimensionality reduction

We have performed molecular dynamics (MD) and replica-exchange (REMD) simulations of folding of the GB1 hairpin peptide in aqueous solution. REMD results were consistent with a cooperative zipper folding model. 120 MD trajectories at 320 K yielded relaxation times of 1.8 and 100 ns, with the slower assigned to global folding. The MD folding/unfolding transitions also followed the cooperative zipper model, specifying nucleation at the central turn followed by consecutive hydrogen bond formation. Formation of hydrogen bonds and hydrophobic contacts were highly correlated. Coarse-grained kinetic models constructed with the Optimal Dimensionality Reduction (ODR) approach found a folding time of 3.3 and unfolding time of 4.0 Additionally, relaxation times in the 130–170 ns range could be assigned to formation of the transition state and off-path intermediates. The unfolded state was the most highly populated and, significantly, most heterogenous, assembling the largest number of microstates, primarily composed of extended and turn structures. The folded state was also heterogenous, but a to a lesser degree, involving the fully folded and partially folded in-register hairpins at early stages of the zipper pathway. The transition state corresponded to the nucleated hairpin, with central turn and first beta-sheet hydrogen bond, while the off-path intermediates were off-register partial hairpins. Our simulation results were in excellent agreement with experimental data on folded fraction, relaxation time and folding mechanism. The new findings from this work suggest a highly cooperative zipper folding mechanism, nascent hairpin transition state and ∼100 ns relaxation related to intermediate formation. Communicated by Ramaswamy H. Sarma