Mapping the mutual information network of enzymatic families in the protein structure to unveil functional features

Amino acids committed to a particular function correlate tightly along evolution and tend to form clusters in the 3D structure of the protein. Consequently, a protein can be seen as a network of co-evolving clusters of residues. The goal of this work is two-fold: first, we have combined mutual information and structural data to describe the amino acid networks within a protein and their interactions. Second, we have investigated how this information can be used to improve methods of prediction of functional residues by reducing the search space. As a main result, we found that clusters of co-evolving residues related to the catalytic site of an enzyme have distinguishable topological properties in the network. We also observed that these clusters usually evolve independently, which could be related to a fail-safe mechanism. Finally, we discovered a significant enrichment of functional residues (e.g. metal binding, susceptibility to detrimental mutations) in the clusters, which could be the foundation of new prediction tools.     

Molecular dynamics of mesophilic-like mutants of a cold-adapted enzyme: insights into distal effects induced by the mutations

Networks and clusters of intramolecular interactions, as well as their "communication" across the three-dimensional architecture have a prominent role in determining protein stability and function. Special attention has been dedicated to their role in thermal adaptation. In the present contribution, seven previously experimentally characterized mutants of a cold-adapted α-amylase, featuring mesophilic-like behavior, have been investigated by multiple molecular dynamics simulations, essential dynamics and analyses of correlated motions and electrostatic interactions. Our data elucidate the molecular mechanisms underlying the ability of single and multiple mutations to globally modulate dynamic properties of the cold-adapted α-amylase, including both local and complex unpredictable distal effects. Our investigation also shows, in agreement with the experimental data, that the conversion of the cold-adapted enzyme in a warm-adapted variant cannot be completely achieved by the introduction of few mutations, also providing the rationale behind these effects. Moreover, pivotal residues, which are likely to mediate the effects induced by the mutations, have been identified from our analyses, as well as a group of suitable candidates for protein engineering. In fact, a subset of residues here identified (as an isoleucine, or networks of mesophilic-like salt bridges in the proximity of the catalytic site) should be considered, in experimental studies, to get a more efficient modification of the features of the cold-adapted enzyme.     

Conservative mutation Met8 --> Leu affects the folding process and structural stability of squash trypsin inhibitor CMTI-I

Protein molecules can accommodate a large number of mutations without noticeable effects on their stability and folding kinetics. On the other hand, some mutations can have quite strong effects on protein conformational properties. Such mutations either destabilize secondary structures, e.g., alpha-helices, are incompatible with close packing of protein hydrophobic cores, or lead to disruption of some specific interactions such as disulfide cross links, salt bridges, hydrogen bonds, or aromatic-aromatic contacts. The Met8 --> Leu mutation in CMTI-I results in significant destabilization of the protein structure. This effect could hardly be expected since the mutation is highly conservative, and the side chain of residue 8 is situated on the protein surface. We show that the protein destabilization is caused by rearrangement of a hydrophobic cluster formed by side chains of residues 8, Ile6, and Leu17 that leads to partial breaking of a hydrogen bond formed by the amide group of Leu17 with water and to a reduction of a hydrophobic surface buried within the cluster. The mutation perturbs also the protein folding. In aerobic conditions the reduced wild-type protein folds effectively into its native structure, whereas more then 75% of the mutant molecules are trapped in various misfolded species. The main conclusion of this work is that conservative mutations of hydrophobic residues can destabilize a protein structure even if these residues are situated on the protein surface and partially accessible to water. Structural rearrangement of small hydrophobic clusters formed by such residues can lead to local changes in protein hydration, and consequently, can affect considerably protein stability and folding process.     

Active-site residues move independently from the rest of the protein in a 200 ns molecular dynamics simulation of cytochrome P450 CYP119

The conformational dynamics of cytochrome P450 enzymes are critical to their catalytic activity. In this study, the correlated motion between residues in a 200 ns molecular dynamics trajectory of the thermophilic CYP119 was analyzed to parse out conformational relationships. Residues that are structurally related, for example residues within a helix, generally have highly correlated motion. In addition, clusters of non-adjacent residues that show correlated motion (“hot spots”) are seen in various regions, including at the base of the F and G helices that make up the most dynamic region of the enzyme. A modified k-means algorithm that clusters residues based on their correlated motion indicates that functionally related residues are in the same cluster (e.g., the catalytic threonines and the heme). Tightly coupled clusters form a solvent-exposed “shell” around the enzyme, whereas less coupling between clusters is seen in regions that are critical to ligand interactions, redox partner interactions, and catalysis. Most notably, we find that residues in the active site move independently from the rest of the enzyme, effectively insulating the catalytic machinery from other regions of the protein.     

Molecular analysis of type I-A (tyrosinase negative) oculocutaneous albinism

Type I oculocutaneous albinism (OCA) is caused by the reduction in or absence of activity of tyrosinase in melanocytes in skin, hair, and the eyes, the result of mutations of the tyrosinase gene. To date, a total of 22 unique mutations in the coding region of tyrosinase have been described in the literature. In this report we present 5 additional mutations of the tyrosinase gene associated with type I-A OCA in four individuals, including 2 missense, 1 frameshift and 2 nonsense mutations, and review the relevant literature on all published mutations. Analysis of the distribution of all identified missense mutations (n = 17) shows that most cluster in three areas of the gene and involve amino acids conserved between humans and the mouse. Two clusters involve the copper A and copper B binding sites and may disrupt the metal ion-protein interaction necessary for enzyme function. The third cluster in exon I could represent a functional domain important in enzyme function such as the tyrosine or the dihydroxyphenylalanine (DOPA) binding site of the enzyme. Small deletions or insertions resulting in frameshift mutations and nonsense mutations are distributed throughout the coding region and do not appear to cluster.     

The key role of a non-active-site residue Met148 on the catalytic efficiency of meta-cleavage product hydrolase BphD

meta-Cleavage product (MCP) hydrolases (EC 3.7.1.9) can catalyze a specific C–C bond fission during the microbial aerobic degradation of aromatics. The previous studies on structure–function relationship of MCP hydrolases mainly focus on the active site residues by site-directed mutagenesis. However, the information about the role of the non-active-site residues is still unclear. In this study, a non-active-site residue Met148 of MCP hydrolase BphD was selected as the mutagenesis site according to the sequence alignments, structure superimpose and the tunnel analysis, which underwent the saturation mutagenesis resulting 19 mutants. The catalytic efficiencies of the mutants on 6-oxo-6-phenylhexa-2,4-dienoic acid (HOPDA) were all decreased compared with the wild-type one except for the M148D mutant. Especially, the M148P mutant exhibited 290-fold lower k _cat/K _m than that of the wild-type BphD. Transient kinetic analyses of M148P showed the reciprocal relaxation time corresponded to C–C bond cleavage and product release steps (9.6 s^?1) was 4.08-fold lower than BphD WT (39.2 s^?1). Tunnel cluster analysis of BphD WT, M148P and M148W demonstrated that only the bulky Trp148 could block tunnel T2 in the BphD WT, but it exhibited slight effects on the catalytic efficiency (0.94-fold of BphD WT). Therefore, product release was not the main reason for the efficiency decrease of M148P. On the other hand, molecular dynamics simulations on the BphD WT and BphD M148P in complex with HOPDA indicated that the dramatic decrease of the catalytic efficiencies of BphD M148P should be due to the unproductive binding of HOPDA. The study demonstrated the catalytic efficiency of MCP hydrolase can be engineered by modification of non-active site residue.     

How conformational changes can affect catalysis,inhibition and drug resistance of enzymes with induced-fit binding mechanism such as the HIV-1 protease

A central question is how the conformational changes of proteins affect their function and the inhibition of this function by drug molecules. Many enzymes change from an open to a closed conformation upon binding of substrate or inhibitor molecules. These conformational changes have been suggested to follow an induced-fit mechanism in which the molecules first bind in the open conformation in those cases where binding in the closed conformation appears to be sterically obstructed such as for the HIV-1 protease. In this article, we present a general model for the catalysis and inhibition of enzymes with induced-fit binding mechanism. We derive general expressions that specify how the overall catalytic rate of the enzymes depends on the rates for binding, for the conformational changes, and for the chemical reaction. Based on these expressions, we analyze the effect of mutations that mainly shift the conformational equilibrium on catalysis and inhibition. If the overall catalytic rate is limited by product unbinding, we find that mutations that destabilize the closed conformation relative to the open conformation increase the catalytic rate in the presence of inhibitors by a factor exp(ΔΔG_C/RT) where ΔΔG_C is the mutation-induced shift of the free-energy difference between the conformations. This increase in the catalytic rate due to changes in the conformational equilibrium is independent of the inhibitor molecule and, thus, may help to understand how non-active-site mutations can contribute to the multi-drug-resistance that has been observed for the HIV-1 protease. A comparison to experimental data for the non-active-site mutation L90M of the HIV-1 protease indicates that the mutation slightly destabilizes the closed conformation of the enzyme. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.     

Structural mapping of Kelch13 mutations associated with artemisinin resistance in malaria

Mutations in Plasmodium falciparum gene kelch13 (pfkelch13) are strongly and causally associated with resistance to anti-malarial drug artemisinin, but their effects on PfKelch13 structure and function remain unclear. Utilizing the publicly available three-dimensional structure of PfKech13 (PDB ID: 4yy8), we find that most of the mutations in its propeller domain occur in two spatial clusters. Of these, one cluster is enriched in surface exposed residues which may drive PfKelch13-centered protein interactions, and the second cluster mostly contains residues which are buried and whose mutations may destabilize PfKelch13 structure. The most prevalent resistant mutations C580Y and Y493H are distal from the above two clusters. The C580Y mutation creates sterically unfavourable contacts while Y493H possibly alters the hydrophobic core of the propeller domain. These analyses will facilitate further experimental studies aimed at understanding how mutations in pfkelch13 lead to artemisinin resistance.     

Alteration of hydrogen bonding in the vicinity of histidine 48 disrupts millisecond motions in RNase A

The motion of amino acid residues on the millisecond (ms) time scale is involved in the tight regulation of catalytic function in numerous enzyme systems. Using a combination of mutational, enzymological, and relaxation-compensated (15)N Carr-Purcell-Meiboom-Gill (CPMG) methods, we have previously established the conformational significance of the distant His48 residue and the neighboring loop 1 in RNase A function. These studies suggested that RNase A relies on an intricate network of hydrogen bonding interactions involved in propagating functionally relevant, long-range ms motions to the catalytic site of the enzyme. To further investigate the dynamic importance of this H-bonding network, this study focuses on the individual replacement of Thr17 and Thr82 with alanine, effectively altering the key H-bonding interactions that connect loop 1 and His48 to the rest of the protein. (15)N CPMG dispersion studies, nuclear magnetic resonance (NMR) chemical shift analysis, and NMR line shape analysis of point mutants T17A and T82A demonstrate that the evolutionarily conserved single H-bond linking His48 to Thr82 is essential for propagating ms motions from His48 to the active site of RNase A on the time scale of catalytic turnover, whereas the T17A mutation increases the off rate and conformational exchange motions in loop 1. Accumulating evidence from our mutational studies indicates that residues experiencing conformational exchange in RNase A can be grouped into two separate clusters displaying distinct dynamical features, which appear to be independently affected by mutation. Overall, this study illuminates how tightly controlled and finely tuned ms motions are in RNase A, suggesting that designed modulation of protein motions may be possible.     

The barley stripe mosaic virus gamma b gene encodes a multifunctional cysteine-rich protein that affects pathogenesis.

Barley stripe mosaic virus contains seven genes, one of which specifies a 17-kD cysteine-rich protein, gamma b, that is known to affect virulence. To further characterize the role of gamma b in pathogenesis, we mutagenized sequences encoding amino acids within two clusters of cysteine and histidine residues in the cysteine-rich domain and a group of basic amino acids located between the clusters and determined the effects of these mutations on the symptom phenotype in barley. Three single amino acid substitutions in cluster 1 and two amino acid exchanges in the basic region caused bleached symptoms associated with pronounced elevations in accumulation of gamma b protein. In contrast, three single amino acid substitutions in cluster 2 and a mutation in the basic motif resulted in attenuated ("null") symptoms typical of those produced when the gamma b gene is deleted. Tissue infected with these "null" mutants accumulated slightly elevated amounts of the gamma b protein but significantly lower levels of coat protein and the putative movement protein beta b. Genetic complementation tests revealed that cluster 1 mutations are dominant over the wild-type gamma b gene, whereas those in cluster 2 are recessive. These results highlight the pivotal role of gamma b in pathogenesis and suggest that the two cysteine-rich clusters are functionally distinct and that they affect different aspects of disease development.     

Structure and function of the mitochondrial bc1 complex. A mutational analysis of the yeast Rieske iron-sulfur protein

Respiratory-defective mutants of Saccharomyces cerevisiae assigned to a single complementation group (G12) have been determined to have lesions in the iron-sulfur protein (Rieske protein) of ubiquinol: cytochrome c reductase. Mutants capable of expressing the protein were chosen for further studies. The genes from 13 independent isolates were cloned and their mutations sequenced. Twelve mutations were ascertained to cause single amino acid substitutions in the carboxyl-terminal regions of the protein between residues 127 and 173. This region is proposed to be part of the catalytic domain with the ligands responsible for co-ordinating the two irons of the 2Fe-2S cluster. Based on the catalytic properties of the ubiquinol: cytochrome c reductase complex and the electron paramagnetic resonance (e.p.r.) signals of the iron-sulfur protein, the mutants describe two different phenotypes. A subset of mutants have no detectable iron-sulfur cluster and are completely deficient in ubiquinol: cytochrome c reductase activity. These strains identify mutations in residues considered to be essential for binding of the iron or for maintaining a proper tertiary structure of the catalytic domain. A second group of mutants have reduced levels of enzymatic activity and exhibit e.p.r. spectra characteristic of the Rieske iron-sulfur cluster. The mutations in the latter strains have been ascribed to residues that influence the redox properties of the cluster by distorting the iron-binding pocket. A secondary and tertiary structure model is presented of the carboxyl-terminal 65 residues constituting the catalytic domain of the iron-sulfur protein. It is postulated that the two irons of the cluster are co-ordinated by three cysteine and a single histidine residue located in a loop structure. The catalytic domain also contains two short alpha-helices and three beta-strands that form a partial beta-barrel. Most of the hydrophilic amino acids are present in turns that map to one pole of the domain. When viewed in the context of the model, mutations that abolish the iron-sulfur cluster are mostly in residues defining the boundaries of the alpha-helices and beta-strands. The notable exception is a cysteine residue that has been assigned to the loop with the iron ligands. This cysteine residue is proposed to co-ordinate one iron of the cluster. Mutations that reduce ubiquinol: cytochrome c reductase activity and alter the redox potential of the cluster occur in residues located in the loop that contains the ligands of the cluster.     

The bacterial ribonuclease P holoenzyme requires specific, conserved residues for efficient catalysis and substrate positioning

RNase P is an RNA-based enzyme primarily responsible for 5′-end pre-tRNA processing. A structure of the bacterial RNase P holoenzyme in complex with tRNA^Phe revealed the structural basis for substrate recognition, identified the active site location, and showed how the protein component increases functionality. The active site includes at least two metal ions, a universal uridine (U52), and P RNA backbone moieties, but it is unclear whether an adjacent, bacterially conserved protein loop (residues 52–57) participates in catalysis. Here, mutagenesis combined with single-turnover reaction kinetics demonstrate that point mutations in this loop have either no or modest effects on catalytic efficiency. Similarly, amino acid changes in the ‘RNR’ region, which represent the most conserved region of bacterial RNase P proteins, exhibit negligible changes in catalytic efficiency. However, U52 and two bacterially conserved protein residues (F17 and R89) are essential for efficient Thermotoga maritima RNase P activity. The U52 nucleotide binds a metal ion at the active site, whereas F17 and R89 are positioned >20 Å from the cleavage site, probably making contacts with N₋₄ and N₋₅ nucleotides of the pre-tRNA 5′-leader. This suggests a synergistic coupling between transition state formation and substrate positioning via interactions with the leader.     

Stepwise improvements in catalytic effectiveness: independence and interdependence in combinations of point mutations of a sluggish triosephosphate isomerase

Second-site suppressor changes that improve the catalytic potency of a sluggish mutant of the enzyme triosephosphate isomerase have been examined both individually and in combination. Each of the second-site mutations increases the specific catalytic activity of a triosephosphate isomerase in which the catalytic base, glutamate-165, has been changed to aspartate. These second-site suppressors are G10S, S96P, S96T, E97D, V167D, and G233R. Not one of these changes enhances the value of kcat/Km for the wild-type enzyme, which is consistent with the knowledge that the reaction catalyzed by the wild-type enzyme is already diffusion-controlled. Indeed, two of the changes, S96P and V167D, are catalytically deleterious to the wild-type isomerase. When pairs of second-site suppressors are combined with the primary lesion E165D, six pairs show additive independence while the effects of eight other pairs are less than additive. The sites fall into two clusters: pairs within a cluster always interfere with one another and do not produce additive improvements in catalytic activity, whereas combinations of changes from different clusters tend to be additive in their effects. No combination of second-site suppressor mutations behaves synergistically, though there seems to be no a priori reason to exclude this possibility. Since the catalytic potency of each of the six second-site suppressor mutants can be further improved by the introduction of (at least) one of the other five changes, it is evident that none of the double mutants lies at a local catalytic maximum. In these cases, therefore, the opportunity exists for at least two "steps" of monotonic catalytic improvement along each of six different "paths" in protein space.     

Aromatic clusters: a determinant of thermal stability of thermophilic proteins

A number of factors have been elucidated as responsible for the thermal stability of thermophilic proteins. However, the contribution of aromatic interactions to thermal stability has not been systematically studied. In the present investigation we used a graph spectral method to identify aromatic clusters in a dataset of 24 protein families for which the crystal structures of both the thermophilic and their mesophilic homologues are known. Our analysis shows a presence of additional aromatic clusters or enlarged aromatic networks in 17 different thermophilic protein families, which are absent in the corresponding mesophilic homologue. The additional aromatic clusters identified in the thermophiles are smaller in size and are largely found on the protein surface. The aromatic clusters are found to be relatively rigid regions of the surface and often the additional aromatic cluster is located close to the active site of the thermophilic enzyme. The residues in the additional aromatic clusters are preferably mutated to Leu, Ser or Ile in the mesophilic homologue. An analysis of the packing geometry of the pairwise aromatic interaction in the additional aromatic clusters shows a preference for a T-shaped orthogonal packing geometry. The present study also provides new insights for protein engineers to design thermostable and thermophilic proteins.     

Long-range interactions in the dimer interface of ornithine decarboxylase are important for enzyme function.

Ornithine decarboxylase (ODC) is a pyridoxal 5'-phosphate dependent enzyme that catalyzes the first committed step in the biosynthesis of polyamines. ODC is a proven drug target for the treatment of African sleeping sickness. The enzyme is an obligate homodimer, and the two identical active sites are formed at the dimer interface. Alanine scanning mutagenesis of dimer interface residues in Trypanosoma brucei ODC was undertaken to determine the energetic contribution of these residues to subunit association. Twenty-three mutant enzymes were analyzed by analytical ultracentrifugation, and none of the mutations were found to cause a greater than 1 kcal/mol decrease in dimer stability. These data suggest that the energetics of the interaction may be distributed across the interface. Most significantly, many of the mutations had large effects (DeltaDeltaG kcat/Km > 2.5 kcal/mol) on the catalytic efficiency of the enzyme. Residues that affected activity included those in or near the substrate binding site but also a number of residues that are distant (15-20 A) from this site. These data provide evidence that long-range energetic coupling of interface residues to the active site is essential for enzyme function, even though structural changes upon ligand binding to wild-type ODC are limited to local conformational changes in the active site. The ODC dimer interface appears to be optimized for catalytic function and not for dimer stability. Thus, small molecules directed to the ODC interfaces could impact biological function without having to overcome the difficult energetic barrier of dissociating the interacting partners.     

Analysis of a nuclear localization signal of simian virus 40 major capsid protein Vp1.

The nuclear localization signal of the major structural protein, Vp1, of simian virus 40 was further defined by mutagenesis. The targeting activity was examined in cells microinjected with SV-Vp1 variant viral DNAs bearing either an initiation codon mutation of the agnoprotein or mutations in the Vp1 coding sequence or microinjected with pSG5-Vp1 and pSG5-Vp1 mutant DNAs in which Vp1 or mutant Vp1 is expressed from simian virus 40 early promoter. The Vp1 nuclear localization signal functioned autonomously without agno-protein once the Vp1 protein was synthesized in the cytoplasm. The targeting activity was localized to the amino-terminal 19 residues. While replacement of cysteine 10 with glycine, alanine, or serine did not affect the activity, replacement of arginine 6 with glycine caused the cytoplasmic phenotype. When multiple mutations were introduced among residue 5, 6, 7, 16, 17, or 19, the targeting activity was found to reside in two clusters of basic residues, a cluster of lysine 5, arginine 6, and lysine 7 and a cluster of lysine 16, lysine 17, and lysine 19. The clusters are independently important for nuclear localization activity.     

Aromatic cluster mutations produce focal modulations of β-sheet structure

Site-directed mutagenesis is a powerful tool for altering the structure and function of proteins in a focused manner. Here, we examined how a model β-sheet protein could be tuned by mutation of numerous surface-exposed residues to aromatic amino acids. We designed these aromatic side chain “clusters” at highly solvent-exposed positions in the flat, single-layer β-sheet of Borrelia outer surface protein A (OspA). This unusual β-sheet scaffold allows us to interrogate the effects of these mutations in the context of well-defined structure but in the absence of the strong scaffolding effects of globular protein architecture. We anticipated that the introduction of a cluster of aromatic amino acid residues on the β-sheet surface would result in large conformational changes and/or stabilization and thereby provide new means of controlling the properties of β-sheets. Surprisingly, X-ray crystal structures revealed that the introduction of aromatic clusters produced only subtle conformational changes in the OspA β-sheet. Additionally, despite burying a large degree of hydrophobic surface area, the aromatic cluster mutants were slightly less stable than the wild-type scaffold. These results thereby demonstrate that the introduction of aromatic cluster mutations can serve as a means for subtly modulating β-sheet conformation in protein design.     

Amino Acid Residues Involved in the Functional Integrity of Escherichia coli Methionine Aminopeptidase

Amino acid residues in the metal-binding and putative substrate-binding sites of Escherichia coli methionine aminopeptidase (MAP) were mutated, and their effects on the function of the enzyme were investigated. Substitution of any amino acid residue at the metal-binding site resulted in complete loss of the two cobalt ions bound to the protein and diminished the enzyme activity. However, only Cys70 and Trp221 at the putative substrate-binding site are involved in the catalytic activity of MAP. Changing either of them caused partial loss of enzyme activity, while mutations at both positions abolished MAP function. Both residues are found to be conserved in type I but not type II MAPs.     

The transmembrane domain 6 of vacuolar H(+)-pyrophosphatase mediates protein targeting and proton transport

Vacuolar H(+)-pyrophosphatase (V-PPase; EC 3.6.1.1) plays a significant role in the maintenance of the pH in cytoplasm and vacuoles via proton translocation from the cytosol to the vacuolar lumen at the expense of PP(i) hydrolysis. The topology of V-PPase as predicted by TopPred II suggests that the catalytic site is putatively located in loop e and exposed to the cytosol. The adjacent transmembrane domain 6 (TM6) is highly conserved and believed to participate in the catalytic function and conformational stability of V-PPase. In this study, alanine-scanning mutagenesis along TM6 of the mung bean V-PPase was carried out to identify its structural and functional role. Mutants Y299A, A306S and L317A exhibited gross impairment in both PP(i) hydrolysis and proton translocation. Meanwhile, mutations at L307 and N318 completely abolished the targeting of the enzyme, causing broad cytosolic localization and implicating a possible role of these residues in protein translocation. The location of these amino acid residues was on the same side of the helix wheel, suggesting their involvement in maintaining the stability of enzyme conformation. G297A, E301A and A305S mutants showed declines in proton translocation but not in PP(i) hydrolysis, consequently resulting in decreases in the coupling efficiency. These amino acid residues cluster at one face of the helix wheel, indicating their direct/indirect participation in proton translocation. Taken together, these data indicate that TM6 is crucial to vacuolar H(+)-pyrophosphatase, probably mediating protein targeting, proton transport, and the maintenance of enzyme structure.     

Role of the electrostatic loop of Cu,Zn superoxide dismutase in the copper uptake process.

Cu,Zn superoxide dismutases are characterized by the presence of four highly conserved charged residues (Lys120, Glu/Asp130, Glu131 and Lys134), which are placed at the edge of the active site channel and have been shown to be individually involved in the electrostatic attraction of the substrate toward the catalytically active copper ion. By genetic engineering we mutated these four residues into neutrally charged ones (Leu120, Gln130, Gln131, Thr134). The effects of these mutations on the rate of superoxide dismutation were not dramatic. In fact, at two different pH and ionic strength values, the mutant enzyme had a catalytic constant even higher with respect to the wild-type protein, showing that electrostatic interaction at these surface sites is not essential for high catalytic efficiency of the enzyme. The mutant and the wild-type enzyme showed the same degree of inhibition by CN(-), and both were not affected by I(-), showing that mutations did not alter the sensitivity of the enzyme to anions. On the other hand, reconstitution of active enzyme from either the wild-type or mutant copper-free enzymes with a copper(I)-glutathione [Cu(I)-GSH] complex showed that metal uptake by the mutant was much slower than by the wild-type enzyme. The demonstration that the 'electrostatic loop' is apparently conserved to assure optimal copper uptake by the enzyme, rather than fast dismutation, may provide further support to the idea that Cu,Zn superoxide dismutase is a bifunctional protein, acting in cellular defense against oxidative stress both as a copper buffer and as a superoxide radical scavenger.