Modelling the polypeptide backbone with 'spare parts' from known protein structures

An automatic procedure for building a protein polyalanine backbonefrom C positions and ‘spare parts’ retrieved froma data base of 66 high-resolution protein structures is described.Protein backbones are constructed from over-lapping fragmentsof variable length, which allows the backbone of regular secondarystructure elements to be built in one block. The procedure isshown to yield backbones which compare very favourably withthose from highly refined X-ray structures (r.m.s. deviationbetween generated and crystal structures <lÅ). Themethod is furthermore quite insensitive to experimental errorsin C positions as well as to the size of the data base, andis seen to yield valuable insight into the relationships betweensequence and 3-D structure: one example on triose phosphateisomerase, a ß-barrel protein, shows that ßloops can be considered as structurally more uncommon than ßloops. The ‘spare parts’ approach is also foundto be useful for general-purpose modelling of local structuralchanges produced by insertion or deletion of residues. It should,however, be used with caution. Crude selection criteria basedsolely on fragment length and geometric fit to the loop baseregions yield realistic backbones in about two-thirds of thetest cases (r.m.s. deviations from refined crystal structure{small tilde}lÅ). In the remaining cases, sequence information,in particular the presence of glycine residues which tend toadopt more unusual backbone conformations, must be consideredto obtain comparable results.     

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.     

Building protein backbones from C{alpha} co-ordinates

An automatic procedure for building polyalanine backbones fromguiding -carbon positions is presented. Polyalanine backbonesare built based on the geometric restraints of angle N-C-C andthe knowledge of main-chain dihedral angle distributions. Abuilding module constructs a list of polyalanine backbones thatfollow exactly the C trace. Then a selection module selectsone backbone with the largest portion of phi-psi pairs in favouredregions. Several test cases on C coordinates from X-ray refinedstructures give acceptable results. Less than 10% of the peptideplanes are incorrectly built, and the result is not sensitiveto random shift up to 0.5 Å of C coordinates.     

Structural interpretation of site-directed mutagenesis and specificity of the catalytic subunit of protein kinase CK2 using comparative modelling

The catalytic subunit of protein kinase casein kinase 2 (CK2),which has specificity for both ATP and GTP, shows significantamino acid sequence similarity to the cyclin-dependent kinase2 (CDK2). We constructed site-directed mutants of CK2 and useda three-dimensional model to investigate the basis for the dualspecificity. Introduction of Phe and Gly at positions 50 and51, in order to restore the pattern of the glycine-rich motif,did not seriously affect the specificity for ATP or GTP. Weshow that the dual specificity probably originates from theloop situated around the position His115 to Asp120 (HVNNTD).The insertion of a residue in this loop in CK2 subunits, comparedwith CDK2 and other kinases, might orient the backbone to interactwith the base A and G; this insertion is conserved in all knownCK2. The mutant N118, the design of which was based on the modelling,showed reduced affinity for GTP as predicted from the model.Other mutants were intended to probe the integrity of the catalyticloop, alter the polarity of a buried residue and explore theimportance of the carboxy terminus. Introduction of Arg to replaceAsn189, which is mapped on the activation loop, results in amutant with decreased k_cat, possibly as a result of disruptionof the interaction between this residue and basic residues inthe vicinity. Truncation at position 331 eliminates the last60 residues of the subunit and this mutant has a reduced catalyticefficiency compared with the wild-type. Catalytic efficiencyis restored in the truncation mutant by the replacement of apotentially buried Glu at position 252 by Lys, probably owingto a higher stability resulting from the formation of a saltbridge between Lys252 and Asp208.     

Model building of disulfide bonds in proteins with known three-dimensional structure

As an aid in the selection of sites in a protein where a disulfidebond might be engineered, a computer program has been developed.The algorithm starts with the generation of Cß positionsfrom the N, C and C atom coordinates available from a three-dimensionalmodel. A first set of residue pairs that might form a disulfidebond is selected on the basis of Cß–Cßdistances between residues. Then, for each residue in this set,S positions are generated, which satisfy the requirement that,with ideal values for the C–Cß and Cß–Sbond lengths and for the bond angle at Cß, the distancebetween S of residue 1 and Cß of residue 2 in a pair(determined by the bond angle at S2) is at, or very close toits ideal value. Usually two acceptable S positions are foundfor each half cystine, resulting in up to four different conformationsfor the disulfide bond. Finally, these conformations are subjectedto an energy minimization procedure to remove large deviationsfrom ideal geometry and their final energies are calculated.User input determines which final conformations are energeticallyacceptable. These conformations are written to a file to allowfurther analysis and e.g. inspection on a computer graphicsdevice.     

A molecular dynamics simulation of bacteriophage T4 lysozyme

An analysis of a 400 ps molecular dynamics simulation of the164 amino acid enzyme T4 lysozyme is presented. The simulationwas carried out with all hydrogen atoms modeled explicitly,the inclusion of all 152 crystallographic waters and at a temperatureof 300 K. Temporal analysis of the trajectory versus energy,hydrogen bond stability, r.m.s. deviation from the startingcrystal structure and radius of gyration, demonstrates thatthe simulation was both stable and representative of the averageexperimental structure. Average structural properties were calculatedfrom the enzyme trajectory and compared with the crystal structure.The mean value of the C displacements of the average simulatedstructure from the X-ray structure was 1.1 ± 0.1 Å;differences of the backbone and angles between the averagesimulated structure and the crystal structure were also examined.Thermal-B factors were calculated from the simulation for heavyand backbone atoms and both were in good agreement with experimentalvalues. Relationships between protein secondary structure elementsand internal motions were studied by examining the positionalfluctuations of individual helix, sheet and turn structures.The structural integrity in the secondary structure units waspreserved throughout the simulation; however, the A helix didshow some unusually high atomic fluctuations. The largest backboneatom r.m.s. fluctuations were found in non-secondary structureregions; similar results were observed for r.m.s. fluctuationsof non-secondary structure and angles. In general, the calculatedvalues of r.m.s. fluctuations were quite small for the secondarystructure elements. In contrast, surface loops and turns exhibitedmuch larger values, being able to sample larger regions of conformationalspace. The C difference distance matrix and super-positioninganalyses comparing the X-ray structure with the average dynamicsstructure suggest that a ‘hinge-bending’ motionoccurs between the N- and C-terminal domains.     

Inhibition of {beta}S-chain dependent polymerization by synergistic complementation of contact site perturbations of {alpha}-chain: application of semisynthetic chimeric {alpha}-chains

Mouse _1–30-horse _31–141 chimeric -chain, a semisyntheticsuper-inhibitory -chain, inhibits ß^S-chain dependent polymerizationbetter than both parent -chains. Although contact site sequencedifferences are absent in the _1–30 region of the chimericchain, the four sequence differences of the region _17-22 couldinduce perturbations of the side chains at ₁₆, ₂₀ and ₂₃, thethree contact sites of the region. A synergistic complementationof such contact site perturbation with that of horse _31–141probably results in the super-inhibitory activity of the chimeric-chain. The inhibitory contact site sequence differences, bythemselves, could also exhibit similar synergistic complementation.Accordingly, the polymerization inhibitory activity of Hb Le-Lamentin(LM) mutation [His20()Gln], a contact site sequence difference,engineered into human–horse chimeric -chain has been investigatedto map such a synergistic complementation. Gln20() has littleeffect on the O₂ affinity of HbS, but in human–horse chimeric-chain it reduces the O₂ affinity slightly. In the chimeric-chain, Gln20() increased sensitivity of the ßßcleft for the DPG influence, reflecting a cross-talk betweenthe ₁ß₁ interface and ßß cleft in this semisyntheticchimeric HbS. In the human -chain frame, the polymerizationinhibitory activity of Gln20() is higher compared with horse_1–30, but lower than mouse _1–30. Gln20() synergisticallycomplements the inhibitory propensity of horse _31–141.However, the inhibitory activity of LM–horse chimeric-chain is still lower than that of mouse–horse chimeric-chain. Therefore, perturbation of multiple contact sites inthe _1–30 region of the mouse–horse chimeric -chainand its linkage with the inhibitory propensity of horse _31–141has been now invoked to explain the super-inhibitory activityof the chimeric -chain. The `linkage-map' of contact sites canserve as a blueprint for designing synergistic complementationof multiple contact sites into -chains as a strategy for generatingsuper-inhibitory antisickling hemoglobins for gene therapy ofsickle cell disease.     

Expression of recombinant {alpha}Ains-crystallin and not {alpha}A-crystallin inhibits bacterial growth

A-Crystallin and A^ins-crystallin are derived from the A-crystallingene via alternative splicing. They are identical except forthe presence of a polypeptide, 23 amino acids long, encodedby the ‘insert’ exon. Evolutionary logic would suggestthat the insertion of a 23 amino acid peptide in the middleof A-crystallin, a protein evolving more slowly than eitherhistone H1, cytochrome c or hemoglobin, would lead to appreciablestructural and functional changes. However, based on physico-chemicalstudies, it is presently believed that A-crystallin and A^ins-crystallinare functionally equivalent and that the presence of the ‘insert’peptide in A^Ins-crystallin is inconsequential. We report herethat the independent expression of recombinant A^Ins-crystallin,and not A-crystallin, inhibits growth of the bacterial host.These observations were confirmed in co-expression experiments,wherein both the proteins were expressed in the same cell. Interestingly,growth inhibition is reversible. Importantly, the data demonstratethat it is catalytic amounts and not the gross accumulationof A^Ins-crystalline which causes growth inhibition. Given theprior knowledge that A-crystallin and A^Ins-crystallin differby a peptide of 23 amino acids, these data suggest that the‘insert peptide’ in A^Ins-crystallin imparts propertieson this protein that are different from A-crystallin.     

Effect of amino acid deletions in the O-glycosylated region of Aspergillus awamori glucoamylase

Aspergillus awamori glucoamylase (GA) contains globular catalyticand starch-binding domains (residues 1–471 and 509–616,respectively). A heavily O-glycosylated sequence comprises twoparts. The first (residues 441–471) in the crystal structurewraps around an /-barrel formed by residues 1–440. Thesecond (residues 472–508) is an extended, semi-rigid linkerbetween the two domains. To investigate the functional roleof this linker, we made internal deletions to remove residues466–512 (GA1), 485–512 (GA2) and 466–483 (GA3).GA2 and GA3 were expressed in Saccharomyces cerevisiae culturesupernatants at 60 and 20% the wild-type level, respectively,while GA1 was almost undetectable. Western blots comparing extracellularand intracellular fractions indicated that the region deletedin GA3 was critical for secretion, while the region deletedin GA2 contributed to the production of a stable enzyme structure.The activities of purified GA2 and GA3 on soluble and insolublestarch were similar to those of wild-type GA, indicating thatfor soluble starch their deletions did not affect the catalyticdomain and for insoluble starch the linker does not coordinatethe activities of the catalytic and starch-binding domains.The deletions had a significant negative effect on GA2 and GA3thermos tabilities.     

Mutational analysis of structure--activity relationships in human tumor necrosis factor-alpha

To determine the region of human tumor necrosis factor-alpha(TNF-), essential for cytotoxic activity against mouse L-M cells,single amino-acid-substituted TNF- mutant proteins (muteins)were produced in Escherichia coli by protein engineering techniques.An expression plasmid for TNF- was mutagenized by passage throughan E.coli mutD5 mutator strain and by oligonucleotide-directedmutagenesis. Approximately 100 single amino-acid-substitutedTNF- muteins were produced and assayed for cytotoxic activity.The cytotoxic activities of purified TNF- muteins, e.g. TNF-31T,-32Y, -82D, -85H, -115L, -141Y, -144K and -146E, were < 1%of that of parent TNF-. These results indicate that the integrityof at least four distinct regions of the TNF- molecule is requiredfor full biological activity. These regions are designated asfollows: region I, from position 30 to 32; region II, from position82 to 89; region III, from position 115 to 117; region FV, fromposition 141 to 146. In addition, TNF-141Y could not completelycompete with parent TNF- for binding to the receptor. This demonstratesthat region IV, and at least aspartk acid at position 141, mustbe involved in the TNF receptor binding site.     

The stability and unfolding of an IgG binding protein based upon the B domain of protein A from Staphylococcus aureus probed by tryptophan substitution and fluorescence spectroscopy

The stability and unfolding of an immunoglobulin (Ig) G bindingprotein based upon the B domain of protein A (SpAB) from Staphylococcusaureus were studied by substituting tryptophan residues at strategiclocations within each of the three a-helical regions (al-a3)of the domain. The role of the C-terminal helix, a3, was investigatedby generating two protein constructs, one corresponding to thecomplete SpAB, the other lacking a part of ct3; the Trp substitutionswere made in both one-and two-domain versions of each of theseconstructs. The fluorescence properties of each of the single-tryptophanmutants were studied in the native state and as a function ofguanidine-HCl-mediated unfolding, and their IgG binding activitieswere determined by a competitive enzyme-linked immunosorbentassay. The free energies of folding and of binding to IgG foreach mutant were compared with those for the native domains.The effect of each substitution upon the overall structure andupon the IgG binding interface was modelled by molecular graphicsand energy minimization. These studies indicate that (i) 3 contributesto the overall stability of the domain and to the formationof the IgG binding site in l and 2, and (ii) al unfolds first,followed by 2 and 3 together.     

Three-dimensional structure of protein C inhibitor predicted from structure of {alpha}1-antitrypsin with computer graphics

The three-dimensional structure of a proteolytically modifiedprotein C inhibitor, a member of the serine protease inhibitorsuperfamily, was constructed with computer graphics based onits amino acid sequence homology with that of the modified 1-antitrypsinwhose structure had been elucidated by X-ray crystallography.The intact form of protein C inhibitor was predicted with an-carbon model based on its hydrophilicity and hydrogen bondpattern. Furthermore, a model of its interaction with activatedprotein C was constructed based on the structure of the complexbetween trypsin and its inhibitor, which had been determinedby X-ray crystallography.     

Semi-empirical simulation of Zn/Cd binding site preference in the metal binding domains of mammalian metallothionein

Metallothionein, a two-domain protein, naturally binds sevengram atoms of divalent ions such as Zn and Cd. Four of the metals(Ml, M5, M6 and M7) are found in the -domain and three (M2,M3 and M4) in the ß-domain. Previous studies haveshown that metals in the -domain are more readily exchangeable,and the level of avidity is site specific. By semi-empiricalMNDO modified neglect of diatomic overlap calculations, we foundthe tendency of binding energy for Cd to be M3 > M2 >M4 in the ß-cluster and M5 > M7 > Ml, M6 inthe -cluster. Thus, the replacement of Zn by Cd can be expectedto follow the order M4 M2 M3 in the ß-domain andMS M7 M1 or M6 in the -domain. This is reflected by energydifferences computed with a series of simulated structures derivedfrom either X-ray crystallography or NMR coordinates.     

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.     

The differences in conformation between {alpha}-human atrial natriuretic polypeptide, {alpha}-hANP, and its derivative, Met(O)-{alpha}-hANP, in solution

The differences in conformation between -human atrial natriureticpolypeptide (-hANP) and its inactive analog, Met(O)--hANP, havebeen analyzed by nuclear magnetic resonance spectroscopy. Allproton resonances for both peptides were assigned by means ofthe sequential assignment procedure. The three-dimensional structureof -hANP in solution had previously been determined by distancegeometry calculation using distance constraints derived fromnuclear Overhauser effects (NOEs). Here, the three-dimensionalstructure of Met(O)--hANP was determined. The conformationaldifferences between these two molecules were as follows: threesegments of -hANP, Serl–Cys7, Arg11–Ala17 and Gln18–Tyr28,have some ordered structures. In Met(O)--hANP the Gln18-Tyr28region has a similar conformation, while the remaining two regionsdo not have the ordered structure found in -hANP. It is suggestedthat the conserved conformation of the Gln18–Tyr28 regionis required for binding to the ANP receptor and that the slightbiological activity of Met(O)-a-hANP is due to loss of the orderedstructures evoked in the Serl–Cys7 and Arg11–Ala17regions of -hANP.     

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.     

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     

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 and functional domains in human tumour necrosis factors

Human tumour necrosis factors (hTNFs) and ß are relatedpleiotropic cytokines which share many activities and competewith each other for binding to two receptor components on manycell types. Although structural and biological data indicatethat the active form of hTNF- may be a symmetrical trimer, themanner in which hTNFs interact with their receptors to triggera myriad of cell type-dependent responses is not clear. A combinationof chemical modification, epitope mapping and site-directedmutagenesis approaches suggest that at least four distinct peptidesequences are Important for the biological activity of hTNF-.In particular, certain peptide sequences between amino acidpositions 11 and 35 in hTNF- appear to be critical for receptorbinding and triggering biological responses. The recent cloningof the two hTNF-/ß receptors opens the way for precisemapping of the functional domains in hTNFs     

Site-directed mutagenesis and X-ray crystallography of two phospholipase A2 mutants: Y52F and Y73F

Tyr52 and Tyr73 are conserved amino acid residues throughoutall vertebrate phospholipases A₂. They are part of an extendedhydrogen bonding system that links the N-terminal -NH⁺₃ -groupto the catalytic residues His48 and Asp99. These tyrosines werereplaced by phenylalanines in a porcine pancreatic phospholipaseA₂ mutant, in which residues 62–66 had been deleted (62–66PLA₂).The mutations did not affect the catalytic properties of theenzyme, nor the folding kinetics. The stability against denaturatlonby guanidine hydrochloride was decreased, however. To analysehow the enzyme compensates for the loss of the tyrosine hydroxylgroup, the X-ray structures of the Y52F and AY73F mutants weredetermined. After crystallographic refinement the final crystallographicR-factors were 18.1% for the %Y52F mutant (data between 7 and2.3 Å resolution) and 19.1% for the Y73F mutant (databetween 7 and 2.4 Å resolution). No conformational changesoccurred in the mutants compared with the 62–66PLA₂, butan empty cavity formed at the site of the hydroxyl group ofthe former tyrosine. In both mutants the Asp99 side chain losesone of its hydrogen bonds and this might explain the observeddestabilization.