Positional dependence of the effects of negatively charged Glu side chains on the stability of two-stranded α-helical coiled-coils

The effects on protein stability of negatively charged Glu side chains at different positions along the length of the α-helix were investigated in the two-stranded α-helical coiled-coil. A native coiled-coil has been designed which consists of two identical 35 residue polypeptide chains with a heptad repeat Q_gV_aG_bA_cL_dQ_eK_f and a Cys residue at position 2 to allow the formation of an interchain 2-2′ disulphide bridge. This coiled-coil contains no intra- or interchain electrostatic interactions and served as a control for peptides in which Glu was substituted for Gln in the e or g heptad positions. The effect of the substitutions on stability was determined by urea denaturation at 20°C with the degree of unfolding monitored by circular dichroism spectroscopy. A Glu substituted for Gln near the N-terminus in each chain of the coiled-coil stabilizes the coiled-coil at pH 7, consistent with the charge–helix dipole interaction model. This stability increase is modulated by pH change and the addition of salt (KCl or guanidine hydrochloride), confirming the electrostatic nature of the effect. In contrast, Glu substitution in the middle of the helix destabilizes the coiled-coil because of the lower helical propensity and hydrophobicity of Glu compared with Gln at pH 7. Taking the intrinsic differences into account, the apparent charge–helix dipole interaction at the N-terminus is approximately 0.35 kcal/mol per Glu substitution. A Glu substitution at the C-terminus destabilizes the coiled-coil more than in the middle owing to the combined effects of intrinsic destabilization and unfavourable charge–helix dipole interaction with the negative pole of the helix dipole. The estimated destabilizing charge–helix dipole interaction of 0.08 kcal/mol is smaller than the stabilizing interaction at the N-terminus. The presence of a 2-2′disulphide bridge appears to have little influence on the magnitude of the charge–helix dipole interactions at either end of the coiled-coil. © 1997 European Peptide Society and John Wiley & Sons, Ltd.     

Surface salt bridges stabilize the GCN4 leucine zipper.

We present a study of the role of salt bridges in stabilizing a simplified tertiary structural motif, the coiled-coil. Changes in GCN4 sequence have been engineered that introduce trial patterns of single and multiple salt bridges at solvent exposed sites. At the same sites, a set of alanine mutants was generated to provide a reference for thermodynamic analysis of the salt bridges. Introduction of three alanines stabilizes the dimer by 1.1 kcal/mol relative to the wild-type. An arrangement corresponding to a complex type of salt bridge involving three groups stabilizes the dimer by 1.7 kcal/ mol, an apparent elevation of the melting temperature relative to wild type of about 22 degrees C. While identifying local from nonlocal contributions to protein stability is difficult, stabilizing interactions can be identified by use of cycles. Introduction of alanines for side chains of lower helix propensity and complex salt bridges both stabilize the coiled-coil, so that combining the two should yield melting temperatures substantially higher than the starting species, approaching those of thermophilic sequences.     

Protein denaturation with guanidine hydrochloride or urea provides a different estimate of stability depending on the contributions of electrostatic interactions.

The objective of this study was to address the question of whether or not urea and guanidine hydrochloride (GdnHCl) give the same estimates of the stability of a particular protein. We previously suspected that the estimates of protein stability from GdnHCl and urea denaturation data might differ depending on the electrostatic interactions stabilizing the proteins. Therefore, 4 coiled-coil analogs were designed, where the number of intrachain and interchain electrostatic attractions (A) were systematically changed to repulsions (R): 20A, 15A5R, 10A10R, and 20R. The GdnHCl denaturation data showed that the 4 coiled-coil analogs, which had electrostatic interactions ranging from 20 attractions to 20 repulsions, had very similar [GdnHCl]1/2 values (average of congruent to 3.5 M) and, as well, their delta delta Gu values were very close to 0 (0.2 kcal/mol). In contrast, urea denaturation showed that the [urea]1/2 values proportionately decreased with the stepwise change from 20 electrostatic attractions to 20 repulsions (20A, 7.4 M; 15A5R, 5.4 M; 10A10R, 3.2 M; and 20R, 1.4 M), and the delta delta Gu values correspondingly increased with the increasing differences in electrostatic interactions (20A-15A5R, 1.5 kcal/mol; 20A-10A10R, 3.7 kcal/mol; and 20A-20R, 5.8 kcal/mol). These results indicate that the ionic nature of GdnHCl masks electrostatic interactions in these model proteins, a phenomenon that was absent when the unchanged urea was used. Thus, GdnHCl and urea denaturations may give vastly different estimates of protein stability, depending on how important electrostatic interactions are to the protein.     

Unique stabilizing interactions identified in the two-stranded alpha-helical coiled-coil: crystal structure of a cortexillin I/GCN4 hybrid coiled-coil peptide

We determined the 1.17 A resolution X-ray crystal structure of a hybrid peptide based on sequences from coiled-coil regions of the proteins GCN4 and cortexillin I. The peptide forms a parallel homodimeric coiled-coil, with C(alpha) backbone geometry similar to GCN4 (rmsd value 0.71 A). Three stabilizing interactions have been identified: a unique hydrogen bonding-electrostatic network not previously observed in coiled-coils, and two other hydrophobic interactions involving leucine residues at positions e and g from both g-a' and d-e' interchain interactions with the hydrophobic core. This is also the first report of the quantitative significance of these interactions. The GCN4/cortexillin hybrid surprisingly has two interchain Glu-Lys' ion pairs that form a hydrogen bonding network with the Asn residues in the core. This network, which was not observed for the reversed Lys-Glu' pair in GCN4, increases the combined stability contribution of each Glu-Lys' salt bridge across the central Asn15-Asn15' core to approximately 0.7 kcal/mole, compared to approximately 0.4 kcal mole(-1) from a Glu-Lys' salt bridge on its own. In addition to electrostatic and hydrogen bonding stabilization of the coiled-coil, individual leucine residues at positions e and g in the hybrid peptide also contribute to stability by 0.7 kcal/mole relative to alanine. These interactions are of critical importance to understanding the stability requirements for coiled-coil folding and in modulating the stability of de novo designed macromolecules containing this motif.     

Defining the minimum size of a hydrophobic cluster in two-stranded alpha-helical coiled-coils: effects on protein stability

The alpha-helical coiled-coil motif is characterized by a heptad repeat pattern (abcdefg)(n) in which residues a and d form the hydrophobic core. Long coiled-coils (e.g., tropomyosin, 284 residues per polypeptide chain) typically do not have a continuous hydrophobic core of stabilizing residues, but rather one that consists of alternating clusters of stabilizing and destabilizing residues. We have arbitrarily defined a cluster as a minimum of three consecutive stabilizing or destabilizing residues in the hydrophobic core. We report here on a series of two-stranded, disulfide-bridged parallel alpha-helical coiled-coils that contain a central cassette of three consecutive hydrophobic core positions (d, a, and d) with a destabilizing cluster of three consecutive Ala residues in the hydrophobic core on each side of the cassette. The effect of adding one to three stabilizing hydrophobes in these positions (Leu or Ile; denoted as [see text]) was investigated. Alanine residues (denoted as [see text]) are used to represent destabilizing residues. The peptide with three Ala residues in the d a d cassette positions ([see text]) was among the least stable coiled-coil (T(m) = 39.3 degrees C and Urea(1/2) = 1.9 M). Surprisingly, the addition of one stabilizing hydrophobe (Leu) to the cassette or two stabilizing hydrophobes (Leu), still interspersed by an Ala in the cassette ([see text]), also did not lead to any gain in stability. However, peptides with two adjacent hydrophobes in the cassette ([see text])([see text]) did show a gain in stability of 0.9 kcal/mole over the peptide with two interspersed hydrophobes ([see text]). Because the latter three peptides have the same inherent hydrophobicity, the juxtaposition of stabilizing hydrophobes leads to a synergistic effect, and thus a clustering effect. The addition of a third stabilizing hydrophobe to the cassette ([see text]) resulted in a further synergistic gain in stability of 1.7 kcal/mole (T(m) = 54.1 degrees C and Urea(1/2) = 3.3M). Therefore, the role of hydrophobicity in the hydrophobic core of coiled-coils is extremely context dependent and clustering is an important aspect of protein folding and stability.     

Design of a leucine zipper coiled coil stabilized 1.4 kcal mol-1 by phosphorylation of a serine in the e position.

Using a dimeric bZIP protein, we have designed a leucine zipper that becomes more stable after a serine in the e position is phosphorylated by protein kinase A (delta delta GP = -1.4 kcal mol-1 dimer-1 or -0.7 kcal mol-1 residue-1). Mutagenesis studies indicate that three arginines form a network of inter-helical (i,i' + 5; i, i' + 2) and intra-helical (i, i + 4) attractive interactions with the phosphorylated serine. When the arginines are replaced with lysines, the stabilizing effect of serine phosphorylation is reduced (delta delta GP = -0.5 kcal mol-1 dimer-1). The hydrophobic interface of the leucine zipper needs a glycine in the d position to obtain an increase in stability after phosphorylation. The phosphorylated protein binds DNA with a 15-fold higher affinity. Using a transient transfection assay, we document a PKA dependent four-fold activation of a reporter gene. Phosphorylation of a threonine in the same e position decreases the stability by delta delta GP = +1.2 kcal mol-1 dimer-1. We present circular dichroism (CD) thermal denaturations of 15 bZIP proteins before and after phosphorylation. These data provide insights into the structural determinants that result in stabilization of a coiled coil by phosphorylation.     

Simple and versatile restraints for the accurate modeling of alpha-helical coiled-coil structures of multiple strandedness, orientation and composition

We present a minimalist approach for the modeling of the three-dimensional structure of multistranded alpha-helical coiled coils. The approach is based on empirical principles introduced by F. H. C. Crick (F. H. C. Crick, Acta Crystallogr, 1953, Vol. 6, pp. 689-697). Crick hypothesized that keeping the distance between the residues at the interacting interface of alpha-helices constant would lead to supercoiling or the formation of a coiled coil through the knobs-into-holes mode of packing. We have implemented the latter hypothesis in a simulating annealing protocol in the simple form of interhelical distance restraints (two per heptad) between Calpha at the interfacial positions and. To demonstrate the authenticity of Crick's hypothesis and the precision and accuracy of our approach, we have modeled the crystal structures of six synthetic coiled coils in dimeric, trimeric, and tetrameric states. The mean root mean square deviations (RMSDs) between the backbone atoms of the ensemble of structures calculated and those of the corresponding geometric averages is always below 0.76 A, indicating that our protocol has an excellent degree of convergence and precision. The RMSDs between the backbone atoms of each of the six geometric average structures and the backbone of the corresponding crystal structures all range between 0.43 and 0.95 A, indicative of excellent accuracy and proving the authenticity of Crick's hypothesis. Moreover, without specifying any dihedral angles, we found that in 81% of the occurrences, the most populated conformer of the side chains at positions and in the ensembles calculated were identical to those observed in the crystal structures. This shows that our simple approach, which is the simplest reported so far, can generate accurate results for the backbone and side chains. Finally, as a test case for a wider application of our approach in the field of structural proteomics, we describe the successful modeling of the overall structure of SNARE and the organization of its interfacial ionic layer known to play an important functional role.     

Packing and hydrophobicity effects on protein folding and stability: effects of beta-branched amino acids, valine and isoleucine, on the formation and stability of two-stranded alpha-helical coiled coils/leucine zippers.

The aim of this study was to examine the differences between hydrophobicity and packing effects in specifying the three-dimensional structure and stability of proteins when mutating hydrophobes in the hydrophobic core. In DNA-binding proteins (leucine zippers), Leu residues are conserved at positions "d," and beta-branched amino acids, Ile and Val, often occur at positions "a" in the hydrophobic core. In order to discern what effect this selective distribution of hydrophobes has on the formation and stability of two-stranded alpha-helical coiled coils/leucine zippers, three Val or three Ile residues were simultaneously substituted for Leu at either positions "a" (9, 16, and 23) or "d" (12, 19, and 26) in both chains of a model coiled coil. The stability of the resulting coiled coils was monitored by CD in the presence of Gdn.HCl. The results of the mutations of Ile to Val at either positions "a" or "d" in the reduced or oxidized coiled coils showed a significant hydrophobic effect with the additional methylene group in Ile stabilizing the coiled coil (delta delta G values range from 0.45 to 0.88 kcal/mol/mutation). The results of mutations of Leu to Ile or Val at positions "a" in the reduced or oxidized coiled coils showed a significant packing effect in stabilizing the coiled coil (delta delta G values range from 0.59 to 1.03 kcal/mol/mutation). Our results also indicate the subtle control hydrophobic packing can have not only on protein stability but on the conformation adopted by the amphipathic alpha-helices. These structural findings correlate with the observation that in DNA-binding proteins, the conserved Leu residues at positions "d" are generally less tolerant of amino acid substitutions than the hydrophobic residues at positions "a."     

Analysis of the role of the leucine zipper motif in regulating the ability of AFAP-110 to alter actin filament integrity

AFAP-110 has an intrinsic ability to alter actin filament integrity as an actin filament crosslinking protein. This capability is regulated by a carboxy terminal leucine zipper (Lzip) motif. The Lzip motif facilitates self-association stabilizing the AFAP-110 multimers. Deletion of the Lzip motif (AFAP-110(Deltalzip)) reduces the stability of the AFAP-110 multimer and concomitantly increases its ability to crosslink actin filaments, in vitro, and to activate cSrc and alter actin filament integrity, in vivo. We sought to determine how the Lzip motif regulates AFAP-110 function. Substitution of the c-Fos Lzip motif in place of the AFAP-110 Lzip motif (AFAP-110(fos)) was predicted to preserve the alpha-helical structure while changing the sequence. To alter the structure of the alpha-helix, a leucine to proline mutation was generated in the AFAP-110 alpha-helical Lzip motif (AFAP-110(581P)), which largely preserved the sequence. The helix mutants, AFAP-110(Deltalzip), AFAP-110(fos), and AFAP-110(581P), demonstrated reduced multimer stability with an increased capacity to crosslink actin filaments, in vitro, relative to AFAP-110. An analysis of opposing binding sites indicated that the carboxy terminus/Lzip motif can contact sequences within the amino terminal pleckstrin homology (PH1) domain indicating an auto-inhibitory mechanism for regulating multimer stability and actin filament crosslinking. In vivo, only AFAP-110(Deltalzip) and AFAP-110(581P) were to activate cSrc and to alter cellular actin filament integrity. These data indicate that the intrinsic ability of AFAP-110 to crosslink actin filaments is dependent upon both the sequence and structure of the Lzip motif, while the ability of the Lzip motif to regulate AFAP-110-directed activation of cSrc and changes in actin filament integrity in vivo is dependent upon the structure or presence of the Lzip motif. We hypothesize that the intrinsic ability of AFAP-110 to crosslink actin filaments or activate cSrc are distinct functions.     

Crystal structure of S. aureus YlaN, an essential leucine rich protein involved in the control of cell shape

The crystal structure of a conserved leucine rich protein, YlaN, from Staphylococcus aureus has been determined by X-ray crystallography to 2.3 A resolution. Whilst the precise function of S. aureus YlaN is unknown its homologue in B. subtilis has been shown to be essential for cell survival and is thought to be involved in controlling cell shape. The structure of S. aureus YlaN provides the first view of its protein family, which reveals that it is a novel homodimer whose subunit architecture is comprised of an antiparallel 3 helix bundle reminiscent of the helical arrangements seen in leucine zipper proteins. Analysis of the pattern of sequence conservation on the structure has led to the identification of two connected solvent exposed patches of conserved residues in each subunit located at one end of but on opposite faces of the molecule. We suggest that YlaN has a binding role in the cell rather than a catalytic function and a search for its ligand is underway to accelerate its exploitation as a target for antibiotic discovery.     

Optimization of the electrostatic interactions between ionized groups and peptide dipoles in proteins.

The three-dimensional optimization of the electrostatic interactions between the charged amino acid residues and the peptide partial charges was studied by Monte-Carlo simulations on a set of 127 nonhomologous protein structures with known atomic coordinates. It was shown that this type of interaction is very well optimized for all proteins in the data set, which suggests that they are a valuable driving force, at least for the native side-chain conformations. Similar to the optimization of the charge-charge interactions (Spassov VZ, Karshikoff AD, Ladenstein R, 1995, Protein Sci 4:1516-1527), the optimization effect was found more pronounced for enzymes than for proteins without enzymatic function. The asymmetry in the interactions of acidic and basic groups with the peptide dipoles was analyzed and a hypothesis was proposed that the properties of peptide dipoles are a factor contributing to the natural selection of the basic amino acids in the chemical composition of proteins.     

The effects of intramolecular and intermolecular electrostatic repulsions on the stability and aggregation of NISTmAb revealed by HDX‐MS,DSC, and nanoDSF

The stability and aggregation of NIST monoclonal antibody (NISTmAb) were investigated by hydrogen/deuterium exchange mass spectrometry (HDX‐MS), differential scanning calorimetry (DSC), and nano‐differential scanning fluorimetry (nanoDSF). NISTmAb was prepared in eight formulations at four different pHs (pH 5, 6, 7, and 8) in the presence and absence of 150 mM NaCl and analyzed by the three methods. The HDX‐MS results showed that NISTmAb is more conformationally stable at a pH near its isoelectric point (pI) in the presence of NaCl than a pH far from its pI in the absence of NaCl. The stabilization effects were global and not localized. The midpoint temperature of protein thermal unfolding transition results also showed the C_H2 domain of the protein is more conformationally stable at a pH near its pI. On the other hand, the onset of aggregation temperature results showed that NISTmAb is less prone to aggregate at a pH far from its pI, particularly in the absence of NaCl. These seemingly contradicting results, higher conformational stability yet higher aggregation propensity near the pI than far away from the pI, can be explained by intramolecular and intermolecular electrostatic repulsion using Lumry‐Eyring model, which separates folding/unfolding equilibrium and aggregation event. The further a pH from the pI, the higher the net charge of the protein. The higher net charge leads to greater intramolecular and intermolecular electrostatic repulsions. The greater intramolecular electrostatic repulsion destabilizes the protein and the greater intermolecular electrostatic repulsion prevents aggregation of the protein molecules at pH far from the pI.     

The Drosophila homolog of human AF10/AF17 leukemia fusion genes (Dalf) encodes a zinc finger/leucine zipper nuclear protein required in the nervous system for maintaining EVE expression and normal growth

Sibling neurons in the embryonic central nervous system (CNS) of Drosophila can adopt distinct states as judged by gene expression and axon projection. In the NB4-2 lineage, two even-skipped (eve)-expressing sibling neuronal cells, RP2 and RP2sib, are formed in each hemineuromere. Throughout embryogenesis, only RP2, but not RP2sib, maintains eve expression. In this report, we describe a P-element induced mutation that alters the expression pattern of EVE in RP2 motoneurons in the Drosophila embryonic CNS. The mutation was mapped to a Drosophila homolog of human AF10/AF17 leukemia fusion genes (alf), and therefore named Dalf. Like its human counterparts, Dalf encodes a zinc finger/leucine zipper nuclear protein that is widely expressed in embryonic and larval tissues including neurons and glia. In Dalf mutant embryos, the RP2 motoneuron no longer maintains EVE expression. The effect of the Dalf mutation on EVE expression is RP2-specific and does not affect other characteristics of the RP2 motoneuron. In addition to the embryonic phenotype, Dalf mutant larvae are retarded in their growth and this defect can be rescued by the ectopic expression of a Dalf transgene under the control of a neuronal GAL4 driver. This indicates a requirement for Dalf function in the nervous system for maintaining gene expression and the facilitation of normal growth.     

Electrostatic interactions in the GCN4 leucine zipper: substantial contributions arise from intramolecular interactions enhanced on binding.

The GCN4 leucine zipper is a peptide homodimer that has been the subject of a number of experimental and theoretical investigations into the determinants of affinity and specificity. Here, we utilize this model system to investigate electrostatic effects in protein binding using continuum calculations. A particularly novel feature of the computations made here is that they provide an interaction-by-interaction breakdown of the electrostatic contributions to the free energy of docking that includes changes in the interaction of each functional group with solvent and changes in interactions between all pairs of functional groups on binding. The results show that (1) electrostatic effects disfavor binding by roughly 15 kcal/mol due to desolvation effects that are incompletely compensated in the bound state, (2) while no groups strongly stabilize binding, the groups that are most destabilizing are charged and polar side chains at the interface that have been implicated in determining binding specificity, and (3) attractive intramolecular interactions (e.g., backbone hydrogen bonds) that are enhanced on binding due to reduced solvent screening in the bound state contribute significantly to affinity and are likely to be a general effect in other complexes. A comparison is made between the results obtained in an electrostatic analysis carried out calculationally and simulated results corresponding to idealized data from a scanning mutagenesis experiment. It is shown that scanning experiments provide incomplete information on interactions and, if overinterpreted, tend to overestimate the energetic effect of individual side chains that make attractive interactions. Finally, a comparison is made between the results available from a continuum electrostatic model and from a simpler surface-area dependent solvation model. In this case, although the simpler model neglects certain interactions, on average it performs rather well.