Cation-pi interactions in protein-protein interfaces

Arginine is an abundant residue in protein-protein interfaces. The importance of this residue relates to the versatility of its side chain in intermolecular interactions. Different classes of protein-protein interfaces were surveyed for cation-pi interactions. Approximately half of the protein complexes and one-third of the homodimers analyzed were found to contain at least one intermolecular cation-pi pair. Interactions between arginine and tyrosine were found to be the most abundant. The electrostatic interaction energy was calculated to be approximately 3 kcal/mol, on average. A distance-based search of guanidinium:aromatic interactions was also performed using the Macromolecular Structure Database (MSD). This search revealed that half of the guanidinium:aromatic pairs pack in a coplanar manner. Furthermore, it was found that the cationic group of the cation-pi pair is frequently involved in intermolecular hydrogen bonds. In this manner the arginine side chain can participate in multiple interactions, providing a mechanism for inter-protein specificity. Thus, the cation-pi interaction is established as an important contributor to protein-protein interfaces.     

Open-and-shut cases in coiled-coil assembly: alpha-sheets and alpha-cylinders

The coiled coil is a ubiquitous protein-folding motif. It generally is accepted that coiled coils are characterized by sequence patterns known as heptad repeats. Such patterns direct the formation and assembly of amphipathic alpha-helices, the hydrophobic faces of which interface in a specific manner first proposed by Crick and termed "knobs-into-holes packing". We developed software, SOCKET, to recognize this packing in protein structures. As expected, in a trawl of the protein data bank, we found examples of canonical coiled coils with a single contiguous heptad repeat. In addition, we identified structures with multiple, overlapping heptad repeats. This observation extends Crick's original postulate: Multiple, offset heptad repeats help explain assemblies with more than two helices. Indeed, we have found that the sequence offset of the multiple heptad repeats is related to the coiled-coil oligomer state. Here we focus on one particular sequence motif in which two heptad repeats are offset by two residues. This offset sets up two hydrophobic faces separated by approximately 150 degrees -160 degrees around the alpha-helix. In turn, two different combinations of these faces are possible. Either similar or opposite faces can interface, which leads to open or closed multihelix assemblies. Accordingly, we refer to these two forms as alpha-sheets and alpha-cylinders. We illustrate these structures with our own predictions and by reference to natural variants on these designs that have recently come to light.     

Accommodation of a highly symmetric core within a symmetric protein superfold

An alternative core packing group, involving a set of five positions, has been introduced into human acidic FGF-1. This alternative group was designed so as to constrain the primary structure within the core region to the same threefold symmetry present in the tertiary structure of the protein fold (the beta-trefoil superfold). The alternative core is essentially indistinguishable from the WT core with regard to structure, stability, and folding kinetics. The results show that the beta-trefoil superfold is compatible with a threefold symmetric constraint on the core region, as might be the case if the superfold arose as a result of gene duplication/fusion events. Furthermore, this new core arrangement can form the basis of a structural "building block" that can greatly simplify the de novo design of beta-trefoil proteins by using symmetric structural complementarity. Remaining asymmetry within the core appears to be related to asymmetry in the tertiary structure associated with receptor and heparin binding functionality of the growth factor.     

Stoichiometry of lipid interactions with transmembrane proteins--Deduced from the 3D structures

The stoichiometry of the first shell of lipids interacting with a transmembrane protein is defined operationally by the population of spin-labeled lipid chains whose motion is restricted directly by the protein. Interaction stoichiometries have been determined experimentally for a wide range of alpha-helical integral membrane proteins by using spin-label ESR spectroscopy. Here, we determine the spatially defined number of first-shell lipids at the hydrophobic perimeter of integral membrane proteins whose 3D structure has been determined by X-ray crystallography and lipid-protein interactions characterized by spin-labeling. Molecular modeling is used to build a single shell of lipids surrounding transmembrane structures derived from the PDB. Constrained energy optimization of the protein-lipid assemblies is performed by molecular mechanics. For relatively small proteins (up to 7-12 transmembrane helices), the geometrical first shell corresponds to that defined experimentally by perturbation of the lipid-chain dynamics. For larger, multi-subunit alpha-helical proteins, the lipids perturbed directly by the protein may either exceed or be less in number than those that can be accommodated at the intramembranous perimeter. In these latter cases, the motionally restricted spin-labeled lipids can be augmented by intercalation, or can correspond to a specific subpopulation at the protein interface, respectively. For monomeric beta-barrel proteins, the geometrical lipid stoichiometry corresponds to that determined from lipid mobility for a 22-stranded barrel, but fewer lipids are motionally restricted than can be accommodated around an eight-stranded barrel. Deviations from the geometrical first shell, in the beta-barrel case, are for the smaller protein with a highly curved barrel.     

Importance of secondary structural specificity determinants in protein folding: insertion of a native beta-sheet sequence into an alpha-helical coiled-coil

To examine how a short secondary structural element derived from a native protein folds when in a different protein environment, we inserted an 11-residue beta-sheet segment (cassette) from human immunoglobulin fold, Fab new, into an alpha-helical coiled-coil host protein (cassette holder). This de novo design protein model, the structural cassette mutagenesis (SCM) model, allows us to study protein folding principles involving both short- and long-range interactions that affect secondary structure stability and conformation. In this study, we address whether the insertion of this beta-sheet cassette into the alpha-helical coiled-coil protein would result in conformational change nucleated by the long-range tertiary stabilization of the coiled-coil, therefore overriding the local propensity of the cassette to form beta-sheet, observed in its native immunoglobulin fold. The results showed that not only did the nucleating helices of the coiled-coil on either end of the cassette fail to nucleate the beta-sheet cassette to fold with an alpha-helical conformation, but also the entire chimeric protein became a random coil. We identified two determinants in this cassette that prevented coiled-coil formation: (1) a tandem dipeptide NN motif at the N-terminal of the beta-sheet cassette, and (2) the hydrophilic Ser residue, which would be buried in the hydrophobic core if the coiled-coil structure were to fold. By amino acid substitution of these helix disruptive residues, that is, either the replacement of the NN motif with high helical propensity Ala residues or the substitution of Ser with Leu to enhance hydrophobicity, we were able to convert the random coil chimeric protein into a fully folded alpha-helical coiled-coil. We hypothesized that this NN motif is a "secondary structural specificity determinant" which is very selective for one type of secondary structure and may prevent neighboring residues from adopting an alternate protein fold. These sequences with secondary structural specificity determinants have very strong local propensity to fold into a specific secondary structure and may affect overall protein folding by acting as a folding initiation site.     

Cation-pi (Na+-Trp) interactions in the crystal structure of tetragonal lysozyme.

Experimental evidence of a cation-pi interaction between a sodium cation (Na+) and the indole ring of residue Trp123 in a structure (2.0 A) of hen egg-white lysozyme is presented. The geometry of the metal ion-pi interaction observed in the protein structure (distance between the aromatic plane and the cation approximately 4 A) is consistent with geometries observed among small molecules crystal structures and quantum chemistry ab initio calculations. The present crystal structure of lysozyme provides unique structural information about the geometry of binding of cations to pi systems in proteins. It shows that the metal ion-pi interaction within proteins is not significantly different from similar bindings found in small molecules and that it can be modeled by theoretical methods.     

Solvent effects on the conformational transition of a model polyalanine peptide

We have investigated the folding of polyalanine by combining discontinuous molecular dynamics simulation with our newly developed off-lattice intermediate-resolution protein model. The thermodynamics of a system containing a single Ac-KA(14)K-NH(2) molecule has been explored by using the replica exchange simulation method to map out the conformational transitions as a function of temperature. We have also explored the influence of solvent type on the folding process by varying the relative strength of the side-chain's hydrophobic interactions and backbone hydrogen bonding interactions. The peptide in our simulations tends to mimic real polyalanine in that it can exist in three distinct structural states: alpha-helix, beta-structures (including beta-hairpin and beta-sheet-like structures), and random coil, depending upon the solvent conditions. At low values of the hydrophobic interaction strength between nonpolar side-chains, the polyalanine peptide undergoes a relatively sharp transition between an alpha-helical conformation at low temperatures and a random-coil conformation at high temperatures. As the hydrophobic interaction strength increases, this transition shifts to higher temperatures. Increasing the hydrophobic interaction strength even further induces a second transition to a beta-hairpin, resulting in an alpha-helical conformation at low temperatures, a beta-hairpin at intermediate temperatures, and a random coil at high temperatures. At very high values of the hydrophobic interaction strength, polyalanines become beta-hairpins and beta-sheet-like structures at low temperatures and random coils at high temperatures. This study of the folding of a single polyalanine-based peptide sets the stage for a study of polyalanine aggregation in a forthcoming paper.     

Long-range side-chain-main-chain interactions play crucial roles in stabilizing the (betaalpha)8 barrel motif of the alpha subunit of tryptophan synthase

The role of hither-to-fore unrecognized long-range hydrogen bonds between main-chain amide hydrogens and polar side chains on the stability of a well-studied (betaalpha)8, TIM barrel protein, the alpha subunit of tryptophan synthase (alphaTS), was probed by mutational analysis. The F19-D46 and I97-D124 hydrogen bonds link the N terminus of a beta-strand with the C terminus of the succeeding antiparallel alpha-helix, and the A103-D130 hydrogen bond links the N terminus of an alpha-helix with the C terminus of the succeeding antiparallel beta-strand, forming clamps for the respective betaalpha or alphabeta hairpins. The individual replacement of these aspartic acid side chains with alanine leads to what appear to be closely related partially folded structures with significantly reduced far-UV CD ellipticity and thermodynamic stability. Comparisons with the effects of eliminating another main-chain-side-chain hydrogen bond, G26-S33, and two electrostatic side-chain-side-chain hydrogen bonds, D38-H92 and D112-H146, all in the same N-terminal folding unit of alphaTS, demonstrated a unique role for the clamp interactions in stabilizing the native barrel conformation. Because neither the asparagine nor glutamic acid variant at position 46 can completely reproduce the spectroscopic, thermodynamic, or kinetic folding properties of aspartic acid, both size and charge are crucial to its unique role in the clamp hydrogen bond. Kinetic studies suggest that the three clamp hydrogen bonds act in concert to stabilize the transition state leading to the fully folded TIM barrel motif.     

A model of dynamic side-chain--side-chain interactions in the alpha-lactalbumin molten globule

Proteins in the molten globule state contain high levels of secondary structure, as well as a rudimentary, nativelike tertiary topology. Thus, the structural similarity between the molten globule and native proteins may have a significant bearing in understanding the protein-folding problem. To explore the nature of side-chain--side-chain interactions in the alpha-lactalbumin (alpha-LA) molten globule, we determined the effective concentration for formation of the 28--111 disulfide bond in 14 double-mutant proteins, each containing two hydrophobic core residues replaced by alanine. We compared our results with those of single-alanine substitutions using the framework of double-mutant cycle analysis and found that, in the majority of cases, the effects of two alanine substitutions are additive. Based on these results, we propose a model of side-chain-side-chain interactions in the alpha-LA molten globule, which takes into consideration the dynamic nature of this partially folded species.     

The role of context on alpha-helix stabilization: host-guest analysis in a mixed background peptide model.

The helix content of a series of peptides containing single substitutions of the 20 natural amino acids in a new designed host sequence, succinyl-YSEEEEKAKKAXAEEAEKKKK-NH2, has been determined using CD spectroscopy. This host is related to one previously studied, in which triple amino acid substitutions were introduced into a background of Glu-Lys blocks completely lacking alanine. The resulting free energies show that only Ala and Glu- prove to be helix stabilizing, while all other side chains are neutral or destabilizing. This agrees with results from studies of alanine-rich peptide modela, but not the previous Glu-Lys block oligomers in which Leu and Met also stabilize helix. The helix propensity scale derived from the previous block oligomers correlated well with the frequencies of occurrence of different side chains in helical sequences of proteins, whereas the values from the present series do not. The role of context in determining scales of helix propensity values is discussed, and the ability of algorithms designed to predict helix structure from sequence is compared.     

The CXXC motif at the N terminus of an alpha-helical peptide

An active site containing a CXXC motif is always found in the thiol-disulphide oxidoreductase superfamily. A survey of crystal structures revealed that the CXXC motif had a very high local propensity (26.3 +/- 6.2) for the N termini of alpha-helices. A helical peptide with the sequence CAAC at the N terminus was synthesized to examine the helix-stabilizing capacity of the CXXC motif. Circular dichroism was used to confirm the helical nature of the peptide and study behavior under titration with various species. With DTT, a redox potential of E(o) = -230 mV was measured, indicating that the isolated peptide is reducing in nature and similar to native human thioredoxin. The pK(a) values of the individual Cys residues could not be separated in the titration of the reduced state, giving a single transition with an apparent pK(a) of 6.74 (+/-0.06). In the oxidized state, the N-terminal pK(a) is 5.96 (+/-0.05). Analysis of results with the modified helix-coil theory indicated that the disulfide bond stabilized the alpha-helical structure by 0.5 kcal/mol. Reducing the disulfide destabilizes the helix by 0.9 kcal/mol.     

Stability of HAMLET--a kinetically trapped alpha-lactalbumin oleic acid complex

The stability toward thermal and urea denaturation was measured for HAMLET (human alpha-lactalbumin made lethal to tumor cells) and alpha-lactalbumin, using circular dichroism and fluorescence spectroscopy as well as differential scanning calorimetry. Under all conditions examined, HAMLET appears to have the same or lower stability than alpha-lactalbumin. The largest difference is seen for thermal denaturation of the calcium free (apo) forms, where the temperature at the transition midpoint is 15 degrees C lower for apo HAMLET than for apo alpha-lactalbumin. The difference becomes progressively smaller as the calcium concentration increases. Denaturation of HAMLET was found to be irreversible. Samples of HAMLET that have been renatured after denaturation have lost the specific biological activity toward tumor cells. Three lines of evidence indicate that HAMLET is a kinetic trap: (1) It has lower stability than alpha-lactalbumin, although it is a complex of alpha-lactalbumin and oleic acid; (2) its denaturation is irreversible and HAMLET is lost after denaturation; (3) formation of HAMLET requires a specific conversion protocol.     

Effect of the N1 residue on the stability of the alpha-helix for all 20 amino acids

N1 is the first residue in an alpha-helix. We have measured the contribution of all 20 amino acids to the stability of a small helical peptide CH(3)CO-XAAAAQAAAAQAAGY-NH(2) at the N1 position. By substituting every residue into the N1 position, we were able to investigate the stabilizing role of each amino acid in an isolated context. The helix content of each of the 20 peptides was measured by circular dichroism (CD) spectroscopy. The data were analyzed by our modified Lifson-Roig helix-coil theory, which includes the n1 parameter, to find free energies for placing a residue into the N1 position. The rank order for free energies is Asp(-), Ala > Glu(-) > Glu(0) > Trp, Leu, Ser > Asp(0), Thr, Gln, Met, Ile > Val, Pro > Lys(+), Arg, His(0) > Cys, Gly > Phe > Asn, Tyr, His(+). N1 preferences are clearly distinct from preferences for the preceding N-cap and alpha-helix interior. pK(a) values were measured for Asp, Glu, and His, and protonation-free energies were calculated for Asp and Glu. The dissociation of the Asp proton is less favorable than that of Glu, and this reflects its involvement in a stronger stabilizing interaction at the N terminus. Proline is not energetically favored at the alpha-helix N terminus despite having a high propensity for this position in crystal structures. The data presented are of value both in rationalizing mutations at N1 alpha-helix sites in proteins and in predicting the helix contents of peptides.     

The Bacillus licheniformis BlaP beta-lactamase as a model protein scaffold to study the insertion of protein fragments

Using genetic engineering technologies, the chitin-binding domain (ChBD) of the human macrophage chitotriosidase has been inserted into the host protein BlaP, a class A beta-lactamase produced by Bacillus licheniformis. The product of this construction behaved as a soluble chimeric protein that conserves both the capacity to bind chitin and to hydrolyze beta-lactam moiety. Here we describe the biochemical and biophysical properties of this protein (BlaPChBD). This work contributes to a better understanding of the reciprocal structural and functional effects of the insertion on the host protein scaffold and the heterologous structured protein fragments. The use of BlaP as a protein carrier represents an efficient approach to the functional study of heterologous protein fragments.     

Elongation of the BH8 β-hairpin peptide: Electrostatic interactions in β-hairpin formation and stability

An elongated version of the de novo designed beta-hairpin peptide, BH8, has allowed us to gain insight into the role of electrostatic interactions in beta-hairpin stability. A Lys-Glu electrostatic pair has been introduced by adding a residue at the beginning and at the end of the N-terminal and C-terminal strands, respectively, of the beta-hairpin structure, in both orientations. The two resulting peptides and controls having Ala residues at these positions and different combinations of Ala with Lys, or Glu residues, have been analyzed by nuclear magnetic resonance (NMR), under different pH and ionic strength conditions. All of the NMR parameters, in particular the conformational shift analysis of Calpha protons and the coupling constants, (3)J(HNalpha), correlate well and the population estimates are in reasonable agreement among the different methods used. In the most structured peptides, we find an extension of the beta-hairpin structure comprising the two extra residues. Analysis of the pH and salt dependence shows that ionic pairs contribute to beta-hairpin stability. The interaction is electrostatic in nature and can be screened by salt. There is also an important salt-independent contribution of negatively charged groups to the stability of this family of beta-hairpin peptides.     

Solution structure of alpha t alpha, a helical hairpin peptide of de novo design.

alpha t alpha is a 38-residue peptide designed to adopt a helical hairpin conformation in solution (Fezoui Y, Weaver DL Osterhout JJ, 1995, Protein Sci 4:286-295). A previous study of the carboxylate form of alpha t alpha by CD and two-dimensional NMR indicated that the peptide was highly helical and that the helices associated in approximately the intended orientation (Fezoui Y, Weaver DL, Osterhout JJ, 1994, Proc Natl Acad Sci USA 91:3675-3679). Here, the solution structure of alpha t alpha as determined by two-dimensional NMR is reported. A total of 266 experimentally derived distance restraints and 20 dihedral angle restraints derived from J-couplings were used. One-hundred initial structures were generated by distance geometry and refined by dynamical simulated annealing. Twenty-three of the lowest-energy structures consistent with the experimental restraints were analyzed. The results presented here show that alpha t alpha is comprised of two associating helices connected by a turn region.     

Design and characterization of the anion-sensitive coiled-coil peptide.

As a model for analyzing the role of charge repulsion in proteins and its shielding by the solvent, we designed a peptide of 27 amino acid residues that formed a homodimeric coiled-coil. The interface between the coils consisted of hydrophobic Leu and Val residues, and 10 Lys residues per monomer were incorporated into the positions exposed to solvent. During the preparation of a disulfide-linked dimer in which the two peptides were linked in parallel by the two disulfide bonds located at the N and C terminals, a cyclic monomer with an intramolecular disulfide bond was also obtained. On the basis of CD and 1H-NMR, the conformational stabilities of these isomers and several reference peptides were examined. Whereas all these peptides were unfolded in the absence of salt at pH 4.7 and 20 degrees C, the addition of NaClO4 cooperatively stabilized the alpha-helical conformation. The crosslinking of the peptides by disulfide bonds significantly decreased the midpoint salt concentration of the transition. The 1H-NMR spectra in the presence of NaClO4 suggested that, whereas the disulfide-bonded dimer assumed a native-like conformation, the cyclic monomer assumed a molten globule-like conformation with disordered side chains. However, the cyclic monomer exhibited cooperative transitions against temperature and Gdn-HCl that were only slightly less cooperative than those of the disulfide-bonded parallel dimer. These results indicate that the charge repulsion critically destabilizes the native-like state as well as the molten globule-like state, and that the solvent-dependent charge repulsion may be useful for controlling the conformation of designed peptides.     

Effect of hexafluoroisopropanol alcohol on the structure of melittin: a molecular dynamics simulation study

The molecular mechanism by which HFIP stabilizes the alpha-helical structure of peptides is not well understood. In the present study, we use melittin as a model to gain insight into the details of the atomistic interactions of HFIP with the peptide. We have performed extensive comparative molecular dynamics simulations (up to 100 nsec) in the absence and in the presence of HFIP. In agreement with recent NMR experiments, the simulations show rapid loss of tertiary structure in water at pH 2 but much higher helicity in 35% HFIP. The MD simulations also indicate that melittin adopts a highly dynamic global structure in 35% HFIP solution with two alpha-helical segments sampling a wide range of angular orientations. The analysis of the HFIP distribution shows the tendency of HFIP to aggregate around the peptide, increasing the local cosolvent concentration to more than two times that in the bulk concentration. The correlation of local peptide structure with HFIP coating suggests that displacement of water at the peptide surface is the main contribution of HFIP in stabilizing the secondary structure of melittin. Finally, a stabilizing effect promoted by the presence of counter-ions was also observed in the simulations.     

Effect of the N2 residue on the stability of the α-helix for all 20 amino acids

N2 is the second position in the alpha-helix. All 20 amino acids were placed in the N2 position of a synthetic helical peptide (CH(3)CO-[AXAAAAKAAAAKAAGY]-NH(2)) and the helix content was measured by circular dichroism spectroscopy at 273K. The dependence of peptide helicity on N2 residue identity has been used to determine a free-energy scale by analysis with a modified Lifson-Roig helix-coil theory that includes a parameter for the N2 energy (n2). The rank order of DeltaDeltaG((relative to Ala)) is Glu(-), Asp(-) > Ala > Glu(0), Leu, Val, Gln, Thr, Ile, Ser, Met, Asp(0), His(0), Arg, Cys, Lys, Phe > Asn, > Gly, His(+), Pro, Tyr. The results correlate very well with N2 propensities in proteins, moderately well with N1 and helix interior preferences, and not at all with N-cap preferences. The strongest energetic effects result from interactions with the helix dipole, which favors negative charges at the helix N terminus. Hydrogen bonds to side chains at N2, such as Gln, Ser, and Thr, are weak, despite occurring frequently in protein crystal structures, in contrast to the N-cap position. This is because N-cap hydrogen bonds are close to linear, whereas N2 hydrogen bonds have poor geometry. These results can be used to modify protein stability rationally, help design helices, and improve prediction of helix location and stability.     

Effect of the N3 residue on the stability of the alpha-helix

N3 is the third position from the N terminus in the alpha-helix with helical backbone dihedral angles. All 20 amino acids have been placed in the N3 position of a synthetic helical peptide (CH(3)CO-[AAX AAAAKAAAAKAGY]-NH(2)) and the helix content measured by circular dichroism spectroscopy at 273 K. The dependence of peptide helicity on N3 residue identity has been used to determine a free energy scale by analysis with a modified Lifson-Roig helix coil theory that includes a parameter for the N3 energy (n3). The most stabilizing residues at N3 in rank order are Ala, Glu, Met/Ile, Leu, Lys, Ser, Gln, Thr, Tyr, Phe, Asp, His, and Trp. Free energies for the most destabilizing residues (Cys, Gly, Asn, Arg, and Pro) could not be fitted. The results correlate with N1, N2, and helix interior energies and not at all with N-cap preferences. This completes our work on studying the structural and energetic preferences of the amino acids for the N-terminal positions of the alpha-helix. These results can be used to rationally modify protein stability, help design helices, and improve prediction of helix location and stability.