Analysis of Interactions of Buried Polar Side Chains in β-Proteins

Qualitative and quantitative analysis of polar side chains inaccessible to water molecules, as well as their interactions in 100 globular β-sheet proteins, was performed. It was shown that completely buried polar side chains are widespread in β-proteins, with their vast majority being involved in side chain-side chain or side chain-main chain interactions. An analysis of frequency of occurrence of different side chain-partner pairs demonstrated that these interactions are selective. The results were compared with similar data obtained earlier for α-helical proteins.     

A Systematic Analysis of Buried Polar Side Chains and Their Interactions in α-Helical Proteins

Examination of 80 -helical proteins and domains demonstrates that they contain from 1 to more than 20 completely buried (water-inaccessible) polar side chains. As a rule the latter have partners for H-bonding but the resulting H-bond system is often not saturating. Basing on statistical analysis, we determined the optimal number of H-bonds for every type of polar side chain, and discuss the structural role of vacant donors and acceptors. About half of the H-bonds formed by buried side chains pertain to interhelix contacts of the (side chain)–(side chain) and (side chain)–(main chain) types. Such interactions appear to be a most important factor determining the mutual arrangement of -helices in proteins. Analysis of the frequency of occurrence of various interacting pairs reveals that these interactions are selective.     

Systematic analysis of buried polar side chains and their interactions in alpha-helical proteins

Examination of 80 alpha-helical proteins and domains demonstrates that they contain from 1 to more than 20 completely buried (water-inaccessible) polar side chains. As a rule the latter have partners for H-bonding but the resulting H-bond system is often not exhaustive. Basing on statistical analysis, we determined the optimal number of H-bonds for every type of polar side chain, and discuss the structural role of vacant donors and acceptors. About half of the H-bonds formed by buried side chains pertain to interhelix contacts of the (side chain)-(side chain) and (side chain)-(main chain) types. Such interactions appear to be a most important factor determining the mutual arrangement of alpha-helices in proteins. Analysis of the frequency of occurrence of various interacting pairs reveals that these interactions are selective.     

Patterns in ionizable side chain interactions in protein structures

In a selected set of 44 high-resolution, non-homologous protein structures, the intramolecular hydrogen bonds or salt bridges formed by ionizable amino acid side chains were identified and analyzed. The analysis was based on the investigation of several properties of the involved residues such as their solvent exposure, their belonging to a certain secondary structural element, and their position relative to the N- and C-termini of their respective structural element. It was observed that two-thirds of the interactions made by basic or acidic side chains are hydrogen bonds to polar uncharged groups. In particular, the majority (78%) of the hydrogen bonds between ionizable side chains and main chain polar groups (sch:mch bonds) involved at least one buried atom, and in 42% of the cases both interacting atoms were buried. In α-helices, the sch:mch bonds observed in the proximity of the C- and N-termini show a clear preference for acidic and basic side chains, respectively. This appears to be due to the partial charges of peptide group atoms at the termini of α-helices, which establish energetically favorable electrostatic interactions with side chain carrying opposite charge, at distances even greater than 4.5 Å. The sch:mch interactions involving ionizable side chains that belong either to β-strands or to the central part of α-helices are based almost exclusively on basic residues. This results from the presence of main chain carbonyl oxygen atoms in the protein core which have unsatisfied hydrogen bonding capabilities.     

Position-dependence of stabilizing polar interactions of asparagine in transmembrane helical bundles

Recent studies with model peptides and statistical analyses of the crystal structures of membrane proteins have shown that buried polar interactions contribute significantly to the stabilization of the three-dimensional structures of membrane proteins. Here, we probe how the location of these polar groups along the transmembrane helices affect their free energies of interaction. Asn residues were placed singly and in pairs at three positions within a model transmembrane helix, which had previously been shown to support the formation of trimers in micelles. The model helix was designed to form a transmembrane coiled coil, with Val side chains at the "a" positions of the heptad repeat. Variants of this peptide were prepared in which an Asn residue was introduced at one or more of the "a" positions, and their free energies of association were determined by analytical ultracentrifugation. When placed near the middle of the transmembrane helix, the formation of trimers was stabilized by at least -2.0 kcal/mol per Asn side chain. When the Asn was placed at the interface between the hydrophobic and polar regions of the peptide, the substitution was neither stabilizing nor destabilizing (0.0 +/- 0.5 kcal/mol of monomer). Finally, it has previously been shown that a Val-for-Asn mutation in a water-soluble coiled coil destabilizes the structure by approximately 1.5 kcal/mol of monomer [Acharya, A., et al. (2002) Biochemistry 41, 14122-14131]. Thus, the headgroup region of a micelle appears to have a conformational impact intermediate between that of bulk water and the apolar region of micelle. A similarly large dependence on the location of the polar residues was found in a statistical survey of helical transmembrane proteins. The tendency of different types of residues to be buried in the interiors versus being exposed to lipids was analyzed. Asn and Gln show a very strong tendency to be buried when they are located near the middle of a transmembrane helix. However, when placed near the ends of transmembrane helices, they show little preference for the surface versus the interior of the protein. These data show that Asn side chains within the apolar region of the transmembrane helix provide a significantly larger driving force for association than Asn residues near the apolar/polar interface. Thus, although polar interactions are able to strongly stabilize the folding of membrane proteins, the energetics of association depend on their location within the hydrophobic region of a transmembrane helix.     

Extension of the fragment method to calculate amino acid zwitterion and side chain partition coefficients

The fragment method of calculating partition coefficients (P) has been extended to include the common amino acids (AAs). The results indicate that polar and charged side chains influence the hydrophobicity of atoms in the side chain in a predictable manner. Field effects, as evidenced through polar proximity factors and bond factors, need to be considered for accurate estimation of transfer phenomena. The calculated log P and delta G degree ' values of the 20 AAs agree well with the observed values. Pro calculates to be more hydrophilic than the observed log P. Hydrophobicity scales for peptide side chain residues are compared and evaluated in terms of suitability. Calculated pi values for nonpolar side chain residues agree well with the observed values; calculated values for uncharged polar side chain residues deviate by about 0.6 log units except for Gln and Cys; and polar side chain residues with charged side chains calculate as too hydrophilic. Reasons for the differences are explored. We also suggest that tightly bound water to polar moieties in amino acids and peptides may be transferred into the octanol phase during partitioning experiments. A quantitative methodology is presented which characterizes the thermodynamic partitioning of groups and individual atoms in amino acids and proteins.     

Structural basis of neurophysin hormone specificity: Geometry, polarity, and polarizability in aromatic ring interactions

The structural origins of the specificity of the neurophysin hormone-binding site for an aromatic residue in peptide position 2 were explored by analyzing the binding of a series of peptides in the context of the crystal structure of liganded neurophysin. A new modeling method for describing the van der Waals surface of binding sites assisted in the analysis. Particular attention was paid to the unusually large (5 kcal/mol) difference in binding free energy between Phe and Leu in position 2, a value representing more than three times the maximum expected based on hydrophobicity alone, and additionally remarkable since modeling indicated that the Leu side chain was readily accommodated by the binding pocket. Although evidence was obtained of a weak thermodynamic linkage between the binding interactions of the residue 2 side chain and of the peptide alpha-amino group, two factors are considered central. (1) The bound Leu side chain can establish only one-third of the van der Waals contacts available to a Phe side chain. (2) The bound Phe side chain appears to be additionally stabilized relative to Leu by more favorable dipole and induced dipole interactions with nonaromatic polar and sulfur ligands in the binding pocket, as evidenced by examination of its interactions in the pocket, analysis of the detailed energetics of transfer of Phe and Leu side chains from water to other phases, and comparison with thermodynamic and structural data for the binding of residue 1 side chains in this system. While such polar interactions of aromatic rings have been previously observed, the present results suggest their potential for significant thermodynamic contributions to protein structure and ligand recognition.     

One- and two-body decomposable Poisson-Boltzmann methods for protein design calculations

Successfully modeling electrostatic interactions is one of the key factors required for the computational design of proteins with desired physical, chemical, and biological properties. In this paper, we present formulations of the finite difference Poisson-Boltzmann (FDPB) model that are pairwise decomposable by side chain. These methods use reduced representations of the protein structure based on the backbone and one or two side chains in order to approximate the dielectric environment in and around the protein. For the desolvation of polar side chains, the two-body model has a 0.64 kcal/mol RMSD compared to FDPB calculations performed using the full representation of the protein structure. Screened Coulombic interaction energies between side chains are approximated with an RMSD of 0.13 kcal/mol. The methods presented here are compatible with the computational demands of protein design calculations and produce energies that are very similar to the results of traditional FDPB calculations.     

Stably folded de novo proteins from a designed combinatorial library

Binary patterning of polar and nonpolar amino acids has been used as the key design feature for constructing large combinatorial libraries of de novo proteins. Each position in a binary patterned sequence is designed explicitly to be either polar or nonpolar; however, the precise identities of these amino acids are varied extensively. The combinatorial underpinnings of the "binary code" strategy preclude explicit design of particular side chains at specified positions. Therefore, packing interactions cannot be specified a priori. To assess whether the binary code strategy can nonetheless produce well-folded de novo proteins, we constructed a second-generation library based upon a new structural scaffold designed to fold into 102-residue four-helix bundles. Characterization of five proteins chosen arbitrarily from this new library revealed that (1) all are alpha-helical and quite stable; (2) four of the five contain an abundance of tertiary interactions indicative of well-ordered structures; and (3) one protein forms a well-folded structure with native-like features. The proteins from this new 102-residue library are substantially more stable and dramatically more native-like than those from an earlier binary patterned library of 74-residue sequences. These findings demonstrate that chain length is a crucial determinant of structural order in libraries of de novo four-helix bundles. Moreover, these results show that the binary code strategy--if applied to an appropriately designed structural scaffold--can generate large collections of stably folded and/or native-like proteins.     

Single-channel studies on linear gramicidins with altered amino acid side chains. Effects of altering the polarity of the side chain at position 1 in gramicidin A.

The modulation of gramicidin A single-channel characteristics by the amino acid side chains was investigated using gramicidin A analogues in which the NH2 terminal valine was chemically replaced by other amino acids. The replacements were chosen such that pairs of analogues would have essentially isosteric side chains of different polarities at position 1 (valine vs. trifluorovaline or hexafluorovaline; norvaline vs. S-methyl-cysteine; and norleucine vs. methionine). Even though the side chains are not in direct contact with the permeating ions, the single-channel conductances for Na+ and Cs+ are markedly affected by the changes in the physico-chemical characteristics of the side chains. The maximum single-channel conductance for Na+ is decreased by as much as 10-fold in channels formed by analogues with polar side chains at position 1 compared with their counterparts with nonpolar side chains, while the Na+ affinity is fairly insensitive to these changes. The relative conductance changes seen with Cs+ were less than those seen with Na+; the ion selectivity of the channels with polar side chains at position 1 was increased. Hybrid channels could form between compounds with a polar side chain at position 1 and either valine gramicidin A or their counterparts with a nonpolar side chain at position 1. The structure of channels formed by the modified gramicidins is thus essentially identical to the structure of channels formed by valine gramicidin A. The polarity of the side chain at position 1 is an important determinant of the permeability characteristics of the gramicidin A channel. We discuss the importance of having structural information when interpreting the functional consequences of site-directed amino acid modifications.     

Alpha helix capping in synthetic model peptides by reciprocal side chain–main chain interactions: Evidence for an N terminal “capping box”

A significant fraction of the amino acids in proteins are alpha helical in conformation. Alpha helices in globular proteins are short, with an average length of about twelve residues, so that residues at the ends of helices make up an important fraction of all helical residues. In the middle of a helix, H-bonds connect the NH and CO groups of each residue to partners four residues along the chain. At the ends of a helix, the H-bond potential of the main chain remains unfulfilled, and helix capping interactions involving bonds from polar side chains to the NH or CO of the backbone have been proposed and detected. In a study of synthetic helical peptides, we have found that the sequence Ser-Glu-Asp-Glu stabilizes the alpha helix in a series of helical peptides with consensus sequences. Following the report by Harper and Rose, which identifies SerXaaXaaGlu as a member of a class of common motifs at the N termini of alpha helices in proteins that they refer to as “capping boxes,” we have reexamined the side chain–main chain interactions in a varient sequence using ¹H NMR, and find that the postulated reciprocal side chain-backbone bonding between the first Ser and last Glu side chains and their peptide NH partners can be resolved: Deletion of two residues N terminal to the Ser-Glu-Asp-Glu sequence in these peptides has no effect on the initiation of helical structure, as defined by two-dimensional (2D) NMR experiments on this variant. Thus the capping box sequence Ser-Glu-Asp-Glu inhibits N terminal fraying of the N terminus of alpha helix in these peptides, and shows the side chain–main chain interactions proposed by Harper and Rose. It thus acts as a helix initiating signal. Since normal a helix cannot propagate beyond the N terminus of this structure, the box acts as a termination signal in this direction as well. © 1994 John Wiley & Sons, Inc.     

High resolution X-ray studies of mammalian intestinal and muscle fatty acid-binding proteins provide an opportunity for defining the chemical nature of fatty acid: protein interactions

The structure ofE. coli-derived rat intestinal fatty acid-binding protein has recently been refined to 1.2 Å without bound fatty acid and to 2.0 Å and 1.75 Å with bound hexadecanoate (palmitate) and 9Z-octadecenoate (oleate), respectively. The structure ofE. coli-derived human muscle fatty acid-binding protein has also been solved to 2.1 Å with a C16 bacterial fatty acid. Both proteins contain 10 anti-parallel -strands in a+1, +1, +1... motif. The strands are arranged in two -pleated sheets that are orthogonally oriented. In each case, the fatty acid is enclosed by the -sheets and is bound to the proteins by feeble forces. These feeble forces consist of (i) a hydrogen bonding network between the fatty acid's carboxylate group, ordered solvent, and side chains of polar/ionizable amino acid residues; (ii) van der Waals contacts between the methylene chain of the fatty acid and the side chain atoms of hydrophobic and aromatic residues; (iii) van der Waals interactions between the methyl and the component methenyls of the phenyl side chain of a Phe which serves as an adjustable terminal sensor situated over a surface opening or portal connecting interior and exterior solvent; and (iv) van der Waals contacts between methylenes of the alkyl chain and oxygens of ordered waters that have been located inside the binding cavity. These waters are positioned over one face of the ligand and are held in place by hydrogen bonding with one another and with the side chains of protein's polar and ionizable residues. Binding of the fatty acid ligand is associated with minimal adjustments of the positions of main chain or side chain atoms. However, acquisition of ligand is associated with removal of ordered interior solvent suggesting that the free energy of dehydration of the binding site may be as important for the energy of the binding reaction as the free energy of stabilization of the fatty acid: protein complex.     

Methodological advancements for characterising protein side chains by NMR spectroscopy

The surface of proteins is covered by side chains of polar amino acids that are imperative for modulating protein functionality through the formation of noncovalent intermolecular interactions. However, despite their tremendous importance, the unique structures of protein side chains require tailored approaches for investigation by nuclear magnetic resonance spectroscopy and so have traditionally been understudied compared with the protein backbone. Here, we review substantial recent methodological advancements within nuclear magnetic resonance spectroscopy to address this issue. Specifically, we consider advancements that provide new insight into methyl-bearing side chains, show the potential of using non-natural amino acids and reveal the actions of charged side chains. Combined, the new methods promise unprecedented characterisations of side chains that will further elucidate protein function.     

Electrostatic interactions in aromatic oligopeptides contribute to protein stability

Recent X-ray crystallographic studies of aromatic oligopeptides have shown that aromatic amino acid side chains participate in enthalpically-favorable, weakly polar interactions that stabilize oligopeptide folds. These interactions are important in peptides used as model therapeutic agents for sickle-cell disease, in vasopressin (antidiuretic hormone) and in [Leu]-enkephalin. The aromatic groups of globular proteins display similar behavior and thereby contribute to the stability of the three-dimensional structure of proteins.     

Structural basis for the unusual specificity of Escherichia coli aminopeptidase N

Aminopeptidase N from Escherichia coli is a M1 class aminopeptidase with the active-site region related to that of thermolysin. The enzyme has unusual specificity, cleaving adjacent to the large, nonpolar amino acids Phe and Tyr but also cleaving next to the polar residues Lys and Arg. To try to understand the structural basis for this pattern of hydrolysis, the structure of the enzyme was determined in complex with the amino acids L-arginine, L-lysine, L-phenylalanine, L-tryptophan, and L-tyrosine. These amino acids all bind with their backbone atoms close to the active-site zinc ion and their side chain occupying the S1 subsite. This subsite is in the form of a cylinder, about 10 A in cross-section and 12 A in length. The bottom of the cylinder includes the zinc ion and a number of polar side chains that make multiple hydrogen-bonding and other interactions with the alpha-amino group and the alpha-carboxylate of the bound amino acid. The walls of the S1 cylinder are hydrophobic and accommodate the nonpolar or largely nonpolar side chains of Phe and Tyr. The top of the cylinder is polar in character and includes bound water molecules. The epsilon-amino group of the bound lysine side chain and the guanidinium group of arginine both make multiple hydrogen bonds to this part of the S1 site. At the same time, the hydrocarbon part of the lysine and arginine side chains is accommodated within the nonpolar walls of the S1 cylinder. This combination of hydrophobic and hydrophilic binding surfaces explains the ability of ePepN to cleave Lys, Arg, Phe, and Tyr. Another favored substrate has Ala at the P1 position. The short, nonpolar side chain of this residue can clearly be bound within the hydrophobic part of the S1 cylinder, but the reason for its facile hydrolysis remains uncertain.     

Calculations of the conformations of Iturin A in relation with NMR studies

Iturin A is an antifungal antibiotic which was isolated from a strain of Bacillus subtilis, and contains a lipophilic beta amino acid closing an heptapeptide cycle with polar L and D residues. Iturin A belongs to a lipopeptide family of which the LDDLLDL sequence is kept constant. NMR spectroscopy and semi-empirical energy calculations are combined to design the conformations of Iturin A in pyridine solution. J coupling constants and nOes (nuclear Overhauser enhancements) are used as guiding line for energy calculations. This preliminary study shows that Iturin A in pyridine appears as rather rigid, especially in the L Pro 5-D Asn 6 region, probably involved in a beta turn. The polar side chains can form different networks of intramolecular hydrogen bonds. The Tyr side chain, relatively mobile, could be involved in interactions with an hydrophobic environment as the beta amino acid side chain found away from the peptide cycle.     

The transmembrane 7-alpha-bundle of rhodopsin: distance geometry calculations with hydrogen bonding constraints.

A 3D model of the transmembrane 7-alpha-bundle of rhodopsin-like G-protein-coupled receptors (GPCRs) was calculated using an iterative distance geometry refinement with an evolving system of hydrogen bonds, formed by intramembrane polar side chains in various proteins of the family and collectively applied as distance constraints. The alpha-bundle structure thus obtained provides H bonding of nearly all buried polar side chains simultaneously in the 410 GPCRs considered. Forty evolutionarily conserved GPCR residues form a single continuous domain, with an aliphatic "core" surrounded by six clusters of polar and aromatic side chains. The 7-alpha-bundle of a specific GPCR can be calculated using its own set of H bonds as distance constraints and the common "average" model to restrain positions of the helices. The bovine rhodopsin model thus determined is closely packed, but has a few small polar cavities, presumably filled by water, and has a binding pocket that is complementary to 11-cis (6-s-cis, 12-s-trans, C = N anti)-retinal or to all-trans-retinal, depending on conformations of the Lys296 and Trp265 side chains. A suggested mechanism of rhodopsin photoactivation, triggered by the cis-trans isomerization of retinal, involves rotations of Glu134, Tyr223, Trp265, Lys296, and Tyr306 side chains and rearrangement of their H bonds. The model is in agreement with published electron cryomicroscopy, mutagenesis, chemical modification, cross-linking, Fourier transform infrared spectroscopy, Raman spectroscopy, electron paramagnetic resonance spectroscopy, NMR, and optical spectroscopy data. The rhodopsin model and the published structure of bacteriorhodopsin have very similar retinal-binding pockets.     

Use of a potential of mean force to analyze free energy contributions in protein folding.

A method for calculation of the free energy of residues as a function of residue burial is proposed. The method is based on the potential of mean force, with a reaction coordinate expressed by residue burial. Residue burials are calculated from high-resolution protein structures. The largest individual contributions to the free energy of a residue are found to be due to the hydrophobic interactions of the nonpolar atoms, interactions of the main chain polar atoms, and interactions of the charged groups of residues Arg and Lys. The contribution to the free energy of folding due to the uncharged side chain polar atoms is small. The contribution to the free energy of folding due to the main chain polar atoms is favorable for partially buried residues and less favorable or unfavorable for fully buried residues. Comparison of the accessible surface areas of proteins and model spheres shows that proteins deviate considerably from a spherical shape and that the deviations increase with the size of a protein. The implications of these results for protein folding are also discussed.     

Structural stability and dynamics of an amyloid-forming peptide GNNQQNY from the yeast prion sup-35

A seven amino acid yeast prion sup-35 fragment (GNNQQNY) forms amyloid fibrils. The availability of its detailed atomic oligomeric structure makes it a good model for studying the early stage of aggregation. Here we perform long all-atom explicit solvent molecular simulations of various sizes and arrangements of oligomer seeds of the wild-type and its mutants to study its stability and dynamics. Previous studies have suggested that the early stage rate-limiting step of oligomer formation occurs in high-order oligomers. Our simulations show that with the increase in the number of strands even from a dimer to a trimer, oligomer stability increases dramatically. This suggests that the minimal nucleus seed for GNNQQNY fibril formation could be small and is likely three or four peptides, in agreement with experiment, and that higher-order oligomers do not dissociate quickly since they have small diffusion coefficients and thus slow kinetics. Further, for the hydrophilic polar GNNQQNY, there are no hydrogen bonds and no hydrophobic interactions between adjacent beta-sheets. Simulations suggest that within the sheet, the driving forces to associate and stabilize are interstrand backbone-backbone and side chain-side chain hydrogen bonds, whereas between the sheets, shape-complementary by the dry polar steric zipper via the side chains of Asn-2, Gln-4, and Asn-6 holds the sheets together, as proposed in an earlier study. Since the polar side chains of Asn-2, Gln-4, and Asn-6 act as a hook to bind two neighboring sheets together, these geometric restraints reduce the conformational search for the correct side chain packing to a two-dimensional problem of intersheet side chain interactions. Mutant simulations show that substitution of Asn-2, Gln-4, or Asn-6 by Ala would disrupt this steric zipper, leading to unstable oligomers.