Catalysis of decarboxylation by an adjacent negative charge: a theoretical approach

We have calculated the rate acceleration in decarboxylation reactions that can be accomplished by charge-charge repulsion between the substrate carboxylate and an adjacent negative charge in media of various dielectric constants. It has been shown that a full negative charge or surrounding partial negative charges will have the same effect. It is concluded that the rate of decarboxylation could be greatly accelerated by the presence of a negative charge nearby. For example, in media with dielectric constants from 5.62 to 20.7, a 10(8)-fold rate acceleration could be achieved by a negative charge placed 3.77 A away from the substrate carboxylate group. However, pK(a) perturbation on two carboxylate groups at close proximity may limit the extent of catalysis. It should also be noted that the extent of catalysis does not change much when the dielectric constant varies from 5.62 to 20.7.     

Formation and stability of the enolates of N-protonated proline methyl ester and proline zwitterion in aqueous solution: a nonenzymatic model for the first step in the racemization of proline catalyzed by proline racemase

Rate constants for the hydrolysis of L-proline methyl ester to form proline and methanol in D(2)O buffered at neutral pD and 25 degrees C and the deuterium enrichment of the proline product determined by electrospray ionization mass spectrometry are reported. The data give k(DO) = 5.3 +/- 0.5 M(-1) s(-1) as the second-order rate constant for carbon deprotonation of N-protonated proline methyl ester by deuterioxide ion in D(2)O at 25 degrees C and I = 1.0 (KCl). The data provide good estimates of carbon acidities of pK(a) = 21 for N-protonated proline methyl ester and pK(a) = 29 for proline zwitterion in water and of the second-order rate constant k(HO) = 4.5 x 10(-5) M(-1) s(-1) for carbon deprotonation of proline zwitterion by hydroxide ion at 25 degrees C. There is no detectable acceleration of the deprotonation of N-protonated proline methyl ester by the Br?nsted base 3-quinuclidinone in water, and it is not clear that such Br?nsted catalysis would make a significant contribution to the rate acceleration for deprotonation of bound proline at proline racemase. A comparison of the first-order rate constants k(HO)[HO(-)] = 4.5 x 10(-11) s(-1) for deprotonation of free proline zwitterion in water at pH 8 and k(cat) = 2600 s(-1) for deprotonation of proline bound to the active site of proline racemase at pH 8 shows that the enzymatic rate acceleration for proline racemase is ca. 10(13)-fold. This corresponds to a 19 kcal/mol stabilization of the transition state for deprotonation of the enzyme-bound carbon acid substrate by interaction with the protein catalyst. It is suggested that (1) much of the rate acceleration of the enzymatic over the nonenzymatic reaction in water may result from transfer of the substrate proline zwitterion from the polar solvent water to a nonpolar enzyme active site and (2) the use of thiol anions rather than oxygen anions as Br?nsted bases at this putative nonpolar enzyme active site may be favored, because of the smaller energetic price for desolvation of thiol anions than for desolvation of the more strongly solvated oxygen anions.     

Electrostatic Exclusion of Neutral Solutes from Condensed DNA and Other Charged Phases

Motivated by experiments on condensed DNA phases in binary mixtures of water and a low-dielectric solute, we develop a theory for the electrostatic contribution to solute exclusion from a highly charged phase, within the continuum approximation of the medium. Because the electric field is maximum at the surface of each ion, the electrostatic energy is dominated by the Born energy; interactions between charges are of secondary importance. Neglecting interactions and considering only the competition between the Born energy and the free energy of mixing, we predict that low dielectric solutes are excluded from condensed DNA phases in water-cosolvent mixtures. This suggests that the traditional continuum electrostatic approach of modeling binary mixtures with a uniform dielectric constant needs to be modified. The linking of solute exclusion to solute dielectric properties also suggests a mechanism for predicting the electrostatic contribution to preferential hydration of polar and charged surfaces.     

Catalytic roles for proton transfer and protonation in ribozymes

Utilization of proton transfer in catalysis, which is well known in the mechanisms of protein enzymes, has been described only relatively recently for RNA enzymes. In this article, we present a current understanding of proton transfer by nucleic acids. Rate enhancement and specificity conferred by general acid-base catalysis are discussed. We also present possibilities for electrostatic catalysis from general acids and bases as well as cationic base pairs. The microenvironments of a large RNA provide the possibility of histidine-like pK(a)s for proton transfer, as well as lysine- and arginine-like pK(a)s for electrostatic catalysis. Discussion on proton transfer focuses on the hepatitis delta virus (HDV) and hairpin ribozymes, with select examples drawn from the protein literature. Discussion on electrostatic catalysis also draws on these two ribozymes, and a postulate for electrostatic catalysis by a cationic base pair in the mechanism of peptidyl transfer in the ribosome is presented. We also provide a perspective on possibilities for phosphoryl transfer mechanisms involving phosphorane intermediates and unusual tautomeric forms of the bases. Lastly, a distinction is made between ground state and "transition state" pK(a)s. We favor a model in which changes in pH lead to changes in the distribution of reactive and nonreactive ionizations of the ribozyme molecules in the ground state, and therefore suggest that "pK(a) changes in the transition state" do not provide an acceptable explanation for observed pH-rate profiles.     

Molecular dynamics effects on protein electrostatics

Electrostatic calculations have been carried out on a number of structural conformers of tuna cytochrome c. Conformers were generated using molecular dynamics simulations with a range of solvent simulating, macroscopic dielectric formalisms, and one solvent model that explicitly included solvent water molecules. Structures generated using the lowest dielectric models were relatively tight, with side chains collapsed on the surface, while those from the higher dielectric models had more internal and external fluidity, with surface side chains exploring a fuller range of conformational space. The average structure generated with the explicitly solvated model corresponded most closely with the crystal structure. Individual pK values, overall titration curves, and electrostatic potential surfaces were calculated for average structures and structures along each simulation. Differences between structural conformers within each simulation give rise to substantial changes in calculated local electrostatic interactions, resulting in pK value fluctuations for individual sites in the protein that vary by 0.3-2.0 pK units from the calculated time average. These variations are due to the thermal side chain reorientations that produce fluctuations in charge site separations. Properties like overall titration curves and pH dependent stability are not as sensitive to side chain fluctuations within a simulation, but there are substantial effects between simulations due to marked differences in average side chain behavior. These findings underscore the importance of proper dielectric formalism in molecular dynamics simulations when used to generate alternate solution structures from a crystal structure, and suggest that conformers significantly removed from the average structure have altered electrostatic properties that may prove important in episodic protein properties such as catalysis.     

On the calculation of electrostatic interactions in proteins

In this paper we present a classical treatment of electrostatic interactions in proteins. The protein is treated as a region of low dielectric constant with spherical charges embedded within it, surrounded by an aqueous solvent of high dielectric constant, which may contain a simple electrolyte. The complete analysis includes the effects of solvent screening, polarization forces, and self energies, which are related to solvation energies. Formulae, and sample calculations of forces and energies, are given for the special case of a spherical protein. Our analysis and model calculations point out that any consistent treatment of electrostatic interactions in proteins should account for the following. Solvent polarization is an important factor in the calculation of pairwise electrostatic interactions. Solvent polarization substantially affects both electrostatic energies and forces acting upon charges. No simple expression for the effective dielectric constant, Deff, can generally be valid, since Deff is a sensitive function of position. Solvent screening of pairwise interactions involving dipolar groups is less effective than the screening of charges. In fact for many interactions involving dipoles, solvent screening can be essentially ignored. The self energy of charges makes a large contribution to the total electrostatic energy of a protein. This must be compensated by specific interactions with other groups in the protein. Strategies for applying our analysis to proteins whose structures are known are discussed.     

A microscopic electrostatic model for the amphotericin B channel.

A microscopic model of an amphotericin B channel is proposed. The structure of the pores is generated using the atomic coordinates of the molecule in the structure determined experimentally by X-ray diffraction. The net charges of the atoms are determined by Mulliken analysis. With these charges the electrostatic energy profiles are calculated for a monovalent ion passing through the channels formed by different number of antibiotic molecules having different radii. The water inside the channel was considered through a continuum medium using the dielectric constant of the bulk, and the membrane contribution was included using the virtual images of the pore in a dielectric slab of epsilon = 3. The model satisfactorily explains the permeability and selectivity characteristics as well as other observations yet unexplained. The electrostatic profiles obtained reinforce the hypothesis of the existence of channels formed by a variable number of units.     

pH Dependence of catalysis by Pseudomonas aeruginosa isochorismate-pyruvate lyase: implications for transition state stabilization and the role of lysine 42

An isochorismate-pyruvate lyase with adventitious chorismate mutase activity from Pseudomonas aerugionsa (PchB) achieves catalysis of both pericyclic reactions in part by the stabilization of reactive conformations and in part by electrostatic transition-state stabilization. When the active site loop Lys42 is mutated to histidine, the enzyme develops a pH dependence corresponding to a loss of catalytic power upon deprotonation of the histidine. Structural data indicate that the change is not due to changes in active site architecture, but due to the difference in charge at this key site. With loss of the positive charge on the K42H side chain at high pH, the enzyme retains lyase activity at ～100-fold lowered catalytic efficiency but loses detectable mutase activity. We propose that both substrate organization and electrostatic transition state stabilization contribute to catalysis. However, the dominant reaction path for catalysis is dependent on reaction conditions, which influence the electrostatic properties of the enzyme active site amino acid side chains.     

The effect of cationic electrolytes on the adhesion of cells

The adhesion rate of cells under charge regulation onto a spherical collector with constant potential is investigated in this paper. Particularly, the effect of the presence of cationic electrolytes in the suspension medium on the adhesion rate is examined. The result reveals that the presence of cationic electrolytes in the suspension medium raises the electrostatic repulsion force between cell and collector surface, when the separation distance between them is small than a critical value. This has the effect of decreasing the adhesion rate of cells. The adhesion rate of cells is quite sensitive to the value of Hamaker constant, especially at a high ionic strength value.     

The effect of cationic electrolytes on the electrostatic behavior of cellular surfaces with ionizable groups

Based on the assumption that the electrostatic charges on the surface of sheep leukocytes arise from the dissociation of ionogenic groups, together with the presence of divalent cation (or trivalent cation) in the suspending medium of low ionic strength (or high ionic strength), the non-linear Poisson-Boltzmann equation for cell interaction with a solid surface with constant potential (or constant charge) is numerically solved in this paper. The cellular surface potential and the repulsive (or attractive) force is expressed as the function of separation distance. Because of shrinking the thickness of the electrostatic double layer at high ionic strength, the presence of cationic electrolyte has a less influential role on both the cellular surface potential and interaction force than at low ionic strength. However, due to the continuous equilibration of the ionogenic groups on the cellular surface as separation distance decreases, the presence of cationic electrolyte will not always reduce the interaction force during the whole adhesion period. The distance at which the cationic electrolyte changes its effect from positive to negative is termed the critical separation distance in this paper.     

Correlations of cationic charges with salt sensitivity and microbial specificity of cystine-stabilized beta -strand antimicrobial peptides

The electrostatic interaction of the charge cluster of an amphipathic peptide antibiotic with microbial membranes is a salt-sensitive step that often determines organism specificity. We have examined the correlation between charge clusters and salt insensitivity and microbial specificity in linear, cyclic, and retro-isomeric cystine-stabilized beta-strand (CSbeta) tachyplesin (TP) in a panel of 10 test organisms. Cyclic tachyplesins consisting of 14 and 18 amino acids are constrained by an end-to-end peptide backbone and two or three disulfide bonds to cross-brace the anti-parallel beta-strand that approximates a "beta-tile" structure. Circular dichroism measurements of beta-tile TPs showed that they displayed ordered structures. Control peptides containing the same number of basic amino acids as TP but lacking disulfide constraints were highly salt sensitive. Cyclic TP analogues with six cationic charges were more broadly active and salt-insensitive than those with fewer cationic charges. Reducing their proximity or number of cationic charges, particularly those with three or fewer basic amino acids, led to a significant decrease in potency and salt insensitivity, but an increased selectivity to certain Gram-positive bacteria. An end-group effect of the dibasic N-terminal Lys of TP in the open-chain TP and its retroisomer was observed in certain Gram-negative bacteria under high-salt conditions, an effect that was not found in the cyclic analogs. These results suggest that a stable folded structure together with three or more basic amino acids closely packed in a charged region in CSbeta peptides is important for salt insensitivity and organism specificity.     

A screened electrostatic potential model of hydration for conformational analysis of biomolecules

A screened electrostatic potential model of hydration which could avoid the time consuming Monte-Carlo procedures is derived from the study of the dielectric effect of a water molecule in the vicinity of two charges. It is shown that the molecular dielectric constant depends on the distance, sizes and signs of the interacting charges. The model is used to evaluate the hydration energy of the tetrapeptide (TYR-GLY-GLY-PHE) conformations. Results are in good agreement with Monte-Carlo calculations.     

Mapping the electrostatic potential within the ribosomal exit tunnel

Electrostatic potentials influence interactions among proteins and nucleic acids, the orientation of dipoles and quadrupoles, and the distribution of mobile charges. Consequently, electrostatic potentials can modulate macromolecular folding and conformational stability, as well as rates of catalysis and substrate binding. The ribosomal exit tunnel, along with its resident nascent peptide, is no less susceptible to these consequences. Yet, the electrostatics inside the tunnel have never been measured. Here we map both the electrostatic potential and accessibilities along the length of the tunnel and determine the electrostatic consequences of introducing a charged amino acid into the nascent peptide. To do this we developed novel probes and strategies. Our findings provide new insights regarding the dielectric of the tunnel and the dynamics of its local electric fields.     

Purification and properties of 4-hydroxybenzoate 1-hydroxylase (decarboxylating), a novel flavin adenine dinucleotide-dependent monooxygenase from Candida parapsilosis CBS604.

A novel flavoprotein monooxygenase, 4-hydroxybenzoate 1-hydroxylase (decarboxylating), from Candida parapsilosis CBS604 was purified to apparent homogeneity. The enzyme is induced when the yeast is grown on either 4-hydroxybenzoate, 2,4-dihydroxybenzoate, or 3,4-dihydroxybenzoate as the sole carbon source. The purified monooxygenase is a monomer of about 50 kDa containing flavin adenine dinucleotide as weakly bound cofactor. 4-Hydroxybenzoate 1-hydroxylase from C. parapsilosis catalyzes the oxidative decarboxylation of a wide range of 4-hydroxybenzoate derivatives with the stoichiometric consumption of NAD(P)H and oxygen. Optimal catalysis is reached at pH 8, with NADH being the preferred electron donor. By using (18)O2, it was confirmed that the oxygen atom inserted into the product 1,4-dihydroxybenzene is derived from molecular oxygen. 19F nuclear magnetic resonance spectroscopy revealed that the enzyme catalyzes the conversion of fluorinated 4-hydroxybenzoates to the corresponding hydroquinones. The activity of the enzyme is strongly inhibited by 3,5-dichloro-4-hydroxybenzoate, 4-hydroxy-3,5-dinitrobenzoate, and 4-hydroxyisophthalate, which are competitors with the aromatic substrate. The same type of inhibition is exhibited by chloride ions. Molecular orbital calculations show that upon deprotonation of the 4-hydroxy group, nucleophilic reactivity is located in all substrates at the C-1 position. This, and the fact that the enzyme is highly active with tetrafluoro-4-hydroxybenzoate and 4-hydroxy-3-nitrobenzoate, suggests that the phenolate forms of the substrates play an important role in catalysis. Based on the substrate specificity, a mechanism is proposed for the flavin-mediated oxidative decarboxylation of 4-hydroxybenzoate.     

Molecular dynamics simulation of α-melanocyte stimulating hormone in a water-membrane model interface

The conformation of the tridecapeptide α-melanocyte stimulating hormone in the presence of a double water-membrane interface was studied by molecular dynamics simulation, using the computational package THOR. In this program the solvent is represented by a continuous medium with dielectric constant ɛ, and the interface between different media is simulated by a surface of discontinuity of the dielectric constant. The electrostatic image method was used to write down the terms, added to the force field, that describe the polarisation effects induced in the interface by the atomic charges. The program was further improved by the introduction of a second surface, parallel to the first one, to mimic the membrane. A conformational search using the software Prelude was employed to find an initial geometry for the peptide in water. The molecular dynamics simulation performed during 10 ns showed that the peptide structure is flexible in water, without stabilisation of any preferential conformation. In the presence of the model membrane, the peptide moved to the medium representing the interior of the membrane. Inside the low dielectric constant medium, the structure of the peptide showed a turn in the central sequence of amino acids and a packed conformation remained stabilised during more than 7.0 ns of simulation. Received: 27 November 1998 / Revised version: 11 March 1999 / Accepted: 8 April 1999     

Experimental evaluation of the effective dielectric constant of proteins

Chemical modifications that alter the net charge of residues in reduction-oxidation proteins influence the redox potential of the protein by changing the electrostatic potential at the redox center. If the locations of the modified charges are known, the shift in redox potential may be used to determine the effective dielectric constant for the interactions between the redox center and modified residues. From the shift in redox potential upon charge neutralization of specific lysines in the hemoprotein cytochrome c, an effective dielectric constant of approximately 50 is calculated for the electrostatic interaction between the modified lysines and heme iron in the native protein.     

Reaction of Oxyhemoglobin with Carbon Monoxide

The reaction of oxyhemoglobin and carbon monoxide was studied kinetically at pH 7.8 in a variety of suspending media. The dielectric constant of the suspending media, as well as the viscosity (and hence the Fick diffusion coefficients), was varied with the use of glycine, glycerol, and sucrose. The results showed that the reaction was unaltered by the various additions to the media, provided that the pO₂ and the concentration of carbon monoxide were held constant. Since the concentration of oxygen varies from medium to medium at constant pO₂ while the pCO varies at constant concentration of carbon monoxide, the differences in the reactions with oxygen and carbon monoxide were emphasized. The lack of variation of the rate constants with changes in dielectric constant can be interpreted as indicating that electrostatic effects are unimportant in this reaction.     

Role of an active site guanine in hairpin ribozyme catalysis probed by exogenous nucleobase rescue

The hairpin ribozyme is a small catalytic RNA with reversible phosphodiester cleavage activity. Biochemical and structural studies exclude a requirement for divalent metal cation cofactors and implicate one active site nucleobase in particular, G8, in the catalytic mechanism. Our previous work demonstrated that the cleavage activity that is lost when G8 is replaced by an abasic residue is restored when certain nucleobases are provided in solution. The specificity and pH dependence of exogenous nucleobase rescue were consistent with several models of the rescue mechanism, including general acid base catalysis, electrostatic stabilization of negative charge in the transition state or a requirement for protonation to facilitate exogenous nucleobase binding. Detailed analyses of exogenous nucleobase rescue for both cleavage and ligation reactions now allow us to refine models of the rescue mechanism. Activity increased with increasing pH for both unmodified ribozyme reactions and unrescued reactions of abasic variants lacking G8. This similarity in pH dependence argues against a role for G8 as a general base catalyst, because G8 deprotonation could not be responsible for the pH-dependent transition in the abasic variant. Exogenous nucleobase rescue of both cleavage and ligation activity increased with decreasing pH, arguing against a role for rescuing nucleobases in general acid catalysis, because a nucleobase that contributes general acid catalysis in the cleavage pathway should provide general base catalysis in ligation. Analysis of the concentration dependence of cytosine rescue at high and low pH demonstrated that protonation promotes catalysis within the nucleobase-bound ribozyme complex but does not stabilize nucleobase binding in the ground state. These results support an electrostatic stabilization mechanism in which exogenous nucleobase binding counters negative charge that develops in the transition state.     

The PLP cofactor: lessons from studies on model reactions

Experimental probes of the acidity of weak carbon acids have been developed and used to determine the carbon acid pK(a)s of glycine, glycine derivatives and iminium ion adducts of glycine to the carbonyl group, including 5'-deoxypyridoxal (DPL). The high reactivity of the DPL-stabilized glycyl carbanion towards nucleophilic addition to both DPL and the glycine-DPL iminium ion favors the formation of Claisen condensation products at enzyme active sites. The formation of the iminium ion between glycine and DPL is accompanied by a 12-unit decrease in the pK(a) of 29 for glycine. The complicated effects of formation of glycine iminium ions to DPL and other aromatic and aliphatic aldehydes and ketones on carbon acid pK(a) are discussed. These data provide insight into the contribution of the individual pyridine ring substituents to the catalytic efficiency of DPL. It is suggested that the 5'-phosphodianion group of PLP may play an important role in enzymatic catalysis of carbon deprotonation by providing up to 12 kcal/mol of binding energy that is utilized to stabilize the transition state for the enzymatic reaction. This article is part of a Special Issue entitled: Pyridoxal Phospate Enzymology.     

Role of aspartate 400, arginine 262, and arginine 401 in the catalytic mechanism of human coproporphyrinogen oxidase

Coproporphyrinogen oxidase (CPO) is the sixth enzyme in the heme biosynthetic pathway, catalyzing two sequential oxidative decarboxylations of propionate moieties on coproporphyrinogen-III forming protoporphyrinogen-IX through a monovinyl intermediate, harderoporphyrinogen. Site-directed mutagenesis studies were carried out on three invariant amino acids, aspartate 400, arginine 262, and arginine 401, to determine residue contribution to substrate binding and/or catalysis by human recombinant CPO. Kinetic analyses were performed on mutant enzymes incubated with three substrates, coproporphyrinogen-III, harderoporphyrinogen, or mesoporphyrinogen-VI, in order to determine catalytic ability to perform the first and/or second oxidative decarboxylation. When Asp400 was mutated to alanine no divinyl product was detected, but the production of a small amount of monovinyl product suggested the K(m) value for coproporphyrinogen-III did not change significantly compared to the wild-type enzyme. Upon mutation of Arg262 to alanine, CPO was again a poor catalyst for the production of a divinyl product, with a catalytic efficiency <0.01% compared to wild-type, including a 15-fold higher K(m) for coproporphyrinogen-III. The efficiency of divinyl product formation for mutant enzyme Arg401Ala was approximately 3% compared to wild-type CPO, with a threefold increase in the K(m) value for coproporphyrinogen-III. These data suggest Asp400, Arg262, and Arg401 are active site amino acids critical for substrate binding and/or catalysis. Possible roles for arginine 262 and 401 include coordination of carboxylate groups of coproporphyrinogen-III, while aspartate 400 may initiate deprotonation of substrate, resulting in an oxidative decarboxylation.