Intrinsic Helical Propensities Research Articles

Crystal Structures and Nuclear Magnetic Resonance Studies of the Apo Form of the c-MYC:MAX bHLHZip Complex Reveal a Helical Basic Region in the Absence of DNA.

The c-MYC transcription factor is a master regulator of cell growth and proliferation and is an established target for cancer therapy. This basic helix–loop–helix Zip protein forms a heterodimer with its obligatory partner MAX, which binds to DNA via the basic region. Considerable research efforts are focused on targeting the heterodimerization interface and the interaction of the complex with DNA. The only available crystal structure is that of a c-MYC:MAX complex artificially tethered by an engineered disulfide linker and prebound to DNA. We have carried out a detailed structural analysis of the apo form of the c-MYC:MAX complex, with no artificial linker, both in solution using nuclear magnetic resonance (NMR) spectroscopy and by X-ray crystallography. We have obtained crystal structures in three different crystal forms, with resolutions between 1.35 and 2.2 Å, that show extensive helical structure in the basic region. Determination of the α-helical propensity using NMR chemical shift analysis shows that the basic region of c-MYC and, to a lesser extent, that of MAX populate helical conformations. We have also assigned the NMR spectra of the c-MYC basic helix–loop–helix Zip motif in the absence of MAX and showed that the basic region has an intrinsic helical propensity even in the absence of its dimerization partner. The presence of helical structure in the basic regions in the absence of DNA suggests that the molecular recognition occurs via a conformational selection rather than an induced fit. Our work provides both insight into the mechanism of DNA binding and structural information to aid in the development of MYC inhibitors.

Impact of γ-Amino Acid Residue Preorganization on α/γ-Peptide Foldamer Helicity in Aqueous Solution.

α/γ-Peptide foldamers containing either γ(4)-amino acid residues or ring-constrained γ-amino acid residues have been reported to adopt 12-helical secondary structure in nonpolar solvents and in the solid state. These observations have engendered speculation that the seemingly flexible γ(4) residues have a high intrinsic helical propensity and that residue-based preorganization may not significantly stabilize the 12-helical conformation. However, the prior studies were conducted in environments that favor intramolecular H-bond formation. Here, we use 2D-NMR to compare the ability of γ(4) residues and cyclic γ residues to support 12-helix formation in more challenging environments, methanol and water. Both γ residue types support 12-helical folding in methanol, but only the cyclically constrained γ residues promote helicity in water. These results demonstrate the importance of residue-based preorganization strategies for achieving stable folding among short foldamers in aqueous solution.

Autoinhibition of ETV6 DNA Binding Is Established by the Stability of Its Inhibitory Helix

The ETS transcriptional repressor ETV6 (or TEL) is autoinhibited by an α-helix that sterically blocks its DNA-binding ETS domain. The inhibitory helix is marginally stable and unfolds when ETV6 binds to either specific or non-specific DNA. Using NMR spectroscopy, we show that folding of the inhibitory helix requires a buried charge–dipole interaction with helix H1 of the ETS domain. This interaction also contributes directly to autoinhibition by precluding a highly conserved dipole-enhanced hydrogen bond between the phosphodiester backbone of bound DNA and the N terminus of helix H1. To probe further the thermodynamic basis of autoinhibition, ETV6 variants were generated with amino acid substitutions introduced along the solvent exposed surface of the inhibitory helix. These changes were designed to increase the intrinsic helical propensity of the inhibitory helix without perturbing its packing interactions with the ETS domain. NMR-monitored amide hydrogen exchange measurements confirmed that the stability of the folded inhibitory helix increases progressively with added helix-promoting substitutions. This also results in progressively reinforced autoinhibition and decreased DNA-binding affinity. Surprisingly, locking the inhibitory helix onto the ETS domain by a disulfide bridge severely impairs, but does not abolish DNA binding. Weak interactions still occur via an interface displaced from the canonical ETS domain DNA-binding surface. Collectively, these studies establish a direct thermodynamic linkage between inhibitory helix stability and ETV6 autoinhibition, and demonstrate that helix unfolding does not strictly precede DNA binding. Modulating inhibitory helix stability provides a potential route for the in vivo regulation of ETV6 activity.

The N-Terminal Peptides of the Three Human Isoforms of the Mitochondrial Voltage-Dependent Anion Channel Have Different Helical Propensities.

The voltage-dependent anion channel (VDAC) is the main mitochondrial porin allowing the exchange of ions and metabolites between the cytosol and the mitochondrion. In addition, VDAC was found to actively interact with proteins playing a fundamental role in the regulation of apoptosis and being of central interest in cancer research. VDAC is a large transmembrane β-barrel channel, whose N-terminal helical fragment adheres to the channel interior, partially closing the pore. This fragment is considered to play a key role in protein stability and function as well as in the interaction with apoptosis-related proteins. Three VDAC isoforms are differently expressed in higher eukaryotes, for which distinct and complementary roles are proposed. In this work, the folding propensity of their N-terminal fragments has been compared. By using multiple spectroscopic techniques, and complementing the experimental results with theoretical computer-assisted approaches, we have characterized their conformational equilibrium. Significant differences were found in the intrinsic helical propensity of the three peptides, decreasing in the following order: hVDAC2 > hVDAC3 > hVDAC1. In light of the models proposed in the literature to explain voltage gating, selectivity, and permeability, as well as interactions with functionally related proteins, our results suggest that the different chemicophysical properties of the N-terminal domain are possibly correlated to different functions for the three isoforms. The overall emerging picture is that a similar transmembrane water accessible conduit has been equipped with not identical domains, whose differences can modulate the functional roles of the three VDAC isoforms.

Optimizing the Photocontrol of bZIP Coiled Coils with Azobenzene Crosslinkers: Role of the Crosslinking Site.

DNA binding by bZIP-type coiled-coil proteins can be inhibited by dominant negative versions of the proteins in which the N-terminal basic region is replaced by an acidic extension. Photocontrol of bZIP function can be achieved by introducing intramolecular azobenzene-based crosslinkers into dominant negatives. We show that the largest degree of photocontrol is achieved when the crosslinker is introduced into the zipper region of the dominant negative between Cys residues placed at f sites in the heptad segment showing the highest intrinsic helical propensity. The overall affinity of the dominant negative can then be tuned by varying the length of the acidic extension.

The intrinsic helical propensities of the helical fragments in prion protein under neutral and low pH conditions: a replica exchange molecular dynamics study

Replica exchange molecular dynamics simulations in neutral and acidic aqueous solutions were employed to study the intrinsic helical propensities of three helices in both Syrian hamster (syPrP) and human (huPrP) prion proteins. The helical propensities of syPrP HA and huPrP HA are very high under both pH conditions, which implies that HA is barely involved in the helix-to-β transition. The SyPrP HB chain has a strong tendency to adopt an extended conformation, which is possibly involved in the mechanism of infectious prion diseases in Syrian hamster. HuPrP HC has more of a preference for the extended conformation than huPrP HA and huPrP HB do, which leads to the conjecture that it is more likely to be the source of β-rich structure for human prion protein. We also noticed that the presence of salt bridges is not correlated with helical propensity, indicating that salt bridges do not stabilize helices.

The effect of azobenzene cross-linker position on the degree of helical peptide photo-control

Photo-switchable α-helical peptides offer a valuable tool to probe protein-biomacromolecule interactions in a spatiotemporally controlled manner. We synthesized a series of 32 residue peptides (AA32 s) with the core structure Ac-W-(E1A2A3A4R5)6-Q-NH2 and introduced the azobenzene based cross-linker BSBCA via reaction with Cys residues spaced at i, i + 7 intervals in different positions along the helix. UV/Vis measurements show that the composition of the photostationary state as well as thermal relaxation rates do not change considerably with changes in cross-linker position. CD analysis shows that photo-control of helix folding/unfolding is most effective when the cross-linker is targeted to the middle of the peptide so long as this segment has a high intrinsic helical propensity. The largest degree of photo-controlled conformational change occurred when a cross-linked central region of high intrinsic helicity was followed on the C-terminal side by a region of lower intrinsic helicity. This indicates the BSBCA cross-linker can act as a nucleation site for N-to-C propagation of a helix. These data help to guide the choice of cross-linking site in larger peptides and proteins where photo-control of conformation is desired.

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.

Folding of the C-Terminal Fragment V111-D143 of Staphylococcal Nuclease in Aqueous Solution

Studies of conformational features of fragments SNase(111-143) and SNase(118-143) and segment E122-K136 in 1-139 fragment (SNase139) suggest that the high intrinsic helical propensity can drive segment E122-K136 fold into a stable helix only when the segments V111-H121 and L137-D143 flanked on segment E122-K136 in staphylococcal nuclease (SNase) have stable folding.

Probing site-specific conformational distributions in protein folding with solid-state NMR

We demonstrate an experimental approach to structural studies of unfolded and partially folded proteins in which conformational distributions are probed at a site-specific level by 2D solid-state 13C NMR spectroscopy of glassy frozen solutions. Experiments on chemical denaturation of the 35-residue villin headpiece subdomain, a model three-helix-bundle protein with a known folded structure, reveal that 13C-labeled residues in the three helical segments of the folded state have markedly different conformational distributions in the unfolded state. Moreover, the 2D solid-state NMR line shapes near the unfolding midpoint do not fit a simple two-state model, in which the conformational distributions of the unfolded component are assumed to be independent of denaturant concentration. Comparison with solid-state NMR spectra of peptides containing the individual helical segments suggests an alternative two-step description of conformational distributions in partially folded states of the helical villin headpiece subdomain, in which chemical denaturation is viewed as a disruption of tertiary contacts followed by equilibration of local secondary structure according to the intrinsic helical propensities of individual segments.

Fast folding of a helical protein initiated by the collision of unstructured chains.

To examine whether helix formation necessarily precedes chain collision, we have measured the folding of a fully helical coiled coil that has been specially engineered to have negligible intrinsic helical propensity but high overall stability. The folding rate approaches the diffusion-limited value and is much faster than possible if folding is contingent on precollision helix formation. Therefore, the collision of two unstructured chains is the initial step of the dominant kinetic pathway, whereas helicity exerts its influence only at a later step. Folding from an unstructured encounter complex may be efficient and robust, which has implications for any biological process that couples folding to binding.

Assemblages of prion fragments: novel model systems for understanding amyloid toxicity

We report the conformational and toxic properties of two novel fibril-forming prion amyloid sequences, GAVVGGLG (PrP(119–126)) and VVGGLGG (PrP(121–127)). The conformational preferences of these fragments were studied in differing microenvironments of TFE/water mixtures and SDS solution. Interestingly, with an increase in TFE concentration, PrP(119–126) showed a helical conformational propensity, whereas PrP(121–127) adopted a more random coil structure. In 5% SDS, PrP(119–126) showed more α-helical content than in TFE solution, and PrP(121–127) exhibited a predominantly random coil conformation. However, both peptides took a random coil conformation in water, and over time the random coil transformed into a β-sheet structure with a significant percentage of helical conformation and β-turn structure in PrP(119–126) and PrP(121–127), respectively, as observed with CD spectroscopy. The aged fibrils of PrP(119–126) were insoluble in SDS, and PrP(121–127) was extractable with SDS solution. These fibrils were characterized by transmission electron microscopy. Both PrP(119–126) and PrP(121–127) formed stable monolayer's consisting of multimeric assemblages at the air–water interface. Monomeric PrP(119–126) was more toxic to astrocytes than the control Aβ peptide; however, the fibrillar form of PrP(119–126) was less toxic to astrocytes. PrP(121–127) elicited moderate toxicity in both soluble and fibrillar forms on astrocytes. Furthermore, quenching experiments using acroyl-labeled PrP(119–126) and PrP(121–127) with eosin-labeled synaptosomal membrane revealed that these prion fragments bind to anion-exchange protein. The binding of PrP(119–126) and PrP(121–127) with a membrane microdomain (lipid raft) was also analyzed using pyrenated derivatives. We conclude that the formation of PrP(119–126) and PrP(121–127) fibrils is a concentration-dependent process that involves coil to sheet conversion with aging. PrP(119–126), the sequence with intrinsic helical propensity, is more toxic in monomer form, and the fibril formation in this case seems to be protective to cells. For PrP(121–127), the SDS-soluble fibrils are more cytotoxic, indicating that a higher order assemblage structure is required for cytotoxic activity of this peptide.

Temperature-induced formation of a non-native intermediate state of the all beta-sheet protein CD2.

Domain 1 of the cell adhesion protein CD2 (CD2-1) has an all beta-structure typical of proteins belonging to the immunoglobin superfamily. It has a remarkable ability to fold as a native monomer or a metastable intertwined dimer. To understand the origin of structural rearrangements of CD2-1, we have studied equilibrium unfolding of the protein using various biophysical spectroscopic techniques. At temperatures above approx 68 degrees C, a partially folded state of CD2-1 (H state) with a distinct secondary structure, involving largely exposed aromatic and hydrophobic residues and a substantially perturbed tertiary structure, is observed. In contrast, an unfolded state (D state) of CD2-1 with random-coil-like secondary and tertiary structures is observed in 6 M GuHCl. This partially folded high-temperature state has increased negative molar ellipticity at 222 nm in far-ultraviolet CD spectra, implying formation of a non-native helical conformation. The existence of this non-native high-temperature intermediate is consistent with relatively high intrinsic helical propensities in the primary sequence of CD2-1. This conformational flexibility may be important in the observed domain swapping of CD2-1.

15N NMR relaxation as a probe for helical intrinsic propensity: the case of the unfolded D2 domain of annexin I.

The isolated D2 domain of annexin I is unable to adopt a tertiary fold but exhibits both native and non-native residual structures. It thus constitutes an attractive model for the investigation of dynamics of partially folded states in the context of protein folding and stability. 15N relaxation parameters of the D2 domain have been acquired at three different magnetic fields, 500, 600 and 800 MHz. This enables the estimation of the contribution of conformational exchange to the relaxation parameters on the micro- to millisecond time scale, thus providing a suitable data set for the description of motions on the pico- and nanosecond time scale. The analysis of the seven spectral densities obtained (J(0), J(50 MHz), J(60 MHz), J(80 MHz), <J(500 MHz)>, <J(600 MHz)>, <J(800 MHz)>) provides complementary and meaningful results on the conformational features of the D2 domain structure previously depicted by chemical shift and NOE data. Especially, residual helix segments exhibit distinct dynamical behaviors that are related to their intrinsic helical propensity. Beside the spectral density analysis, a series of models derived from the Lipari and Szabo model-free approach are investigated. Two models containing three parameters are able to reproduce equally well the experimental data within experimental errors but provide different values of order parameters and correlation times. The inability to find a unique model to describe the data emphasizes the difficulty to use and interpret the model-free parameters in the case of partially or fully unfolded proteins consisting of a wide range of interconverting conformers.

Intrinsic helical propensities and stable secondary structure in a membrane-bound fragment (S4) of the shaker potassium channel

The location and stability of helical secondary structure in a fragment comprising an extended sequence of the S4 transmembrane segment of the Shaker potassium channel was determined in methanol, and when bound to vesicles composed of egg phosphatidylcholine: egg phosphatidylglycerol (4:1; mol:mol) in water. The N-acetylated, C-amidated peptide corresponds to the sequence comprising residues A355-I384 in the Shaker potassium channel. Although NOEs characteristic of helical structure encompass essentially the full peptide sequence in methanol, analysis of amide and CHα chemical shifts, and amide exchange protection factors establish that stable helical structure comprises only around the first 22 amino acids of the 30 residue peptide. This sequence corresponds to that predicted to have the highest helical stability in water, indicating that while helical structure is considerably stabilised in methanol, the relative helical propensities of amino acids in methanol may be similar to those in water.In the presence of vesicles containing negatively charged lipids, helical structure corresponding to a maximum of around 40% of the extended S4 peptide is induced; no helical structure is induced in the presence of vesicles composed only of neutral lipids. The location of stable helical structure in the membrane-bound peptide was determined by amide hydrogen-deuterium exchange trapping, and was shown to encompass the sequence between residues near M2 and I18. This sequence is similar to that having high helix propensity in water and methanol, supporting the idea that intrinsic helical propensities are important in defining the location of stable helical structure in polypeptides bound in the interfacial region of lipid bilayers. The study defines an approach to determining the location of, and contributions to, the stability of helical secondary structure in membrane-reconstituted polypeptides.

Transition state heterogeneity in GCN4 coiled coil folding studied by using multisite mutations and crosslinking.

We have investigated the folding behavior of dimeric and covalently crosslinked versions of the 33-residue alpha-helical GCN4-p1 coiled coil derived from the leucine zipper region of the transcriptional activator GCN4. The effects of multisite substitutions indicate that folding occurs along multiple routes with nucleation sites located throughout the protein. The similarity in activation energies of the different routes together with an analysis of intrinsic helical propensities indicate that minimal helix is present before a productive collision of the two chains. However, approximately one-third to one-half of the total helical structure is formed in the postcollision transition state ensemble. For the crosslinked, monomeric version, folding occurs along a single robust pathway. Here, the region nearest the crosslink, with the least helical propensity, is structured in the transition state whereas the region farthest from the tether, with the most propensity, is completely unstructured. Hence, the existence of transition state heterogeneity and the selection of folding routes critically depend on chain topology.

Position dependence of amino acid intrinsic helical propensities II: non-charged polar residues: Ser, Thr, Asn, and Gln.

The assumption that the intrinsic alpha-helical propensities of the amino acids are position independent was critical in several helix/coil transition theories. In the first paper of these series, we reported that this is not the case for Gly and nonpolar aliphatic amino acids (Val, Leu, Met, and Ile). Here we have analyzed the helical intrinsic propensities of noncharged polar residues (Ser, Thr, Asn, and Gln) at different positions of a model polyalanine-based peptide. We found that Thr is more favorable (by approximately 0.3 kcal/mol) at positions N1 and N2 than in the helix center, although for Ser, Asn, and Gln the differences are smaller (+/-0.2 kcal/mol), and in many cases within the experimental error. There is a reasonable agreement (+/-0.2 kcal/mol) between the calculated free energies, using the ECEPP/2 force field equipped with a hydration potential, and the experimental data, except at position N1.

Elucidating the folding problem of alpha-helices: local motifs, long-range electrostatics, ionic-strength dependence and prediction of NMR parameters.

The information about the conformational behavior of monomeric helical peptides in solution, as well as the alpha-helix stability in proteins, has been previously utilized to derive a database with the energy contributions for various interactions taking place in an alpha-helix: intrinsic helical propensities, side-chain-side-chain interactions, main-chain-main-chain hydrogen bonds, and capping effects. This database was implemented in an algorithm based on the helix/coil transition theory (AGADIR). Here, we have modified this algorithm to include previously described local motifs: hydrophobic staple, Schellman motif and Pro-capping motif, new variants of these, and newly described side-chain-side-chain interactions. Based on recent experimental data we have introduced a position dependence of the helical propensities for some of the 20 amino acid residues. A new electrostatic model that takes into consideration all electrostatic interactions up to 12 residues in distance in the helix and random-coil conformations, as well as the effect of ionic strength, has been implemented. We have synthesized and analyzed several peptides, and used data from peptides already analysed by other groups, to test the validity of our electrostatic model. The modified algorithm predicts, with an overall standard deviation value of 6.6 (maximum helix is 100%), the helical, content of 778 peptides of which 223 correspond to wild-type and modified protein fragments. To improve the prediction potential of the algorithm and to have a direct comparison with nuclear magnetic resonance data, the algorithm now predicts the conformational shift of the CalphaH protons, 13Calpha and 3JalphaN values. We have found that for those peptides correctly predicted from the point of view of circular dichroism, the prediction of the NMR parameters is very good.

Position dependence of non-polar amino acid intrinsic helical propensities

Until now and based on the success of the helix/coil transition theory it has been assumed that the α-helical propensities of the amino acids are position independent. This has been critical to derive the set of theoretical parameters for the 20 natural amino acids. Here, we have analyzed the behavior of several non-polar residues, Val, Ile, Leu, Met and Gly at the N-cap, at each position of the first helical turn and at a central helical position of a 16-residue peptide model system that starts with eight consecutive alanine residues. We have interpreted the results from these experiments with the model of the helix/coil transition (AGADIR), that indicates that the intrinsic helical propensity is position dependent. Gly, Val and Ile are more favorable at the first turn than in the middle of the α-helix, while for Leu and Met we observe the opposite behavior. The differences between the observed helical propensities are as large as 1.0 kcal/mol in some cases. Molecular modeling calculations using the ECEPP/2 force-field equipped with a hydration potential show that this effect can be explained by the combination of three factors: (a) the side-chains in the first helix turn are more solvent-exposed; (b) they have fewer intramolecular van der Waals’ contacts; and (c) they posses higher configurational entropy than that in the central position of an α-helix. The position-dependent results of the calculations are in reasonable agreement with the experimental estimates and with the intrinsic propensities of the amino acids derived from the statistical analysis of the protein structure database.

Tolerance of a protein helix to multiple alanine and valine substitutions

Protein stability is influenced by the intrinsic secondary structure propensities of the amino acids and by tertiary interactions, but which of these factors dominates is not known in most cases. We have used combinatorial mutagenesis to examine the effects of substituting a good helix-forming residue (alanine) and a poor helix-forming residue (valine) at many positions in an alpha helix of a native protein. This has allowed us to average over many molecular environments and assess to what extent the results reflect intrinsic helical propensities or are masked by tertiary effects. Alanine or valine residues were combinatorially substituted at 12 positions in alpha-helix lambda repressor. Functional proteins were selected and sequenced to determine the degree to which each residue type was tolerated. On average, valine substitutions were accommodated slightly less well than alanine substitutions. On a positional basis, however, valine was tolerated as well as alanine at the majority of sites. In fact, alanine was preferred over valine statistically significantly only at four sites. Studies of mutant protein and peptide stabilities suggest that tertiary interactions mask the intrinsic secondary structure propensity differences at most of the remaining residue positions in this alpha helix. At the majority of positions in alpha-helix lambda repressor, tertiary interactions with other parts of the protein can be viewed as an environmental "buffer" that help to diminish the helix destabilizing effects of valine mutations and allow these mutations to be tolerated at frequencies similar to alanine mutations.