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.     

Synthesis and characterization of a recombinant fragment of human {alpha}-fetoprotein with antigenic selectivity versus albumin

A DNA sequence coding for human -fetoprotein amino acid sequence38–119 was synthesized and cloned in a bacterial expressionvector. The -fetoprotein sequence was selected as the leasthomologous to albumin, since the two proteins have an overallamino acid identity of %. A chimeric protein was obtained whichwas purified by preparative electrophoresis and characterizedin its primary structure by fast atom bombardment mass spectrometry.About 70% of the -fetoprotein sequence was physically mappedand found to correspond to the amino acids encoded in the syntheticgene. The use of this recombinant protein allowed the selectionof monoclonal antibodies recognizing both the recombinant fragmentand native -fetoprotein. These antibodies should allow the developmentof an immunoassay for -fetoprotein with absolute selectivityversus albumin. This might result in more sensitive clinicaldeterminations, avoiding the possibility of cross-reactions.     

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     

Coordinated amino acid changes in homologous protein families

In the tobamovirus coat protein family, amino acid residuesat some spatially close positions are found to be substitutedin a coordinated manner [Altschuh et al. (1987) J. Mol. Biol.,193,693]. Therefore, these positions show an identical patternof amino acid substitutions when amino acid sequences of thesehomologous proteins are aligned. Based on this principle, coordinatedsubstitutions have been searched for in three additional proteinfamilies: serine proteases, cysteine proteases and the haemoglobins.Coordinated changes have been found in all three protein familiesmostly within structurally constrained regions. This methodworks with a varying degree of success depending on the functionof the proteins, the range of sequence similarities and thenumber of sequences considered. By relaxing the criteria forresidue selection, the method was adapted to cover a broaderrange of protein families and to study regions of the proteinshaving weaker structural constraints. The information derivedby these methods provides a general guide for engineering ofa large variety of proteins to analyse structure–functionrelationships.     

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.     

Structure of {alpha}-{alpha}-hairpins with short connections

This paper presents the results of detailed stereochemical analysisof structures and sequences of --hairpins with short connections.It is shown that --hairpins of each given type have very similarpatterns of hydrophobic, hydrophilic and glycine residues intheir amino acid sequences. These results can be used in theprediction of --hairpin conformation as well as in protein designand engineering.     

Cooperative deformation of a de novo designed protein

A de novo protein design has been made to understand the uniquepacking of natural proteins that have a ß/-barrelfold. A carefully designed 207 amino acid sequence was synthesizedusing an Escherichia coli expression system and the structuraland thermodynamic characteristics of the purified protein werestudied. At neutral pH the protein is soluble and monomeric,with large amounts of secondary structure and a hydrophobiccore, although the broad resonance peaks of its NMR spectrumsuggest that the designed protein does not have a unique structurewith tightly packed side chains. In an H–D exchange experiment,no amido protons of the designed protein exchanged slowly withdeuterons. At acidic pH, thermal unfolding was observedwitha remarkable change in the excess heat capacity measured directlyby a differential scanning microcalorimeter. The enthalpy andentropy differences at 110°C, extrapolated from analyzedthermodynamic parameters, are 1/3 of the common values for naturalproteins. These measurements indicate that the folding is significantlycooperative as expected, but that the protein is still looselypacked.     

A Pore-forming protein with a protease-activated trigger

Hemolysin (HL) is a 293 amino acid pore-forming toxin, whichis secreted as a water-soluble monomer by Staphylococcus aureus.By forming a hexameric pore, HL damages the plasma membranesof target cells. Previous studies established that HL proteinswith nicks near the midpoint of a central glycine-rich loopare held together by a domain-domain interaction and are hemolyticallyactive. In contrast, HL proteins comprising two HL truncationmutants that overlap in the central loop have no or greatlyreduced pore-forming activity, even though the two chains againform a tight complex. Based on these findings, overlap mutantshave now been designed that are activated when redundant aminoacids in the loop are removed by proteases. Further, the identityof the activating enzyme can be specified by additional mutagenesisof the protease recognition site in the overlap sequence. Mutantsof aHL that are activated by tumor-associated proteases mightbe useful components of immunotoxins     

Sequence requirements for cytochromes P450IIA1 and P450IIA2 catalytic activity: evidence for both specific and non-specific substrate binding interactions through use of chimeric cDNAs and cDNA expression

Cytochrome P450s IIA1 and IIA2, encoded by the CYP2A1 and CYP2A2genes, display 88% amino acid sequence similarities. The dissimilaritiesof sequence between these two enzymes are primarily localizedwithin four discrete regions of the polypeptides that are separatedby regions of absolute sequence identity. IIA1 specificallyhydroxylates the prototype substrate testosterone at the 7 and6 position with a predominance of 7 metabolite. IIA2, on theother hand, hydroxylates this steroid at eight positions onthe molecule, with one of the most abundant metabolites being15hydroxytestosterone. To determine those amino acids responsiblefor the difference in testosterone hydroxylation specificities,chimeras were constructed between IIA1 and IIA2 cDNAs and expressedin cell culture using vaccinia-virus-mediated cDNA expression.Chimeras, in which the first 355 amino acids correspond to asingle enzyme, maintain the specificity associated with thatenzyme. Of six chimeras which have substitutions between aminoacids 161 and 276, two are inactive and the remaining four givesimilar metabolite profiles, in which both 7 and 15 hydroxylationspecificities have been lost. Two of these four chimeras arediametric apposites, suggesting that modification of eitherthe N-terminal or central regions of the enzymes results inconformational changes that prevent the specific binding interactionsresponsible for the narrow regioselectivity associated withIIA1 and 15-hydroxytestosterone formation associated with IIA2.     

A holistic approach to protein structure alignment

A method of protein structure comparison developed previouslyis extended to incorporate other aspects of protein structurein addition to the inter-atomic vectors on which it was originallybased. Each additional aspect, which included hydrogen bonding,solvent exposure, torsional angles and sequence, was introducedseparately and evaluated for its ability to improve alignmentquality. The components were then combined, suitably weighted,to produce a more holistic comparison method. The method wastested on a group of remotely related ß/ type proteinsthat share a common feature in their overall chain fold. Theresults indicated that while the original inter-atomic vectorcomponent was sufficient to give the correct alignment of mostpairs of topologically equivalent proteins, the inclusion ofhydrogen bonds, torsion angles and a measure of solvent exposureled to improvements in the more difficult comparisons. Considerationof amino acid properties, including hydrophobicity, had no beneficialeffect. The failure of the latter component was not unexpectedconsidering the almost total lack of sequence similarity amongthe proteins considered.     

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.     

Detection of secondary structure elements in proteins by hydrophobic cluster analysis

Hydrophobic cluster analysis (HCA) is a protein sequence comparisonmethod based on -helical representations of the sequences wherethe size, shape and orientation of the clusters of hydrophobicresidues are primarily compared. The effectiveness of HCA hasbeen suggested to originate from its potential ability to focuson the residues forming the hydrophobic core of globular proteins.We have addressed the robustness of the bidimensional representationused for HCA in its ability to detect the regular secondarystructure elements of proteins. Various parameters have beenstudied such as those governing cluster size and limits, thehydrophobic residues constituting the clusters as well as thepotential shift of the cluster positions with respect to theposition of the regular secondary structure elements. The followingresults have been found to support the -helical bidimensionalrepresentation used in HCA: (i) there is a positive correlation(clearly above background noise) between the hydrophobic clustersand the regular secondary structure elements in proteins; (ii)the hydrophobic clusters are centred on the regular secondarystructure elements; (iii) the pitch of the helical representationwhich gives the best correspondence is that of an -helix. Thecorrespondence between hydrophobic clusters and regular secondarystructure elements suggests a way to implement variable gappenalties during the automatic alignment of protein sequences.     

Hyperthermostable mutants of Bacillus licheniformis {alpha}-amylase: multiple amino acid replacements and molecular modelling

We have identified previously two critical positions for thethermostability of the highly thermostable -amylase from Bacilluslicheniformis. We have now introduced all 19 possible aminoacid residues to these two positions, His 133 and Ala209. Themost favourable substitutions were to Ile and Val, respectively,which both increased the half-life of the enzyme at 80°Cby a factor of 3. At both positions a stabilizing effect ofhydrophobic residues was observed, although only in the caseof position 133 could a clear correlation be drawn between thehydrophobicity of the inserted amino acid and the gain in proteinstability. The construction of double mutants showed a cumulativeeffect of the most favourable and/or deleterious substitutions.Computer modelling was used to generate a 3-D structure of thewild-type protein and to model substitutions at position 209,which lies in the conserved (/ß)₈ barrel domain of-amylase; Ala209 would be located at the beginning of the thirdhelix of the barrel, in the bottom of a small cavity facingthe fourth helix. The model suggests that replacement by, forexample, a valine could fill this cavity and therefore increaseintra- and interhelical compactness and hydrophobic interactions.     

Sequence divergence analysis for the prediction of seven-helix membrane protein structures: I. Comparison with bacteriorhodopsin

A method using protein sequence divergence to predict the three-dimensionalstructure of the transmembrane domain of seven-helix membraneproteins is described. The key component in the multistep procedureis the calculation of a hydrophilic and lipophilic variabilityindex for each amino acid in an alignment of a family of homologousproteins. The variability profile, a plot of the calculatedvariability index versus alignment position, can be used topredict a tertiary model of the backbone conformation of thetransmembrane domain. This method was applied to bacteriorhodopsin(BR) and the model obtained was compared with the known structureof this protein. Using an alignment of the amino acid sequencesof BR and closely related (20% identity) proteins, the boundariesof the transmembrane regions, their secondary structures andorientations inside the membrane bilayer were predicted basedon the variability profile. Additional information about theshape of the helix bundle was also obtained from the averagevariability of each transmembrane helix with the assumptionthat the helices are packed sequentially and form a closed helixbundle. Correct features of the known structure of BR were foundin the model structure, suggesting that a similar strategy canbe used to predict transmembrane helices and the packing shapeof other membrane proteins with seven transmembrane helices,such as the opsins and other G-protein coupled receptors.     

Protein stability for single substitution mutants and the extent of local compactness in the denatured state

The stability changes caused by single amino acid substitutionsare studied by a simple, empirical method which takes accountof the free energy change in the compact denatured state aswell as in the native state. The conformational free energyis estimated from effective inter-residue contact energies,as evaluated in our previous study. When this method is applied,with a simple assumption about the compactness of the denaturedstate, for single amino acid replacements at Glu49 of the tryptophansynthase subunit and at Ile3 of bacteriophage T4 lysozyme,the estimates of the unfolding Gibbs free energy changes correlatewell with observed values, especially for hydrophobic aminoacids, and it also yields the same magnitudes of energy as theobserved values for both proteins. When it is also applied foramino acid replacements at various positions to estimate theaverage number of contacts at each position in the denaturedstate from the observed value of unfolding free energy change,those values for replacements with Gly and Ala at the same residueposition in staphylococcal nuclease correlate well with eachother. The estimated numbers of contacts indicate that the proteinis not fully expanded in the denatured state and also that thecompact denatured state may have a substantially native-liketopology, like the molten globule state, in that there is aweak correlation between the estimated average number of contactsat each residue position in the denatured state and the numberof contacts in the native structure. These results provide somefurther evidence that the inter-residue contact energies asapplied here (i) properly reflect actual inter-residue interactionsand (ii) can be considered to be a pairwise hydrophobicity scale.Also, the results indicate that characterization of the denaturedstate is critical to understanding the folding process.     

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.     

Aspartic acid 50 and tyrosine 108 are essential for receptor binding and cytotoxic activity of tumour necrosis factor beta (lymphotoxin)

Single amino acid substitutions were generated in predictedhydrophilic loop regions of the human tumour necrosis factorbeta (TNF-ß) molecule, and the mutant proteins wereexpressed in Escherichia coli and purified. Mutants with singleamino acid changes at either of two distinct loop regions, atpositions aspartic acid 50 or tyrosine 108, were found to havegreatly reduced receptor binding and cytotoxic activity. Thesetwo regions in TNF-ß correspond to known loop regionswhere mutations also result in loss of biological activity ofTNF–, a related cytokine which shares the same cellularreceptors with TNF-ß. The two distinct loops at positions31-34 and 84-89 in the known three-dimensional structure ofTNF- (equivalent to positions 46–50 and 105–110respectively in TNF-ß), lie on opposite sides of theTNF- monomer. When the TNF-a monomer forms a trimer, the twoloops, each from a different subunit of the trimer, come togetherand lie in a cleft between adjacent subunits. Together, thesefindings suggest that a TNF receptor binds to a cleft betweensubunits via surface loops at amino acid residues 31–34and 84–89 in TNF–, and similarly via surface loopsincluding amino acids aspartic acid 50 and tyrosine 108 in TNF–ß.     

Carbonyl-carbonyl interactions stabilize the partially allowed Ramachandran conformations of asparagine and aspartic acid

Asparagine and aspartate are known to adopt conformations inthe left-handed -helical region and other partially allowedregions of the Ramachandran plot more readily than any othernon-glycyl amino acids. The reason for this preference has notbeen established. An examination of the local environments ofasparagine and aspartic acid in protein structures with a resolutionbetter than 1.5 Å revealed that their side-chain carbonylsare frequently within 4 Å of their own backbone carbonylor the backbone carbonyl of the previous residue. Calculationsusing protein structures with a resolution better than 1.8 Åreveal that this close contact occurs in more than 80% of cases.This carbonyl–carbonyl interaction offers an energeticsabilization for the partially allowed conformations of asparagineand aspartic acid with respect to all other non-glycyl aminoacids. The non-covalent attractive interactions between thedipoles of two carbonyls has recently been calculated to havean energy comparable to that of a hydrogen bond. The preponderanceof asparagine in the left-handed -helical region, and in generalof aspartic acid and asparagine in the partially allowed regionsof the Ramachandran plot, may be a consequence of this carbonyl–carbonylstacking interaction.     

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.