Thermodynamic studies on the interaction of the third complement component and its inhibitor, compstatin

Compstatin is a 13-residue cyclic peptide that inhibits complement activation by binding to complement component, C3. Although the activity of compstatin has been improved severalfold using combinatorial and rational design approaches, the molecular basis for its interaction with C3 is not yet fully understood. In the present study, isothermal titration calorimetry was employed to dissect the molecular forces that govern the interaction of compstatin with C3 using four different compstatin analogs. Our studies indicate that the C3-compstatin interaction is an enthalpy-driven process. Substitution of the valine and histidine residues at positions 4 and 9 with tryptophan and alanine, respectively, resulted in the increase of enthalpy of the interaction, thereby increasing the binding affinity for C3. The data also suggest that the interaction is mediated by water molecules. These interfacial water molecules could be the source for unfavorable entropy and large negative heat capacity changes observed in the interaction. Although part of the negative heat capacity changes could be accounted for by the water molecules, the rest might be resulting from the conformational changes in C3 and/or compstatin up on binding. Finally, we propose based on the pK(a) values determined from the protonation studies that histidine on compstatin participates in protonation changes and contributes to the specificity of the interaction between compstatin and C3. These protonation changes vary significantly between the binding of different compstatin analogs to C3.     

Structure of compstatin in complex with complement component C3c reveals a new mechanism of complement inhibition

Undesired complement activation is a major cause of tissue injury in various pathological conditions and contributes to several immune complex diseases. Compstatin, a 13-residue peptide, is an effective inhibitor of the activation of complement component C3 and thus blocks a central and crucial step in the complement cascade. The precise binding site on C3, the structure in the bound form, and the exact mode of action of compstatin are unknown. Here we present the crystal structure of compstatin in complex with C3c, a major proteolytic fragment of C3. The structure reveals that the compstatin-binding site is formed by the macroglobulin (MG) domains 4 and 5. This binding site is part of the structurally stable MG-ring formed by domains MG 1-6 and is far away from any other known binding site on C3. Compstatin does not alter the conformation of C3c, whereas compstatin itself undergoes a large conformational change upon binding. We propose a model in which compstatin sterically hinders the access of the substrate C3 to the convertase complexes, thus blocking complement activation and amplification. These insights are instrumental for further development of compstatin as a potential therapeutic.     

Role of the N-terminus of epidermal growth factor in ErbB-2/ErbB-3 binding studied by phage display

Epidermal growth factor (EGF) binds with high affinity to the EGF receptor, also known as ErbB-1, but upon replacement of the N-terminal linear region by neuregulin (NRG) 1 or transforming growth factor (TGF) alpha sequences it gains in addition high affinity for ErbB-2/ErbB-3 heterodimers. However, these chimeras weakly bind to ErbB-3 alone. To further dissect the ligand binding selectivity of the ErbB network, we have applied the phage display technique to examine the role of the linear N-terminal region in EGF for interaction with ErbB-2/ErbB-3 heterodimers. A library of EGF variants was constructed in which residues 2, 3, and 4 were randomly mutated, followed by selection for binding to intact MDA-MB-453 cells that overexpress ErbB-2 and ErbB-3 but lack ErbB-1. Analysis of the selected phage EGF variants revealed clones with high binding affinity to ErbB-2/ErbB-3 while maintaining high affinity to ErbB-1. In these variants, Trp (or alternatively His) was almost exclusively present at position 2, while specific combinations of hydrophobic, basic, and small residues were found at positions 3 and 4. The mitogenic activity of the phage EGF variants corresponded with their relative binding affinity. Two of the selected EGF variants, EGF/WVS and EGF/WRS, were further characterized as recombinant proteins. In contrast to previously characterized chimeras of EGF with NRG-1 or TGF-alpha, these variants did not only show high binding affinity for ErbB-2/ErbB-3 heterodimers but also for ErbB-3 alone. These data show that the linear N-terminal region of EGF-like growth factors is directly involved in binding to ErbB-3.     

Discovery of entry inhibitors for HIV-1 via a new de novo protein design framework

A new (to our knowledge) de novo design framework with a ranking metric based on approximate binding affinity calculations is introduced and applied to the discovery of what we believe are novel HIV-1 entry inhibitors. The framework consists of two stages: a sequence selection stage and a validation stage. The sequence selection stage produces a rank-ordered list of amino-acid sequences by solving an integer programming sequence selection model. The validation stage consists of fold specificity and approximate binding affinity calculations. The designed peptidic inhibitors are 12-amino-acids-long and target the hydrophobic core of gp41. A number of the best-predicted sequences were synthesized and their inhibition of HIV-1 was tested in cell culture. All peptides examined showed inhibitory activity when compared with no drug present, and the novel peptide sequences outperformed the native template sequence used for the design. The best sequence showed micromolar inhibition, which is a 3-15-fold improvement over the native sequence, depending on the donor. In addition, the best sequence equally inhibited wild-type and Enfuvirtide-resistant virus strains.     

Saturation mutagenesis of the Tn10-encoded tet operator O1. Identification of base-pairs involved in Tet repressor recognition

Saturation mutagenesis of Tn10-encoded tet operator O1 was performed by chemical synthesis of 30 sequence variants yielding all possible point mutations of an operator half side. Their effect on Tet repressor binding was scored by an in-vivo repressor titration system. Tet repressor affinities of selected operator mutants were further characterized in vitro by dissociation rate measurements. The O1 sequence spans 19 base-pairs. Out of these, all 18 palindromic base-pairs are involved in Tet repressor recognition. The central base-pair does not contribute to sequence-specific binding of Tet repressor. At position 1 a pyrimidine residue is sufficient for maximal affinity to the repressor. At positions 2, 3 and 4, each mutation reduces repressor binding at least tenfold. Mutations at positions 5, 6, 7, 8 and 9 result in less drastic reductions of Tet repressor binding. Differential effects of mutations at a given position are used to deduce the chemical functions contacted by Tet repressor. The T.A to A.T transversion at position 9 increases Tet repressor affinity slightly, while all other mutations decrease repressor binding. The increased affinity of the wild-type tet operator O2 compared to wild-type O1 results from the addition of two favorable transversions at positions +/- 9 and an unfavorable T.A to C.G transition at position -7. Deletion or palindromic doubling of the central base-pair of the O1 palindrome reveals that the wild-type spacing of both operator half sides is crucial for efficient Tet repressor binding.     

Structure‐kinetic relationship analysis of the therapeutic complement inhibitor compstatin

Compstatin is a 13‐residue peptide that inhibits activation of the complement system by binding to the central component C3 and its fragments C3b and C3c. A combination of theoretical and experimental approaches has previously allowed us to develop analogs of the original compstatin peptide with up to 264‐fold higher activity; one of these analogs is now in clinical trials for the treatment of age‐related macular degeneration (AMD). Here we used functional assays, surface plasmon resonance (SPR), and isothermal titration calorimetry (ITC) to assess the effect of modifications at three key residues (Trp‐4, Asp‐6, Ala‐9) on the affinity and activity of compstatin and its analogs, and we correlated our findings to the recently reported co‐crystal structure of compstatin and C3c. The K_D values for the panel of tested analogs ranged from 10^?6 to 10^?8 M. These differences in binding affinity could be attributed mainly to differences in dissociation rather than association rates, with a >4‐fold range in k_on values (2–10 × 10⁵ M^?1 s^?1) and a k_off variation of >35‐fold (1–37 × 10^?2 s^?1) being observed. The stability of the C3b‐compstatin complex seemed to be highly dependent on hydrophobic effects at position 4, and even small changes at position 6 resulted in a loss of complex formation. Induction of a β‐turn shift by an A9P modification resulted in a more favorable entropy but a loss of binding specificity and stability. The results obtained by the three methods utilized here were highly correlated with regard to the activity/affinity of the analogs. Thus, our analyses have identified essential structural features of compstatin and provided important information to support the development of analogs with improved efficacy. Copyright © 2009 John Wiley & Sons, Ltd.     

Studies of structure-activity relations of complement inhibitor compstatin

Compstatin, a 13-mer cyclic peptide, is a novel and promising inhibitor of the activation of the complement system. In our search for a more active analog and better understanding of structure-functions relations, we designed a phage-displayed random peptide library based on previous knowledge of structure activity relations, in which seven amino acids deemed necessary for structure and activity were kept fixed while the remaining six were optimized. Screening of this library against C3 identified four binding clones. Synthetic peptides corresponding to these clones revealed one analog, called acetylated Ile(1)Leu/His(9)Trp/Thr(13)Gly triple replacement analog of compstatin corresponding to clone 640 (Ac-I1L/H9W/T13G), which was more active than compstatin. This newly identified peptide had 4-fold higher activity when compared with the originally isolated form of compstatin and 1.6-fold higher activity when compared with acetylated compstatin (Ac-compstatin). The structures of Ac-I1L/H9W/T13G and Ac-compstatin were studied by nuclear magnetic resonance, compared with the structure of compstatin, and found to be very similar. The binding of Ac-I1L/H9W/T13G and the equally active acetylated analog with His(9)Ala replacement (Ac-H9A) to C3 was evaluated by surface plasmon resonance, which suggested similarity in their binding mechanism but difference when compared with Ac-compstatin. Compensatory effects of flexibility outside the beta-turn and tryptophan ring stacking may be responsible for the measured activity increase in Ac-I1L/H9W/T13G and acetylated analog with His(9)Ala replacement and the variability in binding mechanism compared with Ac-compstatin. These data demonstrate that tryptophan is a key amino acid for activity. Finally, the significance of the N-terminal acetylation was examined and it was found that the hydrophobic cluster at the linked termini of compstatin is essential for binding to C3 and for activity.     

Selection and design of high affinity DNA ligands for mutant single-chain derivatives of the bacteriophage 434 repressor

Single-chain repressor RRTRES is a derivative of bacteriophage 434 repressor, which contains covalently dimerized DNA-binding domains (amino acids 1-69) of the phage 434 repressor. In this single-chain molecule, the wild type domain R is connected to the mutant domain RTRES by a recombinant linker in a head-to-tail arrangement. The DNA-contacting amino acids of RTRES at the -1, 1, 2, and 5 positions of the a3 helix are T, R, E, S respectively. By using a randomized DNA pool containing the central sequence -CATACAAGAAAGNNNNNNTTT-, a cyclic, in vitro DNA-binding site selection was performed. The selected population was cloned and the individual members were characterized by determining their binding affinities to RRTRES. The results showed that the optimal operators contained the TTAC or TTCC sequences in the underlined positions as above, and that the Kd values were in the 1×10-12 mol/L-1×10-11mol/L concentration range. Since the affinity of the natural 434 repressor to its natural operator sites is in the 1×10-9 mol/L range, the observed binding affinity increase is remarkable. It was also found that binding affinity was strongly affected by the flanking bases of the optimal tetramer binding sites, especially by the base at the 5′ position. We constructed a new homodimeric single-chain repressor RTRESRTRES and its DNA-binding specificity was tested by using a series of new operators designed according to the recog-nition properties previously determined for the RTRES domain. These operators containing the con-sensus sequence GTAAGAAARNTTACN or GGAAGAAARNTTCCN (R is A or G) were recognized by RTRESRTRES specifically, and with high binding affinity. Thus, by using a combination of random selection and rational design principles, we have discovered novel, high affinity protein-DNA inter-actions with new specificity. This method can potentially be used to obtain new binding specificity for other DNA-binding proteins.     

Design of a modified mouse protein with ligand binding properties of its human analog by molecular dynamics simulations: the case of C3 inhibition by compstatin

The peptide compstatin and its derivatives inhibit the complement-component protein C3 in primate mammals and are potential therapeutic agents against the unregulated activation of complement in humans, but are inactive against C3 from lower mammals. Recent molecular dynamics (MD) simulations showed that the most potent compstatin analog comprised entirely of natural amino acids (W4A9) had a smaller affinity for rat C3, due to reproducible changes in the rat protein structure with respect to the human protein, which eliminated or weakened specific protein-ligand interactions seen in the human C3:W4A9 complex. Here, we study by MD simulations three W4A9 complexes with the mouse C3 protein, and two "transgenic" mouse derivatives, containing a small number (6-9) of human C3 substitutions. The mouse complex experiences the conformational changes and affinity reduction of the rat complex. In the "transgenic" complexes, the conformation remains closer to that of the human complex, the protein-ligand interactions are improved, and the affinity for compstatin becomes "human-like." The present work creates new avenues for a compstatin-sensitive animal model. A similar strategy, involving the comparison of a series of complexes by MD simulations, could be used to design "transgenic" sequences in other systems.     

VH-VL orientation prediction for antibody humanization candidate selection: A case study

Antibody humanization describes the procedure of grafting a non-human antibody's complementarity-determining regions, i.e., the variable loop regions that mediate specific interactions with the antigen, onto a β-sheet framework that is representative of the human variable region germline repertoire, thus reducing the number of potentially antigenic epitopes that might trigger an anti-antibody response. The selection criterion for the so-called acceptor frameworks (one for the heavy and one for the light chain variable region) is traditionally based on sequence similarity. Here, we propose a novel approach that selects acceptor frameworks such that the relative orientation of the 2 variable domains in 3D space, and thereby the geometry of the antigen-binding site, is conserved throughout the process of humanization. The methodology relies on a machine learning-based predictor of antibody variable domain orientation that has recently been shown to improve the quality of antibody homology models. Using data from 3 humanization campaigns, we demonstrate that preselecting humanization variants based on the predicted difference in variable domain orientation with regard to the original antibody leads to subsets of variants with a significant improvement in binding affinity.     

Systematic single base-pair substitution analysis of DNA binding by the cAMP receptor protein in cyanobacterium Synechocystis sp. PCC 6803

The cAMP receptor protein SYCRP1 in cyanobacterium Synechocystis sp. PCC 6803 is a regulatory protein that binds to the consensus DNA sequence (5'-AAATGTGATCTAGATCACATTT-3') for the cAMP receptor protein CRP in Escherichia coli. Here we examined the effects of systematic single base-pair substitutions at positions 4-8 (TGTGA) of the consensus sequence on the specific binding of SYCRP1. The consensus sequence exhibited the highest affinity, and the effects of base-pair substitutions at positions 5 and 7 were the most deleterious. The result is similar to that previously reported for CRP, whereas there were differences between SYCRP1 and CRP in the rank order of affinity for each substitution.     

The structural basis of compstatin activity examined by structure-function-based design of peptide analogs and NMR

We have previously identified compstatin, a 13-residue cyclic peptide, that inhibits complement activation by binding to C3 and preventing C3 cleavage to C3a and C3b. The structure of compstatin consists of a disulfide bridge and a type I beta-turn located at opposite sides to each other. The disulfide bridge is part of a hydrophobic cluster, and the beta-turn is part of a polar surface. We present the design of compstatin analogs in which we have introduced a series of perturbations in key structural elements of their parent peptide, compstatin. We have examined the consistency of the structures of the designed analogs compared with compstatin using NMR, and we have used the resulting structural information to make structure-complement inhibitory activity correlations. We propose the following. 1) Even in the absence of the disulfide bridge, a linear analog has a propensity for structure formation consistent with a turn of a 3(10)-helix or a beta-turn. 2) The type I beta-turn is a necessary but not a sufficient condition for activity. 3) Our substitutions outside the type I beta-turn of compstatin have altered the turn population but not the turn structure. 4) Flexibility of the beta-turn is essential for activity. 5) The type I beta-turn introduces reversibility and sufficiently separates the two sides of the peptide, whereas the disulfide bridge prevents the termini from drifting apart, thus aiding in the formation of the hydrophobic cluster. 6) The hydrophobic cluster at the linked termini is involved in binding to C3 and activity but alone is not sufficient for activity. 7) beta-Turn residues Gln(5) (Asn(5))-Asp(6)-Trp(7)(Phe(7))-Gly(8) are specific for the turn formation, but only Gln(5)(Asn(5))-Asp(6)-Trp(7)-Gly(8) residues are specific for activity. 8) Trp(7) is likely to be involved in direct interaction with C3, possibly through the formation of a hydrogen bond. Finally we propose a binding model for the C3-compstatin complex.     

Monitored whole gene in vitro evolution of an anti-hRaf-1 affibody molecule towards increased binding affinity

The use of library technologies for the generation of affinity proteins often includes an affinity maturation step, based on the construction of secondary libraries from which second generation variants with improved affinities are selected. Here, we describe for the first time the affinity maturation of affibody molecules based on step-wise in vitro molecular evolution, involving cycles of error-prone PCR (epPCR) amplification for the introduction of diversity over the entire 58-residue three-helix bundle structure and ribosome display (RD) for the selection of improved variants. The model affibody molecule for the process was Z(RAF322), binding with a 1.9μm equilibrium dissociation constant (K(D)) to human Raf-1 (hRaf-1), a protein kinase of central importance in the MAPK/ERK proliferation pathway. The molecular evolution process was followed on both gene and protein levels via DNA sequencing and a biosensor-based binding analysis of pools of selected variants. After two cycles of diversification and selection, a significant increase in binding response of selected pools was seen. DNA sequencing showed that a dominant alanine to valine substitution had been effectively enriched, and was found in 83% of all selected clones, either alone or in combination with other enriched substitutions. The evolution procedure resulted in variants showing up to 26-fold increases in affinity to the hRaf-1 target. Noteworthy, for the two variants showing the highest affinities, substitutions were also found in affibody framework positions, corresponding to regions of the protein domain not addressed by traditional affibody molecule affinity maturation strategies. Interestingly, thermal melting point (T(m)) analyses showed that an increased affinity could be associated with both higher and lower T(m) values. All investigated variants showed excellent refolding properties and selective binding to hRaf-1, as analysed using a multiplexed bead-based binding assay, making them potentially valuable affinity reagents for cell biology studies.     

Repertoire selection of variant single-chain Cro: toward directed DNA-binding specificity of helix-turn-helix proteins

A single-chain derivative of the lambda Cro repressor (scCro) has been randomly mutated in amino acid residues critical for specific DNA recognition to create libraries of protein variants. Utilizing phage display-afforded affinity selection, scCro variants have been isolated for binding to synthetic DNA ligands. Isolated scCro variants were analyzed functionally, both in fusion with phage particles and after expression of the corresponding free proteins. The binding properties with regard to specificity and affinity in binding to different DNA ligands were investigated by inhibition studies and determination of equilibrium dissociation constants for formed complexes. Variant proteins with altered DNA-sequence specificity were identified, which favored binding of targeted synthetic DNA sequences over a consensus operator sequence, bound with high affinity by wild-type Cro. The specificities were relatively modest (2-3-fold, as calculated from K(D) values), which can be attributed to the inherent properties in the design of the selection system; one half-site of the synthetic DNA sequences maintains the consensus operator sequence, and one "subunit" of the variant single-chain Cro dimers was conserved as wild-type sequence. The anticipated interaction between the wild-type subunit and the consensus DNA half-site of target DNA ligands is, hence, expected to contribute to the overlap in sequence discrimination. The binding affinity for the synthetic DNA sequences, however, was improved 10-30-fold in selected variant proteins as compared to "wild-type" scCro.     

Use of phage display to probe the evolution of binding specificity and affinity in integrins

The specific binding of RGD-containing proteins to integrin is a function of both the conformation of and the local sequence surrounding the RGD motif. To study the effect of these factors on integrin binding affinity and specificity, we obtained RGD-containing ligands specific for different integrins presented on the same protein scaffold. The beta-turn region between two anti-parallel beta-strands on the loop I of tendamistat, an inhibitor of alpha-amylase, was extended by two residues and randomized in a phagemid library. This library and two subsequently constructed RGD-containing loop I libraries were biopanned with purified integrins alphaIIbbeta3, alphaVbeta3 and alphaVbeta5 individually. The sequence analysis of selected tendamistat variants and characterization by phage ELISA revealed that phage adhesion is mediated exclusively by an RGD motif located at only two out of four possible positions on loop I. Further, sequences flanking the RGD motif were specific for different integrin targets. Interestingly, selected tendamistat variants mimic natural integrin ligands, both in sequence similarity and in integrin binding specificity, indicating that various ligand specificity patterns can be generated by driving towards maximum affinity in the integrin-ligand complexes.     

Use of a peptide mimotope to guide the humanization of MRK-16, an anti-P-glycoprotein monoclonal antibody.

A mimotope-guided strategy for engineering antibodies directed against orphan targets or antigens that are difficult to purify was developed and used to humanize the murine MRK-16 monoclonal antibody (mAb). MRK-16 recognizes a conformational epitope of a 170-kDa membrane protein, termed P-glycoprotein (P-gp). Elevated expression of P-gp on tumor cells is associated with resistance to cytotoxic drugs, a major obstacle in chemotherapy. Murine MRK-16 was used to enrich and screen a phage-displayed peptide library to identify reactive mimotopes. One peptide, termed ALR1, was enriched to a greater extent than others and subsequently was expressed as a fusion protein with glutathione S-transferase. ALR1 fusion protein bound MRK-16 specifically and inhibited binding of MRK-16 to cells expressing elevated levels of P-gp. To humanize MRK-16, the murine complementarity determining regions were grafted onto homologous human heavy and light chain variable region frameworks. Framework residues that differed between the murine MRK-16 and the homologous human templates were analyzed and subsequently, five framework positions potentially important for maintaining the specificity and affinity of MRK-16 were identified. A combinatorial library consisting of 32 variants encoding all possible combinations of murine and human residues at the five differing framework positions was expressed in a phage system. In the absence of purified P-gp, ALR1 fusion protein was used as surrogate antigen to screen the antibody library to identify the framework combination that most preserved the binding activity of the mAb. On the basis of the initial screening against the mimotope four antibody variants were selected for further characterization. The binding affinity of these variants for the ALR1 fusion protein correlated with their binding to cells expressing elevated levels of P-gp. Thus, peptide mimotopes which can be identified for virtually any antibody including those that recognize conformational or carbohydrate epitopes, can serve as antigen templates for antibody engineering.     

Structural insights into humanization of anti-tissue factor antibody 10H10

Murine antibody 10H10 raised against human tissue factor is unique in that it blocks the signaling pathway, and thus inhibits angiogenesis and tumor growth without interfering with coagulation. As a potential therapeutic, the antibody was humanized in a two-step procedure. Antigen-binding loops were grafted onto selected human frameworks and the resulting chimeric antibody was subjected to affinity maturation by using phage display libraries. The results of humanization were analyzed from the structural perspective through comparison of the structure of a humanized variant with the parental mouse antibody. This analysis revealed several hot spots in the framework region that appear to affect antigen binding, and therefore should be considered in human germline selection. In addition, some positions in the Vernier zone, e.g., residue 71 in the heavy chain, that are traditionally thought to be crucial appear to tolerate amino acid substitutions without any effect on binding. Several humanized variants were produced using both short and long forms of complementarity-determining region (CDR) H2 following the difference in the Kabat and Martin definitions. Comparison of such pairs indicated consistently higher thermostability of the variants with short CDR H2. Analysis of the binding data in relation to the structures singled out the ImMunoGeneTics information system® germline IGHV1-2*01 as dubious owing to two potentially destabilizing mutations as compared to the other alleles of the same germline and to other human germlines.     

Identifying Primate ACE2 Variants That Confer Resistance to SARS-CoV-2

SARS-CoV-2 infects humans through the binding of viral S-protein (spike protein) to human angiotensin I converting enzyme 2 (ACE2). The structure of the ACE2-S-protein complex has been deciphered and we focused on the 27 ACE2 residues that bind to S-protein. From human sequence databases, we identified nine ACE2 variants at ACE2–S-protein binding sites. We used both experimental assays and protein structure analysis to evaluate the effect of each variant on the binding affinity of ACE2 to S-protein. We found one variant causing complete binding disruption, two and three variants, respectively, strongly and mildly reducing the binding affinity, and two variants strongly enhancing the binding affinity. We then collected the ACE2 gene sequences from 57 nonhuman primates. Among the 6 apes and 20 Old World monkeys (OWMs) studied, we found no new variants. In contrast, all 11 New World monkeys (NWMs) studied share four variants each causing a strong reduction in binding affinity, the Philippine tarsier also possesses three such variants, and 18 of the 19 prosimian species studied share one variant causing a strong reduction in binding affinity. Moreover, one OWM and three prosimian variants increased binding affinity by >50%. Based on these findings, we proposed that the common ancestor of primates was strongly resistant to and that of NWMs was completely resistant to SARS-CoV-2 and so is the Philippine tarsier, whereas apes and OWMs, like most humans, are susceptible. This study increases our understanding of the differences in susceptibility to SARS-CoV-2 infection among primates.     

The heparin-binding domain of extracellular superoxide dismutase C and formation of variants with reduced heparin affinity.

A fundamental property of the secretory tetrameric extracellular superoxide dismutase (EC-SOD) is its affinity for heparin and analogues, in vivo, mediating attachment to heparan sulfate proteoglycans located on cell surfaces and in the connective tissue matrix. EC-SOD is in vivo heterogeneous with regard to heparin affinity and can be divided into subclasses; A which lacks heparin affinity, B with intermediate affinity, and C with strong heparin affinity. The EC-SOD C subunits contain 222 amino acids and among the last 20 carboxyl-terminal amino acids, 10 are positively charged and six of these are located in a cluster in positions 210-215. To analyze if this local accumulation of basic amino acids is responsible for heparin binding we produced three series of recombinant EC-SOD (rEC-SOD) variants, six containing amino acid exchanges in the carboxyl-terminal end, four with truncations, and two with both truncations and substitutions. Exchange of positively or negatively charged amino acids on the carboxyl-terminal side of the cluster results in only minor modifications in heparin affinity, whereas substitution of three of the amino acids in the cluster abrogates the heparin binding. Insertions of stop codons at different positions resulted in either C or A but not B class EC-SOD. In an attempt to produce EC-SODs with intermediate heparin affinities, plasmids defining C and A class EC-SOD were cotransfected into Chinese hamster ovary cells. In addition to the parental A and C class EC-SOD forms, two variants with intermediate heparin affinities were formed. Coincubation of EC-SOD C and A resulted in the appearance of one heterotetramer with intermediate affinity for heparin. We conclude that the cluster of six basic amino acids forms the essential part of the heparin-binding domain and that the composition of the four subunits in the EC-SOD tetramer determines the affinity for heparin. This domain is different from heparin-binding domains of other proteins, and its localization allows the distribution of EC-SOD in vivo to be regulated by proteolytic processing.