Search PubMed⌕ Search

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 379 records · Page 21Linked to original sources
PubMedDate not supplied

The relative contribution of electrostatic interactions to stabilization of collagen fibrils.

Electrostatic energies of interaction between type I collagen molecules were calculated, using models developed by Timasheff and Hill. These energies, along with a contribution from hydrophobic forces, were then incorporated into an equation due to Flory that described phase equilibria of rod-like polymers. The Flory formalism in turn permitted a calculation of the overall free energy of fibril formation (delta Ff), and an assessment of the relative contribution of electrostatic and hydrophobic forces to delta Ff. Lastly, delta Ff was used in a nucleation-growth model relating halftimes of fibril formation (t1/2) to ionic strength (I) and temperature. Because the theory provided no basis for setting absolute levels of the energetic contributions, five parameters in the model had to be derived from experimental data. Based on the fit of theory to experimental results both for intact and pepsinized collagen, it was found that very low electrostatic energies (about -1 kcal/mole per collagen molecule) were sufficient to explain experimental t1/2 vs I relationships. This energy is equivalent to 1 close charge-pair interaction per molecule and appears to be lower than the energy assignable to hydrophobic interactions.

Chemical Phenomena↗
PubMed1990

On the multiple-minima problem in the conformational analysis of polypeptides. IV. Application of the electrostatically driven Monte Carlo method to the 20-residue membrane-bound portion of melittin.

The conformational space of the membrane-bound portion of melittin has been searched using the electrostatically driven Monte Carlo (EDMC) method with the ECEPP/2 (empirical conformational energy program for peptides) algorithm. The former methodology assumes that a polypeptide or protein molecule is driven toward the native structure by the combined action of electrostatic interactions and stochastic conformational changes associated with thermal movements. The algorithm produces a Monte Carlo search in the conformational hyperspace of the polypeptide using electrostatic predictions and a random sampling technique, combined with local minimization of the energy function, to locate low-energy conformations. As a result of 8 test calculations on the 20-residue membrane-bound portion of melittin, starting from six arbitrary and two completely random conformations, the method was able to locate a very low-energy region of the potential with a well-defined structure for the backbone. In all of the cases under study, the method found a cluster of similar low-energy conformations that agree well with the structure deduced from x-ray diffraction experiments and with one computed earlier by the build-up procedure.

Amino Acid Sequence↗
PubMed2006

Zeta potential as a diagnostic tool to evaluate the biomass electrostatic adhesion during ion-exchange expanded bed application.

Expanded bed adsorption is an integrative technology in downstream processing allowing the direct capture of target proteins from biomass (cells or cell debris) containing feedstocks. Potential adhesion of biomass on the surface of adsorbent, however, may hamper the application of this technique. Since the electrostatic forces dominate the interactions between biomass and adsorbent, the concept of zeta potential was introduced to characterize the biomass/adsorbent electrostatic interactions during expanded bed application. The criterion of zeta potential evaluation proposed in the previous paper (Biotechnol Bioeng, 83(2):149-157, 2003) was verified further with the experimental validation. The zeta potential of intact cells and homogenates of four microorganisms (Escherichia coli, Bacillus subtilis, Pichia pastoris, and S. cerevisiae) were measured under varying pH and salt concentration, and two ion-exchange adsorbents (Streamline DEAE and Streamline QXL) were investigated. The biomass transmission index (BTI) from the biomass pulse response experiments was used as the indicator of biomass adhesion in expanded bed. Combining the influences from zeta potential of adsorbent (zeta(a)), zeta potential of biomass (zeta(b)) and biomass size (d), a good relationship was established between the zeta potential parameter (-zeta(a)zeta(b)d) and BTI for all experimental conditions. The threshold value of parameter (-zeta(a)zeta(b)d) can be defined as 120 mV2 microm for BTI above 0.9. This means that the systems with (-zeta(a)zeta(b)d) < 120 show neglectable electrostatic bio-adhesion, and would have a considerable probability of forming stable expanded beds in a biomass suspension under the particular experimental conditions.

Adsorption↗
PubMed2006

Electronic origin for enhanced nonlinear optical response of complexes from tetraalkylammonium halide and carbon tetrabromide: electrostatic potentials of intermolecular donor-acceptor dyads.

Electronic origin for nonresonant enhancement of nonlinear optical response in the complexes formed from tetraalkylammonium halide and carbon tetrabromide is provided in view of electrostatic potentials of intermolecular donor (halide ion)-acceptor (CBr(4)). The calculated electrostatic potentials of donor-acceptor range from -4.83 to -7.70 kcal mol(-1) and show a decreasing order of [Et(4)Cl(-)Br] > [Et(4)Br(-)Br] congruent with [Et(4)I(-)Br] > [Bu(4)Br(-)Br]. The calculated second-order susceptibilities of solid complexes are in an increasing order of [NEt(4)ClCBr(4)] < [NEt(4)BrCBr(4)] congruent with [NEt(4)ICBr(4)] < [NBu(4)BrCBr(4)C(3)H(6)O]. It has been shown that the donor/acceptor dyads make the exclusive contribution to nonlinear optical response. A large size of halide or tetraalkylammonium ion results in a small electrostatic potential and large nonlinear optical response in these charge-transfer complexes. It indicates that a small supermolecular interaction will create a large nonlinear optical response, and it gives a clue to design the molecular complexes with large non-linear optical susceptibility.

Amines↗
PubMed1991

Electrostatic destabilization of the cytochrome b6f complex in the thylakoid membrane.

Three of the membrane-spanning polypeptides of the chloroplast cytochrome (cyt) b6f complex were sequentially released from the thylakoid membrane, in the order cyt b6, suIV and Rieske iron-sulfur protein, as the pH was increased from 10 to 12, a protocol usually employed to remove peripheral proteins from membranes. The fourth polypeptide of the cyt b6f complex, cyt f, which spans the membrane once, was apparently not released. The pH values for half-release at low ionic strength were approximately 10.7, 11.1 and 11.3 respectively. The separation of the polypeptides of the complex and the sequential release is readily seen at pH 11, where the loss from the membrane of cyt b6, suIV and Fe iron-sulfur center is approximately 90%, 50% and 20%, respectively. the release of cyt b6 from the membrane was reflected by the absence of its characteristic reduced minus oxidized absorbance signal. The pH values at which the release occurred increased as the ionic strength was raised, implying that the release of the b6f polypeptides arises from extrusion due to repulsive electrostatic interactions probably caused by deprotonation of tyrosine and lysine residues. The lipid content of the released polypeptides was very low, consistent with the observation of a non-membranous state. It is proposed that the pH-dependent extrusion requires two electrostatic effects at alkaline pH higher than approximately 10.5: (i) increased electrostatic repulsion between neighbouring polypeptides of the complex, arising from increased net negative charge in the peripheral segments of these polypeptides, which can cause separation of the polypeptides from the complex; and (ii) ionization of residues such as tyrosine in the membrane-spanning alpha-helices, and neutralization of residues such as lysine which can bind to the negative membrane surface.(ABSTRACT TRUNCATED AT 250 WORDS)

Blotting, Western↗
PubMed2002

On the irrelevance of electrostatics for the crystal structures and polymorphism of long even n-alkanes.

It is known that the experimental triclinic crystal structures of even n-alkanes are not well reproduced upon energy minimization with current force fields. The inclusion of electrostatics does not solve this, and, moreover, some charge schemes show unphysical features such as positively charged carbon atoms or charge alternation. The effect of the electrostatics on the energies of the crystal structures of the even n-alkanes, and thereby on their polymorphism, has never been established. A new charge scheme is introduced that yields physically sensible charges without constraints. It will also be shown, however, that electrostatics are relevant neither for the structures of the crystals, nor for their energies.

Alkanes↗
PubMed2004

Electrostatic potentials and fields from density expansions of deformed atoms in molecules.

The exact representation of the molecular density by means of atomic expansions, consisting in spherical harmonics times analytical radial factors, is employed for the calculation of electrostatic potentials, fields, and forces. The resulting procedure is equivalent to an atomic multipolar expansion in the long-range regions, but works with similar efficiency and accuracy in the short-range region, where multipolar expansions are not valid. The performances of this procedure are tested on the calculation of the electrostatic potential contour maps and electrostatic field flux lines of water and nitrobenzene, computed from high-quality molecular electron densities obtained with Slater basis sets.

Journal Article↗
PubMed2005

Nonuniform charge scaling (NUCS): a practical approximation of solvent electrostatic screening in proteins.

In molecular mechanics calculations, electrostatic interactions between chemical groups are usually represented by a Coulomb potential between the partial atomic charges of the groups. In aqueous solution these interactions are modified by the polarizable solvent. Although the electrostatic effects of the polarized solvent on the protein are well described by the Poisson--Boltzmann equation, its numerical solution is computationally expensive for large molecules such as proteins. The procedure of nonuniform charge scaling (NUCS) is a pragmatic approach to implicit solvation that approximates the solvent screening effect by individually scaling the partial charges on the explicit atoms of the macromolecule so as to reproduce electrostatic interaction energies obtained from an initial Poisson--Boltzmann analysis. Once the screening factors have been determined for a protein the scaled charges can be easily used in any molecular mechanics program that implements a Coulomb term. The approach is particularly suitable for minimization-based simulations, such as normal mode analysis, certain conformational reaction path or ligand binding techniques for which bulk solvent cannot be included explicitly, and for combined quantum mechanical/molecular mechanical calculations when the interface to more elaborate continuum solvent models is lacking. The method is illustrated using reaction path calculations of the Tyr 35 ring flip in the bovine pancreatic trypsin inhibitor.

Algorithms↗
PubMed2007

Electrostatic energy in the effective fragment potential method: theory and application to benzene dimer.

Evaluation of the electrostatic energy within the effective fragment potential (EFP) method is presented. The performance of two variants of the distributed multipole analysis (DMA) together with two different models for estimating the charge penetration energies was studied using six homonuclear dimers. The importance of damping the higher order multipole terms, i.e. charge dipole, was also investigated. Damping corrections recover more than 70% of the charge penetration energy in all dimers, whereas higher order damping introduces only minor improvement. Electrostatic energies calculated by the numerical DMA are less accurate than those calculated by the analytic DMA. Analysis of bonding in the benzene dimer shows that EFP with inclusion of the electrostatic damping term performs very well compared to the high-level coupled cluster singles, doubles, and perturbative triples method. The largest error of 0.4 kcal/mol occurs for the sandwich dimer configuration. This error is about half the size of the corresponding error in second order perturbation theory. Thus, EFP in the current implementation is an accurate and computationally inexpensive method for calculating interaction energies in weakly bonded molecular complexes.

Algorithms↗
PubMed1984

Calculating the electrostatic potential of molecular models with separate evaluations by conventional, vector, and array processors.

A simple computational scheme for estimating the electrostatic potential about molecular models of moderate size is given. The large amount of calculations required for the evaluation of the hypersurface lends itself to treatment by high speed, unconventional computing machines. The essence of these calculations lies in Coulombic interactions that are computed between hypothetical proton test probes positioned in a gridded region surrounding the model and the partial electrostatic charges (CNDO/2) of each atom in the model. A specific scientific application is discussed which involves the recognition of amino acids and nucleotide bases. Three different evaluations of the potential hypersurface within the context of this approach were made. The first was performed on a VAX 11/780 which is a general purpose machine widely used in the scientific community; the second was performed using a pipelined Vector Processor, the FPS AP-120B; and the third by a processor array, the ILLIAC-IV. A comparison of the architectures and processing speeds of each class of machines is made. The computing power observed is consistent with the design and purpose of each machine. Also discussed are methods for displaying the vast amount of data that result from such calculations. It is determined that computer graphics offers an effective means for extracting information from large amounts of data. Finally, the scientific value of the calculations are briefly discussed. If caution is applied to interpreting the results, then the electrostatic potential (EP) mappings can be useful in identifying sites of potential chemical interactions.

Algorithms↗
PubMed1992

Cumulative atomic multipole moments complement any atomic charge model to obtain more accurate electrostatic properties.

The quality of several atomic charge models based on different definitions has been analyzed using cumulative atomic multipole moments (CAMM). This formalism can generate higher atomic moments starting from any atomic charges, while preserving the corresponding molecular moments. The atomic charge contribution to the higher molecular moments, as well as to the electrostatic potentials, has been examined for CO and HCN molecules at several different levels of theory. The results clearly show that the electrostatic potential obtained from CAMM expansion is convergent up to R-5 term for all atomic charge models used. This illustrates that higher atomic moments can be used to supplement any atomic charge model to obtain more accurate description of electrostatic properties.

Biopolymers↗
PubMedDate not supplied

The Poisson-Boltzmann equation for biomolecular electrostatics: a tool for structural biology.

Electrostatics plays a fundamental role in virtually all processes involving biomolecules in solution. The Poisson-Boltzmann equation constitutes one of the most fundamental approaches to treat electrostatic effects in solution. The theoretical basis of the Poisson-Boltzmann equation is reviewed and a wide range of applications is presented, including the computation of the electrostatic potential at the solvent-accessible molecular surface, the computation of encounter rates between molecules in solution, the computation of the free energy of association and its salt dependence, the study of pKa shifts and the combination with classical molecular mechanics and dynamics. Theoretical results may be used for rationalizing or predicting experimental results, or for suggesting working hypotheses. An ever-increasing body of successful applications proves that the Poisson-Boltzmann equation is a useful tool for structural biology and complementary to other established experimental and theoretical methodologies.

Algorithms↗
PubMedDate not supplied

Molecular electrostatic potential as a factor of drug-receptor recognition.

When a drug molecule approaches a non-specific acceptor both molecules are in electrostatic fields of equal sign which prevents drug-acceptor complex formation. At the same time, the drug-acceptor system does not achieve the thermodynamic global minimum. Otherwise, when a certain drug interacts with its specific receptor, mutual compensation of their molecular electrostatic potentials takes place. Then separate atoms and groups of the drug molecule can bind to the receptor. We show that the fundamental role of molecular electrostatic potential in the process of drug-receptor recognition consists in fast correction of errors.

Chemical Phenomena↗
PubMed1993

Binding of ferredoxin to ferredoxin:NADP+ oxidoreductase: the role of carboxyl groups, electrostatic surface potential, and molecular dipole moment.

The small, soluble, (2Fe-2S)-containing protein ferredoxin (Fd) mediates electron transfer from the chloroplast photosystem I to ferredoxin: NADP+ oxidoreductase (FNR), a flavoenzyme located on the stromal side of the thylakoid membrane. Ferredoxin and FNR form a 1:1 complex, which is stabilized by electrostatic interactions between acidic residues of Fd and basic residues of FNR. We have used differential chemical modification of Fd to locate aspartic and glutamic acid residues at the intermolecular interface of the Fd:FNR complex (both proteins from spinach). Carboxyl groups of free and FNR-bound Fd were amidated with carbodiimide/2-aminoethane sulfonic acid (taurine). The differential reactivity of carboxyl groups was assessed by double isotope labeling. Residues protected in the Fd:FNR complex were D-26, E-29, E-30, D-34, D-65, and D-66. The protected residues belong to two domains of negative electrostatic surface potential on either side of the iron-sulfur cluster. The negative end of the molecular dipole moment vector of Fd (377 Debye) is close to the iron-sulfur cluster, in the center of the area demarcated by the protected carboxyl groups. The molecular dipole moment and the asymmetric surface potential may help to orient Fd in the reaction with FNR. In support, we find complementary domains of positive electrostatic potential on either side of the FAD redox center of FNR. The results allow a binding model for the Fd:FNR complex to be constructed.

Amino Acid Sequence↗
PubMed1994

Do salt bridges stabilize proteins? A continuum electrostatic analysis.

The electrostatic contribution to the free energy of folding was calculated for 21 salt bridges in 9 protein X-ray crystal structures using a continuum electrostatic approach with the DELPHI computer-program package. The majority (17) were found to be electrostatically destabilizing; the average free energy change, which is analogous to mutation of salt bridging side chains to hydrophobic isosteres, was calculated to be 3.5 kcal/mol. This is fundamentally different from stability measurements using pKa shifts, which effectively measure the strength of a salt bridge relative to 1 or more charged hydrogen bonds. The calculated effect was due to a large, unfavorable desolvation contribution that was not fully compensated by favorable interactions within the salt bridge and between salt-bridge partners and other polar and charged groups in the folded protein. Some of the salt bridges were studied in further detail to determine the effect of the choice of values for atomic radii, internal protein dielectric constant, and ionic strength used in the calculations. Increased ionic strength resulted in little or no change in calculated stability for 3 of 4 salt bridges over a range of 0.1-0.9 M. The results suggest that mutation of salt bridges, particularly those that are buried, to "hydrophobic bridges" (that pack at least as well as wild type) can result in proteins with increased stability. Due to the large penalty for burying uncompensated ionizable groups, salt bridges could help to limit the number of low free energy conformations of a molecule or complex and thus play a role in determining specificity (i.e., the uniqueness of a protein fold or protein-ligand binding geometry).

Crystallography, X-Ray↗
PubMed1994

A negative electrostatic determinant mediates the association between the Escherichia coli trp repressor and its operator DNA.

The electrostatic potential surfaces were characterized for trp repressor models that bind to DNA with sequence specificity, without specificity, and not at all. Comparisons among the surfaces were used to isolate protein surface features likely to be important in DNA binding. Models that differ in protein conformation and tryptophan-analogue binding consistently showed positive potential associated with the protein surfaces that interact with the DNA major groove. However, negative potential is associated with the trp repressor surface that contacts the DNA minor groove. This negative potential is significantly neutralized in the protein conformation that is bound to DNA. Positive potential is also associated with the tryptophan binding-site surface, a consequence of the tryptophan- or tryptophan analogue-induced allosteric change. This protein region is complementary to the strongest negative potential associated with the DNA phosphate backbone and is also present in the isolated protein structure from the protein-DNA complex. The effects of charge-change mutation, pH dependence, and salt dependence on the electrostatic potential surfaces were also examined with regard to their effects on protein-DNA binding constants. A consistent model is formed that defines a role for long-range electrostatics early in the protein-DNA association process and complements previous structural, molecular association, and mutagenesis studies.

Allosteric Site↗
PubMed1995

Protein destabilization by electrostatic repulsions in the two-stranded alpha-helical coiled-coil/leucine zipper.

The destabilizing effect of electrostatic repulsions on protein stability has been studied by using synthetic two-stranded alpha-helical coiled-coils as a model system. The native coiled-coil consists of two identical 35-residue polypeptide chains with a heptad repeat QgVaGbAcLdQeKf and a Cys residue at position 2 to allow formation of an interchain disulfide bridge. This peptide, designed to contain no intrahelical or interhelical electrostatic interactions, forms a stable coiled-coil structure at 20 degrees C in benign medium (50 mM KCl, 25 mM PO4, pH 7) with a [urea]1/2 value of 6.1 M. Four mutant coiled-coils were designed to contain one or two Glu substitutions for Gln per polypeptide chain. The resulting coiled-coils contained potential i to i' + 5 Glu-Glu interchain repulsions (denoted as peptide E2(15,20)), i to i' + 2 Glu-Glu interchain repulsions (denoted E2(20,22)), or no interchain ionic interactions (denoted E2(13,22) and E1(20)). The stabilities of the coiled-coils were determined by measuring the ellipticities at 222 nm as a function of urea or guanidine hydrochloride concentration at 20 degrees C in the presence and absence of an interchain disulfide bridge. At pH 7, in the presence of urea, the stabilities of E2(13,22) and E2(20,22) were identical suggesting that the potential i to i' + 2 interchain Glu-Glu repulsion in the E2(20,22) coiled-coil does not occur. In contrast, the mutant E2(15,20) is substantially less stable than E2(13,22) or E2(15,20) by 0.9 kcal/mol due to the presence of two i to i' + 5 interchain Glu-Glu repulsions, which destabilize the coiled-coil by 0.45 kcal/mol each. At pH 3 the coiled-coils were found to increase in stability as the number of Glu substitutions were increased. This, combined with reversed-phase HPLC results at pH 7 and pH 2, supports the conclusion that the protonated Glu side chains present at low pH are significantly more hydrophobic than Gln side chains which are in turn more hydrophobic than the ionized Glu side chains present at neutral pH. The protonated Glu residues increase the hydrophobicity of the coiled-coil interface leading to higher coiled-coil stability. The guanidine hydrochloride results at pH 7 show similar stabilities between the native and mutant coiled-coils indicating that guanidine hydrochloride masks electrostatic repulsions due to its ionic nature and that Glu and Gln in the e and g positions of the heptad repeat have very similar effects on coiled-coil stability in the presence of GdnHCl.

Amino Acid Sequence↗
PubMed1996

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.

Animals↗