Modeling of the electrostatic potential field of plastocyanin

The DelPhi computer program is used to calculate the electrostatic potential field of the photosynthetic electron transport protein plastocyanin. Knowledge of the potential field is important for understanding the mechanisms by which plastocyanin interacts with other charged reagents. The program uses a macroscopic, continuum approach in which the protein and solvent are assigned different dielectric constants, the crystal structure of the protein defines the dielectric boundary, and the ionic strength of the solvent is taken into account. The potential field is determined by numerically solving the Poisson-Boltzmann equation. The field surrounding plastocyanin is characterized by a region of positive potential over the copper center active site, and a region of negative potential over the adjacent association site containing tyrosine 83. The shape and magnitude of the potential field shows a strong dependence on the ionic strength and pH of the solvent. The program is able to accurately predict the effect of the copper center oxidation state on the pKa of a tetranitromethane derivative of tyrosine 83 using an intrinsic protein dielectric constant of 2 to 4. Evidence is also presented that the glutamate 68 side chain is exposed to the solvent to a greater extent in the solution structure of plastocyanin than in the crystal structure.     

Differential electrostatic stabilization of A-, B-, and Z-forms of DNA

The static accessibility discrete charge algorithm for protein charge interactions is extended to the case of linear polyelectrolytes. In this model, the effective dielectric value between surface charge sites depends predominantly on the solvent ionic strength and the solvent accessibilities of the charge sites. This treatment accounts for the phenomena of specific ion binding in the context of a general electrostatic effect [Matthew and Richards (1982) Biochemistry 21 , 4989]. Specific ion sites are determined by locating areas of high electrostatic potential at the solvent interface of the macromolecule. At a given ionic strength the calculated potential at a site is taken to describe a binding constant and therefore the ion site occupancy. For a 20-base-pair fragment of B-DNA, net charge of ?40, 16 ion sites are indicated in the minor groove. The partial occupancy of each site increases from 0.2 to 0.5 as the ionic strength is increased from 0.01 to 0.50. Over the same range of ionic strength, the electrostatic free energy of this charge array is calculated to change from +0.6 to ?0.05 kcal/bp. Parallel behavior is predicted for A- and Z-DNA charge geometries. The most stable configuration, based on electrostatic criteria, at high ionic strength (I = 0.1–0.5) is that of Z-DNA. In this range, the ratio of “bound” sodium to phosphate is predicted to be less than 0.4.     

溶液介电常数对铜锌超氧化物歧化酶活性的影响

溶液介电常数对天然酶和修饰酶的活性影响不同,天然酶随介电常数增加而酶活性下降,修饰酶则反之,这表明静电相互作用在铜锌超氧化物歧化酶(Cu·Zn-SOD)与超氧阴离子(-O_2~(·-))反应过程中起着重要作用,酶分子活性中心附近ε-NH_3~+为O_2~(·-)进入活性中心提供静电吸引力。在有机溶剂中,SOD的构象会发生变化,从而导致酶活性降低。实验还表明,Cl~-对SOD有明显的抑制作用。     

Electrostatic effect upon association of reduced nicotinamide adenine dinucleotide and equine liver alcohol dehydrogenase

The rate of association of equine liver alcohol dehydrogenase and its coenzymes exhibits a large pH dependence with slower rates at basic pH and an observed kinetic pKa value of approximately 9-9.5. This pH dependence has been explained by invoking local active site electrostatic effects which result in repulsion of the negatively charged coenzyme and the ionized hydroxyl anion form of the zinc-bound water molecule. We have examined a simpler hypothesis, namely, that the pH dependence results from the electrostatic interaction of the coenzyme and the enzyme which changes from an attractive interaction of the negatively charged coenzyme and the positively charged enzyme to a repulsive interaction between the two negatively charged species at the isoelectric point for the enzyme (pH 8.7). We have tested this proposal by examining the ionic strength dependence of the association rate constant at various pH values. These data have been interpreted by using the Wherland-Gray equation, which we have shown can be applied to the kinetics of enzyme-coenzyme association. Our results indicate that the shielding of the buffer electrolyte changes from a negative to a positive value as the charge on the protein changes at the isoelectric point. This result is exactly that which is predicted for electrostatic effects that depend on the charge of the protein molecule and is not consistent with predictions based upon the local active site effects. At low ionic strength values of 10 mM or less, approximately 75% of the observed pH dependence results from the enzyme electrostatic effects; the remaining pH dependence may result from active site effects.(ABSTRACT TRUNCATED AT 250 WORDS)     

Nicotinic acetylcholine receptor channel electrostatics determined by diffusion-enhanced luminescence energy transfer

The electrostatic potentials within the pore of the nicotinic acetylcholine receptor (nAChR) were determined using lanthanide-based diffusion-enhanced fluorescence energy transfer experiments. Freely diffusing Tb3+ -chelates of varying charge constituted a set of energy transfer donors to the acceptor, crystal violet, a noncompetitive antagonist of the nAChR. Energy transfer from a neutral Tb3+ -chelate to nAChR-bound crystal violet was reduced 95% relative to the energy transfer to free crystal violet. This result indicated that crystal violet was strongly shielded from solvent when bound to the nAChR. Comparison of energy transfer from positively and negatively charged chelates indicate negative electrostatic potentials of -25 mV in the channel, measured in low ionic strength, and -10 mV measured in physiological ionic strength. Debye-Hückel analyses of potentials determined at various ionic strengths were consistent with 1-2 negative charges within 8 A of the crystal violet binding site. To complement the energy transfer experiments, the influence of pH and ionic strength on the binding of [3H]phencyclidine were determined. The ionic strength dependence of binding affinity was consistent with -3.3 charges within 8 A of the binding site, according to Debye-Hückel analysis. The pH dependence of binding had an apparent pKa of 7.2, a value indicative of a potential near -170 mV if the titratable residues are constituted of aspartates and glutamates. It is concluded that long-range potentials are small and likely contribute little to selectivity or conductance whereas close interactions are more likely to contribute to electrostatic stabilization of ions and binding of noncompetitive antagonists within the channel.     

Point charge distributions and electrostatic steering in enzyme/substrate encounter: Brownian dynamics of modified copper/zinc superoxide dismutases

The electrostatic steering mechanism of bovine erythrocyte Cu/Zn superoxide dismutase (SOD) was investigated through the use of Brownian dynamics. Simulations of enzyme/substrate encounter were carried out on 14 different SOD models defined by simple changes in the enzyme's point charge distribution. The magnitude and ionic strength dependence of reaction rates, rates for collision anywhere on the enzyme surface, and collision efficiency factors were analyzed to elucidate both the general and specific roles for point charges associated with amino acid residues. Collision rates for the general enzyme surface appear to be solely determined by the net charge on the enzyme. At physiological ionic strength this effect is negligible, with only 6% variation in collision rates observed as the net charge ranges from +2e to -10e. With the exception of a few charged residues in the active-site channel of SOD, point charge modifications had modest effects on reaction rates. For a large region within and surrounding the channel, reaction rates increased or decreased by only 10-15% with the addition or subtraction of a protonic unit of charge, respectively. This effect simply disappeared with increasing distance from the active site. More dramatic effects were seen at only three residues: arginine-141, glutamate-131, and lysine-134. Implications for rate enhancement through site-directed mutagenesis are discussed.     

Carbamoylation of Cu,Zn-superoxide dismutase by cyanate: Role of lysines in the enzyme action

Reaction with cyanate leads to a reversible change of the EPR spectrum of Cu,Zn-superoxide dismutase and to time-dependent carbamoylation of the lysine residues of the enzyme, producing a stable covalent derivative with more negative charge. The carbamoylated enzyme is less active than the native enzyme in spite of unaltered EPR spectra. The extent of this inactivation is much less when the enzyme activity is measured at low ionic strength. These results show that integrity of the active site is not the sole factor playing a role in the enzyme mechanism and that the ionic strength effect is related to electrostatic interactions between O⁻₂ and surface charges of the protein.     

Simulation of the diffusion-controlled reaction between superoxide and superoxide dismutase. II. Detailed models

A two-stage Brownian dynamics simulation method is used to study the diffusion-influenced bimolecular reaction between superoxide and superoxide dismutase (SOD). The crystal structure of the dimeric enzyme is used in constructing detailed topographical and electrostatic models. Several electrostatic models are considered. In the most realistic, the excluded volume of the protein, which is impermeable to penetration by mobile ions, is assigned a dielectric constant of 2 and the surrounding “solvent” is assigned a value of 78. A finite difference method is used to solve the linearized Poisson–Boltzmann equation. For native SOD, the simulations reproduce the pronounced salt dependence of the rate constant observed experimentally. This salt dependence is attributed to electrostatic interactions between enzyme and substrate that are inherently attractive and amplified by the low dielectric constant of the protein interior. The simulation method is also applied to a modified enzyme, acylated SOD.     

The electrostatic potential of Escherichia coli dihydrofolate reductase

Escherichia coli dihydrofolate reductase (DHFR) carries a net charge of -10 electrons yet it binds ligands with net charges of -4 (NADPH) and -2 (folate or dihydrofolate). Evaluation and analysis of the electrostatic potential of the enzyme give insight as to how this is accomplished. The results show that the enzyme is covered by an overall negative potential (as expected) except for the ligand binding sites, which are located inside "pockets" of positive potential that enable the enzyme to bind the negatively charged ligands. The electrostatic potential can be related to the asymmetric distribution of charged residues in the enzyme. The asymmetric charge distribution, along with the dielectric boundary that occurs at the solvent-protein interface, is analogous to the situation occurring in superoxide dismutase. Thus DHFR is another case where the shape of the active site focuses electric fields out into solution. The positive electrostatic potential at the entrance of the ligand binding site in E. coli DHFR is shown to be a direct consequence of the presence of three positively charged residues at positions 32, 52, and 57--residues which have also been shown recently to contribute significantly to electronic polarization of the ligand folate. The latter has been postulated to be involved in the catalytic process. A similar structural motif of three positively charged amino acids that gives rise to a positive potential at the entrance to the active site is also found in DHFR from chicken liver, and is suggested to be a common feature in DHFRs from many species. It is noted that, although the net charges of DHFRs from different species vary from +3 to -10, the enzymes are able to bind the same negatively charged ligands, and perform the same catalytic function.     

Glycocalyx electrostatic potential profile analysis: ion, pH, steric, and charge effects

The Poisson-Boltzmann equation is modified to consider charge ionogenicity, steric exclusion, and charge distribution in order to describe the perimembranous electrostatic potential profile in a manner consistent with the known morphology and biochemical composition of the cell's glycocalyx. Exact numerical and approximate analytical solutions are given for various charge distributions and for an extended form of the Donnan potential model. The interrelated effects of ionic conditions, bulk pH, ion binding, local dielectric, steric volume exclusion, and charge distribution on the local potential, pH, and charge density within the glycocalyx are examined. Local charge-induced, potential-mediated pH reductions cause glycocalyx charge neutralization. Under certain conditions, local potentials may be insensitive to ionic strength or may decrease in spite of increasing charge density. The volume exclusion of the glycocalyx reduces the local ion concentration, thereby increasing the local potential. With neutral lipid membranes, the Donnan and surface potential agree if the glycocalyx charge distribution is both uniform and several times thicker than the Debye length (approximately 20 A in thickness under physiological conditions). Model limitations in terms of application to microdomains or protein endo- and ectodomains are discussed.     

Kinetics of reduction of high redox potential ferredoxins by the semiquinones of Clostridium pasteurianum flavodoxin and exogenous flavin mononucleotide. Electrostatic and redox potential effects

We have measured the ionic strength dependence of the rate constants for the electron-transfer reactions of flavin mononucleotide (FMN) and flavodoxin semiquinones with 10 high redox potential ferredoxins (HiPIP's). The rate constants were extrapolated to infinite ionic strength by using a theoretical model of electrostatic interactions developed in our laboratory. In all cases, the sign of the electrostatic interaction was the same as the protein net charge, but the magnitudes were much smaller. The results are consistent with a model in which the electrical charges are approximately uniformly distributed over the HiPIP surface and in which there are both short- and long-range electrostatic interactions. An electrostatic field calculation for Chromatium vinosum HiPIP is consistent with this. The presumed site of electron transfer includes that region of the protein surface to which the iron-sulfur cluster is nearest and appears to be relatively hydrophobic. The principal short-range electrostatic interaction would involve the negative charge on the iron-sulfur cluster. For some net negatively charged proteins, this effect is magnified, and for net positively charged HiPIP's, it is counterbalanced. The rate constants extrapolated to infinite ionic strength can be correlated with redox potential differences between the reactants, as has previously been shown for cytochrome-flavin semiquinone reactions. Both electrostatic and redox potential effects are magnified for the flavodoxin semiquinone as compared to the FMN semiquinone-HiPIP reactions. This was also observed previously for the flavin semiquinone-cytochrome reactions.(ABSTRACT TRUNCATED AT 250 WORDS)     

Energetics of charge-charge interactions in proteins

Electrostatic interactions between pairs of atoms in proteins are calculated with a model based on the linearized Poisson-Boltzmann equation. The equation is solved accurately by a method that takes into account the detailed shape of the protein. This paper presents applications to several systems. Experimental data for the interaction of ionized residues with an active site histidine in subtilisin BPN' allow the model to be tested, using various assumptions for the electrical properties of the protein and solvent. The electrostatic stabilization of the active site thiolate of rhodanese is analyzed, with attention to the influence of alpha-helices. Finally, relationships between electrostatic potential and charge-charge distance are reported for large and small globular proteins. The above results are compared with those of simpler electrostatic models, including Coulomb's law with both a distance-dependent dielectric constant (epsilon = R) and a fixed dielectric constant (epsilon = 2), and Tanford-Kirkwood theory. The primary conclusions are as follows: 1) The Poisson-Boltzmann model agrees with the subtilisin data over a range of ionic strengths; 2) two alpha-helices generate a large potential in the active site of rhodanese; 3) epsilon = R overestimates weak electrostatic interactions but yields relatively good results for strong ones; 4) Tanford-Kirkwood theory is a useful approximation to detailed solutions of the linearized Poisson-Boltzmann equation in globular proteins; and 5) the modified Tanford-Kirkwood theory over-screens the measured electrostatic interactions in subtilisin.     

A new method for computing the macromolecular electric potential

A general methodology is developed for the rigorous computation of the electrostatic potential for a protein of arbitrary shape, assuming the presence of linear dielectric media. The theory proceeds by considering the distribution of induced polarization charge at the dielectric interface, rather than by attempting a direct solution of Poisson's equation (as in the finite-difference approach of Warwicker & Watson). The method is applied to a study of two-dimensional model proteins, where it is shown that the presence of a cleft is associated with a region of relatively high potential in the solvent medium. The results of a preliminary calculation in three dimensions for the protein lysozyme are also discussed; again, a region of enhanced potential is observed near the cleft at the active site. Our computational evidence supports the suggestion of Warwicker & Watson that clefts are associated with important electrostatic effects.     

Kinetic, equilibrium, and spectroscopic studies on dealkylation ("aging") of alkyl organophosphonyl acetylcholinesterase. Electrostatic control of enzyme topography

The mechanism of dealkylation ("aging") of branched-alkyl organophosphonyl conjugates of acetylcholinesterase and the consequence of this reaction on enzyme conformation were examined by employing kinetic, equilibrium, and spectroscopic techniques. Aging of cycloheptyl methylphosphono-acetylcholinesterase proceeded as a unimolecular reaction in which the enzyme became refractory to oxime reactivation and was accelerated with increases in temperature and decreases in pH and ionic strength of the medium. While aging occurred in a manner invariant with the nature of the salt in buffers containing Na+, K+, Rb+, Cs+, Cl-, CH3COO-, SO2-(4), and PO3-(4), the influence of ionic strength on aging was opposite to that predicted for a mechanism requiring charge separation during formation of the polar transition state. Examination of the equilibrium enzyme conformation with decidium, a fluorescent active center-selective ligand, revealed marked alterations in ligand association and a greater ionic strength dependence for binding after aging. The explanation for this behavior focuses on the high net negative surface charge of the enzyme and proposes that acetylcholinesterase topography is governed by the strength of electrostatic interactions between charged, contiguous, mobile protein regions within the subunit. As such, these studies reveal a reciprocal relationship between acetylcholinesterase topography, surface charge, and ionic strength of the medium.     

The electrostatic potential of B-DNA

Electrostatic potentials around DNA are obtained by solving the nonlinear Poisson-Boltzmann (PB) equation. The detailed charge distribution of the DNA and the different polarizabilities of the macromolecule and solvent are included explicitly in the calculations. The PB equation is solved using extensions of a finite difference approach applied previously to proteins. Electrical potentials and ion concentrations are compared to those obtained with simpler models. It is found that the shape of the dielectric boundary between the macromolecule and solvent has significant effects on the calculated potentials near the surface, particularly in the grooves. Sequence-specific patterns are found, the most surprising result being the existence of positive regions of potential near the bases in both the major and minor grooves. The effect of solvent and ionic atmosphere screening of phosphate-phosphate repulsions is studied, and an effective dielectric function, appropriate for molecular mechanics simulations, is derived.     

Asymmetric electrostatic effects on the gating of rat brain sodium channels in planar lipid membranes.

The effects of ionic strength (10-1,000 mM) on the gating of batrachotoxin-activated rat brain sodium channels were studied in neutral and in negatively charged lipid bilayers. In neutral bilayers, increasing the ionic strength of the extracellular solution, shifted the voltage dependence of the open probability (gating curve) of the sodium channel to more positive membrane potentials. On the other hand, increasing the intracellular ionic strength shifted the gating curve to more negative membrane potentials. Ionic strength shifted the voltage dependence of both opening and closing rate constants of the channel in analogous ways to its effects on gating curves. The voltage sensitivities of the rate constants were not affected by ionic strength. The effects of ionic strength on the gating of sodium channels reconstituted in negatively charged bilayers were qualitatively the same as in neutral bilayers. However, important quantitative differences were noticed: in low ionic strength conditions (10-150 mM), the presence of negative charges on the membrane surface induced an extra voltage shift on the gating curve of sodium channels in relation to neutral bilayers. It is concluded that: (a) asymmetric negative surface charge densities in the extracellular (1e-/533A2) and intracellular (1e-/1,231A2) sides of the sodium channel could explain the voltage shifts caused by ionic strength on the gating curve of the channel in neutral bilayers. These surface charges create negative electric fields in both the extracellular and intracellular sides of the channel. Said electric fields interfere with gating charge movements that occur during the opening and closing of sodium channels; (b) the voltage shifts caused by ionic strength on the gating curve of sodium channels can be accounted by voltage shifts in both the opening and closing rate constants; (c) net negative surface charges on the channel's molecule do not affect the intrinsic gating properties of sodium channels but are essential in determining the relative position of the channel's gating curve; (d) provided the ionic strength is below 150 mM, the gating machinery of the sodium channel molecule is able to sense the electric field created by surface changes on the lipid membrane. I propose that during the opening and closing of sodium channels, the gating charges involved in this process are asymmetrically displaced in relation to the plane of the bilayer. Simple electrostatic calculations suggest that gating charge movements are influenced by membrane electrostatic potentials at distances of 48 and 28 A away from the plane of the membrane in the extracellular sides of the channel, respectively.     

Calculation of the electric potential in the active site cleft due to alpha-helix dipoles

A macroscopic dielectric model has been used to set up the electrostatic equation for the protein-solvent system. A numerical method of solution has been applied, enabling calculation of the electric potential outside a protein due to the charges within the protein. The glycolytic enzyme phosphoglycerate mutase, which is an α/β protein binding negatively charged substrates, has been studied. Modelling the helix dipoles with positive and negative charges shows that the α-helical structure could stabilize negatively charged substrates in the active site cleft of an enzyme with an energy of a few kT.     

Electrical potential of transfer RNAs: codon-anticodon recognition

Calculations of the electrostatic potentials were made around yeast elongator phenylalanine, aspartate tRNAs, and yeast initiator methionine tRNA in aqueous solution at physiological ionic strength. The calculations were carried out with a finite difference algorithm for solving the nonlinear Poisson-Boltzmann equation that incorporates the screening effects of the electrolyte, the exclusion of ions by the molecule, the molecular shape, and the different polarizabilities of the solvent and the tRNA. The initiator tRNA is surrounded by uniformly spaced contours of negative potential. The elongator tRNAs are also surrounded by a similar contour pattern except in the anticodon region where there is a pronounced "hole" in the potential surface. This hole is caused by an invagination of the potential contours, which also results in an increase in the local field strength. The effect of this hole is that the anticodon region in the elongator tRNAs is the least negative, or conversely the most positive, region of the molecule. This hole, which is not found when simple Coulombic potentials are used, is due both to the structure of the elongator tRNA anticodon loops and to the different polarizabilities of the solvent and tRNA. The existence of the potential hole in elongator tRNAs may account in part for their ability to associate with other negatively charged macromolecules, in particular mRNA. Moreover, it suggests that the anticodon loop of elongator tRNAs is the energetically most favorable point of approach of mRNA to tRNA.     

Role of charged residues in the catalytic mechanism of hepatitis C virus NS3 protease: electrostatic precollision guidance and transition-state stabilization

Maturational cleavage of the hepatitis C virus polyprotein involves the viral chymotrypsin-like serine protease NS3. The substrate binding site of this enzyme is unusually flat and featureless. We here show that NS3 has a highly asymmetric charge distribution that is characterized by strong positive potentials in the vicinity of its active site and in the S5/S6 region. Using electrostatic potential calculations, we identified determinants of this positive potential, and the role of six different residues was explored by site-directed mutagenesis. Mutation of residues in the vicinity of the active site led to changes in k(cat) values of a peptide substrate indicating that basic amino acids play a role in the stabilization of the transition state. Charge neutralization in the S5/S6 region increased the K(m) values of peptide substrates in a manner that depended on the presence of negatively charged residues in the P5 and P6 positions. K(i) values of hexapeptide acids spanning P6-P1 (product inhibitors) were affected by charge neutralization in both the active site region and the S5/S6 region. Pre-steady-state kinetic data showed that the electrostatic surface potential is used by this enzyme to enhance collision rates between peptidic ligands and the active site. Calculations of the interaction energies of protease-substrate or protease-inhibitor complexes showed that electrostatic interaction energies oppose the formation of a tightly bound complex due to an unfavorable change in the desolvation energy. We propose that desolvation costs are minimized by avoiding the formation of individual ion pair interactions through the use of clusters of positively charged residues in the generation of local electrostatic potentials.     

Enzyme/non-enzyme discrimination and prediction of enzyme active site location using charge-based methods

Calculations of charge interactions complement analysis of a characterised active site, rationalising pH-dependence of activity and transition state stabilisation. Prediction of active site location through large DeltapK(a)s or electrostatic strain is relevant for structural genomics. We report a study of ionisable groups in a set of 20 enzymes, finding that false positives obscure predictive potential. In a larger set of 156 enzymes, peaks in solvent-space electrostatic properties are calculated. Both electric field and potential match well to active site location. The best correlation is found with electrostatic potential calculated from uniform charge density over enzyme volume, rather than from assignment of a standard atom-specific charge set. Studying a shell around each molecule, for 77% of enzymes the potential peak is within that 5% of the shell closest to the active site centre, and 86% within 10%. Active site identification by largest cleft, also with projection onto a shell, gives 58% of enzymes for which the centre of the largest cleft lies within 5% of the active site, and 70% within 10%. Dielectric boundary conditions emphasise clefts in the uniform charge density method, which is suited to recognition of binding pockets embedded within larger clefts. The variation of peak potential with distance from active site, and comparison between enzyme and non-enzyme sets, gives an optimal threshold distinguishing enzyme from non-enzyme. We find that 87% of the enzyme set exceeds the threshold as compared to 29% of the non-enzyme set. Enzyme/non-enzyme homologues, "structural genomics" annotated proteins and catalytic/non-catalytic RNAs are studied in this context.