Maintenance of the hydrophobic face of the diphtheria toxin amphipathic transmembrane helix 1 is essential for the efficient delivery of the catalytic domain to the cytosol of target cells

Abstract The transmembrane (T) domain of diphtheria toxin (DT) comprisesnine -helices and has been shown to play an essential role inthe efficent delivery of the catalytic (C) domain ofDT acrossthe eukaryotic cell membrane and into the cytosol. We have demonstratedrecently thatthe first three amphipathic helixes of the T domain,although not necessary for either channel formation or receptorbinding, are required for the efficient transmembrane deliveryof the Cdomain.In the present study,we have performed a detailedstructure-function analysis of T domainhelix 1 (TH1) of theDT-related fusion protein DAB₃₈₉lL-2. We performed exchangeandsite-directed mutagenesis of TH1 and the resulting mutantfusion toxins were analyzed by gel electrophoresis and testedfor their efficiencies in the delivery of the C domain to thecell cytosol. We demonstrate that the overall charge distributionand hydrophobicity of amino acids in the amphipathic helix TH1,rather than a specific amino acid sequence, are critical forthe function of this helix. The insertion of a charged residuein the hydrophobic face of TH1 abolishes cytotoxic activity,whereas replacement of a hydrophobic residue by a charged aminoacid in the hydrophilic face of the helix has little, if any,effect on cytotoxic activity. In addition,we have identifiedSer220 by site-directed mutagenesis as a residue that appearsto be criticalfor correct folding of the fusion toxin. Mutationsin this position result in fusion proteins that are extremelysensitive to proteolytic attack.     

A simple method for predicting transmembrane {alpha} helices with better accuracy

The prediction of a protein's structure from its amino acidsequence has been a long-standing goal of molecular biology.In this work, a new set of conformational parameters for membranespanning helices was developed using the information from thetopology of 70 membrane proteins. Based on these conformationalparameters, a simple algorithm has been formulated to predictthe transmembrane helices in membrane proteins. A FORTRAN programhas been developed which takes the amino acid sequence as inputand gives the predicted transmembrane -helices as output. Thepresent method correctly identifies 295 transmembrane helicalsegments in 70 membrane proteins with only two overpredictions.Furthermore, this method predicts all 45 transmembrane helicesin the photosynthetic reaction center, bacteriorhodopsin andcytochrome c oxidase to an 86% level of accuracy and so is betterthan all other methods published to date.     

A 3D model of the {delta} opioid receptor and ligand-receptor complexes

A model for the 3D structure of the transmembrane domain ofthe opioid receptor was predicted from the sequence divergenceanalysis of 42 sequences of G-protein coupled peptide hormonereceptors belonging to the opioid, somatostatin and angiotensinreceptor families. No template was used in the prediction steps,which include multiple sequence alignment, calculation of avariability profile of the aligned sequences, use of the variabilityprofile to identify the boundaries of transmembrane regions,prediction of their secondary structure, optimization of thepacking shape in a helix bundle, prediction of side chain conformationsand structural refinement The general shape of the model issimilar to that of the low resolution rhodopsin structure inthat the TM3 and TM7 helices are most buried in the bundle andthe TM1 and TM4 helices are most exposed to the lipid phase.An initial assessment of this model was made by determiningto what extent a binding site identified using four structurallydisparate high affinity opioid ligands was consistent withknown mutational studies. With the assumption that the pro-tonatedamine nitrogen, a feature common to all opioid ligands, interactswith the highly conserved Aspl27 in TM3, a pocket was foundthat satisfied the criteria of complementarity to the requirementsfor receptor recognition for these four diverse ligands, two selective antagonists (the fused ring naltrindole and the peptideTyr-Tic-Phe-Phe-NH₂) and the two agonists lofentanil and BW373U86deduced from previous studies of the ligands alone. These ligandscould be accommodated in a similar region of the receptor. Thereceptor binding site identified in the optimized complexescontained many residues in positions known to affect ligandbinding in G-protein coupled receptors. These results also allowedidentification of key residues as candidates for point mutationsfor further assessment and refinement of this model as wellas preliminary indications of the requirements for recognitionof this receptor.     

Sequence divergence analysis for the prediction of seven-helix membrane protein structures: II. A 3-D model of human rhodopsin

A three-dimensional (3-D) model of the transmembrane domainof human rhodopsin was predicted from the sequence divergenceanalysis of 42 sequences of rhodopsins and visual pigments withouta template. The prediction steps include multiple sequence alignment,calculation of a variability profile of the aligned sequences,use of the variability profile to identify the boundaries oftransmembrane regions, their secondary structure and packingshape in a helix bundle, prediction of side-chain conformationsand structure refinement. The identification of the retinalbinding site was assisted by its known covalent linkage withK296. The structural features of the predicted 3-D model arein good agreement with a low resolution electron density mapof bovine rhodopsin and with residues in contact with retinalas determined experimentally.     

Physicochemical factors for discriminating between soluble and membrane proteins: hydrophobicity of helical segments and protein length

The average hydrophobicity of a polypeptide segment is consideredto be the most important factor in the formation of transmembranehelices, and the partitioning of the most hydrophobic (MH) segmentinto the alternative nonpolar environment, a membrane or hydrophobiccore of a globular protein may determine the type of proteinproduced. In order to elucidate the importance of the MH segmentin determining which of the two types of protein results froma given amino acid sequence, we statistically studied the characteristicsof MH helices, longer than 19 residues in length, in 97 membraneproteins whose three-dimensional structure or topology is known,as well as 397 soluble proteins selected from the Protein DataBank. The average hydrophobicity of MH helices in membrane proteinshad a characteristic relationship with the length of the protein.All MH helices in membrane proteins that were longer than 500residues had a hydrophobicity greater than 1.75 (Kyte and Doolittlescale), while the MH helices in membrane proteins smaller than100 residues could be as hydrophilic as 0.1. The possibilityof developing a method to discriminate membrane proteins fromsoluble ones, based on the effect of size on the type of proteinproduced, is discussed.     

Recognition of transmembrane {alpha}-helical segments with environmental profiles

A method for assessing the environmental properties of membrane-spanning-helical peptides in proteins has beenproposed. The algorithmemploys a set of environmental preference parameters derivedfor amino acid residues based on the analysis of the 3-D structuresof membrane domains in bacteriorhodopsin and photoreaction centersRhodopseudomonasviridis and Rhodobacter sphaeroides. The resulting3-D–1-D scores for transmembrane segments are significantlydifferent from those derived for -helices in globular proteins.The parameters obtained havebeen used to construct environmentalprofiles for membrane -helices in bacteriorhodopsin and photoreactioncenters. The profiles successfully recognize their own sequencesin several specially designed large databases. The method hasbeenapplied to several membrane proteins with unknown spatial structures.Most of their membrane-spanning peptides were efficiently recognizedby the profiles. The predicted environment of the residues inthe membrane segments fits the experimental data well. The approachis independent of any homology data and can be employed to delineatethe membrane segments of a protein with environmental characteristicsclose to those of bacteriorhodopsin and photoreaction centers.The alignment of these segments with the reference profilesprovides a considerable amount of data about their lipid andprotein exposure.     

Computational analysis of alpha-helical membrane protein structure: implications for the prediction of 3D structural models

Relatively little has been known about the structure of alpha-helical membrane proteins, since until recently few structures had been crystallized. These limited data have restricted structural analyses to the prediction of secondary structure, rather than tertiary folds. In order to address this, this paper describes an analysis of the 23 available membrane protein structures. A number of findings are made that are of particular relevance to transmembrane helix packing: (1) on average lipid-tail-accessible transmembrane residues are significantly more hydrophobic, less conserved and contain different residue types to buried residues; (2) charged residues are not always buried and, when accessible to membrane lipid tails, few are paired with another charge and instead they often interact with phospholipid head-groups or with other residue types; (3) a significant proportion of lipid-tail-accessible charged and polar residues form hydrogen bonds only with residues one turn away in the same helix (intra-helix); (4) pore-lining residues are usually hydrophobic and it is difficult to distinguish them from buried residues in terms of either residue type or conservation; and (5) information was gained about the proportion of helices that tend to contribute to lining a pore and the resulting pore diameter. These findings are discussed with relevance to the prediction of membrane protein 3D structure.     

Design and characterization of a soluble fragment of the extracellular ligand-binding domain of the peptide hormone receptor guanylyl cyclase-C

The intestinal guanylyl cyclase-C (GC-C) was originally identifiedas an Escherichia coli heat-stable enterotoxin (STa) receptor.STa stimulates GC-C to much higher activity than the endogenousligands guanylin and uroguanylin, causing severe diarrhea. Toinvestigate the interactions of the endogenous and bacterialligands with GC-C, we designed and characterized a soluble andproperly folded fragment of the extracellular ligand-bindingdomain of GC-C. The membrane-bound guanylyl cyclases exhibita single transmembrane spanning helix and a globularly foldedextracellular ligand-binding domain that comprises about 410of 1050 residues. Based on the crystal structure of the dimerized-bindingdomain of the guanylyl cyclase-coupled atrial natriuretic peptidereceptor and a secondary structure-guided sequence alignment,we generated a model of the extracellular domain of GC-C comprisedof two subdomains. Mapping of mutational and cross-link dataonto this structural model restricts the ligand-binding regionto the membrane proximal subdomain. We thus designed miniGC-C,a 197 amino acid fragment that mimics the ligand-binding membraneproximal subdomain. Cloning, expression and spectroscopic studiesreveal miniGC-C to be a soluble and properly folded proteinwith a distinct secondary and tertiary structure. MiniGC-C bindsSTa with nanomolar affinity.     

The prediction and orientation of {alpha}-helices from sequence alignments: the combined use of environment-dependent substitution tables, Fourier transform methods and helix capping rules

Amino acid substitution tables are used to estimate the extentto which amino acids in families of homologous proteins areexposed to the solvent. The approach depends on the comparisonof difference environment-dependent tables for solvent accessible/inaccessibleresidues with amino acid substitutions at each position in analigned set of sequences. The periodicity in the predicted accessible/inaccessibleresidues is calculated using a Fourier transform procedure modifiedfrom that used to calculate hydrophobic moments. a-Helices areidentified from the characteristic periodicities and the solventaccessible face of the helix is defined. The initial helix predictionsare refined using rules for identifying the N- and C-terminiof helices from sequence alignments. These rules have been definedfrom a study of protein structures. The combined method correctlypredicts 79% of the residues in helices and incorrectly predictsonly 12% of the nonhelical residues as helical. In addition,since the method is reliable at predicting the correct numberof helices in the correct position in the sequence and sinceit also predicts the internal face of each helix, the resultscan be used to postulate 3-D arrangements of the secondary structureelements.     

Composition analysis of {alpha}-helices in thermophilic organisms

We present a statistical comparison of the amino acid compositionin a secondary structure element, the -helix, of proteins stableat high temperatures with those which are less so. This studyhas shown that the temperature-dependent Zimm-Bragg helix propagationvalue s is not a good predictor for the helix-forming tendencyof an amino acid in thermostable proteins. However, we haveshown that s, the change in s from 20 to 60°C, accuratelypredicts the direction of the probability shift for 15 aminoacids in thermostable protein a-helices, although it does notpredict the magnitude of that change. The residues tyrosine,glycine and glutamine show a significant increase in residencyin a-helices for thermostable proteins over their nonthermostablecounterparts. Significant decreases in -helix residency occurfor the residues valine, glutamic acid, histidine, cysteineand aspartic acid in proteins from thermophilic organisms. Aromaticinteractions, hydrogen bonding and a reduction of charge mayexplain the increase observed for tyrosine and glutamine andthe decrease in glutamic acid and aspartic acid, although packingconsiderations cannot be ruled out The only physical explanationfor the increase in glycine would seem to be its positive svalue     

Secondary structure prediction of the H5 pore of potassium channels

The ‘H5’ segment located between the putative fifthand sixth transmembrane helices is the most highly conservedregion in voltage-gated potassium channels and it is believedto constitute a major part of the ion conduction path (pore).Here we present a two-step procedure, comprising secondary structureprediction and hydrophobic moment profiling, to predict thestructure of this important region. Combined results from theapplication of the procedure to the H5 region of four voltage-gatedand five other K⁺ channel sequences lead to the prediction ofa ß-strand-turn-(3-strand structure for H5. The reasonsfor the application of these soluble protein methods to partsof membrane proteins are: (i) that pore-lining residues areaccessible to water and (ii) that a large enough database ofhighresolution membrane protein structures does not yet existThe results are compared with experimental results, in particularspectroscopic studies of two peptides based on the H5 sequenceof SHAKER potassium channel. The procedure developed here maybe applicable to wateraccessible regions of other membrane proteins.     

Inverse protein folding by the residue pair preference profile method: estimating the correctness of alignments of structurally compatible sequences

The residue pair preference profile (R3P) method is an inversefolding method that combines environmental profiles and pairpreference profiles. The method uses statistical preferencesfor residue pairs which score the likelihood of finding a profiledresidue to be paired with a residue within its local environmentAll pairs are characterized by their dihedral angles, secondarystructure and number of neighboring residues as a function ofresidue type. Each residue pair preference is expressed forall 20 amino acids of the profiled residue and is weighted bythe compatibility of the environment residue with its own localenvironment The R3P method produces an initial profile-sequencealignment which is then refined by converting the initial profileinto a profile of a target sequence threaded into the structureof the initial profile. We have tested this method by evaluatingalignments of sequences with known 3-D structures using structuralsuperposition alignments as reference. R3P-sequence alignmentsare 50% correct on average for sequences whose 3-D structurepairs superimpose with an r.m.s. deviation of 1.97 Å.The average improvement in correctness during this iterativerefinement is 14%. The R3P-sequence alignments are comparedwith sequence-sequence and 3-D profile-sequence alignments.When all three methods are combined, on average 50% of the alignmentsare correct for pairs of 3-D structures that superimpose within2.12 Å. A 3-D model of HisA is predicted with the combinedmethod.     

Functional diversity of the phosphoglucomutase superfamily: structural implications

Three-dimensional structural models of three members of thephosphoglucomutase (PGM) superfamily, parafusin, phosphoglucomutase-relatedprotein and sarcoplasmic reticulum phosphoglucomutase, wereconstructed by homology modeling based on the known crystalstructure of rabbit muscle phosphoglucomutase. Parafusin, phosphoglucomutase-relatedprotein and sarcoplasmic reticulum phosphoglucomutase each have50% or more identity with rabbit muscle phosphoglucomutase atthe amino acid level and all are reported to exhibit no or minorphosphoglucomutase activity. There are four major insertionsand two deletions in the parafusin sequence relative to PGM,all of which are located in surface-exposed loops connectingsecondary structural elements. The remaining amino acid substitutionsare distributed throughout the sequence and are not predictedto alter the polypeptide fold. Parafusin contains a putativeprotein kinase C site located on a surface loop in domain IIthat is not present in the homologs. Although the general domainstructure and the active site of rabbit muscle phosphoglucomutaseare preserved in the model of phosphoglucomutase-related protein,a major structural difference is likely to occur in domain 1due to the absence of 55 amino acid residues in PGM-RP. Thisdeletion predicts the loss of three -helices and one ß-strandfrom an anti-parallel ß-sheet in this domain as comparedwith the rabbit muscle phosphoglucomutase.     

Analysis of heregulin symmetry by weighted evolutionary tracing

Heregulins are members of the protein family of EGF-like growthand differentiation factors. The primary cell-surface targetsof heregulins are heterodimers of the EGF-receptor homolog HER2with either HER3 or HER4. We used a weighted evolutionary traceanalysis to identify structural features that distinguish theEGF-like domain (hrg) of heregulins from other members of theEGF family. In this analysis, each amino acid sequence is weightedaccording to its uniqueness and the variability in each positionis assigned by an amino acid substitution matrix. Conservedresidues in heregulin that are variable in other EGF-like domainsare considered possible specificity-conferring residues. Thisanalysis identifies two clusters of residues at the foot ofthe boot-shaped hrg domain. The residues in one cluster arerecruited from the N-terminus; those in the other are from the-loop region and show a weak sequence similarity to the N-terminalresidues at the opposite side of the boot. The remaining residueswith high conservation scores distribute themselves into thesetwo distinct surfaces on hrg. This pseudo-twofold symmetry andthe presence of two distinct interfaces may reflect the preferenceof hrg for heterodimeric versus homodimeric HER complexes.     

The Influence of Fatty Acids on the GpA Dimer Interface by Coarse-Grained Molecular Dynamics Simulation

The hydrophobic thickness of membranes, which is manly defined by fatty acids, influences the packing of transmembrane domains of proteins and thus can modulate the activity of these proteins. We analyzed the dynamics of the dimerization of Glycophorin A (GpA) by molecular dynamics simulations to describe the fatty acid dependence of the transmembrane region assembly. GpA represents a well-established model for dimerization of single transmembrane helices containing a GxxxG motif in vitro and in silico. We performed simulations of the dynamics of the NMR-derived dimer as well as self-assembly simulations of monomers in membranes composed of different fatty acid chains and monitored the formed interfaces and their transitions. The observed dimeric interfaces, which also include the one known from NMR, are highly dynamic and converted into each other. The frequency of interface formation and the preferred transitions between interfaces similar to the interface observed by NMR analysis strongly depend on the fatty acid used to build the membrane. Molecular dynamic simulations after adaptation of the helix topology parameters to better represent NMR derived structures of single transmembrane helices yielded an enhanced occurrence of the interface determined by NMR in molecular dynamics simulations. Taken together we give insights into the influence of fatty acids and helix conformation on the dynamics of the transmembrane domain of GpA.     

Effects of amino acid substitutions in the hydrophobic core of {alpha}-lactalbumin on the stability of the molten globule state

Five mutant –lactalbumins, with one or two amino acidsubstitution(s) in the B helix, were engineered to examine therelation between the stability of the molten globule state andthe hydrophobicity of these amino acids. The mutation sites(Thr29, Ala30 and Thr33) have been chosen on the basis of comparisonof the amino acid sequences of goat, bovine and gunea pig –lactalbumin,in which the guinea pig protein shows a remarkably more stablemolten globule than the other proteins. The recombinant proteinswere expressed Escherichia coli and then purified and refoldedefficiently to produce the active proteins. The stability ofthe molten globule state of these engineered proteins has beeninvestigated by urea–induced unfolding transition underan acidic condition (pH 2.0), where the molten globule stateis stable in the absence of urea. The results show that themolten globule state is stabilized by the amino acid substitutionswhich raise the hydrophobicity of the residues, suggesting thatthe hydrophobic core in a globular protein plays an importantrole in the stability of the molten globule state. The changein stabilization free energy of the molten globule state causedby each amino acid substitution has been evaluated, and molecularmechanisms of stabilization of the molten globule state arediscussed.     

A multivariate analysis method for discriminating protein secondary structural segments

Using discriminant analysis, three types of protein secondarystructure segments—helices, ß-strands and coils—arediscriminated by amino acid sequence information alone. A variablein the discriminant analysis is defined by the amino acid indexused to represent the sequence data and by the calculation methodused to extract a feature in this representation. Thus, thethree types of secondary structure segments derived from a setof non-homologous proteins from the Protein Data Bank are analyzedby 888 variables, which correspond to the mean, standard deviation,3.6-residue periodicity and 2-residue periodicity for the numericalprofiles determined from 222 published amino acid indices. Thesevariables are combined to obtain best discrimination of thethree types of segments. When up to three variables are combined,the best discrimination rate was 75%. The variables selectedconsist of the mean of propensity (or turn propensity), themean of ß propensity, and the 3.6-residue periodicityof hydrophobicity. This variable selection procedure can alsobe applied to other types of discrimination problem, once groupsof sequence data are properly organized.     

A combinatorial library of an {alpha}-helical bacterial receptor domain

The construction and characterization of a combinatorial libraryof a solvent-exposed surface of an -helical domain derived froma bacterial receptor is described. Using a novel solid-phaseapproach, the library was assembled in a directed and successivemanner utilizing single-stranded oligonucleotides containingmultiple random substitutions for the variegated segments ofthe gene fragment The simultaneous substitution of 13 residuesto all 20 possible amino acids was carried out in a region spanning81 nucleotides. The randomization was made in codons for aminoacids that were modelled to be solvent accessible at a surfacemade up from two of the three a-helices of a monovalent Fc-bindingdomain of staphylococcal protein A. After cloning of the PCR-amplifiedlibrary into a phagemid vector adapted for phage display ofthe mutants, DNA sequencing analysis suggested a random distributionof codons in the mutagenized positions. Four members of thelibrary with multiple substitutions were produced in Escherichiacoli as fusions to an albumin-binding affinity tag derived fromstreptococcal protein G. The fusion proteins were purified byhuman serum albumin affinity chromatography and subsequentlycharacterized by SDSelectrophoresis, CD spectroscopy and biosensoranalysis. The analyses showed that the mutant protein A derivativescould all be secreted as soluble full-length proteins. Furthermore,the CD analysis showed that all mutants, except one with a prolineintroduced into helix 2, have secondary structures in closeagreement with the wild-type domain. These results proved thatmembers of this -helical receptor library with multiple substitutionsin the solvent-exposed surface remain stable and soluble inE.coli. The possibility of using this library for a phenotypicselection strategy to obtain artificial antibodies with novelfunctions is discussed.     

Proline residues in transmembrane helices of channel and transport proteins: a molecular modelling study

Proline residues are commonly found in putative transbilayerhelices of many integral membrane proteins which act as transporters,channels and receptors. Intramembranous prolines are often conservedbetween homologous proteins. It has been suggested that suchintrahelical prolines provide liganding sites for cations viaexposure of the backbone carbonyl oxygen atoms of residues i-3and i-4 (relative to the proline). Molecular modelling studieshave been carried out to evaluate this proposal. Bundles ofparallel proline-kinked helices are considered as simplifiedmodels of ion channels. The energetics of K⁺ ion-helix bundleinteractions are explored. It is shown that carbonyl oxygensexposed by the proline-induced kink and at the C-terminus ofthe helices may provide cation-liganding sites. ‘Hybrid’bundles of antiparallel helices, only some of which containproline residues, are considered as models of transport proteins.Again, praline-exposed carbonyi oxygens are shown to be capableof liganding cations. The roles of -helix dipoles and of thegeometry of helix packing are considered in relation to cation-bundleinteractions. Implications with respect to modelling of ionchannel and transport proteins are discussed