The {alpha}/{beta} hydrolase fold

We have identified a new protein fold—the /ßhydrolase fold—that is common to several hydrolytic enzymesof widely differing phylogenetic origin and catalytic function.The core of each enzyme is similar: an /ß sheet, notbarrel, of eight ß-sheets connected by -helices. Theseenzymes have diverged from a common ancestor so as to preservethe arrangement of the catalytic residues, not the binding site.They all have a catalytic triad, the elements of which are borneon loops which are the best-conserved structural features inthe fold. Only the histidine in the nucleophile-histidine-acidcatalytic triad is completely conserved, with the nucleophileand acid loops accommodating more than one type of amino acid.The unique topological and sequence arrangement of the triadresidues produces a catalytic triad which is, in a sense, amirror-image of the serine protease catalytic triad. There arenow four groups of enzymes which contain catalytic triads andwhich are related by convergent evolution towards a stable,useful active site: the eukaryotic serine proteases, the cysteineproteases, subtilisins and the /ß hydrolase fold enzymes.     

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.     

Structural models of the evolutionarily conservative central domain of silk-moth chorion proteins

Silk-moth chorion proteins belong to a small number of families:A, B, C, Hc-A and Hc-B. The central domain is an evolutionarilyconservative region in each family, of variable length and compositionbetween families. This domain shows dear 6-fold periodicitiesfor various amino acid residues, e.g. glycine. The periodicities,together with the well-documented prevalence of ß-sheetand ß-turn secondary structure of chorion proteins,strongly support a structural model in which four-residue ß-strandsalternate with ß-turns, forming a compact antiparallel,probably twisted ß-sheet. Conformational analysis,aided by interactive graphics refinement and recent experimentalfindings, further suggest that this structure consists of ß-strands,alternating with I' and II' ß-turns, and apparentlyforms the basis for the molecular and supramolecular assemblyof chorion.     

Characterization of the DNA binding domain of the mouse IRF-2 protein

The DNA binding domain of the interferon regulatory factor-2protein (IRF-2) has been produced and characterized, -chymotrypsindigestion of the purified IRF-2 protein bound to a syntheticbinding site yields a peptide fragment of 14 K in molecularweight. N-terminal analysis of this peptide fragment showedthat its sequence is the same as that of the intact IRF-2. Apeptide fragment of {small tilde} 14 K, IRF-2(113), which correspondsto the N-terminal 113 amino acids of the intact IRF-2 protein,has been expressed in a functional form in Escherichia coli.The first methionine was processed during the expression andthe purified IRF-2(113) thus contains 112 amino acids. DNaseI footprinting and gel retardation assaying showed that IRF-2(113)binds to a synthetic DNA having the consensus binding site andto the upstream regulatory sequence of the IFN-ß geneas intact IRF-2 does. These results showed that this peptidefragment, IRF-2(113), may be a good material for investigationof the DNA binding domain of IRF-2 and of the DNA–proteininteraction.     

Native-like {beta}-hairpin structure in an isolated fragment from ferredoxin: NMR and CD studies of solvent effects on the N-terminal 20 residues

The conformational properties of protein fragments have beenwidely studied as models of the earliest initiation events inprotein folding. While native-like -helices and ß-turnshave been identified, less is known about the factors that underlyß-sheet formation, in particular ß-hairpins,where considerably greater long-range order is required. TheN-terminal 20 residue sequence of native ferredoxin I (fromthe blue-green alga Aphanothece sacrum ) forms a ß-hairpinin the native structure and has been studied in isolation byNMR and CD spectroscopy. Local native-like interactions aloneare unable to stabilize significantly a folded conformationof the 20-residue fragment in purely aqueous solution. However,we show that the addition of low levels of organic co-solventspromotes formation of native-like ß-hairpin structure.The results suggest an intrinsic propensity of the peptide toform a native-like ß-hairpin structure, and that theorganic co-solvent acts in lieu of the stabilizing influenceof tertiary interactions (probably hydrophobic contacts) whichoccur in the folding of the complete ferredoxin sequence. Thestructure of the isolated hairpin, including the native-likeregister of interstrand hydrogen bonding interactions, appearsto be determined entirely by the amino acid sequence. The solventconditions employed have enabled this intrinsic property tobe established.     

Model structures and action of interleukin 1 and its antagonist

A comparison has been made between the homology and hydrophobkityprofiles of six interleukin amino add sequences and that ofthe human interleukin 1ß (IL-lß) for whicha crystal structure exists. The resulting sequence alignmentwas used to build model structures for the sequences for threeIL-l, two IL-1ß and an interleukin receptor antagonist.Analysis of these structures demonstrates that the interleukinmolecule has a strong electric dipole which is generated bythe topological position of the amino acids in the sequence.Electrostatic surface calculations implicate a particular residues(Lysl45) as being fundamental to interleukin activity and thissupports site-directed mutation evidence that this residue isrequired for activity.     

Permuteins of interleukin 1{beta}--a simplified approach for the construction of permutated proteins having new termini

A technique for the rapid and simple generation of permutatedversions of the interleukin-1ß (IL-1ß) geneis described. In this method, the human IL-1ß cDNAis twice amplified by the polymerase chain reaction (PCR) andthe resulting DNA fragments are ligated in tandem. Between thetwo genes, the DNA sequence encodes a short four amino acidloop to link the native N- and C-terminal ends of the IL-1ßprotein. By using PCR amplification from this starting template,a new version of the IL-1ß cDNA was obtained thatencodes a permutated form of the IL-1ß protein wherethe new N- and C-terminal amino acids correspond to residues65 and 64 of the native IL-1ß sequence, respectively.The name ‘permutein’ is proposed to describe proteinsgenerated by this technology. The molecular profile (IL-1 receptorbinding, biologic activity and solution properties) of the IL-1permutein produced by this technology, permutein 65/64, is shownto be identical to that of native IL-1ß The approachshould be useful to define further the structural features ofthis protein that are important for its function.     

6-Phospho-{beta}-galactosidases of Gram-positive and 6-phospho-{beta}-glucosidase B of Gram-negative bacteria: comparison of structure and function by kinetic and immunological methods and mutagenesis of the lacG gene of staphylococcus aureus

The 6-phospho-ß-galactosidase of Staphylococcus aureus,Lactococcus lactis and Lactobacillus casei and 6-phospho-ßglucosidaseB of Escherichia coli build a subfamily inside a greater enzymefamily, named the glycosal hydrolase family 1, which, hi addition,contains nine ß-glycosidases of different origins.Kinetic and immunological evidence is provided in this reportwhich strengthens the relationship of the four 6-phospho-ß-glycosidases.It is shown that the 6-phospho-ß-galactosidases and6-phospho-ß-glucosidase B are able to split aromaticß-galactoside phosphates and ß-glucosidephosphates. The turnover numbers of hydrolysis of substrateswith different epimerization at C-4 of the glycon vary up to15-fold only. Two polydonal antisera, one derived against thenative 6-phospho-ß-galactosidase from S.aureus andthe other derived against the 6-phospho-ß-glucosidaseB, cross-reacted with both enzymes. Peptides of the proteinswere separated by reverse phase HPLC. The cross-reacting peptideswere sequenced and shown to be localized at almost the sameposition in the aligned primary structures of both enzymes.An insertion of nine amino adds near these antigenic domainsis unique for the 6-phospho-ß-glycosidases and missingwithin the sequences of the ß-glycoside-specific membersof the family. The lacG gene of a 6-phospho-ß-galactosidasenegative S.aureus mutant was doned into E.coli and sequenced.In the totally inactive mutant protein only the glycine at position332 was changed to an arginine. This amino acid is part of thesequence insertion near the antigenic domain reacting with bothantisera. These data support the assumption that the regionis of great importance for the function of the enzymes and thatit is possible it determines the specificity of the phosphorylatedform of the substrates. In addition, the 6-phospho-ß-galactosidaseof S.aureus was modified by sitedirected mutagenesis of thecorresponding lacG gene hi order to replace residues Glul60and Glu375, which were suspected of being involved hi the generalacid catalysis of substrate hydrolysis, with glutamine residues.The mutant protein 160EQ retained some catalytic activity whilethe protein 375EQ was totally inactive. Glu375 is the activesite nudeophile of the 6-phospho-ß-galactosidase ofS.aureus. It is located in the sequence motif ENG where Glu358was identified as the catalytkally active nudeophile hi theß-glucosidase of Agrobacterium.     

Structure-function relationships in the catalytic and starch binding domains of glucoamylase

Sixteen primary sequences from five sub-families of fungal,yeast and bacterial glucoamylases were related to structuralinformation from the model of the catalytic domain of Aspergillusawamori var. X100 glucoamylase obtained by protein crystallography.This domain is composed of thirteen -belices, with five conservedregions defining the active site. Interactions between methyl-maltoside and active site residues were modelled, and the importanceof these residues on the catalytic action of different glucoamylaseswas shown by their presence in each primary sequence. The overallstructure of the starch binding domain of some fungal glucoamylaseswas determined based on homology to the Cterminal domains ofBacillus cyclodextrin glucosyltransferases. Crystallographyindicated that this domain contains 6–8 ß-strandsand homology allowed the attribution of a disulfide bridge inthe glucoamylase starch binding domain. Glucoamylase residuesThr525, Asn530 and Trp560, homologous to Bacillus stearothermophiluscyclodextrin glucosyltransferase residues binding to maltosein the Cterminal domain, could be involved in raw-starch binding.The structure and length of the linker region between the catalyticand starch binding domains in fungal glucoamylases can varysubstantially, a further indication of the functional independenceof the two domains.     

Context dependence of phenotype prediction and diversity in combinatorial mutagenesis

Two different combinatorial mutagenesis experiments on the light-harvestingII (LH2) protein of Rhodobacter capsulatus indicate that heuristicrules relating sequence directly to phenotype are dependenton which sets or groups of residues are mutated simultaneously.Previously reported combinatorial mutagenesis of this chromogenicprotein (based on both phylogenetic and structural models) showedthat substituting amino acids with large molar volumes at Gly^ß31caused the mutated protein to have a spectrum characteristicof light-harvesting I (LH1). The six residues that underwentcombinatorial mutagenesis were modeled to lie on one side ofa transmembrane -helix that binds bacteriochlorophyll. In asecond experiment described here, we have not used structuralmodels or phylogeny in choosing mutagenesis sites. Instead,a set of six contiguous residues was selected for combinatorialmutagenesis. In this latter experiment, the residue substitutedat Gly^ß31 was not a determining factor in whetherLH2 or LH1 spectra were obtained; therefore, we conclude thatthe heuristic rules for phenotype prediction are context dependent.While phenotype prediction is context dependent, the abilityto identify elements of primary structure causing phenotypediversity appears not to be. This strengthens the argument forperforming combinatorial mutagenesis with an arbitrary groupingof residues if structural models are unavailable.     

Binding of the globular domain of linker histones H5/H1 to the nucleosome: a hypothesis

The amino acid sequence of the central globular domain of histoneH1/H5 family members is highly homologous. Twenty-four suchsequences have been compared to establish the conserved andvariable residues. Fitting this to the tertiary structure ofthe H5 globular domain shows which of the conserved and variableresidues are peripheral and which internal. Particular attentionis paid to conserved basic residues on the surface, which wetake to be DNA binding. Variable regions and conserved acidicresidues are assumed not to be sites of contact with DNA. Weconclude that one face of the domain, containing a cluster ofbasic residues, is the principal DNA binding site whilst twoopposing faces, orthogonal to the principal site and also containingconserved basic residues, are subsidiary DNA binding sites.Since the DNA binding surface of the domain covers a full 180°arc, we propose that it contacts a ‘cage’ of threeDNA strands on the 2-fold axis of the chromatosome.     

Cloning of rabbit interleukin-1{beta}: differential evolution of IL-1{alpha} and IL-1{beta} proteins

We have cloned the rabbit IL-1ß cDNA, which encodesa 268 amino acid precursor similar in length to other sequencedIL-1 precursors. Comparison of all published IL-1 and IL-1ßsequences respectively indicates that the IL-1 gene family isevolving faster than the IL-1ß family, and that thetwo genes diverged –270 million years ago. Surprisingly,there are differences in the regions preferentially conservedwithin the two families. The IL-1 family is most conserved atthe amino terminus whereas the IL-1ß family is mostconserved in the carboxy-terminal half. This is despite thefact that the carboxy-terminal half encodes the active portionof both molecules and would be expected to adopt a similar ß-sheetstructure in IL-1 as in the published X-ray structure of matureIL-1ß. These findings suggest that differences inthe function and properties of the IL-1 and IL-1ßprecursor molecules may have been conserved. These differencesmay therefore provide an explanation for the existence of twoIL-1 molecules.     

Structural analysis of the CD2 T lymphocyte antigen by site-directed mutagenesis to introduce a disulphide bond into domain 1

Many proteins have been predicted to contain domains with immunoglobulin-Iikefolds and hence to be members of the immunoglobulin superfamily(IgSF). However, several members lack the Cys residues capableof forming the disulphide bond that forms a characteristic bridgebetween the ß sheets in the Ig fold, e.g. domain 1of the lymphocyte antigen CD2. The assignment of ßstrands in CD2 by sequence analysis was tested by attemptingto introduce a disulphide bridge between the ß sheetsby mutating two residues in the relevant positions to Cys. Mutant,soluble forms of CD2 were expressed in Chinese hamster ovarycell lines and amino add sequencing showed that a disulphidebond was formed as predicted, but not in the control where oneCys residue was misplaced by four residues. Evidence that bothmutated molecules folded correctly is given by the indistinguishablebinding of three monoclonal antibodies recognising differentepftopes on CD2. The 3-D structure of rat CD2 domain 1 has beendetermined by NMR spectroscopy and X-ray crystallography, confirmingthe predictions from the sequence. Applications of this methodof insertion of disulphide residues for probing protein structuresare discussed, together with other structures of IgSF domainslacking the typical inter-sheet disulphide bond.     

Molecular modeling of single polypeptide chain of calcium-binding protein p26olf from dimeric S100B()

P26olf from olfactory tissue of frog, which may be involvedin olfactory transduction or adaptation, is a Ca²⁺-binding proteinwith 217 amino acids. The p26olf molecule contains two homologousparts consisting of the N-terminal half with amino acids 1–109and the C-terminal half with amino acids 110–217. Eachhalf resembles S100 protein with about 100 amino acids and containstwo helix–loop–helix Ca²⁺-binding structural motifsknown as EF-hands: a normal EF-hand at the C-terminus and apseudo EF-hand at the N-terminus. Multiple alignment of thetwo S100-like domains of p26olf with 18 S100 proteins indicatedthat the C-terminal putative EF-hand of each domain containsa four-residue insertion when compared with the typical EF-handmotifs in the S100 protein, while the N-terminal EF-hand ishomologous to its pseudo EF-hand. We constructed a three-dimensionalmodel of the p26olf molecule based on results of the multiplealignment and NMR structures of dimeric S100B(ßß)in the Ca²⁺-free state. The predicted structure of the p26olfsingle polypeptide chain satisfactorily adopts a folding patternremarkably similar to dimeric S100B(ßß). Each domainof p26olf consists of a unicornate-type four-helix bundle andthey interact with each other in an antiparallel manner formingan X-type four-helix bundle between the two domains. The twoS100-like domains of p26olf are linked by a loop with no sterichindrance, suggesting that this loop might play an importantrole in the function of p26olf. The circular dichroism spectraldata support the predicted structure of p26olf and indicatethat Ca²⁺-dependent conformational changes occur. Since theC-terminal putative EF-hand of each domain fully keeps the helix–loop–helixmotif having a longer Ca²⁺-binding loop, regardless of the four-residueinsertion, we propose that it is a new, novel EF-hand, althoughit is unclear whether this EF-hand binds Ca²⁺. P26olf is a newmember of the S100 protein family.     

Functionally essential, invariant glutamate near the C-terminus of strand {beta}5 in various ({alpha}/{beta})8-barrel enzymes as a possible indicator of their evolutionary relatedness

Twelve different (/ß)₈-barrel enzymes belonging tothree structurally distinct families were found to contain,near the C-terminus of their strand ß5, a conservedinvariant glutamic acid residue that plays an important functionalrole in each of these enzymes. The search was based on the ideathat a conserved sequence region of an (/ß)₈-barrelenzyme should be more or less conserved also in the equivalentpart of the structure of the other enzymes with this foldingmotif owing to their mutual evolutionary relatedness. For thispurpose, the sequence region around the well conserved fifthß-strand of a-amylase containing catalytic glutamate(Glu230, Aspergillus oryzae -amylase numbering), was used asthe sequence-structural template. The isolated sequence stretchesof the 12 (/ß)₈-barrels are discussed from both thesequence-structural and the evolutionary point of view, theinvariant glutamate residue being proposed to be a joining featureof the studied group of enzymes remaining from their ancestral(/ß)₈-barrel     

A predicted three-dimensional structure for the human immunodeficiency virus binding domains of CD4 antigen

A predicted three-dimensional structure of the two N-terminalextracellular domains of human CD4 antigen, a cell surface glycoprotein,is reported. This region of CD4, particularly the first domain,has been identified as containing the binding region for theenvelope gp120 protein of the human immuno-deficiency virus.The model was predicted based on the sequence homology of eachdomain with the variable light chain of immunoglobulins. Theframework ß-sheet regions were taken from the crystalcoordinates of REI. For one region in the first domain of CD4there was an ambiguity in the alignment with REI and two alternatemodels are presented. Loops connecting the framework were modeledfrom fragments selected from a database of main chain coordinatesfrom all known protein structures. Residues identified as involvedin binding gp120 have been located in several other studieswithin the first domain of CD4. Epitopes from eight monoclonalantibodies have been mapped onto residues in both domains. Competitionof these antibodies with each other and with gp120 can be interpretedfrom the structural model.     

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–ß.     

Main-chain complementarity in protein-protein recognition

We present a statistical analysis of protein structures basedon interatomic C_a distances. The overall distance distributionsreflect in detail the contents of sequence-specific substructuresmaintained by local interactions (such as -helixes) and longerrange interactions (such as disulfide bridges and ß-sheets).We also show that a volume scaling of the distances makes distancedistributions for protein chains of different length superimposable.Distance distributions were also calculated specifically foramino acids separated by a given number of residues. Specificfeatures in these distributions are visible for sequence separationsof up to 20 amino acid residues. A simple representation, whichpreserves most of the information in the distance distributions,was obtained using six parameters only. The parameters giverise to canonical distance intervals and when predicting coarse-graineddistance constraints by methods such as data-driven artificialneural networks, these should preferably be selected from theseintervals. We discuss the use of the six parameters for determiningor reconstructing 3-D protein structures.     

Family of G protein {alpha} chains: amphipathic analysis and predicted structure of functional domains

The G proteins transduce hormonal and other signals into regulationof enzymes such as adenylyl cyclase and retinal cGMP phosphodiesterase.Each G protein contains an subunit that binds and hydrolyzesguanine nucleotides and interacts with ß subunitsand specific receptor and effector proteins. Amphipathic andsecondary structure analysis of the primary sequences of fivedifferent chains (bovine _s, _t1 and _t2, mouse _i, and rat _o)predicted the secondary structure of a composite chain (_avg).The chains contain four short regions of sequence homologousto regions in the GDP binding domain of bacterial elongationfactor Tu (EF-Tu). Similarities between the predicted secondarystructures of these regions in _avg and the known secondary structureof EF-Tu allowed us to construct a three-dimensional model ofthe GDP binding domain of _avg. Identification of the GDP bindingdomain of _avg defined three additional domains in the compositepolypeptide. The first includes the amino terminal 41 residuesof _avg, with a predicted am phipathic helical structure; thisdomain may control binding of the chains to the ßcomplex. The second domain, containing predicted ßstrands and helices, several of which are strongly amphipathic,probably contains sequences responsible for interaction of chains with effector enzymes. The predicted structure of thethird domain, containing the carhoxy terminal 100 amino acids,is predominantly ß sheet with an amphipathic helixat the carboxy terminus. We propose that this domain is reponsiblefor receptor binding. Our model should help direct further experimentsinto the structure and function of the G protein chain.