Intramolecular cross-linking evaluated as a structural probe of the protein folding transition state

We examine the utility of intramolecular covalent cross-linking to identify the structure present in the folding transition state. In mammalian ubiquitin, cysteine residues located across two beta-strands are cross-linked with dichloroacetone. The kinetic effects of these covalent cross-links in ubiquitin, and engineered disulfide bonds in src SH3 (Grantcharova, V. P., Riddle, D. S., and Baker, D. (2000) Proc. Natl. Acad. Sci. U.S.A. 97, 7084-7089), are compared to the results of psi-analysis where strand association is stabilized by metal ion binding to engineered bihistidine sites (Krantz, B. A., Dothager, R. S., and Sosnick, T. R. (2004) J. Mol. Biol. 337, 463-75) at the same positions. The results for the two methods agree at some of the sites. The cross-linking phi crosslink-values agree with their corresponding psi-values when they have both have values of zero or one, which represent the absence and presence of native structure, respectively. When phi crosslink > psi, the apparent inconsistency is rationalized by the difference between each method's mode of stabilization; cross-linking reduces the configurational entropy of the unfolded state whereas metal binding directly stabilizes the native state. However, when the cross-linking phi-values are smaller than their corresponding psi-values, the apparent underestimation of structure formation is difficult to rationalize while retaining the assumption that the cross-link exclusively affects the entropy of the unfolded state. The interpretation also is problematic for data on cross-links located across strands which are not hairpins, and hence, these sites are likely to be of limited utility in folding studies. We conclude that cross-linking data for sites on hairpins generally report on the amount of structure formed within the enclosed loop while the metal binding data report on the amount structure formed at the site itself.     

Analysis of the high- and low-spin Soret bands of horse-heart metmyoglobin complexes

Interest in the thermodynamics of the iron-binding site in hemoproteins has increased in recent years due to refinements in x-ray crystallographic studies of hemoproteins [see Deathage, J. F., Lee, R. S., Anderson, C. M. & Moffat, K. (1976) J. Mol. Biol. 104 , 687–706; Heidner, E. J., Ladner, R. C. & Perutz, M. F. (1976) J. Mol. Biol. 104 , 707–722; Deathage, J. F., Lee, R. S. & Moffat, K. (1976) J. Mol. Biol. 104 , 723–728; Ladner, R. C., Heidner, E. J. & Perutz, M. F. (1976) J. Mol. Biol. 114 , 385–414; Fermi, G. & Perutz, M. F. (1977) J. Mol. Biol. 114 , 421–431; Takano, T. (1977) J. Mol. Biol. 110 , 537–568 and 569–589], the synthesis and x-ray analysis of model heme compounds [see Scheidt, W. R. (1977) Acc. Chem. Res. 10 , 339–345; Kastner, M. E., Scheidt, W. R., Mashino, T. & Reed, C. A. (1978) J. Am. Chem. Soc. 100 , 666–667; Mashiko, T., Kastner, M. E., Spartalian, K., Scheidt, W. R. & Reed, C. A. (1978) J. Am. Chem. Soc. 100 , 6354–6362; Hill, H. A. O., Skite, P. P., Buchler, J. W., Luchr, H., Tonn, M., Gregson, A. K. & Pellizer, G. (1979) Chem. Commun. 4 , 151–152; and Scheidt, W. R., Cohen, I. A. & Kastner, M. E. (1979) Biochemistry 18 , 3546–3556], and the numerous data on heme–protein interactions that account for the differences observed in ligand binding between the various species of animals. Numerous probes have been used and provide information about the structure and thermodynamics of the binding site, but no single probe can provide the complete picture [see Iizuka, T. & Yonetani, T. (1970) Adv. Biophys. 1 , 157–182; Smith, D. W. & Williams, R. J. P. (1970) Struct. Bond. 7 , 1–45; and Spiro, T. G. (1975) Biochim. Biophys. Acta 416 , 169–189].     

Structural transition of alpha 1-antitrypsin by a peptide sequentially similar to beta-strand s4A

Crystal structure studies have shown that cleaved and intact serpins differ essentially in the topology of beta-sheet A. This is five-stranded in the intact molecules and six-stranded after cleavage by insertion of strand s4A whose C-terminus has become free [L?bermann, H., Tokuoka, R., Deisenhofer, J. & Huber, R. (1984) J. Mol. Biol. 177, 531-556; Wright, T. H., Qian, H. X. & Huber, R. (1990) J. Mol. Biol. 213, 513-528]. The structural transition is accompanied by changes in spectral properties and an increase in thermal stability. We show here that an N alpha-acetyl-tetradecapeptide with the amino acid sequence of strand s4A, residues 345-358 of human alpha 1-antitrypsin, associates with intact alpha 1-antitrypsin and forms a stoichiometric complex with properties very similar to cleaved alpha 1-antitrypsin. Complex generation has the characteristics of a folding process.     

Extending the folding nucleus of ubiquitin with an independently folding beta-hairpin finger: hurdles to rapid folding arising from the stabilisation of local interactions

The N-terminal beta-hairpin sequence of ubiquitin has been implicated as a folding nucleation site. To extend and stabilise the ubiquitin folding nucleus, we have inserted an autonomously folding 14-residue peptide sequence beta4 which in isolation forms a highly populated beta-hairpin (>70%) stabilised by local interactions. NMR structural analysis of the ubiquitin mutant (Ubeta4) shows that the hairpin finger is fully structured and stabilises ubiquitin by approximately 8kJmol(-1). Protein engineering and kinetic (phi(F)-value) analysis of a series of Ubeta4 mutants shows that the hairpin extension of Ubeta4 is also significantly populated in the transition state (phi(F)-values >0.7) and has the effect of templating the formation of native contacts in the folding nucleus of ubiquitin. However, at low denaturant concentrations the chevron plot of Ubeta4 shows a small deviation from linearity (roll-over effect), indicative of the population of a compact collapsed state, which appears to arise from over-stabilisation of local interactions. Destabilising mutations within the native hairpin sequence and within the engineered hairpin extension, but not elsewhere, eliminate this non-linearity and restore apparent two-state behaviour. The pitfall to stabilising local interactions is to present hurdles to the rapid and efficient folding of small proteins down a smooth folding funnel by trapping partially folded or misfolded states that must unfold or rearrange before refolding.     

Folding of homologous proteins: conservation of the folding mechanism of the alpha subunit of tryptophan synthase from Escherichia coli, Salmonella typhimurium, and five interspecies hybrids

The equilibrium and kinetic properties for the urea-induced unfolding of the alpha subunit of tryptophan synthase from Escherichia coli, Salmonella typhimurium, and five interspecies hybrids were compared to determine the role of protein folding in evolution. The parent proteins differ at 40 positions in the sequence of 268 amino acids, and the hybrids differ by up to 15 amino acids from the Escherichia coli alpha subunit. The results show that all the proteins follow the same folding mechanism and are consistent with a previously proposed hypothesis [Hollecker, M., & Creighton, T. E. (1983) J. Mol. Biol. 168, 409; Krebs, H., Schmid, F. X., & Jaenicke, R. (1983) J. Mol. Biol. 169, 619] that the folding mechanisms are conserved in homologous proteins. Analysis of the kinetic data suggests that the 15 positions at which the parent proteins differ in the amino folding unit, residues 1-188, do not play a role in a rate-limiting step in folding that has been previously identified as the association of the amino and carboxyl folding units [Beasty, A. M., Hurle, M. R., Manz, J. T., Stackhouse, T. S., Onuffer, J. J., & Matthews, C. R. (1986) Biochemistry 25, 2965]. One or more of the 25 positions at which the parent proteins differ in the carboxyl folding unit, residues 189-268, do appear to play a role in this same rate-limiting step.     

raf RBD and ubiquitin proteins share similar folds, folding rates and mechanisms despite having unrelated amino acid sequences

Recent experimental and theoretical studies in protein folding suggest that the rates and underlying mechanisms by which proteins attain the native state are largely determined by the topological complexity of a specific fold rather than by the fine details of the amino acid sequences. However, such arguments are based upon the examination of a limited number of protein folds. To test this view, we sought to investigate whether proteins belonging to the ubiquitin superfamily display similar folding behavior. To do so, we compared the folding-unfolding transitions of mammalian ubiquitin (mUbi) with those of its close yeast homologue (yUbi), and to those of the structurally related Ras binding domain (RBD) of the serine/threonine kinase raf that displays no apparent sequence homology with the ubiquitin family members. As demonstrated for mUbi [Krantz, B. A., and Sosnick, T. R. (2000) Biochemistry 39, 11696-11701], we show that a two-state transition model with no burst phase intermediate can describe folding of both yUbi and raf RBD. We further demonstrate that (1) all three proteins refold at rates that are within 1 order of magnitude (1800, 1100, and 370 s(-1) for mUbi, raf RBD, and yUbi, respectively), (2) both mUbi and raf RBD display similar refolding heterogeneity, and (3) the folding free energy barriers of both mUbi and raf RBD display a similar temperature dependence and sensitivity to a stabilizing agent or to mutations of a structurally equivalent central core residue. These findings are consistent with the view that rates and mechanisms for protein folding depend mostly on the complexity of the native structure topology rather than on the fine details of the amino acid sequence.     

Stability and folding kinetics of a ubiquitin mutant with a strong propensity for nonnative beta-hairpin conformation in the unfolded state

A F45W mutant of yeast ubiquitin has been used as a model system to examine the effects of nonnative local interactions on protein folding and stability. Mutating the native TLTGK G-bulged type I turn in the N-terminal beta-hairpin to NPDG stabilizes a nonnative beta-strand alignment in the isolated peptide fragment. However, NMR structural analysis of the native and mutant proteins shows that the NPDG mutant is forced to adopt the native beta-strand alignment and an unfavorable type I NPDG turn. The mutant is significantly less stable (approximately 9 kJ mol(-1)) and folds 30 times slower than the native sequence, demonstrating that local interactions can modulate protein stability and that attainment of a nativelike beta-hairpin conformation in the transition state ensemble is frustrated by the turn mutations. Surprising, alcoholic cosolvents [5-10% (v/v) TFE] are shown to accelerate the folding rate of the NPDG mutant. We conclude, backed-up by NMR data on the peptide fragments, that even though nonnative states in the denatured ensemble are highly populated and their stability further enhanced in the presence of cosolvents, the simultaneous increase in the proportion of nativelike secondary structure (hairpin or helix), in rapid equilibrium with nonnative states, is sufficient to accelerate the folding process. It is evident that modulating local interactions and increasing nonnative secondary structure propensities can change protein stability and folding kinetics. However, nonlocal contacts formed in the global cooperative folding event appear to determine structural specificity.     

Engineering stabilising beta-sheet interactions into a conformationally flexible region of the folding transition state of ubiquitin

Protein engineering studies suggest that the transition state for the folding of ubiquitin is highly polarised towards the N-terminal part of the sequence and involves a nucleus of residues within the beta-hairpin (residues 1-17) and main alpha-helix (residues 23-34). In contrast, the observation of small phi-values for residues in the C-terminal portion of the sequence (residues 35-76), coupled with a folding topology that results in a much higher contact order, suggests that fast folding of ubiquitin is dependent upon configurational flexibility in the C-terminal part of the polypeptide chain to ensure passage down a relatively smooth folding funnel to the native state. We show that the introduction of a small mini-hairpin motif as an extension of the native 43-50 hairpin stabilises local interactions in the C-terminal part of the sequence, resulting largely in a deceleration of the unfolding kinetics without perturbing the apparent two-state folding mechanism. However, a single-point Leu-->Phe substitution within the engineered hairpin sequence leads to the premature collapse of the denatured ensemble through the stabilisation of non-native interactions and the population of a compact intermediate. Non-linear effects in the kinetic data at low concentrations of denaturant suggest that the collapsed state, which is further stabilised in the presence of cosmotropic salts, may subsequently fold directly to the native state through a "triangular" reaction scheme involving internal rearrangement rather than unfolding and refolding.     

Trehalose-6-phosphate synthetase activity in extracts of Dictyostelium discoideum.

Trehalose-6-phosphate (T-6-P) synthetase activity in extracts of Dictyostelium discoideum has been reexamined in an effort to resolve discrepancies between the results of previous studies (R. Roth and M. Sussman (1966). Biochim. Biophys. Acta, 122, 225; K. A. Killick and B. E. Wright (1972). J. Biol. Chem., 247, 2967). We find that T-6-P synthetase is not cold sensitive as reported by Killick and Wright (1972), is not present in bacterial-grown vegetative cells (though subject to some modulation by other nutritional conditions), and is not in our hands unmasked or activated by ammonium sulfate fractionation. We conclude that the pattern of T-6-P synthetase accumulation and disappearance during fruiting body construction in D. discoideum is as originally described by R. Roth and M. Sussman (1968). J. Biol. Chem., 243, 5081) and confirmed elsewhere (P. C. Newell et al. (1972). J. Mol. Biol., 63, 373; R. W. Brackenbury et al. (1974). J. Mol. Biol., 90, 529; B. D. Hames and J. M. Ashworth (1974). Biochem. J., 142, 301).     

Fluorescent lipoprotein probe

The fluorescent cholesterol analog cholesta-5,7,9(11)triene-3-β-ol was used to label high-density and low-density lipoproteins in vivo (rabbit) and in vitro (human). Rabbit feeding experiments demonstrated that in vivo both the esterified and nonesterified forms of the fluorophore were incorporated by these particles. Using in vitro labeling techniques, it was possible to selectively incorporate either the free form of the fluorophore or both the free and the esterified forms depending upon the presence or absence of serum esterases during incubation. Subsequent to labeling, the thermotropic behavior of the low- and high-density lipoproteins was evaluated using temperature-dependent fluorescence intensity measurements. Purified low-density lipoprotein samples (human and rabbit) containing both forms of the fluorophore were observed to undergo a thermotropic transition between 27 and 32°C. However, this transition was not observed in low-density lipoprotein samples containing only the nonesterified form of the probe nor was it observed in any of the high-density lipoprotein samples, even those containing both forms of the fluorophore. These results provide further evidence that the previously reported thermotropic transition in low-density lipoprotein is due to a reordering of the low-density lipoprotein cholesterol ester core (Deckelbaum, R.J., Shipley, G. G., Small, D. M., Lees, R. S., and George, P. K. (1975) Science190, 392–394; Sears, B., Deckelbaum, R. J., Janiak, M. J., Shipley, G. G., and Small, D. M. (1976) Biochemistry15, 4151–4157; Chana, G. S., Sheppard, R. J., Mills, G. L., and Grant, E. H. (1980) Phys. Med. Biol.25, 427–432).     

Insights into the stability of native and partially folded states of ubiquitin: effects of cosolvents and denaturants on the thermodynamics of protein folding

The thermodynamics of the native<-->A state and native<-->unfolded transitions for ubiquitin have been characterized in detail using the denaturants methanol and guanidinium chloride (Gdn.HCl) both separately and in combination. Gdn.HCl destabilizes the partially folded alcohol-induced A state such that the effects of alcoholic solvents on the native<-->unfolded transition can be investigated directly via a two-state model. The combined denaturing effects of methanol and Gdn.HCl appear to conform to a simple additive model. We show that ubiquitin folds and unfolds cooperatively in all cases, forming the same "native" state; however, the thermodynamics of the N<-->U transition change dramatically between alcoholic and Gdn.HCl solutions, with folding in aqueous methanol associated with large negative enthalpy and entropy terms at 298 K with a gradual falloff in DeltaC(p) at higher methanol concentrations, as previously reported for the N<-->A transition (Woolfson, D. N., Cooper, A., Harding, M. M., Williams, D. H., and Evans, P. A. (1993) J. Mol. Biol. 229, 502-511.). Both the N<-->U and the N<-->A transitions are enthalpy driven to a similar extent. We conclude that under these conditions van der Waals interactions in the packing of the nonpolar protein core, which is common to both the N<-->U and the N<-->A transitions, appear to drive folding in the absence of entropic effects associated with release of ordered solvent (hydrophobic effect). Solvent transfer studies of hydrocarbons into alcoholic solvents, with and without Gdn.HCl, are consistent with a large enthalpic driving force for burial of a nonpolar surface, with a linear dependence of protein stability (DeltaG(N)(<-->)(U)) on cosolvent concentration reflected in a similar linear dependence of hydrocarbon solubility. The data demonstrate that the hydrophobic effect is not a prerequisite for specific stabilization of the native state or the A state and that van der Waals packing of the nonpolar core appears to be the dominant factor in stabilization of the native state.     

The folding pathway of ubiquitin from all-atom molecular dynamics simulations

The folding (unfolding) pathway of ubiquitin is probed using all-atom molecular dynamics simulations. We dissect the folding pathway using two techniques: first, we probe the folding pathway of ubiquitin by calculating the evolution of structural properties over time and second, we identify the rate determining transition state for folding. The structural properties that we look at are hydrophobic solvent accessible surface area (SASA) and Calpha-root-mean-square deviation (rmsd). These properties on their own tell us relatively little about the folding pathway of ubiquitin; however, when plotted against each other, they become powerful tools for dissecting ubiquitin's folding mechanism. Plots of Calpha-rmsd against SASA serve as a phase space trajectories for the folding of ubiquitin. In this study, these plots show that ubiquitin folds to the native state via the population of an intermediate state. This is shown by an initial hydrophobic collapse phase followed by a second phase of secondary structure arrangement. Analysis of the structure of the intermediate state shows that it is a collapsed species with very little secondary structure. In reconciling these observations with recent experimental data, the transition that we observe in our simulations from the unfolded state (U) to the intermediate state (I) most likely occurs in the dead-time of the stopped flow instrument. The folding pathway of ubiquitin is probed further by identification of the rate-determining transition state for folding. The method used for this is essential dynamics, which utilizes a principal component analysis (PCA) on the atomic fluctuations throughout the simulation. The five transition state structures identified in silico are in good agreement with the experimentally determined transition state. The calculation of phi-values from the structures generated in the simulations is also carried out and it shows a good correlation with the experimentally measured values. An initial analysis of the denatured state shows that it is compact with fluctuating regions of nonnative secondary structure. It is found that the compactness in the denatured state is due to the burial of some hydrophobic residues. We conclude by looking at a correlation between folding kinetics and residual structure in the denatured state. A hierarchical folding mechanism is then proposed for ubiquitin.     

Targeted Molecular Dynamics Simulation of Conformational Change-Application to the T ↔ R Transition in Insulin

Abstract

A novel method to calculate transition pathways between two known protein conformations is presented. It is based on a molecular dynamics simulation starting from one conformational state as initial structure and using the other for a directing constraint. The method is exemplified with the T ? R transition of insulin. The most striking difference between these conformational states is that in T the 8 N-terminal residues of the B chain are arranged as an extended strand whereas in R they are forming a helix. Both the transition from T to R and from R to T were simulated. The method proves capable of finding a continuous pathway for each direction which are moderately different. The refolding processes are illustrated by a series of transient structures and pairs of Ø, ψ angles selected from the time course of the simulations. In the T → R direction the helix is formed in the →last third of the transition, while in the R → T direction it is preserved during more than half of the simulation period. The results are discussed in comparison with those of an atternative method recently apptied to the T → R transition of insulin which is based on targeted energy minimisation.     

Contribution to the thermodynamics of protein folding from the reduction in water-accessible nonpolar surface area

Protein folding and the transfer of hydrocarbons from a dilute aqueous solution to the pure liquid phase are thermodynamically similar in that both processes remove nonpolar surface from water and both are accompanied by anomalously large negative heat capacity changes. On the basis of a limited set of published surface areas, we previously proposed that heat capacity changes (delta C degrees p) for the transfer of hydrocarbons from water to the pure liquid phase and for the folding of globular proteins exhibit the same proportionality to the reduction in water-accessible nonpolar surface area (delta Anp) [Spolar, R.S., Ha, J.H., & Record, M.T., Jr. (1989) Proc. Natl. Acad. Sci. U.S.A. 86, 8382-8385]. The consequence of this proposal is that the experimental delta C degrees p for protein folding can be used to obtain estimates of delta Anp and of the contribution to the stability of the folded state from removal of a nonpolar surface from water. In this paper, a rigorous molecular surface area algorithm [Richmond, T.J. (1984) J. Mol. Biol. 178, 63-89] is applied to obtain self-consistent values of the water-accessible nonpolar surface areas of the native and completely denatured states of the entire set of globular proteins for which both crystal structures and delta C degrees p of folding have been determined and for the set of liquid and liquefiable hydrocarbons for which delta C degrees p of transfer are known. Both processes (hydrocarbon transfer and protein folding) exhibit the same direct proportionality between delta C degrees p and delta Anp. We conclude that the large negative heat capacity changes observed in protein folding and other self-assembly processes involving proteins provide a quantitative measure of the reduction in the water-accessible nonpolar surface area and of the contribution of the hydrophobic effect to the stability of the native state and to protein assembly.     

Converging experimental and computational views of the knotting mechanism of a small knotted protein

MJ0366 from Methanocaldococcus jannaschii is the smallest topologically knotted protein known to date. 92 residues in length, MJ0366 ties a trefoil (3₁) knot by threading its C-terminal helix through a buttonhole formed by the remainder of the secondary structure elements. By generating a library of point mutations at positions pertinent to the knot formation, we systematically evaluated the contributions of individual residues to the folding stability and kinetics of MJ0366. The experimental Φ-values were used as restraints to computationally generate an ensemble of conformations that correspond to the transition state of MJ0366, which revealed several nonnative contacts. The importance of these nonnative contacts in stabilizing the transition state of MJ0366 was confirmed by a second round of mutagenesis, which also established the pivotal role of F15 in stapling the network of hydrophobic interactions around the threading C-terminal helix. Our converging experimental and computational results show that, despite the small size, the transition state of MJ0366 is formed at a very late stage of the folding reaction coordinate, following a polarized pathway. Eventually, the formation of extensive native contacts, as well as a number of nonnative ones, leads to the threading of the C-terminal helix that defines the topological knot.     

Fast-Folding Pathways of the Thrombin-Binding Aptamer G-Quadruplex Revealed by a Markov State Model

G-quadruplex structures participate in many important cellular processes. For a better understanding of their functions, knowledge of the mechanism by which they fold into the functional native structures is necessary. In this work, we studied the folding process of the thrombin-binding aptamer G-quadruplex. Enabled by a computational paradigm that couples an advanced sampling method and a Markov state model, four folding intermediates were identified, including an antiparallel G-hairpin, two G-triplex structures, and a double-hairpin conformation. Likewise, a misfolded structure with a nonnative distribution of syn/anti guanines was also observed. Based on these states, a transition path analysis revealed three fast-folding pathways, along which the thrombin-binding aptamer would fold to the native state directly, with no evidence of potential nonnative competing conformations. The results also showed that the TGT-loop plays an important role in the folding process. The findings of this research may provide general insight about the folding of other G-quadruplex structures.     

NMR evidence for forming highly populated helical conformations in the partially folded hNck2 SH3 domain

Recent studies of several proteins implied that the folding of β-proteins may follow a nonhierarchical mechanism in which two major transitions are essential, i.e., the collapse of a random coil to form a nonnative helical intermediate, followed by a transformation into the native β-structure. We report that the first hNck2 SH3 domain, assuming an all-β barrel in the native form, can be reversibly transformed into a stable and nonnative helical state by acid-unfolding. We also conducted extensive NMR and mutagenesis studies that led to two striking findings: 1), NMR analysis reveals that in the helical state formed at pH 2.0, the first and last β-strands in the native form become unstructured, whereas the rest is surprisingly converted into two highly populated helices with a significantly limited backbone motion; and 2), a conserved four-residue sequence is identified on the second β-strand, a mutation of which suddenly renders the SH3 domain into a helical state even at pH 6.5, with NMR conformational and dynamic properties highly similar to those of the wild-type at pH 2.0. This observation implies that the region might contribute key interactions to disrupt the helical state, and to facilitate a further transformation into the native SH3 fold in the second transition.     

Loop dependence of the stability and dynamics of nucleic acid hairpins

Hairpin loops are critical to the formation of nucleic acid secondary structure, and to their function. Previous studies revealed a steep dependence of single-stranded DNA (ssDNA) hairpin stability with length of the loop (L) as ~L^{8.5 ± 0.5}, in 100 mM NaCl, which was attributed to intraloop stacking interactions. In this article, the loop-size dependence of RNA hairpin stabilities and their folding/unfolding kinetics were monitored with laser temperature-jump spectroscopy. Our results suggest that similar mechanisms stabilize small ssDNA and RNA loops, and show that salt contributes significantly to the dependence of hairpin stability on loop size. In 2.5 mM MgCl₂, the stabilities of both ssDNA and RNA hairpins scale as ~L^{4 ± 0.5}, indicating that the intraloop interactions are weaker in the presence of Mg²⁺. Interestingly, the folding times for ssDNA hairpins (in 100 mM NaCl) and RNA hairpins (in 2.5 mM MgCl₂) are similar despite differences in the salt conditions and the stem sequence, and increase similarly with loop size, ~L^{2.2 ± 0.5} and ~L^{2.6 ± 0.5}, respectively. These results suggest that hairpins with small loops may be specifically stabilized by interactions of the Na⁺ ions with the loops. The results also reinforce the idea that folding times are dominated by an entropic search for the correct nucleating conformation.     

High temperature unfolding simulations of the TRPZ1 peptide

We report high temperature molecular dynamics simulations of the unfolding of the TRPZ1 peptide using an explicit model for the solvent. The system has been simulated for a total of 6 μs with 100-ns minimal continuous stretches of trajectory. The populated states along the simulations are identified by monitoring multiple observables, probing both the structure and the flexibility of the conformations. Several unfolding and refolding transition pathways are sampled and analyzed. The unfolding process of the peptide occurs in two steps because of the accumulation of a metastable on-pathway intermediate state stabilized by two native backbone hydrogen bonds assisted by nonnative hydrophobic interactions between the tryptophan side chains. Analysis of the un/folding kinetics and classical commitment probability calculations on the conformations extracted from the transition pathways show that the rate-limiting step for unfolding is the disruption of the ordered native hydrophobic packing (Trp-zip motif) leading from the native to the intermediate state. But, the speed of the folding process is mainly determined by the transition from the completely unfolded state to the intermediate and specifically by the closure of the hairpin loop driven by formation of two native backbone hydrogen bonds and hydrophobic contacts between tryptophan residues. The temperature dependence of the unfolding time provides an estimate of the unfolding activation enthalpy that is in agreement with experiments. The unfolding time extrapolated to room temperature is in agreement with the experimental data as well, thus providing a further validation to the analysis reported here.     

The complete genome sequence of Bacillus thuringiensis Al Hakam

Bacillus thuringiensis is an insect pathogen that is widely used as a biopesticide (E. Schnepf, N. Crickmore, J. Van Rie, D. Lereclus, J. Baum, J. Feitelson, D. R. Zeigler, and D. H. Dean, Microbiol. Mol. Biol. Rev. 62:775-806, 1998). Here we report the finished, annotated genome sequence of B. thuringiensis Al Hakam, which was collected in Iraq by the United Nations Special Commission (L. Radnedge, P. Agron, K. Hill, P. Jackson, L. Ticknor, P. Keim, and G. Andersen, Appl. Environ. Microbiol. 69:2755-2764, 2003).