Aspartic acid-66 is the only essential negatively charged residue in the putative hydrophilic loop region of the metal-tetracycline/H+ antiporter encoded by transposon Tn10 of Escherichia coli.

Of the 16 acidic amino acid residues located in the hydrophilic region of the metal-tetracycline/H+ antiporter of transposon Tn10, five glutamic acids and three aspartic acids are conserved among the tetracycline/H+ antiporters of Gram-negative bacteria. When these conserved acidic residues were each replaced by a neutral polar residue, glutamine or asparagine, only the Asp66 substitution mutants completely lost their transport activity. The substitution of Glu274, Asp120, Glu181, or Asp38 caused significant reduction of the transport activity, whereas the substitution of the other three residues had no detectable effect on the activity. These findings led to the conclusion that only Asp66 is essential for the transport function.     

The active site of the junction-resolving enzyme T7 endonuclease I

Endonuclease I is a junction-resolving enzyme encoded by bacteriophage T7, that selectively binds and cleaves four-way DNA junctions. We have recently solved the structure of this dimeric enzyme at atomic resolution, and identified the probable catalytic residues. The putative active site comprises the side-chains of three acidic amino acids (Glu20, Asp55 and Glu65) together with a lysine residue (Lys67), and shares strong similarities with a number of type II restriction enzymes. However, it differs from a typical restriction enzyme as the proposed catalytic residues in both active sites are contributed by both polypeptides of the dimer. Mutagenesis experiments confirm the importance of all the proposed active site residues. We have carried out in vitro complementation experiments using heterodimers formed from mutants in different active site residues, showing that Glu20 is located on a different monomer from the remaining amino acid residues comprising the active site. These experiments confirm that the helix-exchanged architecture of the enzyme creates a mixed active site in solution. Such a composite active site structure should result in unilateral cleavage by the complemented heterodimer; this has been confirmed by the use of a cruciform substrate. Based upon analogy with closely similar restriction enzyme active sites and our mutagenesis experiments, we propose a two-metal ion mechanism for the hydrolytic cleavage of DNA junctions.     

Acidic residues critical for the activity and biological function of yeast DNA polymerase eta

Rad30 is a member of the newly discovered UmuC/DinB/Rad30 family of DNA polymerases. The N-terminal regions of these proteins are highly homologous, and they contain five conserved motifs, I to V, while their C-terminal regions are quite divergent. We examined the contributions of the C-terminal and N-terminal regions of Rad30 to its activity and biological function. Although deletion of the last 54 amino acids has no effect on DNA polymerase or thymine-thymine (T-T) dimer bypass activity, this C-terminal deletion-containing protein is unable to perform its biological function in vivo. The presence of a bipartite nuclear targeting sequence within this region suggests that at least one function of this portion of Rad30 is nuclear targeting. To identify the active-site residues of Rad30 important for catalysis, we generated mutations of nine acidic residues that are invariant or highly conserved among Rad30 proteins from different eukaryotic species. Mutations of the Asp30 and Glu39 residues present in motif I and of the Asp155 residue present in motif III to alanine completely inactivated the DNA polymerase and T-T dimer bypass activities, and these mutations did not complement the UV sensitivity of the rad30Delta mutation. Mutation of Glu156 in motif III to alanine confers a large reduction in the efficiency of nucleotide incorporation, whereas the remaining five Rad30 mutant proteins retain wild-type levels of DNA polymerase and T-T dimer bypass activities. From these observations, we suggest a role for the Asp30, Glu39, and Asp155 residues in the binding of two metal ions required for the reaction of the incoming deoxynucleoside 5'-triphosphate with the 3'-hydroxyl in the primer terminus, while Glu156 may participate in nucleotide binding.     

A mechanism for the evolution of phosphorylation sites

Protein phosphorylation provides a mechanism for the rapid, reversible control of protein function. Phosphorylation adds negative charge to amino acid side chains, and negatively charged amino acids (Asp/Glu) can sometimes mimic the phosphorylated state of a protein. Using a comparative genomics approach, we show that nature also employs this trick in reverse by evolving serine, threonine, and tyrosine phosphorylation sites from Asp/Glu residues. Structures of three proteins where phosphosites evolved from acidic residues (DNA topoisomerase II, enolase, and C-Raf) show that the relevant acidic residues are present in salt bridges with conserved basic residues, and that phosphorylation has the potential to conditionally restore the salt bridges. The evolution of phosphorylation sites from glutamate and aspartate provides a rationale for why phosphorylation sometimes activates proteins, and helps explain the origins of this important and complex process.     

Identification of the active site of poly(A)-specific ribonuclease by site-directed mutagenesis and Fe(2+)-mediated cleavage.

Poly(A)-specific ribonuclease (PARN) is the only mammalian exoribonuclease characterized thus far with high specificity for degrading the mRNA poly(A) tail. PARN belongs to the RNase D family of nucleases, a family characterized by the presence of four conserved acidic amino acid residues. Here, we show by site-directed mutagenesis that these residues of human PARN, i.e. Asp(28), Glu(30), Asp(292), and Asp(382), are essential for catalysis but are not required for stabilization of the PARN x RNA substrate complex. We have used iron(II)-induced hydroxyl radical cleavage to map Fe(2+) binding sites in PARN. Two Fe(2+) binding sites were identified, and three of the conserved acidic amino acid residues were important for Fe(2+) binding at these sites. Furthermore, we show that the apparent dissociation constant ((app)K(d)) values for Fe(2+) binding at both sites were affected in PARN polypeptides in which the conserved acidic amino acid residues were substituted to alanine. This suggests that these residues coordinate divalent metal ions. We conclude that the four conserved acidic amino acids are essential residues of the PARN active site and that the active site of PARN functionally and structurally resembles the active site for 3'-exonuclease domain of Escherichia coli DNA polymerase I.     

Analysis by systematic deletion of amino acids of the action of the ribotoxin restrictocin

A series of contiguous deletions were made in a cDNA encoding the ribonuclease restrictocin with the purpose of identifying the amino acids that are essential for the cleavage of the phosphodiester bond on the 3' side of G4325 in the alpha-sarcin/ricin domain of mammalian (rat) 28S rRNA. In all 93 of 149 amino acids, 62% of the residues in restrictocin, were not essential for the action of the toxin. Of the five residues that have been proposed to constitute the active site, three could be deleted without loss of activity if they were part of a deletion of three or five amino acids but not if they were removed singly. It is likely that the loss of these three residues is compensated for by a neighboring residue that occupies the structural space created by the larger amino acid deletions. This was demonstrated to be the case for the active site residue Glu95 which in the deletion mutant Delta91-95 is replaced by Asp90. Systematic deletion of amino acids is a rapid, cost effective method for identifying the residues in a protein likely to contribute directly to function and, hence, deserving of closer scrutiny. Moreover, a semiquantitative estimate of the contribution of the residue to function can be made. For this reason the method may be useful for functional proteomics.     

Role of acidic residues as substrate determinants for casein kinase I

Sites phosphorylated by casein kinase I have been characterized by the presence of acidic amino acids NH2-terminal to the modified residue. Recently, phosphoserine was shown to be a particularly effective determinant for casein kinase I action when present in the motif -S(P)-X-X-S- (Flotow, H., Graves, P. R., Wang, A., Fiol, C. J., Roeske, R. W., and Roach, P. J. (1990) J. Biol. Chem. 265, 14264-14269). Nonetheless, nonphosphorylated substrates for casein kinase I are well documented. In this study, we examined the efficacy of Asp and Glu residues as determinants of casein kinase I action using synthetic peptide substrates. Peptides with runs of Asp residues in the motif Dn-X-X-S- were substrates for casein kinase I. Peptides with n = 3 or 4 were the most effective substrates, much better than n = 2. The peptide with n = 1, a single Asp residue, was a very poor substrate. A block of 4 Glu residues was a little less effective as a substrate determinant than 4 Asp residues in an otherwise identical peptide. The most effective substrate, with the motif -D-D-D-D-X-X-S-, was specific for casein kinase I and was not detectably phosphorylated by cyclic AMP-dependent protein kinase, casein kinase II, glycogen synthase kinase 3, or phosphorylase kinase and thus will be useful for the specific assay of casein kinase I. This peptide was nonetheless significantly worse as a substrate than peptides in which casein kinase I action was determined by phosphoserine in the -3 position. Still, the fact that Asp or Glu residues can specify a casein kinase I substrate suggests that acidic character has a role in substrate selection by this protein kinase.     

Mutational analysis of human DNase I at the DNA binding interface: implications for DNA recognition,catalysis, and metal ion dependence.

Human deoxyribonuclease I (DNase I), an enzyme used to treat cystic fibrosis patients, has been systematically analyzed by site-directed mutagenesis of residues at the DNA binding interface. Crystal structures of bovine DNase I complexed with two different oligonucleotides have implicated the participation of over 20 amino acids in catalysis or DNA recognition. These residues have been classified into four groups based on the characterization of over 80 human DNase I variants. Mutations at any of the four catalytic amino acids His 134, His 252, Glu 78, and Asp 212 drastically reduced the hydrolytic activity of DNase I. Replacing the three putative divalent metal ion-coordinating residues Glu 39, Asp 168, or Asp 251 led to inactive variants. Amino acids Gln 9, Arg 41, Tyr 76, Arg 111, Asn 170, Tyr 175, and Tyr 211 were also critical for activity, presumably because of their close proximity to the active site, while more peripheral DNA interactions stemming from 13 other positions were of minimal significance. The relative importance of these 27 positions is consistent with evolutionary relationships among DNase I across different species, DNase I-like proteins, and bacterial sphingomyelinases, suggesting a fingerprint for a family of DNase I-like proteins. Furthermore, we found no evidence for a second active site that had been previously implicated in Mn2+-dependent DNA degradation. Finally, we correlated our mutational analysis of human DNase I to that of bovine DNase I with respect to their specific activity and dependence on divalent metal ions.     

Functional analysis of polar amino-acid residues in membrane associated regions of the NHE1 isoform of the mammalian Na+/H+ exchanger.

The NHE1 isoform of the Na+/H+ exchanger is a ubiquitous plasma membrane protein that regulates intracellular pH in mammalian cells. Site-specific mutagenesis was used to examine the functional role of conserved, polar amino-acid residues occurring in segments of the protein associated with the membrane. Seventeen mutant proteins were assessed by characterization of intracellular pH changes in stably transfected cells that lacked an endogenous Na+/H+ exchanger. All of the mutant proteins were targeted correctly to the plasma membrane and were expressed at similar levels. Amino-acid residues Glu262 and Asp267 were critical to Na+/H+ exchanger activity while mutation of Glu391 resulted in only a partial reduction in activity. The Glu262-->Gln mutant was expressed partially as a deglycosylated protein with increased sensitivity to trypsin treatment in presence of Na+. Substitution of mutated Glu262, Asp267 and Glu391 with alternative acidic residues restored Na+/H+ exchanger activity. The Glu262-->Asp mutant had a decreased affinity for Li+, but its activity for Na+ and H+ ions was unaffected. The results support the hypothesis that side-chain oxygen atoms in a few, critically placed amino acids are important in Na+/H+ exchanger activity and the acidic amino-acid residues at positions 262, 267 and 391 are good candidates for being involved in Na+ coordination by the protein.     

PI-PfuI and PI-PfuII, intein-coded homing endonucleases from Pyrococcus furiosus. II. Characterization Of the binding and cleavage abilities by site-directed mutagenesis.

PI- Pfu I and PI- Pfu II from Pyrococcus furiosus are homing endonucleases, as shown in the accompanying paper. These two endonucleases are produced by protein splicing from the precursor protein including ribonucleotide reductase (RNR). We show here that both enzymes specifically interact with their substrate DNA and distort the DNA strands by 73 degrees and 67 degrees, respectively. They have two copies of the amino acid sequence motif LAGLIDADG, which is present in the majority of homing endonucleases and provides some of the catalytic residues necessary for DNA cleavage activity. Site-specific mutagenesis studies showed that two acidic residues in the motifs, Asp149 and Glu250 in PI- Pfu I, and Asp156 and Asp249 in PI- Pfu II, were critical for catalysis. The third residues of the active site triads, as predicted from the structure of PI- Sce I, were Asn225 in PI- Pfu I and Lys224 in PI- Pfu II. Substitution of Asn225 in PI- Pfu I by Ala did not affect catalysis. The cleavage activity of PI- Pfu II was 50-fold decreased by the substitution of Ala for Lys224. The binding affinity of the mutant protein for the substrate DNA also decreased 6-fold. The Lys in PI- Pfu II may play a direct or indirect role in catalysis of the endonuclease activity.     

Mutagenesis of two acidic active site residues in human muscle creatine kinase: implications for the catalytic mechanism

Creatine kinase (CK) catalyzes the reversible phosphorylation of the guanidine substrate, creatine, by MgATP. Although several X-ray crystal structures of various isoforms of creatine kinase have been published, the detailed catalytic mechanism remains unresolved. A crystal structure of the CK homologue, arginine kinase (AK), complexed with the transition-state analogue (arginine-nitrate-ADP), has revealed two carboxylate amino acid residues (Glu225 and Glu314) within 2.8 A of the proposed transphosphorylation site. These two residues are the putative catalytic groups that may promote nucleophilic attack by the guanidine amino group on the gamma-phosphate of ATP. From primary sequence alignments of arginine kinases and creatine kinases, we have identified two homologous creatine kinase acidic amino acid residues (Glu232 and Asp326), and these were targeted for examination of their potential roles in the CK mechanism. Using site-directed mutagenesis, we have made several substitutions at these two positions. The results indicate that of these two residues the Glu232 is the likely catalytic residue while Asp326 likely performs a role in properly aligning substrates for catalysis.     

Transport of taurocholate by mutants of negatively charged amino acids, cysteines, and threonines of the rat liver sodium-dependent taurocholate cotransporting polypeptide Ntcp.

The relevance of functional amino acids for taurocholate transport by the sodium-dependent taurocholate cotransporting polypeptide Ntcp was determined by site-directed mutagenesis. cRNA from 28 single-points mutants of the rat liver Ntcp clone was expressed in Xenopus laevis oocytes. Mutations were generated in five conserved negatively charged amino acids (aspartates and glutamates) which were present in nine members of the SBAT-family, in two nonconserved negatively charged amino acids, in all eight Ntcp-cysteines, and in two threonines from a protein kinase C consensus region of the Ntcp C-terminus. Functional amino acids were Asp115, Glu257, and Cys266, which were found to be essential for the maintenance of taurocholic acid transport. Asp115 is located in the large intracellular loop III, whereas Glu257 and Cys266 are located in the large extracellular loop VI. Four mutations of threonines from the C-terminus of the Ntcp by alanines or tyrosines showed no effects on sodium-dependent taurocholate transport. Introduction of the FLAG(R) motif into several transport negative point mutations demonstrated that all mutated proteins besides one were present within the cell membrane of the oocytes and provided proof that an insertion defect has not caused transport deficiency by these Ntcp mutants. The latter was observed only with the transport negative mutant Asp24Asn. In conclusion, loop amino acids are required for sodium-dependent substrate translocation by the Ntcp.     

Are acidic and basic groups in buried proteins predicted to be ionized?

Ionizable residues play essential roles in proteins, modulating protein stability, fold and function. Asp, Glu, Arg, and Lys make up about a quarter of the residues in an average protein. Multi-conformation continuum electrostatic (MCCE) calculations were used to predict the ionization states of all acidic and basic residues in 490 proteins. Of all 36,192 ionizable residues, 93.5% were predicted to be ionized. Thirty-five percent have lost 4.08 kcal/mol solvation energy (DeltaDeltaG(rxn)) sufficient to shift a pK(a) by three pH units in the absence of other interactions and 17% have DeltaDeltaG(rxn) sufficient to shift pK(a) by five pH units. Overall 85% of these buried residues (DeltaDeltaG(rxn)>5DeltapK units) are ionized, including 92% of the Arg, 86% of the Asp, 77% of the Glu, and 75% of the Lys. Ion-pair interactions stabilize the ionization of both acids and bases. The backbone dipoles stabilize anions more than cations. The interactions with polar side-chains are also different for acids and bases. Asn and Gln stabilize all charges, Ser and Thr stabilize only acids while Tyr rarely stabilize Lys. Thus, hydroxyls are better hydrogen bond donors than acceptors. Buried ionized residues are more likely to be conserved than those on the surface. There are 3.95 residues buried per 100 residues in an average protein.     

Identification of the histidine and aspartic acid residues essential for enzymatic activity of a family I.3 lipase by site-directed mutagenesis

A lipase from Pseudomonas sp. MIS38 (PML) is a member of the lipase family I.3. We analyzed the roles of the five histidine residues (His(30), His(274), His(291), His(313), and His(365)) and five acidic amino acid residues (Glu(253), Asp(255), Asp(262), Asp(275), and Asp(290)), which are fully conserved in the amino acid sequences of family I.3 lipases, by site-directed mutagenesis. We showed that the mutation of His(313) or Asp(255) to Ala almost fully inactivated the enzyme, whereas the mutations of other residues to Ala did not seriously affect the enzymatic activity. Measurement of the far- and near-UV circular dichroism spectra suggests that inactivation by the mutation of His(313) or Asp(255) is not due to marked changes in the tertiary structure. We propose that His(313) and Asp(255), together with Ser(207), form a catalytic triad in PML.     

Identification of structural and catalytic classes of highly conserved amino acid residues in lysine 2,3-aminomutase

Lysine 2,3-aminomutase (LAM) from Clostridium subterminale SB4 catalyzes the interconversion of (S)-lysine and (S)-beta-lysine by a radical mechanism involving coenzymatic actions of S-adenosylmethionine (SAM), a [4Fe-4S] cluster, and pyridoxal 5'-phosphate (PLP). The enzyme contains a number of conserved acidic residues and a cysteine- and arginine-rich motif, which binds iron and sulfide in the [4Fe-4S] cluster. The results of activity and iron, sulfide, and PLP analysis of variants resulting from site-specific mutations of the conserved acidic residues and the arginine residues in the iron-sulfide binding motif indicate two classes of conserved residues of each type. Mutation of the conserved residues Arg134, Asp293, and Asp330 abolishes all enzymatic activity. On the basis of the X-ray crystal structure, these residues bind the epsilon-aminium and alpha-carboxylate groups of (S)-lysine. However, among these residues, only Asp293 appears to be important for stabilizing the [4Fe-4S] cluster. Members of a second group of conserved residues appear to stabilize the structure of LAM. Mutations of arginine 130, 135, and 136 and acidic residues Glu86, Asp165, Glu236, and Asp172 dramatically decrease iron and sulfide contents in the purified variants. Mutation of Asp96 significantly decreases iron and sulfide content. Arg130 or Asp172 variants display no detectable activity, whereas variants mutated at the other positions display low to very low activities. Structural roles are assigned to this latter class of conserved amino acids. In particular, a network of hydrogen bonded interactions of Arg130, Glu86, Arg135, and the main chain carbonyl groups of Cys132 and Leu55 appears to stabilize the [4Fe-4S] cluster.     

Molecular cloning of matrin 3. A 125-kilodalton protein of the nuclear matrix contains an extensive acidic domain

We report here the cloning and sequencing of matrin 3, an acidic internal matrix protein, from a rat insuloma cDNA library. The nucleotide sequence has a single open reading frame encoding a polypeptide of 845 amino acids. The Genbank and National Biomedical Research Foundation databases did not contain any sequences similar to that of matrin 3. The primary structure consists of 33% charged residues and is generally hydrophilic. The amino-terminal region (residues 1-120) is positively charged and contains a large number of amino acids with free hydroxyl groups (26 of the first 100 residues) as in the lamins and several non-lamin intermediate filament proteins. A highly acidic domain (approximately 170 amino acids) near the carboxyl terminus, in which 32% of the amino acid residues are acidic (Glu or Asp), is a characteristic found in other nuclear proteins (Earnshaw, W. C. (1987) J. Cell Biol. 105, 1479-1482). A putative nuclear targeting signal sequence (Ser-Lys-Lys-Lys-Leu-Lys-Lys-Val-Glu) is located in the middle of the highly acidic domain. The corresponding human deduced partial amino acid sequence is 96% identical to the rat sequence, indicating that matrin 3 is a highly conserved protein.     

Cloning and characterization of the mitochondrial phosphate transport protein gene from the yeast Saccharomyces cerevisiae.

We have cloned the gene of the Saccharomyces cerevisiae phosphate transport protein (PTP), a member of the mitochondrial anion transport protein gene family. As PTP has a blocked N-terminus, we prepared three peptides. Oligonucleotides, based on their sequences, were used to screen a Yep24-housed genomic library. A total of 2073 bases of clone Y22 code for a 311 amino acid protein (Mr 32,814), which has similarities to the anion transport proteins: a triplicate gene structure and 6 hydrophobic segments. Typical for PTP, the triplicate gene structure possesses the X-Pro-X-(Asp/Glu)-X-X-(Lys/Arg)-X-(Arg/Lys)-X (X is an unspecified amino acid) motif and the very high homology only between the first and second repeat. The 6 hydrophobic segments harbor most of the 116 amino acids that are conserved between the yeast and the beef proteins. An N-terminal-extended signal sequence, as found in the beef protein, is absent. The yeast protein has about 33% fewer basic and acidic amino acids and five fewer Cys residues than the beef protein. The protein is insensitive to N-ethylmaleimide since Cys-42 (beef) has been replaced with a Thr. Mersalyl sensitivity has been retained and must be due to one of its three cysteines. Among these three cysteines, only Cys-28, located in the first hydrophobic segment, is conserved between the yeast and the beef protein.     

Biochemical characterization and identification of catalytic residues in alpha-glucuronidase from Bacillus stearothermophilus T-6.

Alpha-D-glucuronidases cleave the alpha-1,2-glycosidic bond of the 4-O-methyl-D-glucuronic acid side chain of xylan, as a part of an array of xylan hydrolyzing enzymes. The alpha-D-glucuronidase from Bacillus stearothermophilus T-6 was overexpressed in Escherichia coli using the T7 polymerase expression system. The purification procedure included two steps, heat treatment and gel filtration chromatography, and provided over 0.3 g of pure enzyme from 1 L of overnight culture. Based on gel filtration, the native protein is comprised of two identical subunits. Kinetic constants with aldotetraouronic acid as a substrate, at 55 degrees C, were a Km of 0.2 mM, and a specific activity of 42 U x mg(-1) (kcat = 54.9 s(-1)). The enzyme was most active at 65 degrees C, pH 5.5-6.0, in a 10-min assay, and retained 100% of its activity following incubation at 70 degrees C for 20 min. Based on differential scanning calorimetry, the protein denatured at 73.4 degrees C. Truncated forms of the enzyme, lacking either 126 amino acids from its N-terminus or 81 amino acids from its C-terminus, exhibited low residual activity, indicating that the catalytic site is located in the central region of the protein. To identify the potential catalytic residues, site-directed mutagenesis was applied on highly conserved acidic amino acids in the central region. The replacements Glu392-->Cys and Asp364-->Ala resulted in a decrease in activity of about five orders of magnitude, suggesting that these residues are the catalytic pair.     

Predicting DNA-binding amino acid residues from electrostatic stabilization upon mutation to Asp/Glu and evolutionary conservation

Binding of polyanionic DNA depends on the cluster of electropositive atoms in the binding site of a DNA-binding protein. Such a cluster of electropositive protein atoms would be electrostatically unfavorable without stabilizing interactions from the respective electronegative DNA atoms and would likely be evolutionary conserved due to its critical biological role. Consequently, our strategy for predicting DNA-binding residues is based on detecting a cluster of evolutionary conserved surface residues that are electrostatically stabilized upon mutation to negatively charged Asp/Glu residues. The method requires as input the protein structure and sufficient sequence homologs to define each residue's relative conservation, and it yields as output experimentally testable residues that are predicted to bind DNA. By incorporating characteristic DNA-binding site features (i.e., electrostatic strain and amino acid conservation), the new method yields a prediction accuracy of 83%, which is much higher than methods based on only electrostatic strain (57%) or conservation alone (50%). It is also less sensitive to protein conformational changes upon DNA binding than methods that mainly depend on the 3D protein structure.     

Protonation states of important acidic residues in the central Ca²⁺ ion binding sites of the Ca²⁺-ATPase: a molecular modeling study

The P-type ATPases are responsible for the transport of cations across cell membranes. The sarco(endo)plasmic reticulum Ca2?-ATPase (SERCA) transports two Ca2? ions from the cytoplasm to the lumen of the sarco(endo)plasmic reticulum and countertransports two or three protons per catalytic cycle. Two binding sites for Ca2? ions have been located via protein crystallography, including four acidic amino acid residues that are essential to the ion coordination. In this study, we present molecular dynamics (MD) simulations examining the protonation states of these amino acid residues in a Ca2?-free conformation of SERCA. Such knowledge will be important for an improved understanding of atomistic details of the transport mechanism of protons and Ca2? ions. Eight combinations of the protonation of four central acidic residues, Glu309, Glu771, Asp800, and Glu908, are tested from 10 ns MD simulations with respect to protein stability and ability to maintain a structure similar to the crystal structure. The trajectories for the most prospective combinations of protonation states were elongated to 50 ns and subjected to more detailed analysis, including prediction of pK(a) values of the four acidic residues over the trajectories. From the simulations we find that the combination leaving only Asp800 as charged is most likely. The results are compared to available experimental data and explain the observed destabilization upon full deprotonation, resulting in the entry of cytoplasmic K? ions into the Ca2? binding sites during the simulation in which Ca2? ions are absent. Furthermore, a hypothesis for the exchange of protons from the central binding cavity is proposed.