Stereochemical modeling of disulfide bridges. Criteria for introduction into proteins by site-directed mutagenesis

A computer modeling procedure for assessing the stereochemicalsuitability of pairs of residues in proteins as potential sitesfor introduction of cystine disulfide crosslinks has been developed.Residue pairs with C – C distances of 6.5 Å andC^beta;–C^ß distances of 4.5 Å are chosenfor geometrical fixation of S atoms using the program MODIP.The stereochemistry of the modeled disulfides is evaluated usinglimits for the structural parameters of the various torsionangles and S–S bond length in the disulfide bridge. Theability of the procedure to correctly model disulfides has beenchecked with examples of cystine peptides of known crystal structuresand 103 disulfide bridges from 25 available protein crystalstructures determined at 2 Å resolution. An analysis ofresults on three proteins with engineered disulfides, T4 lysozyme,dihydrofolate reductase and subtilisin, is presented. Two positionsfor the introduction of ‘stereochemically optimal’disulfides are identified in subtilisin.     

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.     

The T {leftrightharpoons} R structural transition of insulin; pathways suggested by targeted energy minimization

The transition of insulin between its crystallographically definedstates T and R is connected with considerable change even ofbackbone structure: the N-terminal B chain (residues B1 –B8) refolds from extended conformation in T into helical inR, and vice versa. Although hitherto observed only in hexamersthe transition of the monomer was adequate for developing andtesting the method of ‘targeted energy minimization’(TEM), capable of coping with conformational changes of suchextent at moderate computational expenditure. The simulationis performed in a predetermined number of steps consisting oftwo atomic displacements each, one by force in the directionof the target structure, the second by energy minimization releasingthe constraint caused in the first. The transition pathway isrepresented by the string of energy minimized transient structures.Due to the directedness of the algorithm the simulated pathwayfor R T is not the reversal of that for T R. It is, therefore,not pretended that the minimum energy pathway was identified.In the T R direction the N-terminal B chain first swivels whileremaining largely stretched and then winds up extending thepre-existing helix B9–B19. The A chain advances into thespace abandoned and withdraws from it in the R T simulation.In the latter the extended helix first kinks at B8/B9, and thenthe B1 B8 segment is unwound and stretched. The helical H-bondsof that segment are formed late in T R and are maintained duringalmost half of R T. The A_N helix is less stable and more involvedin the transitions than helix A_C. The two pathways seem plausiblefrom both the energetic and geometric points of view. Knowledgeof them will be of value to suggest mutations to test them byexperimental evidence.     

Effect of replacing helical glycine residues with alanines on reversible and irreversible stability and production of Aspergillus awamori glucoamylase

To decrease irreversible thermoinactivation of Aspergillus awamoriglucoamylase, five Gly residues causing helix flexibility werereplaced with Ala residues. Mutation of Gly57 did not affectthermostability. Mutation of Gly137 doubled it at pHs 3.5 and4.5 but barely changed it at pH 5.5. The Gly139Ala mutationdid not change thermostability at pH 3.5, improved it at pH4.5 and worsened it at pH 5.5. The Gly137/Gly139Ala/Ala mutationgave 1.5–2-fold increased thermostabilities at pHs 3.5–5.5.Mutations of Gly251 and Gly383 decreased it at all pHs. Gly137Alaand Gly137/Gly139Ala/Ala glucoamylases are the most stable yetproduced by mutation. Guanidine treatment at pH 4.5 decreasedthe reversible stabilities of Gly137Ala, Gly139Ala and Gly137/Gly139Ala/Alaglucoamylases at infinite dilution while not changing thoseof Gly251Ala and Gly 383Ala glucoamylases, which is, in general,opposite to what occurred with thermoinactivation. Mutationof Gly57 greatly improved the extracellular glucoamylase productionby yeast, that of Gly137 barely affected it and those of Gly139and of both Gly137 and Gly139 strongly impeded it. These observationssuggest that -helix rigidity can affect reversible and irreversibleglucoamylase stability differently, that the effects of multiplemutations within one -helix to improve stability are not alwaysadditive and that even single mutations can strongly affectextracellular enzyme production.     

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.     

Folding and stability of the active N-terminal domain of tissue inhibitor of metalloproteinases-1 and -2

The truncated forms of tissue inhibitor of metalloproteinase-1and -2 (TIMP-1 and -2), comprising the N-terminal active domain,are ideal molecules for structural analysis by intrinsic fluorescenceas each contains a single conserved tryptophan residue. In thispaper we describe studies on their conformational stability,unfolding/refolding kinetics and the environment of the uniquetryptophan as judged by its fluorescence properties in the nativestate and exposure to an external quencher, acrylamide. Twoforms of TIMP-2 were studied: TIMP-2 T21 derived from the full-lengthcDNA clone isolated from a mixed-tumour library, and TIMP-2A21 containing the highly conserved V18IRAK22 sequence. In allthree TIMP proteins the tryptophan environments in the nativestate appeared to be similar, but substantial differences wereseen in their conformational stabilities and refolding kinetics.TIMP-1 was approximately twice as stable as TIMP-2 T21 and 1.4-foldmore stable than TIMP-2 A21. This stability difference betweenTIMP-1 and TIMP-2 was shown to be independent of N-linked glycosylation.TTMP-1 and TIMP-2 A21 both showed simple two-state refoldingkinetics, whereas TIMP-2 T21 refolding was more complex andbiphasic in character. These differences between TIMP-2 T21and A21 suggest that residue 21 is a structurally importantsite in the TIMP protein.All three truncated molecules can beconsidered as stable independent folding domains ideally suitedfor further structural analysis     

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.     

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.     

The use of alanine scanning mutagenesis to determine the role of the N-terminus of the regulatory chain in the heterotropic mechanism of Escherichia coli aspartate transcarbamoylase

The location of the first seven residues of the regulatory chainof Escherichia coli aspartate transcarbamoylase has been identifiedby X-ray crystallography to be near the binding site of theregulatory nucleotides. In order to determine the function ofthe N-terminus of the regulatory chain of aspartate transcarbamoylasein heterotropic regulation, alanine scanning mutagenesis wasused. Specifically, Thr2r, His3r, Asp4r, Asn5r, Lys6r and Leu7rwere each replaced with alanine. Analyses of these mutant enzymesindicate that none of these substitutions significantly alterthe catalytic properties of the enzyme. However, three of themutant enzymes, Asp4r Ala, Lys6r Ala and Leu7r Ala, exhibitednotable changes in their response to the regulatory nucleotides,while mutations at Thr2r, His3r and Asn5r exhibited only minorchanges in their heterotropic responses. For the Asp4r Alaenzyme, the responses to ATP and CTP were reduced 30 and 40%respectively, compared with the wild-type enzyme. For the Lys6r Ala enzyme, the response to ATP was reduced 70%, while theCTP response was reduced 50%. In the case of the Leu7r Alaenzyme, a 30 and 20% reduction in response to ATP and CTP respectively,was observed. The synergistk inhibition by UTP in the presenceof CTP for the Lys6r Ala enzyme was reduced 40% compared withthat of the wild type enzyme. For the Leu7r Ala enzyme, thesynergistic inhibition was abolished. In addition, UTP decreasedthe CTP binding affinity of the Leu7r Ala enzyme. Analysisof the kinetic data from these mutant enzymes suggests thatresidues Thr2r, His3r and Asn5r have little effect on the heterotropicmechanism, while residues Asp4r, Lys6r and Leu7r play a moresignificant role in the heterotropic response of the enzymetoward the nucleotides. Furthermore, residue Leu7r appears tobe directly involved in the mechanism for synergistic inhibitionof aspartate transcarbamoylase. In this study alanine scanningmutagenesis has provided a rapid method of identifying thoseresidues in the N-terminal region of the regulatory chain ofaspartate transcarbamoylase important for heterotropic regulation.     

Thermosensitive mutants of Aspergillus awamori glucoamylase by random mutagenesis: inactivation kinetics and structural interpretation

Abstract Seven thermosensitive glucoamylase mutants generated by randommutagenesis and expressed inSaccharomyces cerevisiae were sequencedand their inactivation kinetics were determined. Wild-type glucoamylaseexpressed in S.cerevisiae was more glycosylated and more stablethan the native Aspergillus niger enzyme. All mutants had lowerfree energies of inactivation than wild-type glucoamylase. Inthe Ala39 Val, Ala302 Val and Leu410 Phe mutants, small hydrophobicresidues were replaced by larger ones, showing that increasesin size and hydrophobicity of residues included in hydrophobicclusters were destabilizing. The Gly396 Ser and Gly407 Aspmutants had very flexible residues replaced by more rigid ones,and this probably induced changes in the backbone conformationthat destabilized the protein. The Prol28 Ser mutation changeda rigid residue in an a-helix to a more flexible one, and destabilizedthe protein by increasing the entropy of the unfolded state.The Ala residue in the Ala442 Thr mutation is in the highlyO-glycosylated region surrounded by hydrophilk residues, whereitmay be a hydrophobic anchor Unking the O-glycosylated arm tothe catalytic core. It was replaced by a residue that potentiallyis O-glycosylated. In five of the seven mutations, residuesthat were part of hydrophobic microdomains were changed, confirmingthe importance of the latter in protein stability and structure     

A point mutation that decreases the thermal stability of human interferon {gamma}

We have identified a mutation of human gamma-interferon (IFN)causing a temperature-sensitive phenotype. We used a randomizedoligonucleotide to mutagenize a synthetic human IFN gene, thenscreened the resulting mutants produced in Escherichia colifor proteins with altered biological activity. One mutant proteinselected for detailed characterization exhibited < 0.3% ofthe specific biological activity of native IFN in an antiviralactivity assay performed at 37°C. However, the protein boundthe human IFN receptor with native efficiency at 4°C. Sequencingthe plasmid DNA encoding this protein snowed that the mutationchanged the lysine residue at amino acid 43 to glutamic acid(IFN/K43E). Site-specific mutagenesis at amino acid 43 showedthat this protein's phenotype resulted from positioning a negativecharge at position 43. Structural characterization of IFN/K43Eusing CD demonstrated that the protein had native conformationat 25°C, but assumed an altered conformation at 37°C.IFN/K43E in this altered conformation bound poorly to the IFNreceptor at 37°C, providing a rationale for the mutant'sdecreased antiviral activity.     

Computer simulations of signal transduction mechanism in {alpha}1B-adrenergic and m3-muscarinic receptors

Molecular dynamices simulations of the hamster _1Badrenergicand the rat m3-muscarinic seven-helix bundle receptor modelshave been carried out. The free, agonist-bound and antagonist-boundforms have been considered. Moreover, three mutant forms ofthe m3-muscarinic recep-tor (N507A, N507D and N507S) have alsobeen simulated; among these, the N507S mutant shows a constitutiveactivity. A comparative structural/dynamics analysis has beenperformed to elucidate (i) the perturbations induced by thefunctionally different ligands upon binding to their targetreceptor, (ii) the features of the three single-point mutantswith respect to the receptor wild type and (iii) the propertiesshared by the agonist-boundforms of the _1B-adrenergic receptorand the m3-muscarinic receptor and by the constitutively activemutant N507S. The consistency obtained between the structuralrearrangement of the transmembrane seven-helix bundle modelsconsidered, and the experimental pharmacological efficaciesof the ligands and of the mutants, constitute an important validationof the 3-D models obtained and allow the inference of the mechanismof ligand- or mutation-induced receptor activation at the molecularlevel.     

Effects of amino acid replacements around the reactive site of chicken ovomucoid domain 3 on the inhibitory activity toward chymotrypsin and trypsin

We have previously shown that replacing the P1-site residue(Ala) of chicken ovomucoid domain 3 (OMCHI3) with a Met or Lysresults in the acquisition of inhibitory activity toward chymotrypsinor trypsin, respectively. However, the inhibitory activitiesthus induced are not strong. In the present study, we introducedadditional amino acid replacements around the reactive siteto try to make the P1-site mutants more effective inhibitorsof chymotrypsin or trypsin. The amino acid replacement AspTyrat the P2' site of OMCHI3(P1Met) resulted in conversion to a35000-fold more effective inhibitor of chymotrypsin with aninhibitor constant (K_i) of 1.17x10^–11 M. The K_i valueof OMCHI3(P1Met, P2'Ala) indicated that the effect on the interactionwith chymotrypsin of removing a negative charge from the P2'site was greater than that of introducing an aromatic ring.Similarly, enhanced inhibition of trypsin was observed whenthe AspTyr replacement was introduced into the P2' site of OMCHI3(P1Lys).Two additional replacements, AspAla at the P4 site and ArgAlaat the P3' site, made the mutant a more effective inhibitorof trypsin with a K_i value of 1.44x10^–9 M. By contrast,ArgAla replacement at the P3' site of OMCHI3(P1Met, P2'Tyr)resulted in a greatly reduced inhibition of chymotrypsin, andAspAla replacement at the P4 site produced only a small changewhen compared with a natural variant of OMCHI3. These resultsclearly indicate that not only the P1-site residue but alsothe characteristics, particularly the electrostatic properties,of the amino acid residues around the reactive site of the proteaseinhibitor determine the strength of its interactions with proteases.Furthermore, amino acids with different characteristics arerequired around the reactive site for strong inhibition of chymotrypsinand trypsin.     

Chemical synthesis of {gamma}-echistatin analogues: elucidation of status of C-terminal (45-49)

Three analogues of -echistatin, des(45–49)--echistatin,des(46–49)-y-echistatin and des(47–49)--echistatin,were synthesized by solid-phase methodology and their biologicalactivities were measured and compared. The results reveal thatwithout the C-terminal (45–49) of -echistatin, the foldingof the protein to the final active structure is not interferedwith and Lys-45 influences the inhibition of platelet aggregation.     

Engineering of protein epitopes: a single deletion in a snake toxin generates full binding capacity to a previously unrecognized antibody

Structural features associated with the ability of a monoclonalantibody (mAb) to discriminate between protein variants areidentified and engineered. The variants are the curaremimetictoxin from Naja nigricollis and erabutoxin a or b from Laticaudasemifasciata which differ from each other by 16 substitutionsand one insertion. The neutralizing mAb M1 recognizes with highaffinity a topographical epitope on the surface of toxin , butfails to recognize the erabutoxins although they possess mostof the residues forming the presumed epitope. Examinations ofthe toxin and erabutoxin 3-D structures and molecular dynamicssimulations reveal several differences between the variants.In particular, the region involving the ß-turn 17–24is organized differently. Analysis of the differences foundin this region suggests that the insertion (or deletion) atposition 18 of the variant amino add sequences is particularlyimportant in determining the differential cross-reactivity.To test this proposal, residue 18 was deleted in one erabutoxinusing sitedirected mutagenesis, and the biological propertiesof the resulting mutant were examined. We found that full antigenicitywas restored in the previously unrecognized variant. The implicationsof this finding are discussed.     

Altered specificities of genetically engineered {alpha}1 antitrypsin variants

Seven active site variants of human ₁-antitrypsin (₁AT) wereproduced in Escherichia coli following site-specific mutagenesisof the ₁AT complementary DNA. ₁AT (Ala ³⁵⁸), ₁AT (Ile³⁵⁸ and₁AT (Val³⁵⁸), were efficient inhibitors of both neutrophil andpancreatic elastases, but not of cathepsin G. ₁AT (Ala³⁵⁸, Val³⁵⁸)and ₁AT (Phe³⁵⁸ specifically inhibited pancreatic elastase andcathepsin G respectively. The most potent inhibitor of neutrophilelastase was ₁AT (Leu³⁵⁸), which also proved to be effectiveagainst cathepsin G. The ₁AT (Arg³⁵⁸) variant inactivated thrombinwith kinetics similar to antithrombin III in the presence ofheparin. Electrophoretic analysis showed that SDS-stable highmol. wt complexes were formed between the mutant inhibitorsand the cognate proteases in each case. These data indicatethat effective inhibition occurs when the ₁AT P1 residue (position358) corresponds to the primary specificity of the target protease.Moreover, alteration of the P3 residue (position 356) can furthermodify the reactivity of the inhibitor. Two of the variantshave therapeutic potential: ₁AT (Leu³⁵⁸ may be more useful thanplasma ₁AT in the treatment of destructive lung disorders and₁ (Arg³⁵⁸ could be effective in the control of thrombosis.     

Structural study of mutants of Escherichia coli ribonuclease HI with enhanced thermostability

Systematic replacement of the amino acid residues in Escherichiacoli ribonuclease HI with those in the thermophilic counterparthas revealed that two mutations, His62–Pro (H62P) andLys95Gly (K95G), increased the thermostability of the protein.These single-site mutant proteins, together with the mutantproteins His62Ala (H62A), Lys95Asn (K95N) and Lys95Ala (K95A),were crystallized and their structures were determined at 1.8Å resolution. The crystal structures of these mutant proteinsreveal that only the local structure around each mutation siteis essential for the increase in thermostability. For each mutantprotein, the stabilization mechanism is considered to be asfollows: (i) H62P is stabilized because of a decrease in theentropy of the unfolded state, without a change in the nativebackbone structure; (ii) K95G is stabilized since the straincaused by the left-handed backbone structure in the typical3:5 type loop is eliminated; and (iii) K95N is slightly stabilizedby a hydrogen bond formed between the side-chain N-atom of themutated aspargine residue and the main-chain carbonyl oxygenwithin the same residue.     

Co-enzyme specificity of 3-isopropylmalate dehydrogenase from Thermits thermophilus HB8

The co–enzyme specificity of 3–isopropylmalate dehydrogenasefrom an extreme thermophile, Thermus thermophilus HB8, was changedfrom NAD to NADP by site–directed mutagenesis Based onsequence comparison of 3–isopropylmalate dehydrogenasesfrom various organisms with NAD– and NADP–dependentisocitrate dehydrogenases, Ser226, Ser253 and De279 of 3–isopropylmalatedehydrogenase were suggested as determining the co–enzymespecificity. These residues were replaced with the correspondingresidues of NADP–dependent isocitrate dehydrogenases;Arg, Gly and Tyr respectively. The single–mutated enzymes,S226R and I279Y, enhanced the activities towards NADP 10–and 3–fold respectively, whereas S253G reduced the activity.Among the multiple–mutated enzymes, the double–mutatedS226R/I279Y increased the catalytic efficiency against NADP( fold) and shifted the specificity for NAD towards NADP mostsignificantly ( 173–fold).     

Free energy perturbation study on a Trp-binding mutant (Ser88 ->Cys) of the trp-repressor

The Ser⁸⁸Cys mutant of the trp-repressor showed a lower affinityfor the corepressor than the wild-type repressor [G = 1.7 ±0.3 kcal/mol, Chou and Matthews (1989) J. Biol. Chem., 264,18314–18319].A molecular dynamics/free energy cycle perturbation study wasperformed to understand the origin of the decreased affinity.A value (G = 1.58 ± 0.28 kcal/mol) comparable with theexperimental value was obtained by the simulation. Free energycomponent analysis revealed that destabilization of the vander Waals interaction between Ser⁸⁸ and Trp¹⁰⁹ (corepressor)mainly contributed to the decreased affinity of the mutant.The rotational transition of the hydroxyl (sulfhydryl) groupof Ser⁸⁸ (Cys⁸⁸) during the simulations affected the contributionsof Arg⁸⁴ and water to the free energy change in the aporepressorand those of Arg⁸⁴ and Trp ¹⁰⁹ to that in the holorepressor.However, the contributions from different residues compensatedeach other, and the total free energy changes were almost invariablein the various simulations.     

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.