Folding thermodynamics of β‐hairpins studied by replica‐exchange molecular dynamics simulations

We study the differences in folding stability of β‐hairpin peptides, including GB1 hairpin and a point mutant GB1 K10G, as well as tryptophan zippers (TrpZips): TrpZip1, TrpZip2, TrpZip3‐1, and TrpZip4. By performing replica‐exchange molecular dynamics simulations with Amber03* force field (a modified version of Amber ff03) in explicit solvent, we observe ab initio folding of all the peptides except TrpZip3‐1, which is experimentally known to be the least stable among the peptides studied here. By calculating the free energies of unfolding of the peptides at room temperature and folding midpoint temperatures for thermal unfolding of peptides, we find that TrpZip4 and GB1 K10G peptides are the most stable β‐hairpins followed by TrpZip1, GB1, and TrpZip2 in the given order. Hence, the proposed K10G mutation of GB1 peptide results in enhanced stability compared to wild‐type GB1. An important goal of our study is to test whether simulations with Amber 03* model can reproduce experimentally predicted folding stability differences between these peptides. While the stabilities of GB1 and TrpZip1 yield close agreement with experiment, TrpZip2 is found to be less stable than predicted by experiment. However, as heterogenous folding of TrpZip2 may yield divergent thermodynamic parameters by different spectroscopic methods, mismatching of results with previous experimental values are not conclusive of model shortcomings. For most of the cases, molecular simulations with Amber03* can successfully reproduce experimentally known differences between the mutated peptides, further highlighting the predictive capabilities of current state‐of‐the‐art all‐atom protein force fields. Proteins 2015; 83:1307–1315. © 2015 Wiley Periodicals, Inc.     

Spectroscopic studies of structural changes in two beta-sheet-forming peptides show an ensemble of structures that unfold noncooperatively

A characterization of the conformation and stability of model peptide systems that form beta-sheets in aqueous solutions is considerably important in gaining insights into the mechanism of beta-sheet formation in proteins. We have characterized the conformation and equilibrium folding and unfolding of two 20-residue peptides whose NMR spectra suggest a three-stranded beta-sheet topology in aqueous solution: Betanova [Kortemme, T., Ramirez-Alvarado, M., and Serrano, L. (1998) Science 281, 253-256] and (D)P(D)P with d-Pro-Gly segments at the turns [Schenck, H. L., and Gellman, S. H. (1998) J. Am. Chem. Soc. 120, 4869-4870]. Both circular dichroism (CD) and infrared measurements indicate only 20-26% beta-sheet-like structure at 5 degrees C for Betanova and 42-59% beta-sheet for (D)P(D)P. For both peptides, the CD and infrared spectra change nearly linearly with increasing temperatures (or urea concentrations) and lack a sigmoidal signature characteristic of cooperative unfolding. Fluorescence resonance energy transfer (FRET) measurements between donor and acceptor molecules attached to the two ends confirm that Betanova is largely unstructured even at 10 degrees C; the average end-to-end distance estimated from FRET is closer to that of a random coil than a structured beta-sheet. In (D)P(D)P, the FRET results indicate a more compact structure that remains compact even at high temperatures (approximately 80 degrees C) or high urea concentrations (approximately 8 M). These results indicate that both these peptides access an ensemble of conformations at all temperatures or denaturant concentrations, with no significant free energy barrier separating the "folded" and "unfolded" conformations.     

Melting temperatures of nucleic acids: discrepancies in analysis

Melting temperature, T(m), is an important property of nucleic acid duplexes. It is typically determined from spectroscopic or calorimetric melting experiments. More than one analytical method has been used to extract T(m) values from experimental melting data. Unfortunately, different methods do not give the same results; the same melting data can be assigned different T(m) values depending upon which method is used to process that data. Inconsistencies or systematic errors between T(m)s reported in published data sets can be significant and add confusion to the field. Errors introduced from analysis can be greater than experimental errors, ranging from a fraction of degree to several degrees. Of the various methods, the most consistent and meaningful approach defines melting temperature as the temperature at the transition midpoint where half of the base pairs are melted and standard free energy is zero. Assuming a two-state melting behavior, we present here a set of general equations that can be used to reconcile these analytical T(m) differences and convert results to the correct melting temperatures at the transition midpoint. Melting temperatures collected from published sources, which were analyzed using different methods, can now be corrected for these discrepancies and compared on equal footing. The similar corrections apply to T(m) differences between calorimetric and spectroscopic melting curves. New algorithm for selection of linear sloping baselines, 2nd derivative method, is suggested, which can be used to automate melting curve analysis.     

Studies of the minimum hydrophobicity of alpha-helical peptides required to maintain a stable transmembrane association with phospholipid bilayer membranes

The effects of the hydrophobicity and the distribution of hydrophobic residues on the surfaces of some designed alpha-helical transmembrane peptides (acetyl-K2-L(m)-A(n)-K2-amide, where m + n = 24) on their solution behavior and interactions with phospholipids were examined. We find that although these peptides exhibit strong alpha-helix forming propensities in water, membrane-mimetic media, and lipid model membranes, the stability of the helices decreases as the Leu content decreases. Also, their binding to reversed phase high-performance liquid chromatography columns is largely determined by their hydrophobicity and generally decreases with decreases in the Leu/Ala ratio. However, the retention of these peptides by such columns is also affected by the distribution of hydrophobic residues on their helical surfaces, being further enhanced when peptide helical hydrophobic moments are increased by clustering hydrophobic residues on one side of the helix. This clustering of hydrophobic residues also increases peptide propensity for self-aggregation in aqueous media and enhances partitioning of the peptide into lipid bilayer membranes. We also find that the peptides LA3LA2 [acetyl-K2-(LAAALAA)3LAA-K2-amide] and particularly LA6 [acetyl-K2-(LAAAAAA)3LAA-K2-amide] associate less strongly with and perturb the thermotropic phase behavior of phosphatidylcholine bilayers much less than peptides with higher L/A ratios. These results are consistent with free energies calculated for the partitioning of these peptides between water and phospholipid bilayers, which suggest that LA3LA2 has an equal tendency to partition into water and into the hydrophobic core of phospholipid model membranes, whereas LA6 should strongly prefer the aqueous phase. We conclude that for alpha-helical peptides of this type, Leu/Ala ratios of greater than 7/17 are required for stable transmembrane associations with phospholipid bilayers.     

An implicit membrane generalized born theory for the study of structure, stability, and interactions of membrane proteins

Exploiting recent developments in generalized Born (GB) electrostatics theory, we have reformulated the calculation of the self-electrostatic solvation energy to account for the influence of biological membranes. Consistent with continuum Poisson-Boltzmann (PB) electrostatics, the membrane is approximated as an solvent-inaccessible infinite planar low-dielectric slab. The present membrane GB model closely reproduces the PB electrostatic solvation energy profile across the membrane. The nonpolar contribution to the solvation energy is taken to be proportional to the solvent-exposed surface area (SA) with a phenomenological surface tension coefficient. The proposed membrane GB/SA model requires minor modifications of the pre-existing GB model and appears to be quite efficient. By combining this implicit model for the solvent/bilayer environment with advanced computational sampling methods, like replica-exchange molecular dynamics, we are able to fold and assemble helical membrane peptides. We examine the reliability of this model and approach by applications to three membrane peptides: melittin from bee venom, the transmembrane domain of the M2 protein from Influenza A (M2-TMP), and the transmembrane domain of glycophorin A (GpA). In the context of these proteins, we explore the role of biological membranes (represented as a low-dielectric medium) in affecting the conformational changes in melittin, the tilt of transmembrane peptides with respect to the membrane normal (M2-TMP), helix-to-helix interactions in membranes (GpA), and the prediction of the configuration of transmembrane helical bundles (GpA). The present method is found to perform well in each of these cases and is anticipated to be useful in the study of folding and assembly of membrane proteins as well as in structure refinement and modeling of membrane proteins where a limited number of experimental observables are available.     

Studies on the cellular uptake of substance P and lysine-rich, KLA-derived model peptides

In the last decade many peptides have been shown to be internalized into various cell types by different, poorly characterized mechanisms. This review focuses on uptake studies with substance P (SP) aimed at unravelling the mechanism of peptide-induced mast cell degranulation, and on the characterization of the cellular uptake of designed KLA-derived model peptides.Studies on structure-activity relationships and receptor autoradiography failed to detect specific peptide receptors for the undecapeptide SP on mast cells. In view of these findings, a direct interaction of cationic peptides with heterotrimeric G proteins without the participation of a receptor has been proposed. Such a process would require insertion into and translocation of peptides across the plasma membrane.In order to clarify whether a transport of cationic peptides into rat peritoneal mast cells is possible, transport studies were performed by confocal laser scanning microscopy (CLSM) using fluorescence-labeled Arg(3),Orn(7)-SP and its D-amino acid analog, all-D-Arg(3),Orn(7)-SP, as well as by electron microscopic autoradiography using (3)H-labelled SP and (125)I-labelled all-D-SP. The results obtained by CLSM directly showed translocation of SP peptides into pertussis toxin-treated cells. Kinetic experiments indicated that the translocation process was rapid, occurring within a few seconds. Mast cell degranulation induced by analog of magainin 2 amide, neuropeptide Y and the model peptide acetyl-KLALKLALKALKAALKLA-amide was also found to be very fast, pointing to an extensive translocation of the peptides. In order to learn more about structural requirements for the cellular uptake of peptides, the translocation behavior of a set of systematically modified KLA-based model peptides has been studied in detail. By two different protocols for determining the amount of internalized peptide, evidence was found that the structure of the peptides only marginally affects their uptake, whereas the efflux of cationic, amphipathic peptides is strikingly diminished, thus allowing their enrichment within the cells. Although the mechanism of cellular uptake, consisting of energy-dependent and -independent contributions, is not well understood, KLA-derived peptides have been shown to deliver various cargos (PNAs, peptides) into cells. The results obtained with SP- and KLA-derived peptides are discussed in the context of the current literature.     

Identification of carbonylated peptides by tandem mass spectrometry using a precursor ion-like scan in negative ion mode

Reactive oxygen species (ROS) can oxidize proteins at almost any amino acid residue. Whereas some modifications are reversible within the cells, the higher oxidation states are especially irreversible. These irreversible post translational modifications are widely used as biomarkers of oxidative stress, such as protein carbonylation, which refers to aldehydes, ketones and lactams as 'reactive carbonyl groups'. This study relied on a set of synthetic peptides containing a C-terminal aldehyde (arginal) or modification with pyruvic acid (ketone) or 4-hydroxynonenal (aldehyde) at lysine or histidine residues, as well as peptides containing pyroglutamic acid (oxidation product of proline) and 2-amino-3-butyric acid (oxidation product of threonine). The carbonylation sites were specifically derivatized with 2,4-dinitrophenylhydrazine (DNPH) and the fragmentation behavior of the products investigated in electrospray ionization (ESI-) MS. Importantly, the DNPH-labeled carbonylated peptides showed favorable ionization behaviors in negative ion mode ESI, providing a sensitive detection method. Regular peptides were mostly discriminated under these conditions. Among the fragmentation techniques tested for the negatively charged ions, pulsed Q dissociation provided three diagnostic ions at m/z values 152.0, 163.1 and 179.0, specific for DNPH-modified peptides. These marker ions were successfully applied to detect the carbonylated model peptides in a spiked tryptic digest of bovine serum albumin and a complex protein mixture obtained from HeLa cells.     

Studies of DNA dumbbells. II. Construction and characterization of DNA dumbbells with a 16 base-pair duplex stem and Tn end loops (n = 2, 3, 4, 6, 8, 10, 14).

The preparation and characterization of DNA dumbbells that contain the 16 base-pair duplex sequences 5'G-C-A-T-A-G-A-T-G-A-G-A-A-T-G-C3' (set 1) and 5'G-C-A-T-C-A-T-C-G-A-T-G-A-T-G-C3' (set 2) are reported. The dumbbells of set 1 have the duplex stem nucleated on both ends by Tn (n = 2, 3, 4, 6, 8, 10, and 14) loops. The dumbbells of set 2 have Tn (n = 2, 4, 8, 10) end loops. For the molecules of set 1, effects of end loop size on the electrophoretic mobility, CD and UV absorbance spectra, and cleavage by restriction enzymes, were investigated. Effects of loop size on the CD spectra and restriction enzyme cleavage of the molecules of set 2 were also examined. Optical melting curves of the molecules of set 1 were collected as a function of sodium ion concentration from 30 to 120 mM. These investigations revealed that as loop size decreases, the electrophoretic mobilities, rates of enzyme cleavage, and optical melting temperatures increase. For end loops with at least three T's the observed increases are inversely proportional to loop size. The behavior of the dumbbell with T2 end loops departs from this linear dependence and is anomalous in every experimental context. For molecules with end loops comprised of at least four T's CD spectra were virtually indistinguishable. However, these spectra differed considerably from the CD spectrum of the T2-looped molecule. The CD spectrum of the dumbbell with T3 end loops displayed features common to the dumbbells with larger loops and T2 end loops. Thermodynamic evidence that the terminal G.C base pairs (bps) nucleating the T2 end loops were intact was obtained from a comparison of the melting temperature of this molecule with that of a DNA dumbbell containing the 14 central bps of the set 1 duplex sequence linked instead by end loops comprised of the four base sequence, C-T-T-C. The tm of this latter molecule was determined to be 9 degrees C less than that of the former dumbbell assumed to contain a 16-bp stem and T2 end loops.(ABSTRACT TRUNCATED AT 400 WORDS)     

Synthesis and study of peptides with semirigid i and i + 7 side-chain bridges designed for alpha-helix stabilization

A search for conformational constraints on the peptide alpha-helical conformation indicated that para-substituted amino acid derivatives of a benzene ring might be suitable for linking pairs of side chains that are separated by two turns of the helix. A 14-residue synthetic, amphiphilic alpha-helical peptide model system has been used to study the helix stabilizing effects of a series of four such bridges having constitutionally isomeric structures. These bridges were used to link positions 3 and 10 of the model peptides. The peptides were synthesized in good yield by standard solid-phase methods, including cyclization on the solid support. They were then studied for their solution conformations and melting behavior by circular dichroism (CD) spectropolarimetry, and for their elution behavior on reversed-phase HPLC columns. In aqueous solution and in 50% (v/v) trifluoroethanol, the most effective bridge for helix stabilization consisted of a 4-(aminomethyl)phenylacetic acid residue (AMPA) linked by amide bonds to the side chain functional groups of a (S)-2,3-diaminopropionic acid residue (Dap) in position 3 of the model peptide and an aspartic acid residue in position 10. This Dap3(AMPA), Asp10 bridge was about as effective as two Lys(i), Asp(i+4) lactam bridges incorporated linking residues 3 and 7, and 10 and 14, in the same model peptide sequence. This suggests that it is worth about 1 kcal/mol of helix stabilization energy.     

The role of charge and hydrophobicity in peptide-lipid interaction: a comparative study based on tryptophan fluorescence measurements combined with the use of aqueous and hydrophobic quenchers

The interaction of interrelated model peptides with model membranes has been studied by techniques based on tryptophan fluorescence. The peptides used are derivatives of the sequence H-Ala-Met-Leu-Trp-Ala-OH, which was designed for this purpose. Several modifications yielded a set of 13 penta- and hexapeptides varying in net charge, hydrophobicity, charge distribution, and the intramolecular position of the tryptophan residue with respect to the charge(s). The affinity of these peptides for small unilamellar vesicles (SUV) consisting of zwitterionic egg phosphatidylcholine (eggPC) and negatively charged beef heart cardiolipin (bhCL) has been investigated in a comparative way. The criteria for affinity comprise (1) intrinsic fluorescence changes upon titration of the peptides with the lipid vesicles, (2) reduced accessibility of the peptides to aqueous quenchers of tryptophan fluorescence (I- and acrylamide) in the presence of lipid, and (3) exposure to membrane-incorporated fluorescence quenchers, brominated phosphatidylcholines (BrPC). Application of BrPC brominated at different positions along the acyl chains provided information on the membrane topology of the peptides. With respect to the extent of affinity for zwitterionic membranes, the overall hydrophobicity of the peptides is the main determinant. A comparison of the affinity for PC of equally hydrophobic peptides carrying either a single positive or negative charge reveals preferential interaction of the cationic peptide. Both hydrophobic and electrostatic interactions determine the affinity of positively charged mono- and divalent peptides for CL vesicles. The distribution of the charged moieties in divalent positively charged peptides, either both at one end of the molecule or one at each end, has little influence on the affinity of these peptides for CL but does affect the extent of exposure to BrPC. Upon decreasing the surface charge density of the vesicles by diluting CL with increasing amounts of PC, both types of peptides show different behavior. The position of the tryptophan relative to the charged moiety in the peptide molecule is shown to affect the fluorescent properties upon interaction with vesicles. Concerning the membrane topology, all peptides adopt a localization near the membrane surface, with the neutral peptides inserting slightly deeper into the bilayer than the charged peptides. The results allow a comparative analysis of the factors determining the extents and modes of lipid-model peptide interaction; in addition, the validity of the methods applied is discussed.     

Alpha1-microglobulin chromophores are located to three lysine residues semiburied in the lipocalin pocket and associated with a novel lipophilic compound

Alpha1-microglobulin (alpha1m) is an electrophoretically heterogeneous plasma protein. It belongs to the lipocalin superfamily, a group of proteins with a three-dimensional (3D) structure that forms an internal hydrophobic ligand-binding pocket. Alpha1m carries a covalently linked unidentified chromophore that gives the protein a characteristic brown color and extremely heterogeneous optical properties. Twenty-one different colored tryptic peptides corresponding to residues 88-94, 118-121, and 122-134 of human alpha1m were purified. In these peptides, the side chains of Lys92, Lys118, and Lys130 carried size heterogeneous, covalently attached, unidentified chromophores with molecular masses between 122 and 282 atomic mass units (amu). In addition, a previously unknown uncolored lipophilic 282 amu compound was found strongly, but noncovalently associated with the colored peptides. Uncolored tryptic peptides containing the same Lys residues were also purified. These peptides did not carry any additional mass (i.e., chromophore) suggesting that only a fraction of the Lys92, Lys118, and Lys130 are modified. The results can explain the size, charge, and optical heterogeneity of alpha1m. A 3D model of alpha1m, based on the structure of rat epididymal retinoic acid-binding protein (ERABP), suggests that Lys92, Lys118, and Lys130 are semiburied near the entrance of the lipocalin pocket. This was supported by the fluorescence spectra of alpha1m under native and denatured conditions, which indicated that the chromophores are buried, or semiburied, in the interior of the protein. In human plasma, approximately 50% of alpha1m is complex bound to IgA. Only the free alpha1m carried colored groups, whereas alpha1m linked to IgA was uncolored.     

Lipid-destabilising properties of a peptide with structural plasticity

The Chameleon peptide (Cham) is a peptide designed from two regions of the GB1 protein, one folded as an alpha-helix and the other as a beta structure. Depending on the environment, the Cham peptide adopts an alpha or a beta conformation when inserted in different locations of GB1. This environment dependence is also observed for tilted peptides. These short protein fragments, able to destabilise organised system, are mainly folded in beta structure in water and in alpha helix in a hydrophobic environment, like the lipid bilayer. In this paper, we tested whether the Cham peptide can be qualified as a tilted peptide. For this, we have compared the properties of Cham peptide (hydrophobicity, destabilising properties, conformation) to those of tilted peptides. The results suggest that Cham is a tilted peptide. Our study, together the presence of tilted fragments in transconformational proteins, suggests a relationship between tilted peptides and structural lability.     

Exploring the free energy landscape of a model β‐hairpin peptide and its isoform

Secondary structural transitions from α‐helix to β‐sheet conformations are observed in several misfolding diseases including Alzheimer's and Parkinson's. Determining factors contributing favorably to the formation of each of these secondary structures is therefore essential to better understand these disease states. β‐hairpin peptides form basic components of anti‐parallel β‐sheets and are suitable model systems for characterizing the fundamental forces stabilizing β‐sheets in fibrillar structures. In this study, we explore the free energy landscape of the model β‐hairpin peptide GB1 and its E2 isoform that preferentially adopts α‐helical conformations at ambient conditions. Umbrella sampling simulations using all‐atom models and explicit solvent are performed over a large range of end‐to‐end distances. Our results show the strong preference of GB1 and the E2 isoform for β‐hairpin and α‐helical conformations, respectively, consistent with previous studies. We show that the unfolded states of GB1 are largely populated by misfolded β‐hairpin structures which differ from each other in the position of the β‐turn. We discuss the energetic factors contributing favorably to the formation of α‐helix and β‐hairpin conformations in these peptides and highlight the energetic role of hydrogen bonds and non‐bonded interactions. Proteins 2014; 82:2394–2402. © 2014 Wiley Periodicals, Inc.     

Role of peptide structure in lipid-peptide interactions: high-sensitivity differential scanning calorimetry and electron spin resonance studies of the structural properties of dimyristoylphosphatidylcholine membranes interacting with pentagastrin-related pentapeptides

The effects of amino acid substitutions in the pentapeptide pentagastrin on the nature of its interactions with dimyristoylphosphatidylcholine (DMPC) are assessed by differential scanning calorimetry and electron spin resonance. In two peptide analogues, the Asp at position 4 in pentagastrin (N-t-Boc-beta-Ala-Trp-Met-Asp-Phe-NH2) is replaced by Gly or Phe. These uncharged, more hydrophobic peptides have little effect on the transition temperature of DMPC, but they broaden the transition and lower the transition enthalpy as do integral membrane proteins. These peptides also mimic the behavior of integral membrane proteins in decreasing the order of a 5-doxylstearic acid spin probe below the transition temperature and in exhibiting a second immobilized lipid component using a 16-doxylstearic acid spin probe in DMPC. Three charged peptides were studied: pentagastrin, an analogue with positions 4 and 5 reversed (i.e., ending in Phe-Asp-NH2), and one with Asp replaced by Arg at position 4. All three of these charged peptides altered the phase transition behavior of DMPC to give two components, one above and one below the transition temperature of the pure lipid. With increasing peptide concentration, the higher melting transition became more prominent. The arginine-containing peptide produced the largest shifts in melting temperature followed by pentagastrin and then the "reversed" peptide. The arginine-containing peptide also increased the enthalpy of the transition. These peptides also increased the ordering of DMPC below the phase transition as measured with both 5- and 16-doxylstearic acid. The ordering effect was most pronounced with the arginine-containing peptide using the 5-doxylstearic acid probe. The results demonstrate that even the zwitterionic DMPC can interact more strongly with positively charged peptides than with negatively charged ones. In addition, peptide sequence as well as composition is important in determining the nature of peptide-lipid interactions. The markedly different effects of these pentagastrin peptides on the phase transition and motional properties of DMPC occur despite the similar depth of burial of these peptides with DMPC.     

Ab initio prediction of peptide-MHC binding geometry for diverse class I MHC allotypes

Since determining the crystallographic structure of all peptide-MHC complexes is infeasible, an accurate prediction of the conformation is a critical computational problem. These models can be useful for determining binding energetics, predicting the structures of specific ternary complexes with T-cell receptors, and designing new molecules interacting with these complexes. The main difficulties are (1) adequate sampling of the large number of conformational degrees of freedom for the flexible peptide, (2) predicting subtle changes in the MHC interface geometry upon binding, and (3) building models for numerous MHC allotypes without known structures. Whereas previous studies have approached the sampling problem by dividing the conformational variables into different sets and predicting them separately, we have refined the Biased-Probability Monte Carlo docking protocol in internal coordinates to optimize a physical energy function for all peptide variables simultaneously. We also imitated the induced fit by docking into a more permissive smooth grid representation of the MHC followed by refinement and reranking using an all-atom MHC model. Our method was tested by a comparison of the results of cross-docking 14 peptides into HLA-A*0201 and 9 peptides into H-2K(b) as well as docking peptides into homology models for five different HLA allotypes with a comprehensive set of experimental structures. The surprisingly accurate prediction (0.75 A backbone RMSD) for cross-docking of a highly flexible decapeptide, dissimilar to the original bound peptide, as well as docking predictions using homology models for two allotypes with low average backbone RMSDs of less than 1.0 A illustrate the method's effectiveness. Finally, energy terms calculated using the predicted structures were combined with supervised learning on a large data set to classify peptides as either HLA-A*0201 binders or nonbinders. In contrast with sequence-based prediction methods, this model was also able to predict the binding affinity for peptides to a different MHC allotype (H-2K(b)), not used for training, with comparable prediction accuracy.     

Characterization of large peptide fragments derived from the N-terminal domain of the ribosomal protein L9: definition of the minimum folding motif and characterization of local electrostatic interactions

A set of peptides derived from the N-terminal domain of the ribosomal protein L9 (NTL9) have been characterized in an effort to define the minimum unit of this domain required to fold and to provide model peptides for the analysis of electrostatic interactions in the unfolded state. NTL9 is a 56-residue alpha-beta protein with a beta1-loop-beta2-alpha1-beta3-alpha2 topology. The beta-sheet together with the first helix comprise a simple example of a common supersecondary motif called the split beta-alpha-beta fold. Peptides corresponding to the beta1-loop-beta2 unit are unstructured even when constrained by an introduced disulfide. The pK(a)s of Asp-8 and Glu-17 in these peptides are slightly lower than the values found for shorter peptides but are considerably higher than the values in NTL9. A 34-residue peptide, which represents the beta1-loop-beta2-alpha1 portion of NTL9, is also unstructured. In contrast, a 39-residue peptide corresponding to the entire split beta-alpha-beta motif is folded and monomeric as judged by near- and far-UV CD, two-dimensional NMR, ANS binding experiments, pK(a) measurements, and analytical ultracentrifugation. The fold is very similar to the structure of this region in the intact protein. Thermal and urea unfolding experiments show that it is cooperatively folded with a DeltaG degrees of unfolding of 1.8-2.0 kcal/mol and a T(m) of 58 degrees C. This peptide represents the first demonstration of the independent folding of an isolated split beta-alpha-beta motif, and is one of only four naturally occurring sequences of fewer than 40 residues that has been shown to fold cooperatively in the absence of disulfides or ligand binding.     

Insights into conformation and dynamics of protein GB1 during folding and unfolding by NMR

Understanding protein stability requires characterization of structural determinants of the folded and unfolded states. Many proteins are capable of populating partially folded states under specific solution conditions. Occasionally, coexistence of the folded and an unfolded state under non- or mildly denaturing conditions can be observed by NMR, allowing us to structurally probe these states under identical conditions. Here we report on a destabilized mutant of the B1 domain of protein G (GB1) whose equilibrium unfolding was systematically investigated. Backbone amide residual dipolar couplings (RDCs), the tryptophan Nepsilon-H resonance and the amide nitrogen transverse relaxation rates (R2s) for varying pH values and different temperatures were measured. The backbone amide RDCs indicate that prior to complete unfolding, two melting hot spots are formed at the turn around T11, L12 and K13 and the N terminus of the helix at A24 and T25. The RDCs for the low pH, thermally unfolded state of GB1 are very small and do not indicate the presence of any native-like structure. Amide nitrogen transverse relaxation rates for GB1 in the folded state at different temperatures exhibit large contributions from exchange processes and the associated dynamics display considerable heterogeneity. Our data provide clear evidence for intermediate conformations and multi-state equilibrium un/folding for this GB1 variant.