Core side-chain packing and backbone conformation in Lpp-56 coiled-coil mutants

Native proteins exhibit precise geometric packing of atoms in their hydrophobic interiors. Nonetheless, controversy remains about the role of core side-chain packing in specifying and stabilizing the folded structures of proteins. Here we investigate the role of core packing in determining the conformation and stability of the Lpp-56 trimerization domain. The X-ray crystal structures of Lpp-56 mutants with alanine substitutions at two and four interior core positions reveal trimeric coiled coils in which the twist of individual helices and the helix-helix spacing vary significantly to achieve the most favored superhelical packing arrangement. Introduction of each alanine "layer" into the hydrophobic core destabilizes the superhelix by 1.4 kcal mol(-1). Although the methyl groups of the alanine residues pack at their optimum van der Waals contacts in the coiled-coil trimer, they provide a smaller component of hydrophobic interactions than bulky hydrophobic side-chains to the thermodynamic stability. Thus, specific side-chain packing in the hydrophobic core of coiled coils are important determinants of protein main-chain conformation and stability.     

Design, synthesis and structural characterization of model heterodimeric coiled-coil proteins.

We report the design and synthesis of model heterodimeric coiled-coil proteins and the packing contribution of interchain hetero-hydrophobic side-chains to coiled-coil stability. The heterodimeric coiled-coils are obtained by oxidizing two 35-residue polypeptide chains, each containing a cysteine residue at position 2 and differing in amino acid sequences in the hydrophobic positions ("a" and "d") responsible for the formation and stabilization of the coiled-coil. In each peptide, a single Ala residue was substituted for Leu at position "a" or "d". The formation and stability of heterodimeric coiled-coils were investigated by circular dichroism studies in the presence and absence of guanidine hydrochloride and compared to the corresponding homodimeric coiled-coils. The coiled-coil proteins with an Ala substitution at position "a" were less stable than those with an Ala substitution at position "d" in both the homodimeric (Ala-Ala interchain interactions) and heterodimeric (Leu-Ala interchain interactions ) coiled-coils. The 70-residue disulfide bridged peptides (homo- and heterodimeric coiled-coils) can be readily separated by reversed-phase chromatography (RPC) even though they have identical amino acid compositions as well as in the hydrophobic "a" and "d" positions. The elution of the 70-residue peptides prior to their corresponding 35-residue monomers suggests that these proteins are retaining a large portion of their coiled-coil structure during RPC at pH2 and their retention behavior correlates with protein stability.     

A potential molecular switch in an alpha-helical coiled coil

The yeast DNA-binding protein GCN4 forms a homo-dimer through a self-complementary coiled-coil interface. In this article, we describe how such coiled-coils might be bistable and, through Molecular Dynamics computations on the GCN4 coiled coil, we show that the coiled coil can indeed switch between the two states by a pathway in which there is a progressive "flipping" of consecutive steps along the interface. We discuss the general implications of potentially bistable coiled-coil interfaces for allosteric signal-transmission mechanisms along homo-dimeric coiled coils and for the packing of helices in globular proteins.     

Helix packing moments reveal diversity and conservation in membrane protein structure

Helical membrane proteins are more tightly packed and the packing interactions are more diverse than those found in helical soluble proteins. Based on a linear correlation between amino acid packing values and interhelical propensity, we propose the concept of a helix packing moment to predict the orientation of helices in helical membrane proteins and membrane protein complexes. We show that the helix packing moment correlates with the helix interfaces of helix dimers of single pass membrane proteins of known structure. Helix packing moments are also shown to help identify the packing interfaces in membrane proteins with multiple transmembrane helices, where a single helix can have multiple contact surfaces. Analyses are described on class A G protein-coupled receptors (GPCRs) with seven transmembrane helices. We show that the helix packing moments are conserved across the class A family of GPCRs and correspond to key structural contacts in rhodopsin. These contacts are distinct from the highly conserved signature motifs of GPCRs and have not previously been recognized. The specific amino acid types involved in these contacts, however, are not necessarily conserved between subfamilies of GPCRs, indicating that the same protein architecture can be supported by a diverse set of interactions. In GPCRs, as well as membrane channels and transporters, amino acid residues with small side-chains (Gly, Ala, Ser, Cys) allow tight helix packing by mediating strong van der Waals interactions between helices. Closely packed helices, in turn, facilitate interhelical hydrogen bonding of both weakly polar (Ser, Thr, Cys) and strongly polar (Asn, Gln, Glu, Asp, His, Arg, Lys) amino acid residues. We propose the use of the helix packing moment as a complementary tool to the helical hydrophobic moment in the analysis of transmembrane sequences.     

Crystal structure of glucoamylase from Aspergillus awamori var. X100 to 2.2-A resolution.

The crystal structure of a catalytically active fragment of glucoamylase-I from Aspergillus awamori var. X100 has been determined to a resolution of 2.2 A. Twelve of its 13 alpha-helices are arranged into an "alpha/alpha-barrel." An inner core of six mutually parallel alpha-helices are connected to each other through a peripheral set of six alpha-helices. The peripheral helices are parallel to each other, but approximately antiparallel to the inner core of alpha-helices. The putative active site lies in the packing void of the inner set of helices. The last 30 residues of the enzyme comprise a separate domain containing 10 sites of O-glycosylation. Each instance of O-glycosylation involves a serine or threonine side chain linked to the alpha-anomer of a single mannosyl residue. The O-glycosylated domain is in an extended conformation, wrapping around the "waist" of the alpha/alpha-barrel. Two additional sites of N-glycosylation contribute well ordered glycosyl chains that lie in proximity to the belt of O-glycosylation. The model developed for glucoamylase is a rare and valuable structural example of a glycoprotein and an exo-acting amylolytic enzyme.     

X-ray structure of a maquette scaffold

Maquettes are de novo designed mimicries of nature used to test the construction and engineering criteria of oxidoreductases. One type of scaffold used in maquette construction is a four-alpha-helical bundle. The sequence of the four-alpha-helix bundle maquettes follows a heptad repeat pattern typical of left-handed coiled-coils. Initial designs were molten globular due partly to the minimalist approach taken by the designers. Subsequent iterative redesign generated several structured scaffolds with similar heme binding properties. Variant [I(6)F(13)](2), a structured scaffold, was partially resolved with NMR spectroscopy and found to have a set of mobile inter-helical packing interfaces. Here, the X-ray structure of a similar peptide ([I(6)F(13)M(31)](2) i.e. ([CGGG EIWKL HEEFLKK FEELLKL HEERLKKM](2))(2) which we call L31M), has been solved using MAD phasing and refined to 2.8A resolution. The structure shows that the maquette scaffold is an anti-parallel four-helix bundle with "up-up-down-down" topology. No pre-formed heme-binding pocket exists in the protein scaffold. We report unexpected inter-helical crossing angles, residue positions and translations between the helices. The crossing angles between the parallel helices are -5 degrees rather than the expected +20 degrees for typical left-handed coiled-coils. Deviation of the scaffold from the design is likely due to the distribution and size of hydrophobic residues. The structure of L31M points out that four identical helices may interact differently in a bundle and heptad repeats with an alternating [HPPHHPP]/[HPPHHPH] (H: hydrophobic, P: polar) pattern are not a sufficient design criterion to generate left-hand coiled-coils.     

Periodicity of amide proton exchange rates in a coiled-coil leucine zipper peptide.

The two-stranded coiled-coil motif, which includes leucine zippers, is a simple protein structure that is well suited for studies of helix-helix interactions. The interaction between helices in a coiled coil involves packing of "knobs" into "holes", as predicted by Crick in 1953 and confirmed recently by X-ray crystallography for the GCN4 leucine zipper [O'Shea, E.K., Klemm, J.D., Kim, P.S., & Alber, T. (1991) Science 254, 539]. A striking periodicity, extending over six helical turns, is observed in the rates of hydrogen-deuterium exchange for amide protons in a peptide corresponding to the leucine zipper of GCN4. Protons at the hydrophobic interface show the most protection from exchange. The NMR chemical shifts of amide protons in the helices also show a pronounced periodicity which predicts a short H-bond followed by a long H-bond every seven residues. This variation was anticipated in 1953 by Pauling and is sufficient to give rise to a local left-handed superhelical twist characteristic of coiled coils. The amide protons that lie at the base of the "hole" in the "knobs-into-holes" packing show slow amide proton exchange rates and are predicted to have short H-bond lengths. These results suggest that tertiary interactions can lead to highly localized, but substantial, differences in stability and dynamics within a secondary structure element and emphasize the dominant nature of packing interactions in determining protein structure.     

An amino acid code to define a protein's tertiary packing surface

One difficult aspect of the protein‐folding problem is characterizing the nonspecific interactions that define packing in protein tertiary structure. To better understand tertiary structure, this work extends the knob‐socket model by classifying the interactions of a single knob residue packed into a set of contiguous sockets, or a pocket made up of 4 or more residues. The knob‐socket construct allows for a symbolic two‐dimensional mapping of pockets. The two‐dimensional mapping of pockets provides a simple method to investigate the variety of pocket shapes to understand the geometry of protein tertiary surfaces. The diversity of pocket geometries can be organized into groups of pockets that share a common core, which suggests that some interactions in pockets are ancillary to packing. Further analysis of pocket geometries displays a preferred configuration that is right‐handed in α‐helices and left‐handed in β‐sheets. The amino acid composition of pockets illustrates the importance of nonpolar amino acids in packing as well as position specificity. As expected, all pocket shapes prefer to pack with hydrophobic knobs; however, knobs are not selective for the pockets they pack. Investigating side‐chain rotamer preferences for certain pocket shapes uncovers no strong correlations. These findings allow a simple vocabulary based on knobs and sockets to describe protein tertiary packing that supports improved analysis, design, and prediction of protein structure. Proteins 2016; 84:201–216. © 2015 Wiley Periodicals, Inc.     

Refined structure of spinach glycolate oxidase at 2 A resolution

The amino acid sequence of glycolate oxidase from spinach has been fitted to an electron density map of 2.0 A nominal resolution and the structure has been refined using the restrained parameter least-squares refinement of Hendrickson and Konnert. A final crystallographic R-factor of 18.9% was obtained for 32,888 independent reflections from 5.5 to 2 A resolution. The geometry of the model, consisting of 350 amino acid residues, the cofactor flavin mononucleotide and 298 solvent molecules, is close to ideal with root-mean-square deviations of 0.015 A in bond lengths and 2.6 degrees in bond angles. The expected trimodal distribution with preference for staggered conformation is obtained for the side-chain chi 1-angles. The core of the subunit is built up from the eight beta-strands in the beta/alpha-barrel. This core consists of two hydrophobic layers. One in the center is made up of residues pointing in from the beta-strands towards the barrel axis and the second, consisting of two segments of residues, pointing out from the beta-strands towards the eight alpha-helices of the barrel and pointing from the helices towards the strands. The hydrogen bond pattern for the beta-strands in the beta/alpha-barrel is described. There are a number of residues with 3(10)-helix conformation, in particular there is one left-handed helix. The ordered solvent molecules are organized mainly in clusters. The average isotropic temperature factor is quite high, 27.1 A2, perhaps a reflection of the high solvent content in the crystal. The octameric glycolate oxidase molecule, which has 422 symmetry, makes strong interactions around the 4-fold axis forming a tight tetramer, but only weak interactions between the two tetramers forming the octamer.     

The structure of yeast enolase at 2.25-A resolution. An 8-fold beta + alpha-barrel with a novel beta beta alpha alpha (beta alpha)6 topology

The three-dimensional structure of yeast enolase has been determined by the multiple isomorphous replacement method followed by the solvent flattening technique. A polypeptide model, corresponding with the known amino acid sequence, has been fitted to the electron density map. Crystallographic restrained least-squares refinement of the model without solvent gave R = 20.0% for 6-2.25-A resolution with good geometry. A model with 182 water molecules and 1 sulfate which is still being refined has presently R = 17.0%. The molecule is a dimer with subunits related by 2-fold crystallographic symmetry. The subunit has dimensions 60 X 55 X 45 A and is built from two domains. The smaller N-terminal domain has an alpha + beta structure based on a three-stranded antiparallel meander and four helices. The main domain is an 8-fold beta + alpha-barrel. The enolase barrel is, however, different from the triose phosphate isomerase barrel; its topology is beta beta alpha alpha (beta alpha)6 rather than (beta alpha)8 as found in triose phosphate isomerase. The inner beta-barrel is not entirely parallel, the second strand is antiparallel to the other strands, and the direction of the first helix is also reversed with respect to the other helices. This supports the hypothesis that some enzymes evolved independently producing the stable structure of beta alpha barrels with either enolase or triose phosphate isomerase topology. The active site of enolase is located at the carboxylic end of the barrel. A fragment of the N-terminal domain and two long loops protruding from the barrel domain form a wide crevice leading to the active site region. Asp246, Glu295, and Asp320 are the ligands of the conformational cation. Other residues in the active site region are Glu168, Asp321, Lys345, and Lys396.     

The jigsaw puzzle model: search for conformational specificity in protein interiors

The jigsaw puzzle model postulates that the predominant factor relating primary sequence to three-dimensional fold lies in the stereospecific packing of interdigitating side-chains within densely packed protein interiors. An attempt has been made to check the validity of the model by means of a surface complementarity function. Out of a database of 100 highly resolved protein structures the contacts between buried hydrophobic residues (Leu, Ile, Val, Phe) and their neighbours have been categorized in terms of the extent of side-chain surface area involved in a contact (overlap) and their steric fit (Sm). The results show that the majority of contacts between a buried residue and its immediate neighbours (side-chains) are of high steric fit and in the case of extended overlap at least one of the angular parameters characterizing interresidue geometry to have pronounced deviation from a random distribution, estimated by chi(2). The calculations thus tend to support the "jigsaw puzzle" model in that 75-85% of the contacts involving hydrophobic residues are of high surface complementarity, which, coupled to high overlap, exercise fairly stringent constraints over the possible geometrical orientations between interacting residues. These constraints manifest in simple patterns in the distributions of orientational angles. Approximately 60-80% of the buried side-chain surface packs against neighbouring side-chains, the rest interacting with main-chain atoms. The latter partition of the surface maintains an equally high steric fit (relative to side-chain contacts) emphasizing a non-trivial though secondary role played by main-chain atoms in interior packing. The majority of this class of contacts, though of high complementarity, is of reduced overlap. All residues whether hydrophobic or polar/charged show similar surface complementarity measures upon burial, indicating comparable competence of all amino acids in packing effectively with their atomic environments. The specificity thus appears to be distributed over the entire network of contacts within proteins. The study concludes with a proposal to classify contacts as specific and non-specific (based on overlap and fit), with the former perhaps contributing more to the specificity between sequence and fold than the latter.     

Characterization of the backbone geometry of protein native state structures

We present a simple method for analyzing the geometry of noncovalent residue-residue interactions stabilizing the protein structure, which takes into account the constraints on the local backbone geometry. We find that the principal geometrical constraints are amino acid aspecific and are associated with hydrogen bond formation in helices and sheets. In contrast, amino acid residues in nonhelical and nonextended conformations, which make noncovalent interactions stabilizing the protein tertiary structure, display greater flexibility. We apply the method to an analysis of the packing of helices in helical bundle proteins requiring an efficient packing of amino acid side-chains of the interacting helices.     

Peptides related to the active fragment of "proline rich polypeptide", an immunoregulatory protein of the ovine colostrum. Spectroscopic and computer modeling studies

The preferred solution conformation of the PRP-hexapeptide (Tyr-Val-Pro-Leu-Phe-Pro) and of some of its structural analogues was investigated by NMR-spectroscopy, spectrofluorimetry and computer simulation technic. It was found that the preferred conformation is characterized by cis'-conformation of Pro3 and the gamma-turn on the Leu4-residue: for Val2 and Phe5 a beta-structure seems to be privileged. In such a conformation Val2 and Leu4 residues occupy exactly the same positions in space as residues i and i + 3 in an alpha-helix. It suggests that the PRP-hexapeptide can interact with receptor protein inducing or stabilizing its helical conformation by "knobs into holes" packing.     

Exploring conformational states and helical packings in the P2X receptor transmembrane domain by molecular dynamics simulation

The P2X receptor is a trimeric transmembrane protein that acts as an ATP-gated ion channel. Its transmembrane domain (TMD) contains only six helices and three of them, the M2 helices, line the ion conduction pathway. Here, using molecular dynamics simulation, I identify four conformational states of the TMD that are associated with four types of packing between M2 helices. Packing in the extracellular half of the M2 helix produces closed conformations, while packing in the intracellular half produces both open and closed conformations. State transition is observed and supports a mechanism where iris-like twisting of the M2 helices switches the location of helical packing between the extracellular and the intracellular halves of the helices. In addition, this twisting motion alters the position and orientation of residue side-chains relative to the pore and therefore influences the pore geometry and possibly ion permeation. Helical packing, on the other hand, may restrict the twisting motion and generate discrete conformational states.     

Structural basis for molecular recognition in an affibody:affibody complex

Affibody molecules constitute a class of engineered binding proteins based on the 58-residue three-helix bundle Z domain derived from staphylococcal protein A (SPA). Affibody proteins are selected as binders to target proteins by phage display of combinatorial libraries in which typically 13 side-chains on the surface of helices 1 and 2 in the Z domain have been randomized. The Z(Taq):anti-Z(Taq) affibody-affibody complex, consisting of Z(Taq), originally selected as a binder to Taq DNA polymerase, and anti-Z(Taq), selected as binder to Z(Taq), is formed with a dissociation constant K(d) approximately 100 nM. We have determined high-precision solution structures of free Z(Taq) and anti-Z(Taq), and the Z(Taq):anti-Z(Taq) complex under identical experimental conditions (25 degrees C in 50 mM NaCl with 20 mM potassium phosphate buffer at pH 6.4). The complex is formed with helices 1 and 2 of anti-Z(Taq) in perpendicular contact with helices 1 and 2 of Z(Taq). The interaction surface is large ( approximately 1670 A(2)) and unusually non-polar (70 %) compared to other protein-protein complexes. It involves all varied residues on anti-Z(Taq), most corresponding (Taq DNA polymerase binding) side-chains on Z(Taq), and several additional side-chain and backbone contacts. Other notable features include a substantial rearrangement (induced fit) of aromatic side-chains in Z(Taq) upon binding, a close contact between glycine residues in the two subunits that might involve aliphatic glycine Halpha to backbone carbonyl hydrogen bonds, and four hydrogen bonds made by the two guanidinium N(eta)H(2) groups of an arginine side-chain. Comparisons of the present structure with other data for affibody proteins and the Z domain suggest that intrinsic binding properties of the originating SPA surface might be inherited by the affibody binders. A thermodynamic characterization of Z(Taq) and anti-Z(Taq) is presented in an accompanying paper.     

Theoretical study of the packing of alpha-helices into possible transmembrane bundles. Sequences including alanines, leucines and serines

Energy optimizations are carried out on packages of two and five alpha-helices containing leucines on their faces of contact and made otherwise of alanines. The effect of these bulky side-chains on the optimal arrangements is analysed and compared to the results previously obtained for pure poly(L-alanine) packages; the essential pairing properties are conserved (near antiparallelism, preponderous role of the non-bonded interactions, possibility of existence of parallel pairs); five alpha-helices made of 8 alanines and 6 leucines (three on each interface) can pack in different stable P5L bundles including various holes, according to the tilt and relative sliding of the helices. Substitution of serines to the alanines lying on the inner wall affects very little the interhelix packing. The seryl side chains adapt their conformation at best to their surroundings. The P5L packages can be used to represent individual subunits arranged in 'superbundles' around a central pit in a channel-forming protein.     

Mechanisms of domain closure in proteins

Certain enzymes respond to the binding of substrates and coenzymes by the closure of an active site that lies in a cleft between two domains. We have examined the mechanism of the domain closure in citrate synthase, for which atomic co-ordinates are available for "open" and "closed" forms. We show that the mechanism of domain closure involves small shifts and rotations of packed helices within the two domains and at their interface. Large motions of distant segments of the structure are the cumulative effect of the small relative shifts in intervening pairs of packed segments. These shifts are accommodated not by changes in packing but rather by small conformational changes in side-chains. We call this the helix interface shear mechanism of domain closure. The relative movements of packed helices follow the principles suggested by our recent study of insulin. This mechanism of domain closure is quite different from the hinge mechanisms that allow the rigid body movements of domains in immunoglobulins. The large interface between the domains of citrate synthase precludes a simple hinge mechanism for its conformational change. The helix interface shear mechanism of conformational change occurs in other enzymes that contain extensive domain-domain interfaces.     

The Alacoil: A very tight,antiparallel coiled-coil of helices

The Alacoil is an antiparallel (rather than the usual parallel) coiled-coil of α-helices with Ala or another small residue in every seventh position, allowing a very close spacing of the helices (7.5–8.5 Å between local helix axes), often over four or five helical turns. It occurs in two distinct types that differ by which position of the heptad repeat is occupied by Ala and by whether the closest points on the backbone of the two helices are aligned or are offset by half a turn. The aligned, or ROP, type has Ala in position “d” of the heptad repeat, which occupies the “tip-to-tip” side of the helix contact where the Cα–Cβ bonds point toward each other. The more common offset, or ferritin, type of Alacoil has Ala in position “a” of the heptad repeat (where the Cα-Cβ bonds lie back-to-back, on the “knuckle-touch” side of the helix contact), and the backbones of the two helices are offset vertically by half a turn. In both forms, successive layers of contact have the Ala first on one and then on the other helix. The Alacoil structure has much in common with the coiled-coils of fibrous proteins or leucine zippers: both are α-helical coiled-coils, with a critical amino acid repeated every seven residues (the Leu or the Ala) and a secondary contact position in between. However, Leu zippers are between aligned, parallel helices (often identical, in dimers), whereas Alacoils are between antiparallel helices, usually offset, and much closer together. The Alacoil, then, could be considered as an “Ala anti-zipper.” Leu zippers have a classic “knobs-into-holes” packing of the Leu side chain into a diamond of four residues on the opposite helix; for Alacoils, the helices are so close together that the Ala methyl group must choose one side of the diamond and pack inside a triangle of residues on the other helix. We have used the ferritin-type Alacoil as the basis for the de novo design of a 66-residue, coiled helix hairpin called “Alacoilin.” Its sequence is: cmSP DQWDKE A AQYDAHA QE FEKKS HRNng TPEA DQYRHM A SQY QAMA QK LKAIA NQLKK Gseter (with “a” heptad positions underlined and nonhelical parts in lowercase), which we will produce and test for both stability and uniqueness of structure.     

Designation of sequences involved in the “coiled-coil” interdomainal connections in fibrinogen: Construction of an atomic scale model

Primary structure studies on human fibrinogen (α₂β₂γ₂) have revealed certain unusual features which are compatible with the existence of a three-stranded set of supercoiled α-helices thought to be characteristic of the keratin family of fibrous proteins. In particular, each of the three non-identical chains has two characteristic braces of cysteines separated by 111 residues. The three chains are apparently bound together at these two junctures in unique six-cysteine rings. The amino acid sequences between these unusual cysteine pairs (themselves separated in all six cases by three residues) are helix-permissive over significant portions of their lengths. Moreover, the non-polar residues tend to vary rhythmically. To test the proposition that these sequences do indeed correspond to the “coiled-coils” long ago predicted on the basis of wide-angle X-ray diffraction studies, we constructed a detailed, atomic scale model of a part of this region. To this end, we fashioned three α-helical segments, each 29 residues long and corresponding to the designated sequences of the α, β and γ-chains, respectively. In each case we incorporated a pitch of approximately 200 Å. We were then able to fit the three helices together in the two possible combinations which yield a pseudo-3-fold axis. In either case all polar residues extend away from the parallel three-stranded rope, and almost all the non-polar side-chains are directed toward the interior. We also constructed a separate model showing how the six cysteines at each end of the proposed rod-like segment are best arranged. The co-ordinates from both models were collected and utilized in a computer-graphic Molecular Modeling System which can display features of the models selectively. Various projections were plotted automatically, some of which are reproduced here.