Primary folding dynamics of sperm whale apomyoglobin: core formation

The structure, thermodynamics, and kinetics of heat-induced unfolding of sperm whale apomyoglobin core formation have been studied. The most rudimentary core is formed at pH(*) 3.0 and up to 60 mM NaCl. Steady state for ultraviolet circular dichroism and fluorescence melting studies indicate that the core in this acid-destabilized state consists of a heterogeneous composition of structures of approximately 26 residues, two-thirds of the number involved for horse heart apomyoglobin under these conditions. Fluorescence temperature-jump relaxation studies show that there is only one process involved in Trp burial. This occurs in 20 micro s for a 7 degrees jump to 52 degrees C, which is close to the limits placed by diffusion on folding reactions. However, infrared temperature jump studies monitoring native helix burial are biexponential with times of 5 micro s and 56 micro s for a similar temperature jump. Both fluorescence and infrared fast phases are energetically favorable but the slow infrared absorbance phase is highly temperature-dependent, indicating a substantial enthalpic barrier for this process. The kinetics are best understood by a multiple-pathway kinetics model. The rapid phases likely represent direct burial of one or both of the Trp residues and parts of the G- and H-helices. We attribute the slow phase to burial and subsequent rearrangement of a misformed core or to a collapse having a high energy barrier wherein both Trps are solvent-exposed.     

The folding process of apomyoglobin

Apomyoglobin (apoMb) folds through at least two partially folded forms that are detected both as transient intermediates during folding/unfolding kinetics or as stable intermediates at equilibrium. Here, I summarize the results of recent kinetic studies, which combined with detailed characterizations of equilibrium forms of the protein, provide a very detailed picture of apoMb folding process. The data are consistent with a linear U<->Ia<->Ib<->N model where compaction and structure are progressively acquired.     

Molecular hinges in protein folding: the urea-denatured state of apomyoglobin

Unfolded apomyoglobin in 8 M urea at pH 2.3 displays distinct regions with different backbone mobility, as monitored by NMR relaxation. These variations in backbone mobility can be correlated with intrinsic properties of the amino acids in the sequence. Clusters of small amino acids such as glycine and alanine show increased backbone mobility compared to the average. In contrast, local hydrophobic interactions that persist in urea denaturant cause some restriction of backbone motions on a picosecond to nanosecond time scale. The model derived from the behavior of apoMb in urea depends only on the most fundamental properties of the local amino acid sequence, and thus provides a feasible paradigm for the initiation of folding.     

Morphogenesis of a protein: folding pathways and the energy landscape

Current knowledge on the reaction whereby a protein acquires its native three-dimensional structure was obtained by and large through characterization of the folding mechanism of simple systems. Given the multiplicity of amino acid sequences and unique folds, it is not so easy, however, to draw general rules by comparing folding pathways of different proteins. In fact, quantitative comparison may be jeopardized not only because of the vast repertoire of sequences but also in view of a multiplicity of structures of the native and denatured states. We have tackled the problem of the relationships between the sequence information and the folding pathway of a protein, using a combination of kinetics, protein engineering and computational methods, applied to relatively simple systems. Our strategy has been to investigate the folding mechanism determinants using two complementary approaches, i.e. (i) the study of members of the same family characterized by a common fold, but substantial differences in amino acid sequence, or (ii) heteromorphic pairs characterized by largely identical sequences but with different folds. We discuss some recent data on protein-folding mechanisms by presenting experiments on different members of the PDZ domain family and their circularly permuted variants. Characterization of the energetics and structures of intermediates and TSs (transition states), obtained by Φ-value analysis and restrained MD (molecular dynamics) simulations, provides a glimpse of the malleability of the dynamic states and of the role of the topology of the native states and of the denatured states in dictating folding and misfolding pathways.     

Fast protein folding on downhill energy landscape

Proteins fold in a time range of microseconds to minutes despite the large amount of possible conformers. Molecular dynamics simulations of a three-stranded antiparallel beta-sheet peptide (for a total of 12.6 microsec and 72 folding events) show that at the melting temperature the unfolded state ensemble contains many more conformers than those sampled during a folding event.     

岷江上游森林景观边界效应

应用地理信息系统和遥感技术对岷江上游森林景观在景观尺度上的边界效应进行了分析.结果表明:人工形成的林-农边界较清晰、其植被过渡明显,林地边缘的生物量低于林地内部,影响域为60m,农田边缘的生物量高于农田内部,边界效应范围为60～90m;自然形成的林-草边界过渡缓和,林地边缘的生物量低于林地内部,而草地边缘的生物量高于草地内部,边界效应对林地的影响范围为60m、对草地的影响范围为45～75m;林-灌边界的边界效应类似林-草边界,对林地影响域为60m,对灌木林地的影响域在45～75m之间.     

A new folding intermediate of apomyoglobin from Aplysia limacina: stepwise formation of a molten globule

Apomyoglobin from Aplysia limacina (al-apoMb), despite having only 20 % sequence identity with the more commonly studied mammalian globins (m-apoMbs), properties which result in an increased number of hydrophobic contacts and a loss of most internal salt bridges, shares a number of features of their folding profiles. We show here that it contains an unusually stable core which resists unfolding even at 70 degrees C. The equilibrium intermediate (I(T)) at this high temperature is distinct from the acid unfolded state I(A) which has many properties in common with the acid intermediate observed for the mammalian apoproteins (I(AGH)). It contains a smaller amount of secondary structure (27 % alpha-helical instead of 35 %) and is more highly solvated as evidenced from its fluorescence spectrum (lambda(max)=344 nm instead of 338 nm). Its stability is greatly increased (DeltaDeltaG(w)=-6.75 kcal mol(-1)) in the presence of high salt (2 M KCl), lending support to the view that hydrophobic interactions are responsible for its stability. Kinetic data show classical two-state kinetics between I(A) and the folded state both in the presence and absence of salt. Both I(A) and I(T) can be populated within the dead time of the stopped-flow apparatus, since initiating the refolding reaction from I(T) or I(A) rather than the completely unfolded state does not affect the observed refolding time-course. Our conclusion is that al-apoMb, as other "apo" proteins (including for example alpha-lactalbumin in the absence of Ca(2+)), may be described as "uncoupled" with an unusually high and exploitable tendency to populate partially folded states.     

The early folding kinetics of apomyoglobin.

The folding pathway of apomyoglobin has been experimentally shown to have early kinetic intermediates involving the A, B, G, and H helices. The earliest detected kinetic events occur on a ns to micros time scale. We show that the early folding kinetics of apomyoglobin may be understood as the association of nascent helices through a network of diffusion-collision-coalescence steps G + H <--> GH + A <--> AGH + B <--> ABGH obtained by solving the diffusion-collision model in a chemical kinetics approximation. Our reproduction of the experimental results indicates that the model is a useful way to analyze folding data. One prediction from our fit is that the nascent A and H helices should be relatively more helix-like before coalescence than the other apomyoglobin helices.     

Reshaping the folding energy landscape by chloride salt: impact on molten-globule formation and aggregation behavior of carbonic anhydrase

During chemical denaturation different intermediate states are populated or suppressed due to the nature of the denaturant used. Chemical denaturation by guanidine-HCl (GuHCl) of human carbonic anhydrase II (HCA II) leads to a three-state unfolding process (Cm,NI=1.0 and Cm,IU=1.9 M GuHCl) with formation of an equilibrium molten-globule intermediate that is stable at moderate concentrations of the denaturant (1-2 M) with a maximum at 1.5 M GuHCl. On the contrary, urea denaturation gives rise to an apparent two-state unfolding transition (Cm=4.4 M urea). However, 8-anilino-1-naphthalene sulfonate (ANS) binding and decreased refolding capacity revealed the presence of the molten globule in the middle of the unfolding transition zone, although to a lesser extent than in GuHCl. Cross-linking studies showed the formation of moderate oligomer sized (300 kDa) and large soluble aggregates (>1000 kDa). Inclusion of 1.5 M NaCl to the urea denaturant to mimic the ionic character of GuHCl leads to a three-state unfolding behavior (Cm,NI=3.0 and Cm,IU=6.4 M urea) with a significantly stabilized molten-globule intermediate by the chloride salt. Comparisons between NaCl and LiCl of the impact on the stability of the various states of HCA II in urea showed that the effects followed what could be expected from the Hofmeister series, where Li+ is a chaotropic ion leading to decreased stability of the native state. Salt addition to the completely urea unfolded HCA II also led to an aggregation prone unfolded state, that has not been observed before for carbonic anhydrase. Refolding from this state only provided low recoveries of native enzyme.     

Denatured collapsed states in protein folding: example of apomyoglobin.

Experimental approaches, including circular dichroism, small angle X-ray scattering, steady-state fluorescence, and fluorescence energy transfer, were applied to study the 3D-structure of apomyolgobin in different conformational states. These included the native and molten globules, along with either less ordered conformations induced by the addition of anions or completely unfolded states. The results show that the partially folded forms of apomyoglobin stabilized by KCl and/or Na(2)SO(4) under unfolding conditions (pH 2) exhibit a significant amount of secondary structure (circular dichroism), low packing density of protein molecules (SAXS), and native-like dimensions of the AGH core (fluorescence energy transfer). This finding indicates that a native-like tertiary fold of the polypeptide chain, i.e., the spatial organization of secondary structure elements, most likely emerges prior to the formation of the molten globule state.     

湘江流域中上游景观格局及其变化

以湘江流域中上游地区1980和2000年景观图为基础,运用ArcGIS9.0与Fragstats,对湘江流域中上游20年间的景观格局及其变化进行了分析。结果表明:20年间,受人类活动影响,水田与林地景观面积减少,斑块数目增加,破碎度增大,旱地和草地的斑块数减少,破碎度下降,研究区景观整体破碎度降低;优势斑块的连通性下降,景观优势类别对景观整体的控制作用减弱;分析研究区不同高程分带的景观格局变化发现,在上游山地区,有较大面积的林地被开垦为旱地,从而加大了该区水土流失的风险。     

Predicting repeat protein folding kinetics from an experimentally determined folding energy landscape

The Notch ankyrin domain is a repeat protein whose folding has been characterized through equilibrium and kinetic measurements. In previous work, equilibrium folding free energies of truncated constructs were used to generate an experimentally determined folding energy landscape (Mello and Barrick, Proc Natl Acad Sci USA 2004;101:14102–14107). Here, this folding energy landscape is used to parameterize a kinetic model in which local transition probabilities between partly folded states are based on energy values from the landscape. The landscape‐based model correctly predicts highly diverse experimentally determined folding kinetics of the Notch ankyrin domain and sequence variants. These predictions include monophasic folding and biphasic unfolding, curvature in the unfolding limb of the chevron plot, population of a transient unfolding intermediate, relative folding rates of 19 variants spanning three orders of magnitude, and a change in the folding pathway that results from C‐terminal stabilization. These findings indicate that the folding pathway(s) of the Notch ankyrin domain are thermodynamically selected: the primary determinants of kinetic behavior can be simply deduced from the local stability of individual repeats.     

Effects of confinement in chaperonin assisted protein folding: rate enhancement by decreasing the roughness of the folding energy landscape

Chaperonins, such as the GroE complex of the bacteria Escherichia coli, assist the folding of proteins under non-permissive folding conditions by providing a cavity in which the newly translated or translocated protein can be encapsulated. Whether the chaperonin cage plays a passive role in protecting the protein from aggregation, or an active role in accelerating folding rates, remains a matter of debate. Here, we investigate the role of confinement in chaperonin mediated folding through molecular dynamics simulations. We designed a substrate protein with an alpha/beta sandwich fold, a common structural motif found in GroE substrate proteins and confined it to a spherical hydrophilic cage which mimicked the interior of the GroEL/ES cavity. The thermodynamics and kinetics of folding were studied over a wide range of temperature and cage radii. Confinement was seen to significantly raise the collapse temperature, T(c), as a result of the associated entropy loss of the unfolded state. The folding temperature, T(f), on the other hand, remained unaffected by encapsulation, a consequence of the folding mechanism of this protein that involves an initial collapse to a compact misfolded state prior to rearranging to the native state. Folding rates were observed to be either accelerated or retarded compared to bulk folding rates, depending on the temperature of the simulation. Rate enhancements due to confinement were observed only at temperatures above the temperature T(m), which corresponds to the temperature at which the protein folds fastest. For this protein, T(m) lies above the folding temperature, T(f), implying that encapsulation alone will not lead to a rate enhancement under conditions where the native state is stable (T相似文献   

岷江上游景观边界网络格局分析

格局与过程是景观生态学研究的核心内容,与基于斑块的研究相比,基于边界的研究是格局与过程研究的一个新的切入点.以岷江上游地区为例,基于TM影像数据,应用RS、GIS、Fragstats等软件,选用边界长度、边界密度、边界海拔、结点数量、网眼大小、网络连接度等指数,研究了景观边界网络格局的变化,并分析了与之相关的生态过程.结果表明,在1974～2000年期间,由于人为干扰强度的加剧,岷江上游景观边界网络结构变得更加复杂,在早期以边界长度增加为主,网络连接度变大,在后期以边界数量和结点数量增加为主,网络连接度变小;森林景观与低坡位景观类型间的边界减少,森林下线上移;农田与林地的边界在早期增加,1986年后减少,而与灌木林地的边界持续增加;随着林地的蚀退,森林结构变得简单化,而农田、灌木林地、草地等景观类型的结构变得更加复杂.     

The role of protein homochirality in shaping the energy landscape of folding

The homochirality, or isotacticity, of the natural amino acids facilitates the formation of regular secondary structures such as alpha-helices and beta-sheets. However, many examples exist in nature where novel polypeptide topologies use both l- and d-amino acids. In this study, we explore how stereochemistry of the polypeptide backbone influences basic properties such as compactness and the size of fold space by simulating both lattice and all-atom polypeptide chains. We formulate a rectangular lattice chain model in both two and three dimensions, where monomers are chiral, having the effect of restricting local conformation. Syndiotactic chains with alternating chirality of adjacent monomers have a very large ensemble of accessible conformations characterized predominantly by extended structures. Isotactic chains on the other hand, have far fewer possible conformations and a significant fraction of these are compact. Syndiotactic chains are often unable to access maximally compact states available to their isotactic counterparts of the same length. Similar features are observed in all-atom models of isotactic versus syndiotactic polyalanine. Our results suggest that protein isotacticity has evolved to increase the enthalpy of chain collapse by facilitating compact helical states and to reduce the entropic cost of folding by restricting the size of the unfolded ensemble of competing states.     

Investigation of folding/unfolding kinetics of apomyoglobin

Apomyoglobin kinetic and equilibrium unfolding and folding processes were studied at pH 6.2, 11 degrees C by stopped-flow tryptophan fluorescence. There are two distinct consecutive processes in apomyoglobin folding process, namely, the protein fast transition between the unfolded (U) and an intermediate (I) states (U <----> I) and slow transition between the intermediate and the native (N) states (I <----> N). Accumulation of the intermediate state was observed in the wide range of urea concentrations. The presence of the intermediate state was shown even beyond the middle transition on the unfolding limb. The dependence of observed folding/unfolding rates on urea concentration (chevron plot) was obtained. The shape of this dependence was compared with that of two-state proteins, folding from the U to N state.     

Probing the energy landscape of protein folding/unfolding transition states

Previous molecular dynamics (MD) simulations of the thermal denaturation of chymotrypsin inhibitor 2 (CI2) have provided atomic-resolution models of the transition state ensemble that is well supported by experimental studies. Here, we use simulations to further investigate the energy landscape around the transition state region. Nine structures within approximately 35 ps and 3 A C(alpha) RMSD of the transition state ensemble identified in a previous 498 K thermal denaturation simulation were quenched under the quasi-native conditions of 335 K and neutral pH. All of the structures underwent hydrophobically driven collapse in response to the drop in temperature. Structures less denatured than the transition state became structurally more native-like, while structures that were more denatured than the transition state tended to show additional loss of native structure. The structures in the immediate region of the transition state fluctuated between becoming more and less native-like. All of the starting structures had the same native-like topology and were quite similar (within 3.5 A C(alpha) RMSD). That the structures all shared native-like topology, yet diverged into either more or less native-like structures depending on which side of the transition state they occupied on the unfolding trajectory, indicates that topology alone does not dictate protein folding. Instead, our results suggest that a detailed interplay of packing interactions and interactions with water determine whether a partially denatured protein will become more native-like under refolding conditions.     

Outer membrane protein folding from an energy landscape perspective

The cell envelope is essential for the survival of Gram-negative bacteria. This specialised membrane is densely packed with outer membrane proteins (OMPs), which perform a variety of functions. How OMPs fold into this crowded environment remains an open question. Here, we review current knowledge about OMP folding mechanisms in vitro and discuss how the need to fold to a stable native state has shaped their folding energy landscapes. We also highlight the role of chaperones and the β-barrel assembly machinery (BAM) in assisting OMP folding in vivo and discuss proposed mechanisms by which this fascinating machinery may catalyse OMP folding.     

Free energy landscape of protein folding in water: explicit vs. implicit solvent

The Generalized Born (GB) continuum solvent model is arguably the most widely used implicit solvent model in protein folding and protein structure prediction simulations; however, it still remains an open question on how well the model behaves in these large-scale simulations. The current study uses the beta-hairpin from C-terminus of protein G as an example to explore the folding free energy landscape with various GB models, and the results are compared to the explicit solvent simulations and experiments. All free energy landscapes are obtained from extensive conformation space sampling with a highly parallel replica exchange method. Because solvation model parameters are strongly coupled with force fields, five different force field/solvation model combinations are examined and compared in this study, namely the explicit solvent model: OPLSAA/SPC model, and the implicit solvent models: OPLSAA/SGB (Surface GB), AMBER94/GBSA (GB with Solvent Accessible Surface Area), AMBER96/GBSA, and AMBER99/GBSA. Surprisingly, we find that the free energy landscapes from implicit solvent models are quite different from that of the explicit solvent model. Except for AMBER96/GBSA, all other implicit solvent models find the lowest free energy state not the native state. All implicit solvent models show erroneous salt-bridge effects between charged residues, particularly in OPLSAA/SGB model, where the overly strong salt-bridge effect results in an overweighting of a non-native structure with one hydrophobic residue F52 expelled from the hydrophobic core in order to make better salt bridges. On the other hand, both AMBER94/GBSA and AMBER99/GBSA models turn the beta-hairpin in to an alpha-helix, and the alpha-helical content is much higher than the previously reported alpha-helices in an explicit solvent simulation with AMBER94 (AMBER94/TIP3P). Only AMBER96/GBSA shows a reasonable free energy landscape with the lowest free energy structure the native one despite an erroneous salt-bridge between D47 and K50. Detailed results on free energy contour maps, lowest free energy structures, distribution of native contacts, alpha-helical content during the folding process, NOE comparison with NMR, and temperature dependences are reported and discussed for all five models.