Development of the multiple sequence approximation within the AGADIR model of alpha-helix formation: comparison with Zimm-Bragg and Lifson-Roig formalisms

What determines the activity of antimicrobial and cytolytic peptides in model membranes

We previously proposed three hypotheses relating the mechanism of antimicrobial and cytolytic peptides in model membranes to the Gibbs free energies of binding and insertion into the membrane [Almeida, P. F., and Pokorny, A. (2009) Biochemistry 48, 8083-8093]. Two sets of peptides were designed to test those hypotheses, by mutating of the sequences of δ-lysin, cecropin A, and magainin 2. Peptide binding and activity were measured on phosphatidylcholine membranes. In the first set, the peptide charge was changed by mutating basic to acidic residues or vice versa, but the amino acid sequence was not altered much otherwise. The type of dye release changed from graded to all-or-none according to prediction. However, location of charged residues in the sequence with the correct spacing to form salt bridges failed to improve binding. In the second set, the charged and other key residues were kept in the same positions, whereas most of the sequence was significantly but conservatively simplified, maintaining the same hydrophobicity and amphipathicity. This set behaved completely different from predicted. The type of release, which was expected to be maintained, changed dramatically from all-or-none to graded in the mutants of cecropin and magainin. Finally, contrary to the hypotheses, the results indicate that the Gibbs energy of binding to the membrane, not the Gibbs energy of insertion, is the primary determinant of peptide activity.

Identification of lipopolysaccharide-binding peptide regions within HMGB1 and their effects on subclinical endotoxemia in a mouse model

Lipopolysaccharide (LPS) triggers deleterious systemic inflammatory responses when released into the circulation. LPS-binding protein (LBP) in the serum plays an important role in modifying LPS toxicity by facilitating its interaction with LPS signaling receptors, which are expressed on the surface of LPS-responsive cells. We have previously demonstrated that high mobility group box 1 (HMGB1) can bind to and transfer LPS, consequently increasing LPS-induced TNF-α production in human peripheral blood mononuclear cells (PBMCs). We report here on the identification of two LPS-binding domains within HMGB1. Furthermore, using 12 synthetic HMGB1 peptides, we define the LPS-binding regions within each domain. Among them, synthetic peptides HPep1 and HPep6, which are located in the A and B box domains of HMGB1, bind to the polysaccharide and lipid A moieties of LPS respectively. Both HPep1 and HPep6 peptides inhibited binding of LPS to LBP and HMGB1, LBP-mediated LPS transfer to CD14, and cellular uptake of LPS in RAW264.7 cells. These peptides also inhibited LPS-induced TNF-α release in human PBMCs and induced lower levels of TNF-α in the serum in a subclinical endotoxemia mouse model. These results indicate that HMGB1 has two LPS-binding peptide regions that can be utilized to design anti-sepsis or LPS-neutralizing therapeutics.

Inducing toxicity by introducing a leucine-zipper-like motif in frog antimicrobial peptide, magainin 2

Cytotoxicity, a major obstacle in therapeutic application of antimicrobial peptides, is controlled by leucine-zipper-like sequences in melittin and other naturally occurring antimicrobial peptides. Magainin 2 shows significantly lower cytotoxicity than many naturally occurring antimicrobial peptides and lacks this structural element. To investigate the consequences of introducing a leucine zipper sequence in magainin 2, a novel analogue (Mag-mut) was designed by rearranging only the positions of its hydrophobic amino acids to include this structural element. Both magainin 2 and Mag-mut showed appreciable similarities in their secondary structures in the presence of negatively charged lipid vesicles, in localizing and permeabilizing the selected bacteria and exhibiting bactericidal activities. However, Mag-mut bound and localized strongly on to the mammalian cells tested and exhibited significantly higher cytotoxicity than magainin 2. Only Mag-mut, but not magainin 2, permeabilized human red blood cells and zwitterionic lipid vesicles. In contrast with magainin 2, Mag-mut self-assembled in an aqueous environment and bound co-operatively on to zwitterionic lipid vesicles. The peptides formed pores of different sizes on to a selected mammalian cell. The results of the present study indicate an important role of the leucine zipper sequence in the cytotoxicity of Mag-mut and demonstrate that its introduction into a non-toxic peptide, without altering the amino acid composition, can render cytotoxicity.

Modulating repeat protein stability: the effect of individual helix stability on the collective behavior of the ensemble

Repeat proteins are tandem arrays of a small structural motif, in which tertiary structure is stabilized by interactions within a repeat and between neighboring repeats. Several studies have shown that this modular structure is manifest in modular thermodynamic properties. Specifically, the global stability of a repeat protein can be described by simple linear models, considering only two parameters: the stability of the individual repeated units (H) and the coupling interaction between the units (J). If the repeat units are identical, single values of H and J, together with the number of repeated units, is sufficient to completely describe the thermodynamic behavior of any protein within a series. In this work, we demonstrate how the global stability of a repeat protein can be changed, in a predictable fashion, by modifying only the H parameter. Taking a previously characterized series of consensus tetratricopeptide repeats (TPR) (CTPRa) proteins, we introduced mutations into the basic repeating unit, such that the stability of the individual repeat unit was increased, but its interaction with neighboring units was unchanged. In other words, we increased H but kept J constant. We demonstrated that the denaturation curves for a series of such repeat proteins can be fit and additional curves can be predicted by the one-dimensional Ising model in which only H has changed from the original fit for the CTPRa series. Our results show that we can significantly increase the stability of a repeat protein by rationally increasing the stability of the units (H), whereas the interaction between repeats (J) remains unchanged.

Malleability of folding intermediates in the homeodomain superfamily

Members of the homeodomain superfamily are three-helix bundle proteins whose second and third helices form a helix-turn-helix motif (HTH). Their folding mechanism slides from the ultrafast, three-state framework mechanism for the engrailed homeodomain (EnHD), in which the HTH motif is independently stable, to an apparent two-state nucleation-condensation model for family members with an unstable HTH motif. The folding intermediate of EnHD has nearly native HTH structure, but it is not docked with helix1. The determinant of whether two- or three-state folding was hypothesized to be the stability of the HTH substructure. Here, we describe a detailed Φ-value analysis of the folding of the Pit1 homeodomain, which has similar ultrafast kinetics to that of EnHD. Formation of helix1 was strongly coupled with formation of HTH, which was initially surprising because they are uncoupled in the EnHD folding intermediate. However, we found a key difference between Pit1 and EnHD: The isolated peptide corresponding to the HTH motif in Pit1 was not folded in the absence of H1. Independent molecular dynamics simulations of Pit1 unfolding found an intermediate with H1 misfolded onto the HTH motif. The Pit1 folding pathway is the connection between that of EnHD and the slower folding homeodomains and provides a link in the transition of mechanisms from two- to three-state folding in this superfamily. The malleability of folding intermediates can lead to unstable substructures being stabilized by a variety of nonnative interactions, adding to the continuum of folding mechanisms.

Structural and mechanistic basis of Parl activity and regulation

The mitochondrial rhomboid protease Parl governs apoptosis, morphology, metabolism and might be implicated in Parkinson's disease, but the structural basis of its activity and complex regulation remain unknown. We report the discovery of γ-cleavage, a proteolytic event on the loop connecting the first transmembrane helix (TMH) of Parl to the 6-TMH catalytic rhomboid domain of the protease. This cleavage disrupts the '1+6' structure that defines every mitochondrial rhomboid and generates a new form of Parl, PROD (Parl-rhomboid-domain). Structure-function analysis of Parl suggests that γ-cleavage could be implicated in eliminating Parl proteolytic activity, and structural modeling of PROD reveals structural conservation with the bacterial rhomboid GlpG. However, unlike bacterial rhomboids, which employ a diad-based mechanism of catalysis, Parl appears to use a conserved mitochondrial rhomboid-specific Asp residue on TMH-5 in a triad-based mechanism of catalysis. This work provides unexpected insights into the structural determinants regulating Parl stability and activity in vivo, and reveals a complex cascade of proteolytic events controlling the function of the protease in the mitochondrion.

The denatured state dictates the topology of two proteins with almost identical sequence but different native structure and function

The protein folding problem is often studied by comparing the mechanisms of proteins sharing the same structure but different sequence. The recent design of the two proteins G(A)88 and G(B)88, displaying different structures and functions while sharing 88% sequence identity (49 out of 56 amino acids), allows the unique opportunity for a complementary approach. At which stage of its folding pathway does a protein commit to a given topology? Which residues are crucial in directing folding mechanisms to a given structure? By using a combination of biophysical and computational techniques, we have characterized the folding of both G(A)88 and G(B)88. We show that, contrary to expectation, G(B)88, characterized by a native α+β fold, displays in the denatured state a content of native-like helical structure greater than G(A)88, which is all-α in its native state. Both experiments and simulations indicate that such residual structure may be tuned by changing pH. Thus, despite the high sequence identity, the folding pathways for these two proteins appear to diverge as early as in the denatured state. Our results suggest a mechanism whereby protein topology is committed very early along the folding pathway, being imprinted in the residual structure of the denatured state.

A thermodynamic approach to the mechanism of cell-penetrating peptides in model membranes

We report a first test of the hypothesis that the mechanism of antimicrobial, cytolytic, and amphipathic cell-penetrating peptides in model membranes is determined by the thermodynamics of insertion of the peptide into the lipid bilayer from the surface-associated state. Three peptides were designed with minimal mutations relative to the sequence of TP10W, the Y3W variant of transportan 10, which is a helical, amphipathic cell-penetrating peptide previously studied. Binding to 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) membranes and release of dye from those vesicles were assessed by stopped-flow fluorescence, and the secondary structure of the peptides on the membrane was determined by circular dichroism. The Gibbs energy of binding determined experimentally was in excellent agreement with that calculated using the Wimley-White interfacial hydrophobicity scale, taking into account the helical content of the membrane-associated peptide. Release of dye from POPC vesicles remained graded, as predicted by the hypothesis. More significantly, as the Gibbs energy of insertion into the bilayer became more unfavorable, which was estimated using the Wimley-White octanol hydrophobicity scale, dye release became slower, in quantitative agreement with the prediction.

Structural and antimicrobial properties of human pre-elafin/trappin-2 and derived peptides against Pseudomonas aeruginosa

Background

Pre-elafin/trappin-2 is a human innate defense molecule initially described as a potent inhibitor of neutrophil elastase. The full-length protein as well as the N-terminal "cementoin" and C-terminal "elafin" domains were also shown to possess broad antimicrobial activity, namely against the opportunistic pathogen P. aeruginosa. The mode of action of these peptides has, however, yet to be fully elucidated. Both domains of pre-elafin/trappin-2 are polycationic, but only the structure of the elafin domain is currently known. The aim of the present study was to determine the secondary structures of the cementoin domain and to characterize the antibacterial properties of these peptides against P. aeruginosa.

Results

We show here that the cementoin domain adopts an α-helical conformation both by circular dichroism and nuclear magnetic resonance analyses in the presence of membrane mimetics, a characteristic shared with a large number of linear polycationic antimicrobial peptides. However, pre-elafin/trappin-2 and its domains display only weak lytic properties, as assessed by scanning electron micrography, outer and inner membrane depolarization studies with P. aeruginosa and leakage of liposome-entrapped calcein. Confocal microscopy of fluorescein-labeled pre-elafin/trappin-2 suggests that this protein possesses the ability to translocate across membranes. This correlates with the finding that pre-elafin/trappin-2 and elafin bind to DNA in vitro and attenuate the expression of some P. aeruginosa virulence factors, namely the biofilm formation and the secretion of pyoverdine.

Conclusions

The N-terminal cementoin domain adopts α-helical secondary structures in a membrane mimetic environment, which is common in antimicrobial peptides. However, unlike numerous linear polycationic antimicrobial peptides, membrane disruption does not appear to be the main function of either cementoin, elafin or full-length pre-elafin/trappin-2 against P. aeruginosa. Our results rather suggest that pre-elafin/trappin-2 and elafin, but not cementoin, possess the ability to modulate the expression of some P.aeruginosa virulence factors, possibly through acting on intracellular targets.

Unique features of the folding landscape of a repeat protein revealed by pressure perturbation

The volumetric properties of proteins yield information about the changes in packing and hydration between various states along the folding reaction coordinate and are also intimately linked to the energetics and dynamics of these conformations. These volumetric characteristics can be accessed via pressure perturbation methods. In this work, we report high-pressure unfolding studies of the ankyrin domain of the Notch receptor (Nank1-7) using fluorescence, small-angle x-ray scattering, and Fourier transform infrared spectroscopy. Both equilibrium and pressure-jump kinetic fluorescence experiments were consistent with a simple two-state folding/unfolding transition under pressure, with a rather small volume change for unfolding compared to proteins of similar molecular weight. High-pressure fluorescence, Fourier transform infrared spectroscopy, and small-angle x-ray scattering measurements revealed that increasing urea over a very small range leads to a more expanded pressure unfolded state with a significant decrease in helical content. These observations underscore the conformational diversity of the unfolded-state basin. The temperature dependence of pressure-jump fluorescence relaxation measurements demonstrated that at low temperatures, the folding transition state ensemble (TSE) lies close in volume to the folded state, consistent with significant dehydration at the barrier. In contrast, the thermal expansivity of the TSE was found to be equivalent to that of the unfolded state, indicating that the interactions that constrain the folded-state thermal expansivity have not been established at the folding barrier. This behavior reveals a high degree of plasticity of the TSE of Nank1-7.

What can solid state NMR contribute to our understanding of protein folding?

Complete understanding of the folding process that connects a structurally disordered state of a protein to an ordered, biochemically functional state requires detailed characterization of intermediate structural states with high resolution and site specificity. While the intrinsically inhomogeneous and dynamic nature of unfolded and partially folded states limits the efficacy of traditional X-ray diffraction and solution NMR in structural studies, solid state NMR methods applied to frozen solutions can circumvent the complications due to molecular motions and conformational exchange encountered in unfolded and partially folded states. Moreover, solid state NMR methods can provide both qualitative and quantitative structural information at the site-specific level, even in the presence of structural inhomogeneity. This article reviews relevant solid state NMR methods and their initial applications to protein folding studies. Using either chemical denaturation to prepare unfolded states at equilibrium or a rapid freezing apparatus to trap non-equilibrium, transient structural states on a sub-millisecond time scale, recent results demonstrate that solid state NMR can contribute essential information about folding processes that is not available from more familiar biophysical methods.

Intrinsically disordered regions may lower the hydration free energy in proteins: a case study of nudix hydrolase in the bacterium Deinococcus radiodurans

The proteome of the radiation- and desiccation-resistant bacterium D. radiodurans features a group of proteins that contain significant intrinsically disordered regions that are not present in non-extremophile homologues. Interestingly, this group includes a number of housekeeping and repair proteins such as DNA polymerase III, nudix hydrolase and rotamase. Here, we focus on a member of the nudix hydrolase family from D. radiodurans possessing low-complexity N- and C-terminal tails, which exhibit sequence signatures of intrinsic disorder and have unknown function. The enzyme catalyzes the hydrolysis of oxidatively damaged and mutagenic nucleotides, and it is thought to play an important role in D. radiodurans during the recovery phase after exposure to ionizing radiation or desiccation. We use molecular dynamics simulations to study the dynamics of the protein, and study its hydration free energy using the GB/SA formalism. We show that the presence of disordered tails significantly decreases the hydration free energy of the whole protein. We hypothesize that the tails increase the chances of the protein to be located in the remaining water patches in the desiccated cell, where it is protected from the desiccation effects and can function normally. We extrapolate this to other intrinsically disordered regions in proteins, and propose a novel function for them: intrinsically disordered regions increase the "surface-properties" of the folded domains they are attached to, making them on the whole more hydrophilic and potentially influencing, in this way, their localization and cellular activity.

Prion fibrillization is mediated by a native structural element that comprises helices H2 and H3

Aggregation and misfolding of the prion protein (PrP) are thought to be the cause of a family of lethal neurodegenerative diseases affecting humans and other animals. Although the structures of PrP from several species have been solved, still little is known about the mechanisms that lead to the misfolded species. Here, we show that the region of PrP comprising the hairpin formed by the helices H2 and H3 is a stable independently folded unit able to retain its secondary and tertiary structure also in the absence of the rest of the sequence. We also prove that the isolated H2H3 is highly fibrillogenic and forms amyloid fibers morphologically similar to those obtained for the full-length protein. Fibrillization of H2H3 but not of full-length PrP is concomitant with formation of aggregates. These observations suggest a "banana-peeling" mechanism for misfolding of PrP in which H2H3 is the aggregation seed that needs to be first exposed to promote conversion from a helical to a beta-rich structure.

Mutational studies uncover non-native structure in the dimeric kinetic intermediate of the H2A-H2B heterodimer

The folding pathway of the histone H2A-H2B heterodimer minimally includes an on-pathway, dimeric, burst-phase intermediate, I(2). The partially folded H2A and H2B monomers populated at equilibrium were characterized as potential monomeric kinetic intermediates. Folding kinetics were compared for initiation from isolated, folded monomers and the heterodimer unfolded in 4 M urea. The observed rates were virtually identical above 0.4 M urea, exhibiting a log-linear relationship on the final denaturant concentration. Below approximately 0.4 M urea (concentrations inaccessible from the 4-M urea unfolded state), a rollover in the rates was observed; this suggests that a component of the I(2) ensemble contains non-native structure that rearranges/isomerizes to a more native-like species. The contribution of helix propensity to the stability of the I(2) ensemble was assessed with a set of H2A-H2B mutants containing Ala and Gly replacements at nine sites, focusing mainly on the long, central alpha2 helix. Equilibrium and kinetic folding/unfolding data were collected to determine the effects of the mutations on the stability of I(2) and the transition state between I(2) and N(2). This limited mutational study indicated that residues in the alpha2 helices of H2A and H2B as well as alpha1 of H2B and both the C-terminus of alpha3 and the short alphaC helix of H2A contribute to the stability of the I(2) burst-phase species. Interestingly, at least eight of the nine targeted residues stabilize I(2) by interactions that are non-native to some extent. Given that destabilizing I(2) and these non-native interactions does not accelerate folding, it is concluded that the native and non-native structures present in the I(2) ensemble enable efficient folding of H2A-H2B.

Designed proteins to modulate cellular networks

A major challenge of protein design is to create useful new proteins that interact specifically with biological targets in living cells. Such binding modules have many potential applications, including the targeted perturbation of protein networks. As a general approach to create such modules, we designed a library with approximately 10(9) different binding specificities based on a small 3-tetratricopeptide repeat (TPR) motif framework. We employed a novel strategy, based on split GFP reassembly, to screen the library for modules with the desired binding specificity. Using this approach, we identified modules that bind tightly and specifically to Dss1, a small human protein that interacts with the tumor suppressor protein BRCA2. We showed that these modules also bind the yeast homologue of Dss1, Sem1. Furthermore, we demonstrated that these modules inhibit Sem1 activity in yeast. This strategy will be generally applicable to make novel genetically encoded tools for systems/synthetic biology applications.

Extensive and modular intrinsically disordered segments in C. elegans TTN-1 and implications in filament binding, elasticity and oblique striation

TTN-1, a titin like protein in Caenorhabditis elegans, is encoded by a single gene and consists of multiple Ig and fibronectin 3 domains, a protein kinase domain and several regions containing tandem short repeat sequences. We have characterized TTN-1's sarcomere distribution, protein interaction with key myofibrillar proteins as well as the conformation malleability of representative motifs of five classes of short repeats. We report that two antibodies developed to portions of TTN-1 detect an approximately 2-MDa polypeptide on Western blots. In addition, by immunofluorescence staining, both of these antibodies localize to the I-band and may extend into the outer edge of the A-band in the obliquely striated muscle of the nematode. Six different 300-residue segments of TTN-1 were shown to variously interact with actin and/or myosin in vitro. Conformations of synthetic peptides of representative copies of each of the five classes of repeats--39-mer PEVT, 51-mer CEEEI, 42-mer AAPLE, 32-mer BLUE and 30-mer DispRep--were investigated by circular dichroism at different temperatures, ionic strengths and solvent polarities. The PEVT, CEEEI, DispRep and AAPLE peptides display a combination of a polyproline II helix and an unordered structure in aqueous solution and convert in trifluoroethanol to alpha-helix (PEVT, CEEEI, DispRep) and beta-turn (AAPLE) structures, respectively. The octads in BLUE motifs form unstable alpha-helix-like structures coils in aqueous solution and negligible heptad-based, alpha-helical coiled-coils. The alpha-helical structure, as modeled by threading and molecular dynamics simulations, tends to form helical bundles and crosses based on its 8-4-2-2 hydrophobic helical patterns and charge arrays on its surface. Our finding indicates that APPLE, PEVT, CEEEI and DispRep regions are all intrinsically disordered and highly reminiscent of the conformational malleability and elasticity of vertebrate titin PEVK segments. The proposed presence of long, modular and unstable alpha-helical oligomerization domains in the BLUE region of TTN-1 could bundle TTN-1 and stabilize oblique striation of the sarcomere.

α-Helix peptides designed from EBV-gH protein display higher antigenicity and induction of monocyte apoptosis than the native peptide

We tested the hypothesis that stabilizing α-helix of Epstein-Barr virus gH-derived peptide 11438 used for binding human cells will increase its biological activity. Non-stable α-helix of peptide 11438 was unfolded in an entropy-driven process, despite the opposing effect of the enthalpy factor. Adding and/or changing amino acids in peptide 11438 allowed the designing of peptides 33207, 33208 and 33210; peptides 33208 and 33210 displayed higher helical content due to a decreased unfolding entropy change as was determined by AGADIR, molecular dynamics and circular dichroism analysis. Peptides 33207, 33208 and 33210 inhibited EBV invasion of peripheral blood mononuclear cells and displayed epitopes more similar to native protein than peptide 11438; these peptides could be useful for detecting antibodies induced by native gH protein since they displayed high reactivity with anti-EBV antibodies. Anti-peptide 33207 antibodies showed higher reactivity with EBV than anti-peptide 11438 antibodies being useful for inducing antibodies against EBV. Anti-peptide 33210 antibodies inhibit EBV invasion of epithelial cells better than anti-peptide 11438 antibodies. Peptide 33210 bound to normal T lymphocytes and Raji cells stronger than peptide 11438 and also induced apoptosis of monocytes and Raji cells but not of normal T cells in a similar way to EBV-gH. Peptide 33210 inhibited the monocytes' development toward dendritic cells better than EBV and peptide 11438. In conclusion, stabilizing the α-helix in peptides 33208 and 33210 designed from peptide 11438 increased the antigenicity and the ability of the antibodies induced by peptides of inhibiting EBV invasion of host cells.

Structural and pharmacological characteristics of chimeric peptides derived from peptide E and beta-endorphin reveal the crucial role of the C-terminal YGGFL and YKKGE motifs in their analgesic properties

Peptide E (a 25-amino acid peptide derived from proenkephalin A) and beta-endorphin (a 31-amino acid peptide derived from proopiomelanocortin) bind with high affinity to opioid receptors and share structural similarities but induce analgesic effects of very different intensity. Indeed, whereas they possess the same N-terminus Met-enkephalin message sequence linked to a helix by a flexible spacer and a C-terminal part in random coil conformation, in contrast with peptide E, beta-endorphin produces a profound analgesia. To determine the key structural elements explaining this very divergent opioid activity, we have compared the structural and pharmacological characteristics of several chimeric peptides derived from peptide E and beta-endorphin. Structures were obtained under the same experimental conditions using circular dichroism, computational estimation of helical content and/or nuclear magnetic resonance spectroscopy (NMR) and NMR-restrained molecular modeling. The hot-plate and writhing tests were used in mice to evaluate the antinociceptive effects of the peptides. Our results indicate that neither the length nor the physicochemical profile of the spacer plays a fundamental role in analgesia. On the other hand, while the functional importance of the helix cannot be excluded, the last 5 residues in the C-terminal part seem to be crucial for the expression or absence of the analgesic activity of these peptides. These data raise the question of the true function of peptides E in opioidergic systems.

Residue-specific description of non-native transient structures in the ensemble of acid-denatured structures of the all-beta protein c-src SH3

Secondary chemical shift analysis has been used to characterize the unfolded state of acid-denatured c-src SH3. Even though native c-src SH3 adopts an all-beta fold, we found evidence of transient helicity in regions corresponding to native loops. In particular, residues 40-46, connecting the n-src loop to the third beta-strand, exhibited an apparent helicity of nearly 45%. Furthermore, the RT loop and the diverging turn appeared to adopt non-native-like helical conformations. Interestingly, none of the residues found in transient helical conformations exhibited significant varphi-values [Riddle, D. S., et al. (1999) Nat. Struct. Biol. 6, 1016-1024]. This indicated that the transient helicity has no influence or only a weak influence on the actual protein folding reaction. The residual structural propensities were compared to those of other SH3 domains, revealing heterogeneity in the unfolded ensemble that clearly contrasts with the conserved character of the topology of native state and transition state ensembles typical for SH3 domains.

An Achilles' heel in an amyloidogenic protein and its repair: insulin fibrillation and therapeutic design

Insulin fibrillation provides a model for a broad class of amyloidogenic diseases. Conformational distortion of the native monomer leads to aggregation-coupled misfolding. Whereas beta-cells are protected from proteotoxicity by hexamer assembly, fibrillation limits the storage and use of insulin at elevated temperatures. Here, we have investigated conformational distortions of an engineered insulin monomer in relation to the structure of an insulin fibril. Anomalous (13)C NMR chemical shifts and rapid (15)N-detected (1)H-(2)H amide-proton exchange were observed in one of the three classical alpha-helices (residues A1-A8) of the hormone, suggesting a conformational equilibrium between locally folded and unfolded A-chain segments. Whereas hexamer assembly resolves these anomalies in accordance with its protective role, solid-state (13)C NMR studies suggest that the A-chain segment participates in a fibril-specific beta-sheet. Accordingly, we investigated whether helicogenic substitutions in the A1-A8 segment might delay fibrillation. Simultaneous substitution of three beta-branched residues (Ile(A2) --> Leu, Val(A3) --> Leu, and Thr(A8) --> His) yielded an analog with reduced thermodynamic stability but marked resistance to fibrillation. Whereas amide-proton exchange in the A1-A8 segment remained rapid, (13)Calpha chemical shifts exhibited a more helical pattern. This analog is essentially without activity, however, as Ile(A2) and Val(A3) define conserved receptor contacts. To obtain active analogs, substitutions were restricted to A8. These analogs exhibit high receptor-binding affinity; representative potency in a rodent model of diabetes mellitus was similar to wild-type insulin. Although (13)Calpha chemical shifts remain anomalous, significant protection from fibrillation is retained. Together, our studies define an "Achilles' heel" in a globular protein whose repair may enhance the stability of pharmaceutical formulations and broaden their therapeutic deployment in the developing world.

Low-level expression of a folding-incompetent protein in Escherichia coli: search for the molecular determinants of protein aggregation in vivo

Aggregation of peptides and proteins into insoluble amyloid fibrils or related intracellular inclusions is the hallmark of many degenerative diseases, including Alzheimer's disease, Parkinson's disease, and various forms of amyloidosis. In spite of the considerable progress carried out in vitro in elucidating the molecular determinants of the conversion of purified and isolated proteins into amyloid fibrils, very little is known on factors governing this process in the complex environment of living organisms. Taking advantage of increasing evidence that bacterial inclusion bodies consist of amyloid-like aggregates, we have expressed in Escherichia coli both wild type and 21 single-point mutants of the N-terminal domain of the E. coli protein HypF. All variants were expressed as folding-incompetent units in a controlled manner, at low and comparable levels. Their solubilities were measured by quantifying the protein amount contained in the soluble and insoluble fractions by Western blot analysis. A significant negative correlation was found between the solubility of the variants in E. coli and their intrinsic propensity to form amyloid fibrils, predicted using an algorithm previously validated experimentally in vitro on a number of unfolded peptides and proteins, and considering hydrophobicity, beta-sheet propensity, and charge as major sequence determinants of the aggregation process. These findings show that the physicochemical parameters previously recognized to govern amyloid formation by fully or partially unfolded proteins are largely applicable in vivo and pave the way for the molecular exploration of a process as complex as protein aggregation in living organisms.

Microscopic formulation of the Zimm-Bragg model for the helix-coil transition

A microscopic spin model is proposed for the phenomenological Zimm-Bragg model for the helix-coil transition in biopolymers. This model is shown to provide the same thermophysical properties of the original Zimm-Bragg model and it allows a very convenient framework to compute statistical quantities. Physical origins of this spin model are made transparent by an exact mapping into a one-dimensional Ising model with an external field. However, the dependence on temperature of the reduced external field turns out to differ from the standard one-dimensional Ising model and hence it gives rise to different thermophysical properties, despite the exact mapping connecting them. We discuss how this point has been frequently overlooked in the recent literature.

Factors that influence helical preferences for singly charged gas-phase peptide ions: the effects of multiple potential charge-carrying sites

Ion mobility-mass spectrometry is used to investigate the structure(s) of a series of model peptide [M + H](+) ions to better understand how intrinsic properties affect structure in low dielectric environments. The influence of peptide length, amino acid sequence, and composition on gas-phase structure is examined for a series of model peptides that have been previously studied in solution. Collision cross sections for the [M + H](+) ions of Ac-(AAKAA)(n)Y-NH(2) (n = 3-6) and Ac-Y(AEAAKA)(n)F-NH(2) (n = 2-5) are reported and correlated with candidate structures generated using molecular modeling techniques. The [M + H](+) ions of the AAKAA peptide series each exhibit a single, dominant ion mobility arrival time distribution (ATD) which correlates to partial helical structures, whereas the [M + H](+) ions of the AEAAKA ion series are composed of ATDs which correlate to charge-solvated globules (i.e., the charge is coordinated or solvated by polar peptide functional groups). These data raise numerous questions concerning intrinsic properties (amino acid sequence and composition as well as charge location) that dictate gas-phase peptide ion structure, which may reflect trends for peptide ion structure in low dielectric environments, such as transmembrane segments.

Context-independent, temperature-dependent helical propensities for amino acid residues

Assigned from data sets measured in water at 2, 25, and 60 degrees C containing (13)C=O NMR chemical shifts and [theta](222) ellipticities, helical propensities are reported for the 20 genetically coded amino acids, as well as for norvaline and norleucine. These have been introduced by chemical synthesis at central sites within length-optimized, spaced, solubilized Ala(19) hosts. The resulting polyalanine-derived, quantitative propensity sets express for each residue its temperature-dependent but context-independent tendency to forego a coil state and join a preexisting helical conformation. At 2 degrees C their rank ordering is: P << G < H < C, T, N < S < Y, F, W < V, D < K < Q < I < R, M < L < E < A; at 60 degrees C the rank becomes: H, P < G < C < R, K < T, Y, F < N, V < S < Q < W, D < I, M < E < A < L. The DeltaDeltaG values, kcal/mol, relative to alanine, for the cluster T, N, S, Y, F, W, V, D, Q, imply that at 2 degrees C all are strong breakers: DeltaDeltaG(mean) = +0.63 +/- 0.11, but at 60 degrees C their breaking tendencies are dramatically attenuated and converge toward the mean: DeltaDeltaG(mean) = +0.25 +/- 0.07. Accurate modeling of helix-rich proteins found in thermophiles, mesophiles, and organisms that flourish near 0 degrees C thus requires appropriately matched propensity sets. Comparisons are offered between the temperature-dependent propensity assignments of this study and those previously assigned by the Scheraga group; the special problems that attend propensity assignments for charged residues are illustrated by lysine guest data; and comparisons of errors in helicity assignments from shifts and ellipticity data show that the former provide superior precision and accuracy.

Beta2-strand of salivary S cystatins: a "chemeleon sequence"

Secondary structure prediction of salivary cystatins S, SA, and SN carried out by several methods label the 39-58 sequence (beta2-strand) as predominantly alpha-helical. The helical propensity of a peptide corresponding to beta2-strand of salivary SA cystatin analyzed by CD display high helical propensity in aqueous solution, whereas peptides matching the beta2-strand amino acid sequence of cystatins S and SN, display random coil conformation in aqueous solution but acquire alpha-helical conformation in the presence of trifluoroethanol (TFE). Moreover molecular dynamics simulation performed on the homology modeling of cystatin SA constructed on the basis of recently determined three-dimensional structure of salivary cystatin D, suggests that cystatin SA does not significantly deviate from the starting structure over the course of the simulation. The results obtained indicate that the beta2-strand of salivary S cystatins has high helical propensity when isolated from native protein and acquire the final beta structure by interaction with the rest of the polypeptide chain.

Conformational stability of neuroglobin helix F--possible effects on the folding pathway within the globin family

Neuroglobin is a recently discovered member of the globin family, mainly observed in neurons and retina. Despite the low sequence identity (less than 20% over the whole sequence for the human proteins), the general fold of neuroglobin closely resembles that of myoglobin. The latter is a paradigmatic protein for folding studies, whereas much less is known about the neuroglobin folding pathway. In this work, we show how the structural features of helix F in neuroglobin and myoglobin could represent a pivotal difference in their folding pathways. Former studies widely documented that myoglobin lacks helix F in the apo form. In this study, limited proteolysis experiments on aponeuroglobin showed that helix F does not undergo proteolytic cleavage, suggesting that, also in the apo form, this helix maintains a rigid and structured conformation. To understand better the structural properties of helices F in the two proteins, we analyzed peptides encompassing helix F of neuroglobin and myoglobin in the wild-type and mutant forms. NMR and CD experiments revealed a helical conformation for neuroglobin helix F peptide, at both pH 7 and pH 2, absent in the myoglobin peptide. In particular, NMR data suggest a secondary structure stabilization effect caused by hydrophobic interactions involving Tyr88, Leu89 and Leu92. Molecular dynamics simulations performed on the apo and holo forms of the two proteins reveal the persistence of helix F in neuroglobin even in the absence of heme. Conversely myoglobin shows a higher mobility of the N-terminus of helix F on heme removal, which leads to the loss of secondary structure.

Amino acid insertion reveals a necessary three-helical intermediate in the folding pathway of the colicin E7 immunity protein Im7

The small (87-residue) α-helical protein Im7 (an inhibitor protein for colicin E7 that provides immunity to cells producing colicin E7) folds via a three-state mechanism involving an on-pathway intermediate. This kinetic intermediate contains three of four native helices that are oriented in a non-native manner so as to minimise exposed hydrophobic surface area at this point in folding. The short (6-residue) helix III has been shown to be unstructured in the intermediate ensemble and does not dock onto the developing hydrophobic core until after the rate-limiting transition state has been traversed. After helix III has docked, it adopts an α-helical secondary structure, and the side chains of residues within this region provide contacts that are crucial to native-state stability. In order to probe further the role of helix III in the folding mechanism of Im7, we created a variant that contains an eight-amino-acid polyalanine-like helix stabilised by a Glu-Arg salt bridge and an Asn-Pro-Gly capping motif, juxtaposed C-terminal to the natural 6-residue helix III. The effect of this insertion on the structure of the native protein and its folding mechanism were studied using NMR and ϕ-value analysis, respectively. The results reveal a robust native structure that is not perturbed by the presence of the extended helix III. Mutational analysis performed to probe the folding mechanism of the redesigned protein revealed a conserved mechanism involving the canonical three-helical intermediate. The results suggest that folding via a three-helical species stabilised by both native and non-native interactions is an essential feature of Im7 folding, independent of the helical propensity of helix III.

Helix stabilization precedes aqueous and bilayer-catalyzed fiber formation in islet amyloid polypeptide

Islet amyloid polypeptide (IAPP) is an unstructured polypeptide hormone that is cosecreted with insulin. In patients with type 2 diabetes, IAPP undergoes a transition from its natively disordered state to a highly ordered, all-beta-strand amyloid fiber. Although predominantly disordered, IAPP transiently samples alpha-helical structure in solution. IAPP adopts a fully helical structure when bound to membrane surfaces in a process associated with catalysis of amyloid formation. Here, we use spectroscopic techniques to study the structure of full-length, monomeric IAPP under amyloidogenic conditions. We observe that the residues with helical propensity in solution (1-22) also form the membrane-associated helix. Additionally, reduction of the N-terminal disulfide bond (Cys2-Cys7) decreases the extent of helix formed throughout this region. Through manipulation of sample conditions to increase or decrease the amount of helix, we show that the degree of helix formed affects the rate of amyloid assembly. Formation of helical structure is directly correlated with enhanced amyloid formation both on the membrane surface and in solution. These observations support suggested mechanisms in which parallel helix associations bring together regions of the peptide that could nucleate beta-strand structure. Remarkably, stabilization of non-amyloid structure appears to be a key intermediate in assembly of IAPP amyloid.

Optimized molecular dynamics force fields applied to the helix-coil transition of polypeptides

Obtaining the correct balance of secondary structure propensities is a central priority in protein force-field development. Given that current force fields differ significantly in their alpha-helical propensities, a correction to match experimental results would be highly desirable. We have determined simple backbone energy corrections for two force fields to reproduce the fraction of helix measured in short peptides at 300 K. As validation, we show that the optimized force fields produce results in excellent agreement with nuclear magnetic resonance experiments for folded proteins and short peptides not used in the optimization. However, despite the agreement at ambient conditions, the dependence of the helix content on temperature is too weak, a problem shared with other force fields. A fit of the Lifson-Roig helix-coil theory shows that both the enthalpy and entropy of helix formation are too small: the helix extension parameter w agrees well with experiment, but its entropic and enthalpic components are both only about half the respective experimental estimates. Our structural and thermodynamic analyses point toward the physical origins of these shortcomings in current force fields, and suggest ways to address them in future force-field development.

All-atom model for stabilization of alpha-helical structure in peptides by hydrocarbon staples

Recent work has shown that the incorporation of an all-hydrocarbon "staple" into peptides can greatly increase their alpha-helix propensity, leading to an improvement in pharmaceutical properties such as proteolytic stability, receptor affinity, and cell permeability. Stapled peptides thus show promise as a new class of drugs capable of accessing intractable targets such as those that engage in intracellular protein-protein interactions. The extent of alpha-helix stabilization provided by stapling has proven to be substantially context dependent, requiring cumbersome screening to identify the optimal site for staple incorporation. In certain cases, a staple encompassing one turn of the helix (attached at residues i and i+4) furnishes greater helix stabilization than one encompassing two turns (i,i+7 staple), which runs counter to expectation based on polymer theory. These findings highlight the need for a more thorough understanding of the forces that underlie helix stabilization by hydrocarbon staples. Here we report all-atom Monte Carlo folding simulations comparing unmodified peptides derived from RNase A and BID BH3 with various i,i+4 and i,i+7 stapled versions thereof. The results of these simulations were found to be in quantitative agreement with experimentally determined helix propensities. We also discovered that staples can stabilize quasi-stable decoy conformations, and that the removal of these states plays a major role in determining the helix stability of stapled peptides. Finally, we critically investigate why our method works, exposing the underlying physical forces that stabilize stapled peptides.

Different members of a simple three-helix bundle protein family have very different folding rate constants and fold by different mechanisms

The 15th, 16th, and 17th repeats of chicken brain alpha-spectrin (R15, R16, and R17, respectively) are very similar in terms of structure and stability. However, R15 folds and unfolds 3 orders of magnitude faster than R16 and R17. This is unexpected. The rate-limiting transition state for R15 folding is investigated using protein engineering methods (Phi-value analysis) and compared with previously completed analyses of R16 and R17. Characterisation of many mutants suggests that all three proteins have similar complexity in the folding landscape. The early rate-limiting transition states of the three domains are similar in terms of overall structure, but there are significant differences in the patterns of Phi-values. R15 apparently folds via a nucleation-condensation mechanism, which involves concomitant folding and packing of the A- and C-helices, establishing the correct topology. R16 and R17 fold via a more framework-like mechanism, which may impede the search to find the correct packing of the helices, providing a possible explanation for the fast folding of R15.

The folding mechanism of BBL: Plasticity of transition-state structure observed within an ultrafast folding protein family

Studies on members of protein families with similar structures but divergent sequences provide insights into the effects of sequence composition on the mechanism of folding. Members of the peripheral subunit-binding domain (PSBD) family fold ultrafast and approach the smallest size for cooperatively folding proteins. Phi-Value analysis of the PSBDs E3BD and POB reveals folding via nucleation-condensation through structurally very similar, polarized transition states. Here, we present a Phi-value analysis of the family member BBL and found that it also folds by a nucleation-condensation mechanism. The mean Phi values of BBL, E3BD, and POB were near identical, indicating similar fractions of non-covalent interactions being formed in the transition state. Despite the overall conservation of folding mechanism in this protein family, however, the pattern of Phi values determined for BBL revealed a larger dispersion of the folding nucleus across the entire structure, and the transition state was less polarized. The observed plasticity of transition-state structure can be rationalized by the different helix-forming propensities of PSBD sequences. The very strong helix propensity in the first helix of BBL, relative to E3BD and POB, appears to recruit more structure formation in that helix in the transition state at the expense of weaker interactions in the second helix. Differences in sequence composition can modulate transition-state structure of even the smallest natural protein domains.

Topology is the principal determinant in the folding of a complex all-alpha Greek key death domain from human FADD

In order to elucidate the relative importance of secondary structure and topology in determining folding mechanism, we have carried out a phi-value analysis of the death domain (DD) from human FADD. FADD DD is a 100 amino acid domain consisting of six anti-parallel alpha helices arranged in a Greek key structure. We asked how does the folding of this domain compare with that of (a) other all-alpha-helical proteins and (b) other Greek key proteins? Is the folding pathway determined mainly by secondary structure or is topology the principal determinant? Our Φ-value analysis reveals a striking resemblance to the all-beta Greek key immunoglobulin-like domains. Both fold via diffuse transition states and, importantly, long-range interactions between the four central elements of secondary structure are established in the transition state. The elements of secondary structure that are less tightly associated with the central core are less well packed in both cases. Topology appears to be the dominant factor in determining the pathway of folding in all Greek key domains.

The mechanism of folding of Im7 reveals competition between functional and kinetic evolutionary constraints

Many proteins reach their native state through pathways involving the presence of folding intermediates. It is not clear whether this type of folding landscape results from insufficient evolutionary pressure to optimize folding efficiency, or arises from a conflict between functional and folding constraints. Here, using protein-engineering, ultra-rapid mixing and stopped-flow experiments combined with restrained molecular dynamics simulations, we characterize the transition state for the formation of the intermediate populated during the folding of the bacterial immunity protein, Im7, and the subsequent molecular steps leading to the native state. The results provide a comprehensive view of the folding process of this small protein. An analysis of the contributions of native and non-native interactions at different stages of folding reveals how the complexity of the folding landscape arises from concomitant evolutionary pressures for function and folding efficiency.

Solution NMR studies of amphibian antimicrobial peptides: linking structure to function?

The high-resolution three-dimensional structure of an antimicrobial peptide has implications for the mechanism of its antimicrobial activity, as the conformation of the peptide provides insights into the intermolecular interactions that govern the binding to its biological target. For many cationic antimicrobial peptides the negatively charged membranes surrounding the bacterial cell appear to be a main target. In contrast to what has been found for other classes of antimicrobial peptides, solution NMR studies have revealed that in spite of the wide diversity in the amino acid sequences of amphibian antimicrobial peptides (AAMPs), they all adopt amphipathic alpha-helical structures in the presence of membrane-mimetic micelles, bicelles or organic solvent mixtures. In some cases the amphipathic AAMP structures are directly membrane-perturbing (e.g. magainin, aurein and the rana-box peptides), in other instances the peptide spontaneously passes through the membrane and acts on intracellular targets (e.g. buforin). Armed with a high-resolution structure, it is possible to relate the peptide structure to other relevant biophysical and biological data to elucidate a mechanism of action. While many linear AAMPs have significant antimicrobial activity of their own, mixtures of peptides sometimes have vastly improved antibiotic effects. Thus, synergy among antimicrobial peptides is an avenue of research that has recently attracted considerable attention. While synergistic relationships between AAMPs are well described, it is becoming increasingly evident that analyzing the intermolecular interactions between these peptides will be essential for understanding the increased antimicrobial effect. NMR structure determination of hybrid peptides composed of known antimicrobial peptides can shed light on these intricate synergistic relationships. In this work, we present the first NMR solution structure of a hybrid peptide composed of magainin 2 and PGLa bound to SDS and DPC micelles. The hybrid peptide adopts a largely helical conformation and some information regarding the inter-helix organization of this molecule is reported. The solution structure of the micelle associated MG2-PGLa hybrid peptide highlights the importance of examining structural contributions to the synergistic relationships but it also demonstrates the limitations in the resolution of the currently used solution NMR techniques for probing such interactions. Future studies of antimicrobial peptide synergy will likely require stable isotope-labeling strategies, similar to those used in NMR studies of proteins.

An unusual helix turn helix motif in the catalytic core of HIV-1 integrase binds viral DNA and LEDGF

Background

Integrase (IN) of the type 1 human immunodeficiency virus (HIV-1) catalyzes the integration of viral DNA into host cellular DNA. We identified a bi-helix motif (residues 149–186) in the crystal structure of the catalytic core (CC) of the IN-Phe185Lys variant that consists of the α₄ and α₅ helices connected by a 3 to 5-residue turn. The motif is embedded in a large array of interactions that stabilize the monomer and the dimer.

Principal Findings

We describe the conformational and binding properties of the corresponding synthetic peptide. This displays features of the protein motif structure thanks to the mutual intramolecular interactions of the α₄ and α₅ helices that maintain the fold. The main properties are the binding to: 1- the processing-attachment site at the LTR (long terminal repeat) ends of virus DNA with a K_d (dissociation constant) in the sub-micromolar range; 2- the whole IN enzyme; and 3- the IN binding domain (IBD) but not the IBD-Asp366Asn variant of LEDGF (lens epidermal derived growth factor) lacking the essential Asp366 residue. In our motif, in contrast to the conventional HTH (helix-turn-helix), it is the N terminal helix (α₄) which has the role of DNA recognition helix, while the C terminal helix (α₅) would rather contribute to the motif stabilization by interactions with the α₄ helix.

Conclusion

The motif, termed HTHi (i, for inverted) emerges as a central piece of the IN structure and function. It could therefore represent an attractive target in the search for inhibitors working at the DNA-IN, IN-IN and IN-LEDGF interfaces.

Mutational analysis of the stability of the H2A and H2B histone monomers

The eukaryotic histone heterodimer H2A-H2B folds through an obligatory dimeric intermediate that forms in a nearly diffusion-limited association reaction in the stopped-flow dead time. It is unclear whether there is partial folding of the isolated monomers before association. To address the possible contributions of structure in the monomers to the rapid association, we characterized H2A and H2B monomers in the absence of their heterodimeric partner. By far-UV circular dichroism, the H2A and H2B monomers are 15% and 31% helical, respectively--significantly less than observed in X-ray crystal structures. Acrylamide quenching of the intrinsic Tyr fluorescence was indicative of tertiary structure. The H2A and H2B monomers exhibit free energies of unfolding of 2.5 and 2.9 kcal mol(-1), respectively; at 10 microM, the sum of the stability of the monomers is approximately 60% of the stability of the native dimer. The helical content, stability, and m values indicate that H2B has a more stable, compact structure than H2A. The monomer m values are larger than expected for the extended histone fold motif, suggesting that the monomers adopt an overly collapsed structure. Stopped-flow refolding-initiated from urea-denatured monomers or the partially folded monomers populated at low denaturant concentrations-yielded essentially identical rates, indicating that monomer folding is productive in the rapid association and folding of the heterodimer. A series of Ala and Gly mutations were introduced into H2A and H2B to probe the importance of helix propensity on the structure and stability of the monomers. The mutational studies show that the central alpha-helix of the histone fold, which makes extensive intermonomer contacts, is structured in H2B but only partially folded in H2A.

Electrostatics in the ribosomal tunnel modulate chain elongation rates

Electrostatic potentials along the ribosomal exit tunnel are nonuniform and negative. The significance of electrostatics in the tunnel remains relatively uninvestigated, yet they are likely to play a role in translation and secondary folding of nascent peptides. To probe the role of nascent peptide charges in ribosome function, we used a molecular tape measure that was engineered to contain different numbers of charged amino acids localized to known regions of the tunnel and measured chain elongation rates. Positively charged arginine or lysine sequences produce transient arrest (pausing) before the nascent peptide is fully elongated. The rate of conversion from transiently arrested to full-length nascent peptide is faster for peptides containing neutral or negatively charged residues than for those containing positively charged residues. We provide experimental evidence that extraribosomal mechanisms do not account for this charge-specific pausing. We conclude that pausing is due to charge-specific interactions between the tunnel and the nascent peptide.

Structural insights for designed alanine-rich helices: comparing NMR helicity measures and conformational ensembles from molecular dynamics simulation

The temperature dependence of helical propensities for the peptides Ac-ZGG-(KAAAA)(3)X-NH(2) (Z = Y or G, X = A, K, and D-Arg) were studied both experimentally and by MD simulations. Good agreement is observed in both the absolute helical propensities as well as relative helical content along the sequence; the global minimum on the calculated free energy landscape corresponds to a single alpha-helical conformation running from K4 to A18 with some terminal fraying, particularly at the C-terminus. Energy component analysis shows that the single helix state has favorable intramolecular electrostatic energy due to hydrogen bonds, and that less-favorable two-helix globular states have favorable solvation energy. The central lysine residues do not appear to increase helicity; however, both experimental and simulation studies show increasing helicity in the series X = Ala --> Lys --> D-Arg. This C-capping preference was also experimentally confirmed in Ac-(KAAAA)(3)X-GY-NH(2) and (KAAAA)(3)X-GY-NH(2) sequences. The roles of the C-capping groups, and of lysines throughout the sequence, in the MD-derived ensembles are analyzed in detail.

Helix mutations stabilize a late productive intermediate on the folding pathway of ubiquitin

We have investigated the relative placement of rate-limiting energy barriers and the role of productive or obstructive intermediates on the folding pathway of yeast wild-type ubiquitin ( wt-Ub) containing the F45W mutation. To manipulate the folding barriers, we have designed a family of mutants in which stabilizing substitutions have been introduced incrementally on the solvent-exposed surface of the main alpha-helix (residues 23-34), which has a low intrinsic helical propensity in the native sequence. Although the U --> I and I --> N transitions are not clearly delineated in the kinetics of wt-Ub, we show that an intermediate becomes highly populated and more clearly resolved as the predicted stability of the helix increases. The observed acceleration in the rate of folding correlates with helix stability and is consistent with the I-state representing a productive rather than misfolded state. A Leffler analysis of the effects on kinetics of changes in stability within the family of helix mutants results in a biphasic correlation in both the refolding and unfolding rates that suggest a shift from a nucleation-condensation mechanism (weakly stabilized helix) toward a diffusion-collision model (highly stabilized helix). Through the introduction of helix-stabilizing mutations, we are able to engineer a well-resolved I-state on the folding pathway of ubiquitin which is likely to be structurally distinct from that which is only weakly populated on the folding pathway of wild-type ubiquitin.

Helix forming tendency of valine substituted poly-alanine: a molecular dynamics investigation

In this study, classical molecular dynamics simulations have been carried out on the valine (guest) substituted poly alanine (host) using the host-guest peptide approach to understand the role of valine in the formation and stabilization of helix. Valine has been substituted in the host peptide starting from N terminal to C terminal. Various structural parameters have been obtained from the molecular dynamics simulation to understand the tolerance of helical motif to valine. Depending on the position of valine in the host peptide, it stabilizes (or destabilizes) the formation of the helical structure. The substitution of valine in the poly alanine at some positions has no effect on the helix formation (deformation). It is interesting to observe the coexistence of 3 10 and alpha-helix in the peptides due to the dynamical nature of the hydrogen bonding interaction and sterical interactions.

Statistical estimation of statistical mechanical models: helix-coil theory and peptide helicity prediction

Analysis of biopolymer sequences and structures generally adopts one of two approaches: use of detailed biophysical theoretical models of the system with experimentally-determined parameters, or largely empirical statistical models obtained by extracting parameters from large datasets. In this work, we demonstrate a merger of these two approaches using Bayesian statistics. We adopt a common biophysical model for local protein folding and peptide configuration, the helix-coil model. The parameters of this model are estimated by statistical fitting to a large dataset, using prior distributions based on experimental data. L(1)-norm shrinkage priors are applied to induce sparsity among the estimated parameters, resulting in a significantly simplified model. Formal statistical procedures for evaluating support in the data for previously proposed model extensions are presented. We demonstrate the advantages of this approach including improved prediction accuracy and quantification of prediction uncertainty, and discuss opportunities for statistical design of experiments. Our approach yields a 39% improvement in mean-squared predictive error over the current best algorithm for this problem. In the process we also provide an efficient recursive algorithm for exact calculation of ensemble helicity including sidechain interactions, and derive an explicit relation between homo- and heteropolymer helix-coil theories and Markov chains and (non-standard) hidden Markov models respectively, which has not appeared in the literature previously.

Multiple roles for the C-terminal tail of the chemokine scavenger D6

D6 is a heptahelical receptor that suppresses inflammation and tumorigenesis by scavenging extracellular pro-inflammatory CC chemokines. Previous studies suggested this is dependent on constitutive trafficking of stable D6 protein to and from the cell surface via recycling endosomes. By internalizing chemokine each time it transits the cell surface, D6 can, over time, remove large quantities of these inflammatory mediators. We have investigated the role of the conserved 58-amino acid C terminus of human D6, which, unlike the rest of the protein, shows no clear homology to other heptahelical receptors. We show that, in human HEK293 cells, a serine cluster in this region controls the constitutive phosphorylation, high stability, and intracellular trafficking itinerary of the receptor and drives green fluorescent protein-tagged beta-arrestins to membranes at, and near, the cell surface. Unexpectedly, however, these properties, and the last 44 amino acids of the C terminus, are dispensable for D6 internalization and effective scavenging of the chemokine CCL3. Even in the absence of the last 58 amino acids, D6 still initially internalizes CCL3 but, surprisingly, exposure to ligand inhibits subsequent CCL3 uptake by this mutant. Progressive scavenging is therefore abrogated. We conclude that the heptahelical body of D6 on its own can engage the endocytotic machinery of HEK293 cells but that the C terminus is indispensable for scavenging because it prevents initial chemokine engagement of D6 from inhibiting subsequent chemokine uptake.