Measurement of interhelical electrostatic interactions in the GCN4 leucine zipper

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

The effects on protein stability of negatively charged Glu side chains at different positions along the length of the alpha-helix were investigated in the two-stranded alpha-helical coiled-coil. A native coiled-coil has been designed which consists of two identical 35 residue polypeptide chains with a heptad repeat QgVaGbAcLdQeKf 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 degrees 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 pH7. 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.

Buried asparagines determine the dimerization specificities of leucine zipper mutants

Regulation of gene expression by many transcription factors is controlled by specific combinations of homo- and heterodimers through a short alpha-helical coiled-coil known as a leucine zipper. The dimer interface of a leucine zipper involves side chains of the residues at the a, d, e, and g positions of the (abcdefg)n heptad repeat. To understand the basis for the specificity of dimer formation, we characterized GCN4 leucine zipper mutants with all 16 possible permutations and combinations of isoleucines and asparagines at four a positions in the dimer interface, using a genetic test for the specificity of dimer formation by lambda repressor-leucine zipper fusions. Heterodimers were detected by loss of repressor activity in the presence of a fusion to a dominant-negative mutant form of the DNA-binding domain of repressor. Reconstruction experiments using leucine zippers from GCN4, Jun, Fos, and C/EBP showed that this assay distinguishes pairs that form heterodimers from those that do not. We found that the mutants have novel dimerization specificities determined by the positioning of buried asparagine residues at the a positions. The pattern of buried polar residues could also explain the dimerization specificities of some naturally occurring leucine zippers. The altered specificity mutants described here should be useful for the construction of artificial regulatory circuitry.

Salt effects on protein stability: two-stranded alpha-helical coiled-coils containing inter- or intrahelical ion pairs

An investigation into the role of surface-accessible ion pairs in protein stability was carried out by determining the effects of added salt (KCl, MgCl2 and LaCl3) at neutral and acidic pH on the stability of de novo designed two-stranded alpha-helical coiled-coils. The effects of salt on the stability of coiled-coils containing interhelical i to i' + 5 or intrahelical i to i + 4 and i to i + 3 Lys-Glu ion pairs were compared to the effects on the stability of a control coiled-coil, which contained no intra- or interhelical ion pair. These studies indicate that ionic interactions contribute to coiled-coil stability. The results show that added salt can have complex effects on protein stability, involving stabilizing and destabilizing contributions with the net effect depending on the nature of the charged residues and ionic interactions present in the protein.

Interaction between proteins containing homeodomains associated to leucine zippers from sunflower

A strategy based on the use of PCR with one degenerate oligonucleotide deduced from conserved sequences and lambda gt10 primers was used to isolate homeobox containing sequences from sunflower stem and root cDNA libraries. Six different partial cDNAs coding for the first 48 amino acids of homeodomains and amino terminal sequences were analyzed and found to be members of the HD-Zip superfamily, which contain a homeobox linked to a leucine zipper coding region. A full-length cDNA clone, Hahb-10, was isolated and characterized. The leucine zipper portions of Hahb-10 and of the previously reported Hahb-1 have been utilized to construct fusions with the N-terminal domain of the lambda repressor. These fusions were tested for their ability to bind to lambda promoters in vivo. The expression of a protein containing an active dimerization domain, but not capable of DNA binding, exerts a dominant negative effect on the ability of repressor-zipper fusions to bind to its target DNA. From these experiments, it was concluded that Hahb-1 and -10, when co-expressed, form preferentially homodimers. Exchange of conserved threonines and leucines at positions a1 and d1 of both zippers reduces dimerization efficiency and allows the formation of heterodimers, suggesting that these residues are, among others, determinants of the specificity of interaction, most likely through changes in hydrophobic packing interactions at the dimer interface. The results imply that a great number of interacting molecular entities compose this protein superfamily which is presumably involved in regulating plant developmental responses.

Derivation and testing of pair potentials for protein folding. When is the quasichemical approximation correct?

Many existing derivations of knowledge-based statistical pair potentials invoke the quasichemical approximation to estimate the expected side-chain contact frequency if there were no amino acid pair-specific interactions. At first glance, the quasichemical approximation that treats the residues in a protein as being disconnected and expresses the side-chain contact probability as being proportional to the product of the mole fractions of the pair of residues would appear to be rather severe. To investigate the validity of this approximation, we introduce two new reference states in which no specific pair interactions between amino acids are allowed, but in which the connectivity of the protein chain is retained. The first estimates the expected number of side-chain contracts by treating the protein as a Gaussian random coil polymer. The second, more realistic reference state includes the effects of chain connectivity, secondary structure, and chain compactness by estimating the expected side-chain contrast probability by placing the sequence of interest in each member of a library of structures of comparable compactness to the native conformation. The side-chain contact maps are not allowed to readjust to the sequence of interest, i.e., the side chains cannot repack. This situation would hold rigorously if all amino acids were the same size. Both reference states effectively permit the factorization of the side-chain contact probability into sequence-dependent and structure-dependent terms. Then, because the sequence distribution of amino acids in proteins is random, the quasichemical approximation to each of these reference states is shown to be excellent. Thus, the range of validity of the quasichemical approximation is determined by the magnitude of the side-chain repacking term, which is, at present, unknown. Finally, the performance of these two sets of pair interaction potentials as well as side-chain contact fraction-based interaction scales is assessed by inverse folding tests both without and with allowing for gaps.

Cell transformation mediated by homodimeric E2A-HLF transcription factors

The E2A-HLF fusion gene, created by the t(17;19)(q22;p13) chromosomal translocation in pro-B lymphocytes, encodes an oncogenic protein in which the E2A trans-activation domain is linked to the DNA-binding and protein dimerization domain of hepatic leukemia factor (HLF), a member of the proline- and acidic amino acid-rich (PAR) subfamily of bZIP transcription factors. This fusion product binds to its DNA recognition site not only as a homodimer but also as a heterodimer with HLF and two other members of the PAR bZIP subfamily, thyrotroph embryonic factor (TEF) and albumin promoter D-box binding protein (DBP). Thus, E2A-HLF could transform cells by direct regulation of downstream target genes, acting through homodimeric or heterodimeric complexes, or by sequestering normal PAR proteins into nonfunctional heterocomplexes (dominant-negative interference). To distinguish among these models, we constructed mutant E2A-HLF proteins in which the leucine zipper domain of HLF was extended by one helical turn or altered in critical charged amino acids, enabling the chimera to bind to DNA as a homodimer but not as a heterodimer with HLF or other PAR proteins. When introduced into NIH 3T3 cells in a zinc-inducible vector, each of these mutants induced anchorage-independent growth as efficiently as unaltered E2A-HLF, indicating that the chimeric oncoprotein can transform cells in its homodimeric form. Transformation also depended on an intact E2A activator region, providing further support for a gain-of-function contribution to oncogenesis rather than one based on a dominant-interfering or dominant-negative mechanism. Thus, the tumorigenic effects of E2A-HLF and its mutant forms in NIH 3T3 cells favor a straightforward model in which E2A-HLF homodimers bind directly to promoter/enhancer elements of downstream target genes and alter their patterns of expression in early B-cell progenitors.

A recipe for designing water-soluble, beta-sheet-forming peptides

Based on observations of solubility and folding properties of peptide 33-mers derived from the beta-sheet domains of platelet factor-4 (PF4), interleukin-8 (IL-8), and growth related protein (Gro-alpha), as well as other beta-sheet-forming peptides, general guidelines have been developed to aid in the design of water soluble, self-association-induced beta-sheet-forming peptides. CD, 1H-NMR, and pulsed field gradient NMR self-diffusion measurements have been used to assess the degree of folding and state of aggregation. PF4 peptide forms native-like beta-sheet tetramers and is sparingly soluble above pH 6. IL-8 peptide is insoluble between pH 4.5 and pH 7.5, yet forms stable, native-like beta-sheet dimers at higher pH. Gro-alpha peptide is soluble at all pH values, yet displays no discernable beta-sheet structure even when diffusion data indicate dimer-tetramer aggregation. A recipe used in the de novo design of water-soluble beta-sheet-forming peptides calls for the peptide to contain 40-50% hydrophobic residues, usually aliphatic ones (I, L, V, A, M) (appropriately paired and mostly but not always alternating with polar residues in the sheet sequence), a positively charged (K, R) to negatively charged (E, D) residue ratio between 4/2 and 6/2, and a noncharged polar residue (N, Q, T, S) composition of about 20% or less. Results on four de novo designed, 33-residue peptides are presented supporting this approach. Under near physiologic conditions, all four peptides are soluble, form beta-sheet structures to varying degrees, and self-associate. One peptide folds as a stable, compact beta-sheet tetramer, whereas the others are transient beta-sheet-containing aggregates.

High resolution NMR solution structure of the leucine zipper domain of the c-Jun homodimer

The solution structure of the c-Jun leucine zipper domain has been determined to high resolution using a new calculation protocol designed to handle highly ambiguous sets of interproton distance restraints. The domain comprises a coiled coil of parallel alpha-helices in which most of the hydrophobic residues are buried at the highly symmetrical dimer interface; this interface extends over 10 helical turns and is the most elongated protein domain solved to date using NMR methods. The backbone fold is very similar to that seen in crystal structures of the GCN4 and Jun-Fos leucine zippers; however, in contrast with these crystal structures, the Jun leucine zipper dimer appears to be devoid of favorable intermolecular electrostatic interactions. A polar asparagine residue, located at the dimer interface, forms the sole point of asymmetry in the structure; furthermore, the side chain of this residue is disordered due to motional averaging. This residue, which is highly conserved in the leucine zipper family of transcription factors, provides a destabilizing influence that is likely to facilitate the rapid exchange of zipper strands in vivo.

Sequence-specific DNA-binding dominated by dehydration

Fluorescence spectroscopy and isothermal titration calorimetry were used to study the thermodynamics of binding of the glucocorticoid receptor DNA-binding domain to four different, but similar, DNA-binding sites. The binding sites are two naturally occurring sites that differ in the composition of one base pair, i.e., an A-T to G-C mutation, and two sites containing chemical intermediates of these base pairs. The calorimetrically determined heat capacity change (Delta C(p)o(obs)) for glucocorticoid receptor DNA-binding domain binding agrees with that calculated for dehydration of solvent-accessible surface areas. A dominating effect of dehydration or solvent reorganization on the thermodynamics is also consistent with an observed linear relationship between observed enthalpy change (Delta Ho(obs)) and observed entropy change (Delta So(obs)) with a slope close to the experimental temperature. Comparisons with structural data allow us to rationalize individual differences between Delta Ho(obs) (and Delta So(obs)) for the four complexes. For instance, we find that the removal of a methyl group at the DNA-protein interface is enthalpically favorable but entropically unfavorable, which is consistent with a replacement by an ordered water molecule.

Evaluating electrostatic contributions to binding with the use of protein charge ladders

Electrostatic interactions between charges on ligands and charges on proteins that are remote from the binding interface can influence the free energy of binding (delta Gb). The binding affinities between charged ligands and the members of a charge ladder of bovine carbonic anhydrase (CAII) constructed by random acetylation of the amino groups on its surface were measured by affinity capillary electrophoresis (ACE). The values of delta Gb derived from this analysis correlated approximately linearly with the charge. Opposite charges on the ligand and the members of the charge ladder of CAII were stabilizing; like charges were destabilizing. The combination of ACE and protein charge ladders provides a tool for quantitatively examining the contributions of electrostatics to free energies of molecular recognition in biology.

Investigation of electrostatic interactions in two-stranded coiled-coils through residue shuffling

The effects of electrostatic interactions on the stability of coiled-coils were investigated using the strategy of shuffling the sequence without changing the overall content of amino acid residues in the peptides. Shuffling the sequence provides peptides with thermodynamically similar unfolded states. Therefore, the unfolded state can be used as a universal reference state in comparing the thermodynamic properties of the folded coiled-coil structure of the peptides, while varying the configuration of ionized groups, that is, changing the types and number of potential electrostatic interactions. The relative stabilities of these states were determined by monitoring the temperature-induced folding/unfolding of the peptides in solutions with different pH and ionic strength by circular dichroism spectroscopy and scanning microcalorimetry. It was found that, in solutions with low ionic strength, ionic pairs contribute significantly to the stability of the coiled-coil conformation. The stability increases with an increase in the number of ionized groups in the peptide upon changing pH from acidic to neutral. In contrast, in the solutions with high ionic strength, the coiled-coil becomes less stable at neutral pH than at acidic pH. Most surprisingly, the increase in Gibbs energy of stabilization of the coiled-coil state with increasing pH at low ionic strength proceeds with a decrease in the enthalpy and entropy of unfolding. This observation can be explained only by hydration of ionized groups upon unfolding of coiled-coils which is associated with significant negative enthalpy and entropy effects.

Roles of electrostatic interaction in proteins

More than 60 years after the analyses by Linderstrom-Lang and Kirkwood of their hypothetical 'protein' structures, we have now a plethora of experimental evidence and computational estimates of the electrostatic forces in proteins, with very many protein 3D structures at atomic resolution. In the mean time, there were in the beginning, many arguments and suggestions about the roles of electrostatics, mainly from empirical findings and tendencies. A few experimental results indicated that the electrostatic contribution is of the order of several kcal/mol, which was theoretically difficult to reproduce correctly, because a large opposing reaction field should be subtracted from a large, direct Coulombic field. Although the importance of the reaction field was recognized even 70 years ago, appropriate applications to protein molecules were made only in this decade, with the development of numerical computation. Now, an electrostatic molecular surface is one of the most popular pictures in journals of structural biology, indicating that the electrostatic force is one of the important components contributing to molecular recognition, which is a major focus of current biology and biochemistry. The development of NMR techniques has made it possible to observe the individual ionizations of ionizable groups in a protein, in addition to the determination of the 3D structure. Since it does not require any additional probe, each charge state can report the very local and heterogeneous electrostatic potentials working in the protein, without disturbing the original field. From the pKa values, the contributions of electrostatic interactions, ion pairs, charge-dipole interactions, and hydrogen bonds to protein stability have been correctly evaluated. Protein engineering also provides much more information than that obtainable from the native proteins, as the residues concerned can now be easily substituted with other amino acid residues having electrostatically different characteristics. Those experimental results have revealed smaller contributions than previously expected, probably because we underestimated the reaction field effects. Especially, a single ion pair stabilizes a protein only slightly, although a cooperative salt-bridge network can contribute significantly to protein stability. Marginal stabilities of proteins arise from small difference between many factors with driving and opposing forces. In spite of the small contribution of each single electrostatic interaction to the protein stability, the sum of their actions works to maintain the specific 3D structure of the protein. The 'negative' roles of electrostatics, which might destabilize protein conformation, should be pointed out. Unpaired buried charges are energetically too expensive to exit in the hydrophobic core. Isolated hydrogen bond donors and acceptors also exert negative effects, but they are not as expensive as the unpaired buried charges, with costs of a few kcal/mol. Therefore, statistical analyses of protein 3D structures reveal only rare instances of isolated hydrogen bond donors and acceptors. This must be the main reason why alpha-helices and beta-sheets are only observed in protein cores as the backbone structures. Such secondary structures do not leave any backbone hydrogen donors or acceptors unpaired, because of their intrinsically regular packing. Otherwise, it might be very difficult to construct a backbone structure, in which all the backbone amide and carbonyl groups had their own hydrogen bond partners in the protein core. There are two theoretical approaches to protein electrostatics, the macroscopic or continuum model, and the microscopic or molecular model. As described in this article, the macroscopic model has inherent problems because the protein-solvent system is very hetergeneous from the physical point of view...

Double-mutant cycles: a powerful tool for analyzing protein structure and function

A double-mutant cycle involves wild-type protein, two single mutants and the corresponding double mutant protein. If the change in free energy associated with a structural or functional property of the protein upon a double mutation differs from the sum of changes in free energy due to the single mutations, then the residues at the two positions are coupled. Such coupling reflects either direct or indirect interactions between these residues. Double-mutant cycle analysis can be used to measure the strength of intramolecular and intermolecular pairwise interactions in proteins or protein-ligand complexes with known structure. Double-mutant cycles can also be employed to characterize structures that are inaccessible to NMR and X-ray crystallography, such as those of transition states for protein folding, ligand binding and enzyme catalysis, or of membrane proteins. Multidimensional mutant cycle analysis can be used to measure higher-order cooperativity between intramolecular or intermolecular interactions. In the absence of coupling between residues, prediction of mutational effects is possible by assuming their additivity.

Genetic strategy for analyzing specificity of dimer formation: Escherichia coli cyclic AMP receptor protein mutant altered in its dimerization specificity

Many transcriptional regulators function in homo- or heterodimeric combinations. The same protein can carry out distinct regulatory functions depending on the partner with which it associates. Here, we describe a mutant of the Escherichia coli cAMP receptor protein (CRP) that has an altered dimerization specificity; that is, mutant/mutant homodimers form preferentially over wild-type/mutant heterodimers. CRP dimerization involves the formation of a parallel coiled-coil structure, and our CRP mutant bears an amino acid substitution affecting the first "d" position residue within the alpha-helix that mediates CRP dimerization. The genetic strategy we used to isolate this CRP altered dimerization specificity (ADS) mutant is generalizable and could be utilized to isolate ADS mutants of other dimeric transcriptional regulators.

The structure of Pyrococcus furiosus glutamate dehydrogenase reveals a key role for ion-pair networks in maintaining enzyme stability at extreme temperatures

BACKGROUND

The hyperthermophile Pyrococcus furiosus is one of the most thermostable organisms known, with an optimum growth temperature of 100 degrees C. The proteins from this organism display extreme thermostability. We have undertaken the structure determination of glutamate dehydrogenase from P. furiosus in order to gain further insights into the relationship between molecular structure and thermal stability.

RESULTS

The structure of P. furiosus glutamate dehydrogenase, a homohexameric enzyme, has been determined at 2.2 A resolution and compared with the structure of glutamate dehydrogenase from the mesophile Clostridium symbiosum.

CONCLUSIONS

Comparison of the structures of these two enzymes has revealed one major difference: the structure of the hyperthermophilic enzyme contains a striking series of ion-pair networks on the surface of the protein subunits and buried at both interdomain and intersubunit interfaces. We propose that the formation of such extended networks may represent a major stabilizing feature associated with the adaptation of enzymes to extreme temperatures.

The effects of interhelical electrostatic repulsions between glutamic acid residues in controlling the dimerization and stability of two-stranded alpha-helical coiled-coils

The effects of interhelical electrostatic repulsions in controlling the dimerization and stability of two-stranded alpha-helical coiled-coils have been studied using de novo designed synthetic coiled-coils. A native coiled-coil was snythesized, which consisted of two identical 35-residue polypeptide chains with a heptad repeat QgVaGbAcLdQeKf and a Cys residue at position 2 to allow formation of an interchain 2-2' disulfide bridge. This peptide, designed to contain no intrachain or interchain electrostatic interactions, forms a stable coiled-coil structure at 20 degrees C in benign medium (50 mM KCl, 25 mM PO4, pH 7) with a [urea]1/2 value of 6.1 M. Five mutant coiled-coils were designed in which Gln residues at the e and g positions of the heptad repeat were substituted with Glu systematically from the N terminus toward the C terminus, resulting in each polypeptide chain having 2, 4, 6, 8, or 10 Glu residues. These substituted Glu residues are able to form interchain i to i' +5 electrostatic repulsions across the dimer interface. As the number of interchain repulsions increases, a steady loss of helical content is observed by circular dichroism spectroscopy. The effects of the interchain Glu-Glu repulsions on the coiled-coil structure are partly overcome by the presence of an interchain disulfide bridge; the peptide with six Glu substitutions is only 15% helical in the reduced form but 85% helical in the oxidized form. The stabilities of the coiled-coils were determined by urea and guanidine hydrochloride (GdnHCl) denaturation studies at 20 degrees C. The stabilities of the coiled-coils determined by urea denaturation indicate a decrease in stability, which correlates with an increasing number of interchain repulsions ([urea]1/2 values ranging from 8.4 to 3.7 M in the presence of M KCl). In contrast, all coiled-coils had similar stabilities when determined by GdnHCl denaturation (approximately 2.9 M). KCl could not effectively screen the effects of interchain repulsions on coiled-coil stability as compared to GdnHCl.

Specificity of dimer formation in tropomyosins: influence of alternatively spliced exons on homodimer and heterodimer assembly

Tropomyosins consist of nearly 100% alpha-helix and assemble into parallel and in-register coiled-coil dimers. In vitro it has been established that nonmuscle as well as native muscle tropomyosins can form homodimers. However, a mixture of muscle alpha and beta tropomyosin subunits results in the formation of the thermodynamically more stable alpha/beta heterodimer. Although the assembly preference of the muscle tropomyosin heterodimer can be understood thermodynamically, the presence of multiple tropomyosin isoforms expressed in nonmuscle cells points toward a more complex principle for determining dimer formation. We have investigated the dimerization of rat tropomyosins in living cells by the use of epitope tagging with a 16-aa sequence of the influenza hemagglutinin. Employing transfection and immunoprecipitation techniques, we have analyzed the dimers formed by muscle and nonmuscle tropomyosins in rat fibroblasts. We demonstrate that the information for homo- versus heterodimerization is contained within the tropomyosin molecule itself and that the information for the selectivity is conferred by the alternatively spliced exons. These results have important implications for models of the regulation of cytoskeletal dynamics.

A calorimetric characterization of the salt dependence of the stability of the GCN4 leucine zipper

The effects of different salts (LiCl, NaCl, ChoCl, KF, KCl, and KBr) on the structural stability of a 33-residue peptide corresponding to the leucine zipper region of GCN4 have been studied by high-sensitivity differential scanning calorimetry. These experiments have allowed an estimation of the salt dependence of the thermodynamic parameters that define the stability of the coiled coil. Independent of the nature of the salt, a destabilization of the coiled coil is always observed upon increasing salt concentration up to a maximum of approximately 0.5 M, depending on the specific cation or anion. At higher salt concentrations, this effect is reversed and a stabilization of the leucine zipper is observed. The effect of salt concentration is primarily entropic, judging from the lack of a significant salt dependence of the transition enthalpy. The salt dependence of the stability of the peptide is complex, suggesting the presence of specific salt effects at high salt concentrations in addition to the nonspecific electrostatic effects that are prevalent at lower salt concentrations. The data is consistent with the existence of specific interactions between anions and peptide with an affinity that follows a reverse size order (F- > Cl- > Br-). Under all conditions studied, the coiled coil undergoes reversible thermal unfolding that can be well represented by a reaction of the form N2<==>2U, indicating that the unfolding is a two-state process in which the helices are only stable when they are in the coiled coil conformation.