共查询到20条相似文献，搜索用时 0 毫秒
1.
The free energy landscape and the folding mechanism of the C-terminal beta-hairpin of protein G is studied by extensive replica exchange molecular dynamics simulations (40 replicas and 340 ns total simulation time), using the GROMOS96 force field and the SPC explicit water solvent. The study reveals that the system preferentially adopts a beta-hairpin structure at biologically important temperatures, and that the helix content is low at all temperatures studied. Representing the free energy landscape as a function of several types of reaction coordinates, four local minima corresponding to the folded, partially folded, molten globule, and unfolded states are identified. The findings suggest that the folding of the beta-hairpin occurs as the sequence: collapse of hydrophobic core --> formation of H-bond --> formation of the turn. Identifying the folded and molten globule states as the main conformations, the free energy landscape of the beta-hairpin is consistent with a two-state behavior with a broad transition state. The temperature dependence of the folding-unfolding transition is investigated in some detail. The enthalpy and entropy jumps at the folding transition temperature are found to be about three times lower than the experimental estimates, indicating that the folding-unfolding transition in silico is less cooperative than its in vitro counterpart. 相似文献
2.
We have studied the potential of mean force of two peptides, one known to adopt a beta-hairpin and the other an alpha-helical conformation in solution. These peptides are, respectively, residues 41-56 of the C-terminus (GEWTYDDATKTFTVTE) of the B1 domain of protein G and the 13 residue C-peptide (KETAAAKFERQHM) of ribonuclease A. Extensive canonical ensemble sampling has been performed using a parallel replica exchange method. The effective potential employed in this work consists of the OPLS all-atom force field (OPLS-AA) and an analytical generalized Born (AGB) implicit solvent model including a novel nonpolar solvation free energy estimator (NP). An additional dielectric screening parameter has been incorporated into the AGBNP model. In the case of the beta-hairpin, the nonpolar solvation free energy estimator provides the necessary effective interactions for the collapse of the hydrophobic core (W43, Y45, F52, and V54), which the more commonly used surface-area-dependent nonpolar model does not provide. For both the beta-hairpin and the alpha-helix, increased dielectric screening reduces the stability of incorrectly formed salt bridges, which tend to disrupt the formation of the hairpin and helix, respectively. The fraction of beta-hairpin and alpha-helix content we obtained using the AGBNP model agrees well with experimental results. The thermodynamic stability of the beta-hairpin from protein G and the alpha-helical C-peptide from ribonuclease A as modeled with the OPLS-AA/AGBNP effective potential reflects the balance between the nonpolar effective potential terms, which drive compaction, and the polar and hydrogen bonding terms, which promote secondary structure formation. 相似文献
3.
We study the dependence of folding time, nucleation site, and stability of a model beta-hairpin on the location of a cross-strand hydrophobic pair, using a coarse-grained off-lattice model with the aid of Monte Carlo simulations. Our simulations have produced 6500 independent folding trajectories dynamically, forming the basis for extensive statistical analysis. Four folding pathways, zipping-out, middle-out, zipping-in, and reptation, have been closely monitored and discussed in all seven sequences studied. A hydrophobic pair placed near the beta-turn or in the middle section effectively speed up folding; a hydrophobic pair placed close to the terminal ends or next to the beta-turn encourages stability of the entire chain. 相似文献
4.
Chignolin is a 10-residue peptide (GYDPETGTWG) that forms a stable beta-hairpin structure in water. However, its design template, GPM12 (GYDDATKTFG), does not have a specific structure. To clarify which amino acids give it the ability to form the beta-hairpin structure, we calculated the folding free-energy landscapes of chignolin, GPM12, and their chimeric peptides using multicanonical molecular dynamics (MD) simulation. Cluster analysis of the conformational ensembles revealed that the native structure of chignolin was the lowest in terms of free energy while shallow local minima were widely distributed in the free energy landscape of GPM12, in agreement with experimental observations. Among the chimeric peptides, GPM12(D4P/K7G) stably formed the same beta-hairpin structure as that of chignolin in the MD simulation. This was confirmed by nuclear magnetic resonance (NMR) spectroscopy. A comparison of the free-energy landscapes showed that the conformational distribution of the Asp3-Pro4 sequence was inherently biased in a way that is advantageous both to forming hydrogen bonds with another beta-strand and to initiating loop structure. In addition, Gly7 helps stabilize the loop structure by having a left-handed alpha-helical conformation. Such a conformation is necessary to complete the loop structure, although it is not preferred by other amino acids. Our results suggest that the consistency between the short-range interactions that determine the local geometries and the long-range interactions that determine the global structure is important for stable tertiary structure formation. 相似文献
5.
The stability and structure of several beta-hairpin peptide variants derived from the C-terminus of the B1 domain of protein G were investigated by a number of experimental and computational techniques. Our analysis shows that the structure and stability of this hairpin can be greatly affected by one or a few simple mutations. For example, removing an unfavorable charge near the N-terminus of the peptide (Glu42 to Gln or Thr) or optimization of the N-terminal charge-charge interactions (Gly41 to Lys) both stabilize the peptide, even in water. Furthermore, a simple replacement of a charged residue in the turn (Asp47 to Ala) changes the beta-turn conformation. Finally, we show that the effects of combining these single mutations are additive, suggesting that independent stabilizing interactions can be isolated and evaluated in a simple model system. Our results indicate that the structure and stability of this beta-hairpin peptide can be modulated in numerous ways and thus contributes toward a more complete understanding of this important model beta-hairpin as well as to the folding and stability of larger peptides and proteins. 相似文献
6.
Most secreted bacterial proteases, including alpha-lytic protease (alphaLP), are synthesized with covalently attached pro regions necessary for their folding. The alphaLP folding landscape revealed that its pro region, a potent folding catalyst, is required to circumvent an extremely large folding free energy of activation that appears to be a consequence of its unique unfolding transition. Remarkably, the alphaLP native state is thermodynamically unstable; a large unfolding free energy barrier is solely responsible for the persistence of its native state. Although alphaLP folding is well characterized, the structural origins of its remarkable folding mechanism remain unclear. A conserved beta-hairpin in the C-terminal domain was identified as a structural element whose formation and positioning may contribute to the large folding free energy barrier. In this article, we characterize the folding of an alphaLP variant with a more favorable beta-hairpin turn conformation (alphaLP(beta-turn)). Indeed, alphaLP(beta-turn) pro region-catalyzed folding is faster than that for alphaLP. However, instead of accelerating spontaneous folding, alphaLP(beta-turn) actually unfolds more slowly than alphaLP. Our data support a model where the beta-hairpin is formed early, but its packing with a loop in the N-terminal domain happens late in the folding reaction. This tight packing at the domain interface enhances the kinetic stability of alphaLP(beta-turn), to nearly the same degree as the change between alphaLP and a faster folding homolog. However, alphaLP(beta-turn) has impaired proteolytic activity that negates the beneficial folding properties of this variant. This study demonstrates the evolutionary limitations imposed by the simultaneous optimization of folding and functional properties. 相似文献
7.
We study the folding thermodynamics of a beta-hairpin and two three-stranded beta-sheet peptides using a simplified sequence-based all-atom model, in which folding is driven mainly by backbone hydrogen bonding and effective hydrophobic attraction. The native populations obtained for these three sequences are in good agreement with experimental data. We also show that the apparent native population depends on which observable is studied; the hydrophobicity energy and the number of native hydrogen bonds give different results. The magnitude of this dependence matches well with the results obtained in two different experiments on the beta-hairpin. 相似文献
8.
In this work, we apply a detailed all‐atom model with a transferable knowledge‐based potential to study the folding kinetics of Formin‐Binding protein, FBP28, which is a canonical three‐stranded β‐sheet WW domain. Replica exchange Monte Carlo simulations starting from random coils find native‐like (Cα RMSD of 2.68 Å) lowest energy structure. We also study the folding kinetics of FBP28 WW domain by performing a large number of ab initio Monte Carlo folding simulations. Using these trajectories, we examine the order of formation of two β‐hairpins, the folding mechanism of each individual β‐hairpin, and transition state ensemble (TSE) of FBP28 WW domain and compare our results with experimental data and previous computational studies. To obtain detailed structural information on the folding dynamics viewed as an ensemble process, we perform a clustering analysis procedure based on graph theory. Further, a rigorous Pfold analysis is used to obtain representative samples of the TSEs showing good quantitative agreement between experimental and simulated Φ values. Our analysis shows that the turn structure between first and second β strands is a partially stable structural motif that gets formed before entering the TSE in FBP28 WW domain and there exist two major pathways for the folding of FBP28 WW domain, which differ in the order and mechanism of hairpin formation. Proteins 2011. © 2011 Wiley‐Liss, Inc. 相似文献
9.
Human acidic fibroblast growth factor (FGF-1) is a member of the beta-trefoil superfold, a protein architecture that exhibits a characteristic threefold axis of structural symmetry. FGF-1 contains 11 beta-turns, the majority being type I 3:5; however, a type I 4:6 turn is also found at three symmetry-related locations. The relative uniqueness of the type I 4:6 turn in the FGF-1 structure suggests it may play a key role in the stability, folding, or function of the protein. To test this hypothesis a series of deletion mutations were constructed, the aim of which was to convert existing type I 4:6 turns at two locations into type I 3:5 turns. The results show it is possible to successfully substitute the type I 4:6 turn by a type I 3:5 turn with minimal impact upon protein stability or folding. Thus, these different turn structures, even though they differ in length, exhibit similar energetic properties. Additional sequence swapping mutations within the introduced type I 3:5 turns suggests that the turn sequence primarily affects stability but not turn structure (which appears dictated primarily by the local environment). Although the results suggest that a stable, foldable beta-trefoil protein may be designed utilizing a single turn type (type I 3:5), a type I 4:6 turn at turn 1 of FGF-1 appears essential for efficient mitogenic function. 相似文献
10.
The determination of the folding mechanisms of proteins is critical to understand the topological change that can propagate Alzheimer and Creutzfeld-Jakobs diseases, among others. The computational community has paid considerable attention to this problem; however, the associated time scale, typically on the order of milliseconds or more, represents a formidable challenge. Ab initio protein folding from long molecular dynamics simulations or ensemble dynamics is not feasible with ordinary computing facilities and new techniques must be introduced. Here we present a detailed study of the folding of a 16-residue beta-hairpin, described by a generic energy model and using the activation-relaxation technique. From a total of 90 trajectories at 300 K, three folding pathways emerge. All involve a simultaneous optimization of the complete hydrophobic and hydrogen bonding interactions. The first two pathways follow closely those observed by previous theoretical studies (folding starting at the turn or by interactions between the termini). The third pathway, never observed by previous all-atom folding, unfolding, and equilibrium simulations, can be described as a reptation move of one strand of the beta-sheet with respect to the other. This reptation move indicates that non-native interactions can play a dominant role in the folding of secondary structures. Furthermore, such a mechanism mediated by non-native hydrogen bonds is not available for study by unfolding and Gō model simulations. The exact folding path followed by a given beta-hairpin is likely to be influenced by its sequence and the solvent conditions. Taken together, these results point to a more complex folding picture than expected for a simple beta-hairpin. 相似文献
11.
Mechanical stretching of secondary structures is studied through molecular dynamics simulations of a Go-like model. Force versus displacement curves are studied as a function of the stiffness and velocity of the pulling device. The succession of stretching events, as measured by the order in which contacts are ruptured, is compared to the sequencing of events during thermal folding and unfolding. Opposite cross-correlations are found for an alpha-helix and a beta-hairpin structure. In a tandem of two alpha-helices, the two constituent helices unravel nearly simultaneously. A simple condition for simultaneous versus sequential unraveling of repeat units is presented. 相似文献
12.
Plasticins (23 long-residue glycine-leucine-rich dermaseptin-related peptides produced by the skin of South American hylids) have very similar amino acid sequences, hydrophobicities, and amphipathicities, but differ in their membrane-damaging properties and structurations (i.e. destabilized helix states, beta-hairpin, beta-sheet, and disordered states) at anionic and zwitterionic membrane interfaces. Structural malleability of plasticins in aqueous solutions together with parameters that may govern their ability to fold within beta-hairpin like structures were analyzed through circular dichroism and FTIR spectroscopic studies completed by molecular dynamics simulations in polar mimetic media. The goal of this study was to probe to which extent pre-existent peptide conformations, i.e. intrinsic \"conformational landscape\", may be responsible for variability in bioactive conformation and antimicrobial/hemolytic mechanisms of action of these peptides in relation with their various membrane disturbing properties. All plasticins present a turn region that does not always result in folding into a beta-hairpin shaped conformation. Residue at position 8 plays a major role in initiating the folding, while position 12 is not critical. Conformational stability has no major impact on antimicrobial efficacy. However, preformed beta-hairpin in solution may act as a conformational lock that prevents switch to alpha-helical structure. This lock lowers the antimicrobial efficiency and explains subtle differences in potencies of the most active antimicrobial plasticins. 相似文献
13.
The folding of the amyloidogenic H1 peptide MKHMAGAAAAGAVV taken from the syrian hamster prion protein is explored in explicit aqueous solution at 300 K using long time scale all-atom molecular dynamics simulations for a total simulation time of 1.1 mus. The system, initially modeled as an alpha-helix, preferentially adopts a beta-hairpin structure and several unfolding/refolding events are observed, yielding a very short average beta-hairpin folding time of approximately 200 ns. The long time scale accessed by our simulations and the reversibility of the folding allow to properly explore the configurational space of the peptide in solution. The free energy profile, as a function of the principal components (essential eigenvectors) of motion, describing the main conformational transitions, shows the characteristic features of a funneled landscape, with a downhill surface toward the beta-hairpin folded basin. However, the analysis of the peptide thermodynamic stability, reveals that the beta-hairpin in solution is rather unstable. These results are in good agreement with several experimental evidences, according to which the isolated H1 peptide adopts very rapidly in water beta-sheet structure, leading to amyloid fibril precipitates [Nguyen et al., Biochemistry 1995;34:4186-4192; Inouye et al., J Struct Biol 1998;122:247-255]. Moreover, in this article we also characterize the diffusion behavior in conformational space, investigating its relations with folding/unfolding conditions. 相似文献
14.
In order to check our current knowledge on the principles involved in beta-hairpin formation, we have modified the sequence of a 3:5 beta-hairpin forming peptide with two different purposes, first to increase the stability of the formed 3:5 beta-hairpin, and second to convert the 3:5 beta-hairpin into a 2:2 beta-hairpin. The conformational behavior of the designed peptides was investigated in aqueous solution and in 30% trifluoroethanol (TFE) by analysis of the following nuclear magnetic resonance (NMR) parameters: nuclear Overhauser effect (NOE) data, and C(alpha)H, (13)C(alpha), and (13)C(beta) conformational shifts. From the differences in the ability to adopt beta-hairpin structures in these peptides, we have arrived to the following conclusions: (i) beta-Hairpin population increases with the statistical propensity of residues to occupy each turn position. (ii) The loop length, and in turn, the beta-hairpin type, can be modified as a function of the type of turn favored by the loop sequence. These two conclusions reinforce previous results about the importance of beta-turn sequence in beta-hairpin folding. (iii) Side-chain packing on each face of the beta-sheet may play a major role in beta-hairpin stability; hence simplified analysis in terms of isolated pair interactions and intrinsic beta-sheet propensities is insufficient. (iv) Contributions to beta-hairpin stability of turn and strand sequences are not completely independent. (v) The burial of hydrophobic surface upon beta-hairpin formation that, in turn, depends on side-chain packing also contributes to beta-hairpin stability. (vi) As previously observed, TFE stabilizes beta-hairpin structures, but the extent of the contribution of different factors to beta-hairpin formation is sometimes different in aqueous solution and in 30% TFE. 相似文献
15.
The (β/α)8‐barrel is one of the most abundant folds found in enzymes. To identify the independent folding units and the segment(s) that correspond to a minimum core structure within a (β/α)8‐barrel protein, fragmentation experiments were performed with Escherichia coli phosphoribosylanthranilate isomerase, which has a single (β/α)8‐barrel domain. Our previous studies indicated that the central four β/α segments comprise an independent folding unit; whereas, the role(s) of the first two β/α segments in folding had not been clarified prior to this report. Herein, we report the design and synthesis of a series of N‐terminally deleted fragments starting with (β/α)1–5β6 as the parent construct. Analytical gel filtration and urea‐induced equilibrium unfolding experiments indicated that deletions within the N‐terminal region, that is, within the first two β/α modules, resulted in reduced stability or aggregation of the remaining segments. The (β/α)3–5β6 segment appeared to fold into a stable structure and deletion of β6 from (β/α)3–5β6 yielded (β/α)3–5, which did not form native‐like secondary structures. However, urea‐induced unfolding of (β/α)3–5, monitored by reduction of tryptophan fluorescence, indicated that the fragment contained a loosely packed hydrophobic core. Taken together, the results of our previous and present fragmentation experiments suggest the importance of the central (β/α)3–4β5 module in folding, which is a finding that is compatible with our simulated unfolding study performed previously. Proteins 2010. © 2010 Wiley‐Liss, Inc. 相似文献
16.
Kim KH Aulakh S Paetzel M 《Protein science : a publication of the Protein Society》2012,21(6):751-768
β-Barrel proteins found in the outer membrane of Gram-negative bacteria serve a variety of cellular functions. Proper folding and assembly of these proteins are essential for the viability of bacteria and can also play an important role in virulence. The β-barrel assembly machinery (BAM) complex, which is responsible for the proper assembly of β-barrels into the outer membrane of Gram-negative bacteria, has been the focus of many recent studies. This review summarizes the significant progress that has been made toward understanding the structure and function of the bacterial BAM complex. 相似文献
17.
George Flouret Tadeusz Majewski Lajos Balaspiri William Brieher Kevin Mahan Olivier Chaloin Laird Wilson Jirina Slaninová 《Journal of peptide science》2002,8(7):314-326
Twenty analogues were synthesized of [Pmp1, D-Trp2, Arg8]oxytocin, PA, (Pmp = beta,beta-pentamethylene-beta-mercaptopropionic acid), a potent antagonist of the uterotonic effect of oxytocin in the rat (uterotonic test in vitro, pA2 = 7.77) and in the baboon. Systematic substitution of Pmp1 was made with beta-mercaptopropionic acids featuring replacement of the 4-methylene group of the cyclohexyl ring of Pmp with isosteric O, S, NH or with C=O. Since the more hydrophilic NH and C=O substitutions showed a sharply decreased antagonistic potency (rat uterotonic test in vitro), additional modifications were made to reduce their hydrophilicity. Acylation of the NH group with various acyl groups, and ketalization or thioketalization of C=O with more or less bulky substituents led to a partial restoration of potency, the N-carbamyl- and the 2-mercapto-2-adamantaneacetyl analogues being equipotent with PA. Internal cyclization by amidation of the NH-group with Gly-9, resulted in a bicyclic analogue, (cyclo 1-9)[(HN)Pmp1, Gly9]PA which was equipotent with PA. When Pen-6 was introduced into the bicyclic derivative instead of Cys-6, to reduce the flexibility of the rings, the resulting (cyclo 1-9)[(HN)Pmp1, Pen6, Gly9]PA had somewhat better potency (pA2 = 8.17) in the uterotonic test and no detectable activity in the antidiuretic assay. In the case of substitution of PA with beta,beta-(3-thiapentamethylene)-beta-mercaptopropionic acid, (S)Pmp, there was also an increase in inhibitory potency in the uterotonic test (pA2 = 8.08): the analogue had extremely weak antidiuretic activity. To establish the importance of the steric effects of the Pen-6 substitution, analogues [Pen6]PA and [(S)Pmp1, Pen6]PA were made and found to be very potent, with a pA2 of 8.72 and 8.86, respectively. The high potency of the latter analogue and its extremely weak action in the diuretic assay makes it an attractive candidate for studies on the inhibition of the biological effects of oxytocin and for the prevention of preterm labour. 相似文献
18.
Chung IM Ali M Ahmad A Chun SC Kim JT Sultana S Kim JS Min SK Seo BR 《Phytochemical analysis : PCA》2007,18(2):133-145
Two new compounds, 14-methyl stigmast-9(11)-en-3alpha-ol-3beta-D-glucopyranoside (1) and cholest-11-en-3beta, 6beta, 7alpha, 22beta-tetraol-24-one-3beta-palmitoleate (2), along with the known compound beta-sitosteryl-3beta-D-glucopyranosyl-6'-linoleiate (3), were isolated from the methanolic extract of rice (Oryza sativa) hulls. The structures of the two new compounds were elucidated using one- and two-dimensional NMR in combination with IR, EI/MS, FAB/MS, HR-EI/MS and HR-FAB/MS. In bioassays with blue-green algae, Microcystis aeruginosa UTEX 2388 and duckweed, Lemna paucicostata Hegelm 381, the efficacy of bioactivity of the two new compounds linearly increased as the concentration increased from 0.3 to 300 IgM. Compared with momilactone A, compounds 1 and 2 showed similar and higher inhibitory activities against the growth of M. aeruginosa at a concentration of 300 microM. However, compound 2 was similar to momilactone A in inhibiting L. paucicostata growth at a concentration of 300 microM. As a result, compound 2 appears to have a strong potential for the environmentally friendly control of weed and algae that are harmful to water-logged rice. 相似文献
19.
Our recently developed off-lattice bead model capable of simulating protein structures with mixed alpha/beta content has been extended to model the folding of a ubiquitin-like protein and provides a means for examining the more complex kinetics involved in the folding of larger proteins. Using trajectories generated from constant-temperature Langevin dynamics simulations and sampling with the multiple multi-histogram method over five-order parameters, we are able to characterize the free energy landscape for folding and find evidence for folding through compact intermediates. Our model reproduces the observation that the C-terminus loop structure in ubiquitin is the last to fold in the folding process and most likely plays a spectator role in the folding kinetics. The possibility of a productive metastable intermediate along the folding pathway consisting of collapsed states with no secondary structure, and of intermediates or transition structures involving secondary structural elements occurring early in the sequence, is also supported by our model. The kinetics of folding remain multi-exponential below the folding temperature, with glass-like kinetics appearing at T/T(f) approximately 0.86. This new physicochemical model, designed to be predictive, helps validate the value of modeling protein folding at this level of detail for genomic-scale studies, and motivates further studies of other protein topologies and the impact of more complex energy functions, such as the addition of solvation forces. 相似文献
20.
We show that residues at the interfaces of protein-protein complexes have higher side-chain energy than other surface residues. Eight different sets of protein complexes were analyzed. For each protein pair, the complex structure was used to identify the interface residues in the unbound monomer structures. Side-chain energy was calculated for each surface residue in the unbound monomer using our previously developed scoring function.1 The mean energy was calculated for the interface residues and the other surface residues. In 15 of the 16 monomers, the mean energy of the interface residues was higher than that of other surface residues. By decomposing the scoring function, we found that the energy term of the buried surface area of non-hydrogen-bonded hydrophilic atoms is the most important factor contributing to the high energy of the interface regions. In spite of lacking hydrophilic residues, the interface regions were found to be rich in buried non-hydrogen-bonded hydrophilic atoms. Although the calculation results could be affected by the inaccuracy of the scoring function, patch analysis of side-chain energy on the surface of an isolated protein may be helpful in identifying the possible protein-protein interface. A patch was defined as 20 residues surrounding the central residue on the protein surface, and patch energy was calculated as the mean value of the side-chain energy of all residues in the patch. In 12 of the studied monomers, the patch with the highest energy overlaps with the observed interface. The results are more remarkable when only three residues with the highest energy in a patch are averaged to derive the patch energy. All three highest-energy residues of the top energy patch belong to interfacial residues in four of the eight small protomers. We also found that the residue with the highest energy score on the surface of a small protomer is very possibly the key interaction residue. 相似文献