pH-Dependent structural changes at Ca(2+)-binding sites of coagulation factor IX-binding protein

Coagulation factor IX-binding protein, isolated from Trimeresurus flavoviridis (IX-bp), is a C-type lectin-like protein. It is an anticoagulant consisting of homologous subunits, A and B. Each subunit has a Ca(2+)-binding site with a unique affinity (K(d) values of 14muM and 130muM at pH 7.5). These binding characteristics are pH-dependent and, under acidic conditions, the Ca(2+) binding of the low-affinity site was reduced considerably. In order to identify which site has high affinity and to investigate the pH-dependent Ca(2+) release mechanism, we have determined the crystal structures of IX-bp at pH 6.5 and pH 4.6 (apo form), and compared the Ca(2+)-binding sites with each other and with those of the solved structures under alkaline conditions; pH 7.8 and pH 8.0 (complexed form). At pH 6.5, Glu43 in the Ca(2+)-binding site of subunit A displayed two conformations. One (minor) is that in the alkaline state, and the other (major) is that at pH 4.6. However, the corresponding Gln43 residue of subunit B is in only a single conformation, which is almost identical with that in the alkaline state. At pH 4.6, Glu43 of subunit A adopts a conformation similar to that of the major conformer observed at pH 6.5, while Gln43 of subunit B assumes a new conformation, and both Ca(2+) positions are occupied by water molecules. These results showed that Glu43 of subunit A is much more sensitive to protonation than Gln43 of subunit B, and the conformational change of Glu43 occurs around pH6.5, which may correspond to the step of Ca(2+) release.     

Dissecting electrostatic interactions in Bacillus circulans xylanase through NMR-monitored pH titrations

NMR-monitored pH titration curves of proteins provide a rich source of structural and electrostatic information. Although relatively straightforward to measure, interpreting pH-dependent chemical shift changes to obtain site-specific acid dissociation constants (pK (A) values) is challenging. In order to analyze the biphasic titrations exhibited by the side chain (13)C(γ) nuclei of the nucleophilic Glu78 and general acid/base Glu172 in Bacillus circulans xylanase, we have revisited the formalism for the ionization equilibria of two coupled acidic residues. In general, fitting NMR-monitored pH titration curves for such a system will only yield the two macroscopic pK (A) values that reflect the combined effects of both deprotonation reactions. However, through the use of mutations complemented with ionic strength-dependent measurements, we are able to extract the four microscopic pK (Ai) values governing the branched acid/base equilibria of Glu78 and Glu172 in BcX. These data, confirmed through theoretical calculations, help explain the pH-dependent mechanism of this model GH11 xylanase by demonstrating that the kinetically determined pK (A) values and hence catalytic roles of these two residues result from their electrostatic coupling.     

How do lipases and esterases work: the electrostatic contribution

This work explores the role of one of the factors explaining lipase/esterase activity: the contribution of electrostatic interactions to lipase/esterase activity. The electrostatic potential distribution on the molecular surface of an enzyme as a function of pH determines, to a large extent, the enzyme's pH activity profile. Other important factors include the presence and distribution of polar and hydrophobic residues in the active cleft. We have mapped the electrostatic potential distribution as a function of pH on the molecular surface of nine lipases/esterases for which the 3D structure is experimentally known. A comparison of these potential maps at different pH values with the corresponding pH-activity profile, pH optimum or pH range where the activity displayed by the enzyme is maximum, has revealed a considerable correlation. A negative potential in the active site appears correlated with maximum activity towards triglycerides, which has prompted us to propose a model for product release ('The electrostatic catapult model') after cleavage of an ester bond. At the same time as the bottom of the active site cleft becomes negatively charged, other nearby regions also titrate and become negatively charged when pH becomes more alkaline, for some of the studied lipases. If such lipases also show phospholipase activity (such as guinea pig lipase-related proteins 2 chimera) we raise the hypothesis that such other titratable regions after becoming negatively charged might stabilise the positive charge present in the polar head of phospholipids, such as phosphatidylcholine, phosphatidylethanolamine and phosphatidylserine. The distribution of polar, weak polar and non-polar residues on the molecular surface of each studied lipase, in particular the active site region, was compared for all the lipases studied. The combination of graphical visualisation of the electrostatic potential maps and the polarity maps combined with knowledge about the location of key residues on the protein surface allows us to envision atomic models for lipolytic activity.     

Transfer RNAs: electrostatic patterns and an early stage of recognition by synthetases and elongation factor EF-Tu

Distributions of phosphate backbone-produced electrostatic potentials around several tRNAs were calculated by solving the nonlinear Poisson-Boltzmann equation. The tRNAs were either free or bound to the proteins involved in translation: aminoacyl-tRNA and elongation factor EF-Tu. We identified several regions of strong negative potential related to typical structural patterns of tRNA and invariant throughout the tRNAs. The patterns are conserved upon binding of tRNAs to the synthetase and the EF-Tu. Variation of tRNA charge in our theoretical calculations of electrostatic potential-mediated pK shifts of pH-dependent labels attached to tRNA, compared to experimentally observed pK shifts for those labels, shows that the total charge of tRNA is large, within the interval of -40 to -70 proton charges. The electrostatic field of tRNA is sufficient to cause ionization of histidine residues of ARSase, causing additional free energy of ARSase-tRNA interaction of at least several kcal/mol. This may discriminate proteins with respect to the particular tRNA at large distances. Two types of tRNA-protein electrostatic recognition mechanisms are discussed. One, more specific, involves charges induced on protein by the large electrostatic potential of tRNA, while the other, less specific, does not involve induced charges.     

Electrostatic interactions in the association of proteins: an analysis of the thrombin-hirudin complex.

The role of electrostatic interactions in stabilization of the thrombin-hirudin complex has been investigated by means of two macroscopic approaches: the modified Tanford-Kirkwood model and the finite-difference method for numerical solution of the Poisson-Boltzmann equations. The electrostatic potentials around the thrombin and hirudin molecules were asymmetric and complementary, and it is suggested that these fields influence the initial orientation in the process of the complex formation. The change of the electrostatic binding energy due to mutation of acidic residues in hirudin has been calculated and compared with experimentally determined changes in binding energy. In general, the change in electrostatic binding energy for a particular mutation calculated by the modified Tanford-Kirkwood approach agreed well with the experimentally observed change. The finite-difference approach tended to overestimate changes in binding energy when the mutated residues were involved in short-range electrostatic interactions. Decreases in binding energy caused by mutations of amino acids that do not make any direct ionic interactions (e.g., Glu 61 and Glu 62 of hirudin) can be explained in terms of the interaction of these charges with the positive electrostatic potential of thrombin. Differences between the calculated and observed changes in binding energy are discussed in terms of the crystal structure of the thrombin-hirudin complex.     

Electrostatic analysis and Brownian dynamics simulation of the association of plastocyanin and cytochrome f.

The oxidation of cytochrome f by the soluble cupredoxin plastocyanin is a central reaction in the photosynthetic electron transfer chain of all oxygenic organisms. Here, two different computational approaches are used to gain new insights into the role of molecular recognition and protein-protein association processes in this redox reaction. First, a comparative analysis of the computed molecular electrostatic potentials of seven single and multiple point mutants of spinach plastocyanin (D42N, E43K, E43N, E43Q/D44N, E59K/E60Q, E59K/E60Q/E43N, Q88E) and the wt protein was carried out. The experimentally determined relative rates (k(2)) for the set of plastocyanin mutants are found to correlate well (r(2) = 0.90 - 0.97) with the computed measure of the similarity of the plastocyanin electrostatic potentials. Second, the effects on the plastocyanin/cytochrome f association rate of these mutations in the plastocyanin "eastern site" were evaluated by simulating the association of the wild type and mutant plastocyanins with cytochrome f by Brownian dynamics. Good agreement between the computed and experimental relative rates (k(2)) (r(2) = 0.89 - 0.92) was achieved for the plastocyanin mutants. The results obtained by applying both computational techniques provide support for the fundamental role of the acidic residues at the plastocyanin eastern site in the association with cytochrome f and in the overall electron-transfer process.     

pKa calculations for class A beta-lactamases: methodological and mechanistic implications.

Beta-lactamases are responsible for resistance to penicillins and related beta-lactam compounds. Despite numerous studies, the identity of the general base involved in the acylation step is still unclear. It has been proposed, on the basis of a previous pKa calculation and analysis of structural data, that the unprotonated Lys73 in the active site could act as the general base. Using a continuum electrostatic model with an improved treatment of the multiple titration site problem, we calculated the pKa values of all titratable residues in the substrate-free TEM-1 and Bacillus licheniformis class A beta-lactamases. The pKa of Lys73 in both enzymes was computed to be above 10, in good agreement with recent experimental data on the TEM-1 beta-lactamase, but inconsistent with the proposal that Lys73 acts as the general base. Even when the closest titratable residue, Glu166, is mutated to a neutral residue, the predicted downward shift of the pKa of Lys73 shows that it is unlikely to act as a proton abstractor in either enzyme. These results support a mechanism in which the proton of the active Ser70 is transferred to the carboxylate group of Glu166.     

Computed pore potentials of the nicotinic acetylcholine receptor

Electrostatic surface potentials in the vestibule of the nicotinic acetylcholine receptor (nAChR) were computed from structural models using the University of Houston Brownian Dynamics program to determine their effect on ion conduction and ionic selectivity. To further determine whether computed potentials accurately reflect the electrostatic environment of the channel, the potentials were used to predict the rate constants for diffusion-enhanced fluorescence energy transfer; the calculated energy transfer rates are directly comparable with those determined experimentally (see companion article by Meltzer et al. in this issue). To include any effects on the local potentials by the bound acceptor fluorophore crystal violet, its binding site was first localized within the pore by fluorescence energy transfer measurements from dansyl-C6-choline bound to the agonist sites and also by simulations of binding using Autodock. To compare the computed potentials with those determined experimentally, we used the predicted energy transfer rates from Tb3+ chelates of varying charge to calculate an expected potential using the Boltzmann relationship. This expected potential (from -20 to -40 mV) overestimates the values determined experimentally (from -10 to -25 mV) by two- to fourfold at similar conditions of ionic strength. Although the results indicate a basic discrepancy between experimental and computed surface potentials, both methods demonstrate that the vestibular potential has a relatively small effect on conduction and selectivity.     

Dissecting the electrostatic interactions and pH-dependent activity of a family 11 glycosidase

Previous studies of the low molecular mass family 11 xylanase from Bacillus circulans show that the ionization state of the nucleophile (Glu78, pK(a) 4.6) and the acid/base catalyst (Glu172, pK(a) 6.7) gives rise to its pH-dependent activity profile. Inspection of the crystal structure of BCX reveals that Glu78 and Glu172 are in very similar environments and are surrounded by several chemically equivalent and highly conserved active site residues. Hence, there are no obvious reasons why their apparent pK(a) values are different. To address this question, a mutagenic approach was implemented to determine what features establish the pK(a) values (measured directly by (13)C NMR and indirectly by pH-dependent activity profiles) of these two catalytic carboxylic acids. Analysis of several BCX variants indicates that the ionized form of Glu78 is preferentially stabilized over that of Glu172 in part by stronger hydrogen bonds contributed by two well-ordered residues, namely, Tyr69 and Gln127. In addition, theoretical pK(a) calculations show that Glu78 has a lower pK(a) value than Glu172 due to a smaller desolvation energy and more favorable background interactions with permanent partial charges and ionizable groups within the protein. The pK(a) value of Glu172 is in turn elevated due to electrostatic repulsion from the negatively charged glutamate at position 78. The results also indicate that all of the conserved active site residues act concertedly in establishing the pK(a) values of Glu78 and Glu172, with no particular residue being singly more important than any of the others. In general, residues that contribute positive charges and hydrogen bonds serve to lower the pK(a) values of Glu78 and Glu172. The degree to which a hydrogen bond lowers a pK(a) value is largely dependent on the length of the hydrogen bond (shorter bonds lower pK(a) values more) and the chemical nature of the donor (COOH > OH > CONH(2)). In contrast, neighboring carboxyl groups can either lower or raise the pK(a) values of the catalytic glutamic acids depending upon the electrostatic linkage of the ionization constants of the residues involved in the interaction. While the pH optimum of BCX can be shifted from -1.1 to +0.6 pH units by mutating neighboring residues within the active site, activity is usually compromised due to the loss of important ground and/or transition state interactions. These results suggest that the pH optima of an enzyme might be best engineered by making strategic amino acid substitutions, at positions outside of the "core" active site, that electrostatically influence catalytic residues without perturbing their immediate structural environment.     

Examination of an active-site electrostatic node in the cAMP-dependent protein kinase catalytic subunit.

The active site of the cAMP-dependent protein kinase catalytic subunit harbors a cluster of acidic residues-Asp 127, Glu 170, Glu 203, Glu 230, and Asp 241-that are not conserved throughout the protein kinase family. Based on crystal structures of the catalytic subunit, these amino acids are removed from the site of phosphoryl transfer and are implicated in substrate recognition. Glu 230, the most buried of these acidic residues, was mutated to Ala (rC[E230A]) and Gln (rC[E230Q]) and overexpressed in Escherichia coli. In contrast to the mostly insoluble and destabilized rC[E230A], rC[E230Q] is largely soluble, purifies like wild-type enzyme, and displays wild-type-like thermal stability. The mutation in rC[E230Q] causes an order of magnitude decrease in the affinity for a heptapeptide substrate, Kemptide. In addition, two independent kinetic techniques were used to dissect phosphoryl transfer and product release steps in the reaction pathway. Viscosometric and pre-steady-state quench-flow analyses revealed that the phosphoryl transfer rate constant decreases by an order of magnitude, whereas the product release rate constant remains unperturbed. Electrostatic alterations in the rC[E230Q] active site were assessed using modeling techniques that provide molecular interpretations for the substrate affinity and phosphoryl transfer rate decreases observed experimentally. These observations indicate that subsite recognition elements in the catalytic subunit make electrostatic contributions that are important not only for peptide affinity, but also for catalysis. Protein kinases may, therefore, discriminate substrates by not only binding them tightly, but also by only turning over ones that complement the electrostatic character of the active site.     

pH-dependent stability of sperm whale myoglobin in water-guanidine hydrochloride solutions

An experimental-theoretical approach for the elucidation of protein stability is proposed. The theoretical prediction of pH-dependent protein stability is based on the macroscopic electrostatic model for calculation of the pH-dependent electrostatic free energy of proteins. As a test of the method we have considered the pH-dependent stability of sperm whale metmyoglobin. Two theoretical methods for evaluation of the electrostatic free energy and p K values are applied: the finite-difference Poisson-Boltzmann method and the semiempirical approach based on the modified Tanford-Kirkwood theory. The theoretical results for electrostatic free energy of unfolding are compared with the experimental data for guanidine hydrochloride unfolding under equilibrium conditions over a wide pH range. Using the optical parameters of the Soret absorbance to monitor conformational equilibrium and Tanford's method to estimate the resulting data, it was found that the conformational free energy of unfolding of metmyoglobin is 16.3 kcal mol(-1) at neutral pH values. The total unfolding free energies were calculated on the basis of the theoretically predicted electrostatic unfolding free energies and the experimentally measured midpoints (pH(1/2)) of acidic and alkaline denaturation transitions. Experimental data for alkaline denaturation were used for the first time in theoretical analysis of the pH-dependent unfolding of myoglobin. The present results demonstrate that the simultaneous application of appropriate theoretical and experimental methods permits a more complete analysis of the pH-dependent and pH-independent properties and stability of globular proteins.     

Multiple substrate binding states and chiral recognition in cofactor-independent glutamate racemase: a molecular dynamics study

Glutamate racemase (MurI) catalyzes the racemization of glutamate; two cysteine residues serve as catalytic acid and base. On the basis of the crystal structure of MurI from the hyperthermophilic bacterium Aquifex pyrophilus, we performed molecular dynamics (MD) simulations of six different systems to investigate stereochemistry, substrate ligation, and active site protonation state. The catalytic competence of individual systems was assessed by the abundance of reactive conformers. Only systems in which Cys70 is poised to deprotonate d-Glu were found to be catalytically competent (idem Cys178/l-Glu), in agreement with the experimentally observed stereochemistry of Lactobacillus fermentii MurI [Tanner, M. E. et al. (1993) Biochemistry 32, 3998-4006]. Only systems in which the alpha-amino group of l/d-Glu and the imidazole moiety of His are deprotonated are catalytically competent. The active site of MurI displays an unusual flexibility in substrate ligation, and several transitions between stable binding patterns were observed. In catalytically competent binding states, the conserved threonine residues 72, 114, and 117 ligate the alpha-carboxylate of Glu and the Asn71 amides ligate the alpha-amino group of Glu, whereas the delta-carboxylate of Glu is steered by electrostatic repulsion from the Asp7 and Glu147 side chain carboxylates. A network of hydrogen bonds controls the positioning of each thiol/thiolate. In what we term substrate flipping, Glu suddenly rotates into a binding pattern that resembles the post-racemization state of the other enantiomer, i.e., each enantiomer can be bound in two distinct states. Substrate flipping and unfavorable substrate binding successively trigger dissociation of the substrate, accompanied by an opening of the active site channel. We explain how the weak binding of Glu contributes to catalysis and suggest a mechanism by which binding mismatches are propagated into an opening of the active site.     

Investigations of Sso7d catalytic residues by NMR titration shifts and electrostatic calculations

Sso7d is a small basic protein consisting of 62 amino acids isolated from the thermoacidophilic archeobacterium Sulfolobus solfataricus. The protein is endowed with DNA binding properties, RNase activity, and the capability of rescuing aggregated proteins in the presence of ATP. In this study, the electrostatic properties of Sso7d are investigated by using the Poisson-Boltzmann calculation of the surface potential distribution and following by NMR spectroscopy the proton chemical shift pH titration of acidic residues. Although the details of the catalytic mechanism still have to be defined, the results from NMR experiments confirm the possible involvement of Glu35 as the proton acceptor in the catalytic reaction, as seen by its abnormally high pK(a) value. Poisson-Boltzmann calculations and NMR titration shifts suggest the presence of a possible hydrogen bond between Glu35 and Tyr33, with a consequent rather rigid arrangement at these positions. Comparison with RNase T1 suggests that Tyr7 may be a good candidate for acting as a proton donor in the active site of Sso7d as shown by its low phenolic pK(a) of approximately 9.3. Titration experiments performed with the UpA, a RNA dinucleotide model, showed that the protein residues affected by the interaction are mainly located in a different region with respect to the surface affected by DNA recognition, in good agreement with the surface potential distribution found with electrostatic calculations.     

Electrostatic steering at acetylcholine binding sites

The electrostatic environments near the acetylcholine binding sites on the nicotinic acetylcholine receptor (nAChR) and acetylcholinesterase were measured by diffusion-enhanced fluorescence energy transfer (DEFET) to determine the influence of long-range electrostatic interactions on ligand binding kinetics and net binding energy. Changes in DEFET from variously charged Tb3+ -chelates revealed net potentials of -20 mV at the nAChR agonist sites and -14 mV at the entrance to the AChE active site, in physiological ionic strength conditions. The potential at the alphadelta-binding site of the nAChR was determined independently in the presence of d-tubocurarine to be -14 mV; the calculated potential at the alphagamma-site was approximately threefold stronger than at the alphadelta-site. By determining the local potential in increasing ionic strength, Debye-Hückel theory predicted that the potentials near the nAChR agonist binding sites are constituted by one to three charges in close proximity to the binding site. Examination of the binding kinetics of the fluorescent acetylcholine analog dansyl-C6-choline at ionic strengths from 12.5 to 400 mM revealed a twofold decrease in association rate. Debye-Hückel analysis of the kinetics revealed a similar charge distribution as seen by changes in the potentials. To determine whether the experimentally determined potentials are reflected by continuum electrostatics calculations, solutions to the nonlinear Poisson-Boltzmann equation were used to compute the potentials expected from DEFET measurements from high-resolution models of the nAChR and AChE. These calculations are in good agreement with the DEFET measurements for AChE and for the alphagamma-site of the nAChR. We conclude that long-range electrostatic interactions contribute -0.3 and -1 kcal/mol to the binding energy at the nAChR alphadelta- and alphagamma-sites due to an increase in association rates.     

Non-empirical study of the phosphorylation reaction catalyzed by 4-methyl-5-β-hydroxyethylthiazole kinase: relevance of the theory of intermolecular interactions

The subject of this study was an analysis of the role of active site residues in the phosphoryl transfer reaction catalyzed by 4-methyl-5-β-hydroxyethylthiazole kinase (ThiK). The ThiK-catalyzed reaction is of special interest due to the lack of a highly conserved aspartate residue serving as a catalytic base. ONIOM(B3LYP:PM3) models of stationary points along the reaction pathway consisted of reactants, two magnesium ions and several highly conserved ThiK active site residues. The results indicate that an S_N2-like mechanism of ThiK, with γ-phosphate acting as an alcohol-activating base is reasonable. Geometries of substrates, transition state and products were utilized in the non-empirical analysis of the physical nature of catalytic interactions taking place in the ThiK active site. The role of particular residues was investigated in terms of their ability to preferentially stabilize the transition state relative to substrates (differential transition state stabilization, DTSS) or products (differential product stabilization, DPS). It seems that Mg2, Glu126 and Cys198 play a major catalytic role, whereas Mg1 and the same Cys198 are responsible for product release. It is remarkable that no dominant role of an electrostatic term in the interactions involved in catalytic activity is observed for product release. Determination of catalytic fields expressing differential electrostatic potential of the transition state with respect to substrates revealed the optimal electrostatic features of an ideal catalyst for the studied reaction. The predicted catalytic environment is in agreement with experimental data showing increased catalytic activity of ThiK upon mutation of Cys198 to aspartate. Figure Catalytic fields for ThiK-catalyzed reaction juxtaposed with the positions of active site residues of a model system. Magnesium ions are considered part of the transition state/reactants. The surface of constant electronic density is colored according to differential electrostatic potential of transition state with respect to reactants. The sign of the differential potential reflects the electrostatic properties of a complementary molecular environment. Red (green) color denotes regions where a negative (positive) charge would be optimal for catalytic activity     

Topological role of cytochrome P450 2D6 active site residues

Recent reports have identified Phe120, Asp301, Thr309, and Glu216 as important residues in cytochrome P450 2D6 (CYP2D6) substrate binding and catalysis. Complementary homology models have located these amino acids within the binding pocket of CYP2D6 and in the present study we have used aryldiazenes to test these models and gain further insight in the role these amino acids have in maintaining the integrity of the active site cavity. When Phe120 was replaced to alanine, there was a significant increase in probe migration to pyrrole nitrogens C and D, in agreement with homology models which have located the phenyl side-chain of Phe120 above these two pyrrole rings. No changes in topology were observed with the D301Q mutant, supporting claims that in this mutant the electrostatic interactions with the B/C-loop are largely maintained and the loop retains its native orientation. The T309V mutation resulted in significant topological alteration suggesting that, in addition to its potential role in dioxygen activation, Thr309 plays an important structural role within the active site crevice. Replacement of Ile106 with Glu, engineered to cause electrostatic repulsion with Glu216, had a profound topological effect in the higher region within the active site cavity and impaired the catalytic activity towards CYP2D6 probe substrates.     

Electrostatic attraction by surface charge does not contribute to the catalytic efficiency of acetylcholinesterase.

Acetylcholinesterases (AChEs) are characterized by a high net negative charge and by an uneven surface charge distribution, giving rise to a negative electrostatic potential extending over most of the molecular surface. To evaluate the contribution of these electrostatic properties to the catalytic efficiency, 20 single- and multiple-site mutants of human AChE were generated by replacing up to seven acidic residues, vicinal to the rim of the active-center gorge (Glu84, Glu285, Glu292, Asp349, Glu358, Glu389 and Asp390), by neutral amino acids. Progressive simulated replacement of these charged residues results in a gradual decrease of the negative electrostatic potential which is essentially eliminated by neutralizing six or seven charges. In marked contrast to the shrinking of the electrostatic potential, the corresponding mutations had no significant effect on the apparent bimolecular rate constants of hydrolysis for charged and non-charged substrates, or on the Ki value for a charged active center inhibitor. Moreover, the kcat values for all 20 mutants are essentially identical to that of the wild type enzyme, and the apparent bimolecular rate constants show a moderate dependence on the ionic strength, which is invariant for all the enzymes examined. These findings suggest that the surface electrostatic properties of AChE do not contribute to the catalytic rate, that this rate is probably not diffusion-controlled and that long-range electrostatic interactions play no role in stabilization of the transition states of the catalytic process.     

The effects of modifying the surface charge on the catalytic activity of a thermolysin-like protease

The impact of long range electrostatic interactions on catalysis in the thermolysin-like protease from Bacillus stearothermophilus was studied by analyzing the effects of inserting or removing charges on the protein surface. Various mutations were introduced at six different positions, and double-mutant cycle analysis was used to study the extent to which mutational effects were interdependent. The effects of single point mutations on the k(cat)/K(m) were non-additive, even in cases where the point mutations were located 10 A or more from the active site Zn(2+) and separated from each other by up to 25 A. This shows that catalysis is affected by large electrostatic networks that involve major parts of the enzyme. The interdependence of mutations at positions as much as 25 A apart in space also indicates that other effects, such as active site dynamics, play an important role in determining active site electrostatics. Several mutations yielded a significant increase in the activity, the most active (quadruple) mutant being almost four times as active as the wild type. In some cases the shape of the pH-activity profile was changed significantly. Remarkably, large changes in the pH-optimum were not observed.     

Contributions of counter-charge in a potassium channel voltage-sensor domain

Voltage-sensor domains couple membrane potential to conformational changes in voltage-gated ion channels and phosphatases. Highly coevolved acidic and aromatic side chains assist the transfer of cationic side chains across the transmembrane electric field during voltage sensing. We investigated the functional contribution of negative electrostatic potentials from these residues to channel gating and voltage sensing with unnatural amino acid mutagenesis, electrophysiology, voltage-clamp fluorometry and ab initio calculations. The data show that neutralization of two conserved acidic side chains in transmembrane segments S2 and S3, namely Glu293 and Asp316 in Shaker potassium channels, has little functional effect on conductance-voltage relationships, although Glu293 appears to catalyze S4 movement. Our results suggest that neither Glu293 nor Asp316 engages in electrostatic state-dependent charge-charge interactions with S4, likely because they occupy, and possibly help create, a water-filled vestibule.     

The effect of pH on the transpeptidation and hydrolytic reactions of rat kidney gamma-glutamyltransferase

The effect of pH upon the transpeptidation and hydrolytic reactions of gamma-glutamyltransferase [5-glutamyl)-peptide:amino-acid 5-glutamyltransferase, EC 2.3.2.2) have been investigated. It was found that the enzyme was irreversibly inactivated below pH 7.5 or above pH 9.4. Transpeptidation was markedly pH-dependent, while hydrolysis was pH-independent. The pH optimum for transpeptidation was found to vary for different acceptors. The ascending limb of the pH-optimum curve is attributed to the pK of the alpha-amino group of the acceptor, while the descending limb of the pH-optimum curve is attributed to an ionisable group in the active site of the enzyme. These observations provide much information about the interaction of the enzyme with the acceptor: (1) the true acceptor for gamma-glutamyltransferase is the deprotonated form of the amino acid; (2) glycylglycine has a similar acceptor activity to methionine, its apparent higher activity being due to the low pK of the alpha-amino group; (3) the enzyme is reversibly inactivated at higher pH by the deprotonation of a group in the active site which is involved in both binding of acceptor and catalysis of transpeptidation (this group is not involved in the hydrolysis reaction); (4) at pH 8.5, the normal pH for assay, only 47% of the enzyme is active, while at pH 7.4 gamma-glutamyltransferase is 93% in the active form.