Search PubMedSearch

SEARCH · Search PubMed

Results for “electrostatics”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

Electrostatic effects in a large enzyme complex: subunit interactions and electrostatic potential field of the icosahedral beta 60 capsid of heavy riboflavin synthase.

The role of the electrostatic interactions in the stability of the icosahedral beta 60 capsid of heavy riboflavin synthase from Bacillus subtilis has been investigated using an approach based on the theory of Kirkwood and Tanford. The pH dependence of the electrostatic subunit interactions agrees well with experimental data. The electrostatic subunit interaction energy has a pronounced minimum at pH 8.2 for both the ligated and ligand-free capsid. The latter is characterized by a reduction of the magnitude and the pH range of the electrostatic attraction. It is found that only 8 charged groups, which form one cluster and two ion pairs, provide a significant contribution to the capsid stability. The analysis has shown that the aggregation/disaggregation equilibrium seems to be regulated by electrostatic interactions between beta-subunits forming dimers, which connect the relatively stable pentamers in the beta-60 capsid. The release of the ligand causes a reduction of the electrostatic attraction of the dimers, which may induce disaggregation of the capsid. The electrostatic potential field due to the titratable groups and alpha-helix macrodipoles has been calculated on the basis of the Coulomb relation. Two different values of the dielectric constant have been used for the protein and the surrounding solvent, respectively. The electrostatic potential shows a radially polar distribution with a positive pole at the inner capsid wall and a negative pole outside the capsid. An interesting feature of the electrostatic field is the formation of positive potential "channels" that coincide with the channels constituted by the pentameric and trimeric beta-subunit aggregates. It is supposed that the electrostatic potential field plays a role in enzyme-substrate recognition.

Bacillus subtilis

Theoretical studies of the mechanism of the action of the neurohypophyseal hormones. I. Molecular electrostatic potential (MEP) and molecular electrostatic field (MEF) maps of some vasopressin analogues.

Continuing our theoretical studies of the oxytocin and vasopressin analogues, we have analysed the molecular electrostatic potential (MEP) and the norm of the molecular electrostatic field (MEF) of [1-beta-mercaptopropionic acid]-arginine-vasopressin ([ Mpa1]-AVP), [1-(beta-mercapto-beta,beta-cyclopentamethylene)propionic acid]-arginine-vasopressin ([Cpp']-AVP), and [1-thiosalicylic acid]-arginine-vasopressin ([Ths']-AVP) whose low-energy conformations were calculated in our previous work. These compounds are known from experiment to exhibit different biological activity. The scalar fields mentioned determine the energy of interaction with either charged (MEP) or polar (MEF) species, the energy being in the second case either optimal or Boltzmann-averaged over all the possible orientations of the dipole moment versus the electrostatic field. The electrostatic interactions slowly vanish with distance and can therefore be considered to be the factor determining the molecular shape at greater distances, which can help in both predicting the interactions with the receptor at the stage of remote recognition and in finding the preferred directions of solvation by a polar solvent. In the analysis of the fields three techniques have been used: (i) the construction of maps in certain planes; (ii) the construction of maps on spheres centered in the charge center of the molecule under study and of poles chosen according to the main axes of the quadrupole moment; and (iii) the construction of surfaces corresponding to a given value of potential. The results obtained show that the shapes of both MEP and MEF are similar in the case of [Mpa1]-AVP and [Cpp1]-AVP (biologically active), while some differences emerge when comparing these compounds with [Ths1]-AVP (inactive). It has also been found that both MEP and MEF depend even more strongly on conformation.

Electrochemistry

On the multiple-minima problem in the conformational analysis of polypeptides. V. Application of the self-consistent electrostatic field and the electrostatically driven Monte Carlo methods to bovine pancreatic trypsin inhibitor.

In connection with the accompanying paper to test various models for the hydration of polypeptides, we have explored a limited portion of the conformational energy hyperspace of the small protein bovine pancreatic trypsin inhibitor (BPTI) with the aid of two search methods developed in this laboratory. A series of low-energy conformations was obtained as a result of this study. These conformations constitute a set of local minima in the conformational energy space of the molecule as described by the ECEPP/2 (Empirical Conformational Energy Program for Peptides) potential energy function, without the inclusion of hydration. Five different initial conformations were used in this exploration: the first corresponds to an energy-refined structure based on the crystallographic coordinates (4PTI) provided by Deisenhofer and Steigemann and reported previously by Meirovitch and Scheraga. The remaining four initial conformations were obtained by using a Variable-Target-Function procedure, applied to the experimental Cartesian coordinates (5PTI) reported by Wlodawer et al. The self-consistent electrostatic field (SCEF) and the electrostatically driven Monte Carlo (EDMC) methods were used to search the conformational space. The SCEF and EDMC methodologies assume that a polypeptide or protein molecule is driven toward the native structure mainly by the action of the electrostatic interactions. Application of these methodologies led to a set of conformations (up to 50 kcal/mol lower than the starting ones) with ECEPP/2 energies lower than any of those that we had previously found. Application of both methods to the initial conformation generated from 4PTI led to a series of low-energy conformations exhibiting similar rms deviations with respect to the experimental data (4PTI) as did the starting conformation. However, statistical analysis of the runs that had started from the conformations generated by using the variable-target-function procedure (and applying the EDMC method) indicated that the rms deviations of the atomic positions of the new low-energy conformations tended to increase as the energy improved, when compared with the X-ray data from which the starting conformations had been generated. The structures with the lowest energies also had radii of gyration smaller than the experimentally observed one. These results indicated a need to include hydration in the potential function, and provided the conformations used in the accompanying paper to test various hydration models.

Animals

Electrostatic interactions in sperm whale myoglobin. Site specificity, roles in structural elements, and external electrostatic potential distributions.

The electrostatic free energy contribution to the stability of sperm whale ferrimyoglobin was evaluated according to the static accessibility modified Tanford-Kirkwood model. The electrostatic free energy contribution of each distinct structural element was divided into one term arising from interactions between it and other elements (interelemental) and another from interactions within the particular element itself (intraelemental). At pH 7 the majority of the terms were found to be stabilizing. The interelemental terms are the dominant ones for most structural elements. The small interelemental terms of the C and D helices are compensated by large intraelemental interactions which stabilize these short helices. Perturbations in pH can be accommodated by the structural elements through a redistribution of stabilizing and destabilizing interactions. The electrostatic potentials calculated at the surface of the protein indicate that the internal compensation of local potentials achieved during folding results in a generally neutral protein-solvent interface save for two distinct areas of nonzero potential. The accessibility of each charged atom to solvent was analyzed in terms of the surface area lost to charged, polar and nonpolar atoms separately. The net solvent accessibility lost parallels closely that lost to nonpolar atoms alone, indicating a specific role for nonpolar atoms in defining dielectric shielding of charged atoms, aside from their participation in the well-known hydrophobic interactions.

Amino Acids

Electrostatic interactions at charged lipid membranes. Electrostatically induced tilt.

The changes in bilayer structure induced by surface charges in the case of an ionizable lipid were studied by X-ray diffraction, Raman spectroscopy, and film-balance measurements. With increasing surface charge in the ordered phase, the X-ray results show a decrease in bilayer thickness, whereas the hydrocarbon chain packing stays essentially constant, the Raman data signify that the internal chain ordering does not change, and the monolayer studies show a lateral expansion of the bilayer. These results are interpreted in terms of a tilt of the chains caused by the surface charges on the polar heads. The tilt angle between the direction of the chains and the bilayer normal is obtained by a detailed theoretical evaluation. The tilt allows for a better understanding of the electrostatically induced shift of the phase transition temperature and of the shift induced by the binding of water in the case of lecithin in contrast ethanolamine.

Hydrogen-Ion Concentration

Electrostatic effects and the dynamics of enzyme reactions at the surface of plant cells. 2. The role of pectin methyl esterase in the modulation of electrostatic effects in soybean cell walls.

The pectin methyl esterase from soybean cell walls has been isolated and purified to homogeneity. It is a protein with a relative molecular mass close to 33 000. The enzyme is maximally active at a pH close to 8 and its pH dependence may be explained by a classical Dixon model, where the two interconvertible enzyme ionization states coexist. The outflux of protons from cell walls, upon raising the ionic strength, may be taken as an indirect estimate of the fixed charge density. If the cell-wall fragments are pre-incubated at pH values between 5 and 9, the outflux of protons rises with the pH of pre-incubation. This implies, as postulated from the theory developed in the preceding paper, that alkaline pH favours the activity of pectin methyl esterase and that this enzyme effectively generates the fixed negative charges of the cell wall. Therefore the pectin methyl esterase reaction builds up the Donnan potential, delta psi, at the cell surface. The cell-wall charge density, estimated from the proton outflux, as well as from the titration of methyl groups on the cell wall, reaches a maximum between the third and the fourth day of growth. While the cell-wall volume increases and reaches a plateau, the fixed charge density increases at first and then declines. This is understandable if one assumes that the building up of a high charge density is a co-operative phenomenon and that the local pH inside the wall rises during cell growth. When both the cell-wall volume and the charge density increase together, this suggests that the local pH inside the wall lies within the critical pH range associated with the steep response of the system. When the cell-wall volume increases together with a decrease of the fixed charge density, the local pH should have dropped below this critical pH range. Under these conditions the pectin methyl esterase remains inactive, or poorly active. As the number of fixed negative charges increases, calcium becomes tightly bound to cell walls. This binding is so tight that the net charge density is minimum when the calcium concentration is maximum. The experimental results, presented above, offer experimental support to two important ideas discussed in the preceding paper, namely that pectin methyl esterase reaction builds up the Donnan potential at the cell surface, and that this response may be co-operative with respect to pH.

Calcium

Electrostatic effects and the dynamics of enzyme reactions at the surface of plant cells. 3. Interplay between limited cell-wall autolysis, pectin methyl esterase activity and electrostatic effects in soybean cell walls.

Soybean cell walls display a process of autolysis which results in the release of reducing sugars from the walls. Loosening and autolysis of cell wall are involved in the cell-wall growth process, for autolysis is maximum during both cell extension and cell-wall synthesis. Autolysis goes to completion within about 50 h and is an enzymatic process that results from the activity of cell wall exo- and endo-glycosyltransferases. The optimum pH of autolysis is about 5. Increasing the ionic strength of the bulk phase where cell-wall fragments are suspended, results in a shift of the pH profile towards low pH. This is consistent with the view that at 'low' ionic strength, the local pH in the cell wall is lower than in the bulk phase. One of the main ideas of the model proposed in a preceding paper, is that pectin methyl esterase reaction, by building up a high fixed charge density, results in proton attraction in the wall. Low pH must then activate the wall loosening enzymes involved in autolysis and cell growth. This view may be directly confirmed experimentally. The pH of a cell-wall suspension, initially equal to 5, was brought to 8 for 20 min, then back to 5. Under these conditions, the rate of cell-wall autolysis was enhanced with respect to the rate of autolysis obtained with cell-wall fragments kept at pH 5. The pH response of the multienzyme plant cell-wall system basically relies on opposite pH sensitivities of the two types of enzymes involved in the growth process. Pectin methyl esterase, which generates the cell-wall Donnan potential, is inhibited by protons, whereas the wall-loosening enzymes involved in cell growth are activated by protons.

Autolysis

Role of electrostatic interaction in the stability of the hexamer of constitutive phycocyanin from Fremyella diplosiphon.

The role of electrostatic interactions between the charges carried by the titratable groups in the different aggregates of constitutive phycocyanin from Fremyella diplosiphon has been studied by using a simple theoretical approach based on the modified Tanford-Kirkwood model. The electrostatic potential has also been calculated by means of a numerical solution of the linearized Poisson-Boltzmann equations using the finite-differences technique. The pH dependence of the electrostatic contribution to free energy suggests an electrostatic stabilization of the alpha- and beta-subunits as well as of the (alpha beta)-monomer over a broad pH interval. The charge distributions in the individual alpha- and beta-subunits produce electrostatic complementarity and promote the assembly of the subunits to the (alpha beta)-monomer, as well as of the monomers to the larger trimeric and hexameric aggregates. Trimer-trimer electrostatic interactions exhibit strong pH dependence, predicting an association/dissociation equilibrium with a midpoint at pH 6. The electrostatic trimer-trimer interactions correspond to the steric fit, suggesting that electrostatic interactions may initially help to orient the trimers during aggregation. The distribution of the electrostatic potential of the monomers and of the higher aggregates suggests that it plays an important role also in phycocyanin-linker protein binding.

Computer Simulation

Electrostatic orientation of the electron-transfer complex between plastocyanin and cytochrome c.

To understand the specificity and efficiency of protein-protein interactions promoting electron transfer, we evaluated the role of electrostatic forces in precollision orientation by the development of two new methods, computer graphics alignment of protein electrostatic fields and a systematic orientational search of intermolecular electrostatic energies for two proteins at present separation distances. We applied these methods to the plastocyanin/cytochrome c interaction, which is faster than random collision, but too slow for study by molecular dynamics techniques. Significant electrostatic potentials were concentrated on one-fourth (969 A2) of the plastocyanin surface, with the greatest negative potential centered on the Tyr-83 hydroxyl within the acidic patch, and on one-eighth (632 A2) of the cytochrome c surface, with the greatest positive potential centered near the exposed heme edge. Coherent electrostatic fields occurred only over these regions, suggesting that local, rather than global, charge complementarity controls productive recognition. The three energetically favored families of pre-collision orientations all directed the positive region surrounding the heme edge of cytochrome c toward the acidic patch of plastocyanin but differed in heme plane orientation. Analysis of electrostatic fields, electrostatic energies of precollision orientations with 12 and 6 A separation distances, and surface topographies suggested that the favored orientations should converge to productive complexes promoting a single electron-transfer pathway from the cytochrome c heme edge to Tyr-83 of plastocyanin. Direct interactions of the exposed Cu ligand in plastocyanin with the cytochrome c heme edge are not unfavorable sterically or electrostatically but should occur no faster than randomly, indicating that this is not the primary pathway for electron transfer.

Amino Acid Sequence

A practical procedure for the determination of electrostatic charges of large molecules.

A practical procedure for the precise determination of electrostatic charges, which are evaluated by fitting the rigorous quantum mechanical molecular electrostatic potential to a monopole-monopole expression, is presented. The proposal of this procedure arises from the study of the minimum requirements necessary to obtain reliable electrostatic charges. Such a study is focused on: (i) the dependence of the electrostatic charges on the set of points where the quantum mechanical and the monopole-monopole molecular electrostatic potentials are fitted; thus, both the influence of the number of points and their distribution in layers located out of the van der Waals radii of the atoms are examined, and (ii) the reliability of the use of fractional models for the evaluation of electrostatic charges of large molecules. Results point out that the optimum number of points is defined by a density of points ranging from 0.45 to 0.60 points/A2 when four layers (separated by 0.2 A) are considered. Nevertheless, the use of only two layers (separated by 0.4 A) for large molecules is recommended, thus enabling one to obtain reliable charges at a reduced computational cost. Moreover, results justify the use of fractional models for the determination of electrostatic charges of extremely large molecules, even when aromatic structures exist.

Adenosine

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.

Amino Acid Sequence

Ligand atom partial charges assignment for complementary electrostatic potentials.

The design of molecules to fit into the active site of receptors is a rapidly developing area of pharmacology and medicinal chemistry. A good ligand needs a suitable geometry and also appropriate electrostatic properties. The electrostatic properties of the ligand should complement those of the receptor. We present a method for the assignment of atom-centred point charges for a ligand, based on the electrostatic potential of the receptor. These point charges are chosen to give the best possible complementarity to the receptor electrostatic potential over the van der Waals surface of the ligand. We demonstrate that point charges can be chosen to give good electrostatic complementarity, and suggest that a molecule with similar electrostatic properties should bind well to the receptor.

Binding Sites

On the representation of electrostatic fields around ab initio charge distributions.

We compare two methods (Mulliken charges and a distributed multipole analysis, DMA) of representing an ab initio charge distribution for calculating the electrostatic field and potential outside the molecule, using pyrimidine and the RNA base uracil as examples. This is done using a 3-D graphical display of the electrostatic fields, which, when used with real-time rotation, zooming and clipping, has many advantages for qualitatively assessing the electrostatic interactions of a molecule. The errors involved in using Mulliken point charges may be of similar magnitude to the total electrostatic field in regions which are important in recognition processes. The DMA representation automatically includes the anisotropic electrostatic effects of non-spherical features in the charge distribution of each atom, and yet the displayed electrostatic fields around the atoms which have lone-pair density do not show marked anisotropy.

Chemical Phenomena

Visualization of molecular flexibility and its effects on electrostatic recognition.

To study the effect of protein flexibility on electrostatic recognition, we have devised two novel computer graphic representations of the changes in the electrostatic field of a protein resulting from its internal motions. The atomic structure of Cu, Zn superoxide dismutase was minimized, and the 200 lowest frequency normal modes of the enzyme were determined. Individual and combined normal-mode vibrations were visualized interactively with the program Flex. Normal-mode motions are fast enough (approximately 10(-11) s cycle-1) to evade solvent damping, thus allowing long-range electrostatic interactions to dominate. The changing electrostatic environment of the protein was examined by animating precalculated frames of electrostatic field vectors with GRAMPS. With Vu, changes in electrostatic potential were displayed as variations in the color-coding of dots lying on a consensus surface that maintains the protein's shape. The consensus surface was calculated with the program Sphinx, and was derived from spherical harmonic approximations of expanded molecular surfaces. The ability to view the effects of molecular motions interactively should be useful in understanding the relationships of protein structure to function.

Computer Graphics

Is the activity-linked electrostatic gradient of bovine Cu, Zn superoxide dismutases conserved in homologous enzymes irrespective of the number and distribution of charges?

Electrostatic potential calculations have been performed on three different Cu, Zn superoxide dismutases (superoxide: superoxide oxidoreductase, EC 1.15. 1.1), in order to evaluate the degree of conservation of the pattern of electrostatic interactions between O2- and the active site recently pointed out in bovine Cu Zn SOD. The three Cu, Zn SODs that have been selected for this study, namely the bovine, ovine, and porcine enzymes, are highly homologous as to reasonably assume identical three-dimensional structure but display large differences in their net charge, as shown by their pI's, which span over a wide pH range: 8.0 (sheep), 6.5 (pig), 5.2 (ox). Despite such a large difference in the net protein charge and in the spatial arrangement of electrostatic charges, electrostatic potential calculations show that the electrostatic channel directing the negatively charged substrate toward the positive catalytic site is strictly preserved with the same features for the three proteins. This suggests that the electrostatic funnel for conducting small anions into the active site is a highly conservative property in the evolution of Cu, Zn SOD.

Animals

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)

Benzoquinones