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 721 records · Page 40Linked to original sources

Comparison of alpha-helix stability in peptides having a negatively or positively charged residue block attached either to the N- or C-terminus of an alpha-helix: the electrostatic contribution and anisotropic stability of the alpha-helix.

An estimation of the thermodynamic effects of a charged random coil, which is attached either to the N- or C-terminus of polyalanine, upon alpha-helix stability is attempted. A temperature-induced helix-coil transition of Ala20Lys20Phe and Lys20Ala20Phe was studied under various conditions of salt concentration and pH. By combining the results with previous ones for Ala20Glu20Phe and Glu20Ala20Phe, which have opposite electric charges to the present system [S. Ihara et al. (1982) Biopolymers 21, 131-145], the free energy of the coil to helix transition of the polyalanine block could be separated into two terms--one term for the electrostatic interaction of electric charges in the random-coil block with the alpha-helix dipole, and a second term for the intrinsic stability of the helix. The first term indicates the significance of the helix dipole-charge interactions, which affects the helix stability depending on the attaching side of the charged block and on the sign of the charges. This clearly shows the anisotropic stability of the alpha-helix. Furthermore, analysis of the dependence of these thermodynamic quantities on salt concentrations showed, assuming that the effect of the attached electric charges was symmetric (in other words, the absolute values of the electrostatic interaction terms were independent of the sign of electric charges), that the intrinsic stability of the alpha-helix was dependent on which side of the helix was attached to the random coil: a random coil attached to the N-terminus of the alpha-helix had little effect while that attached to a C-terminal significantly destabilized the helix.

Circular Dichroism↗

The electrostatic contribution to DNA base-stacking interactions.

Base-stacking and phosphate-phosphate interactions in B-DNA are studied using the finite difference Poisson-Boltzmann equation. Interaction energies and dielectric constants are calculated and compared to the predictions of simple dielectric models. No extant simple dielectric model adequately describes phosphate-phosphate interactions. Electrostatic effects contribute negligibly to the sequence and conformational dependence of base-stacking interactions. Electrostatic base-stacking interactions can be adequately modeled using the Hingerty screening function. The repulsive and dispersive Lennard-Jones interactions dominate the dependence of the stacking interactions on roll, tilt, twist, and propellor. The Lennard-Jones stacking energy in ideal B-DNA is found to be essentially independent of sequence.

Base Sequence↗

The electrostatic potential for the phosphodiester group determined from X-ray diffraction.

The charge density distribution in the crystal structure of ammonium dimethylphosphate at 123 K has been determined from x-ray diffraction data (MoK alpha) using 8437 reflections with sin theta/lambda less than 1.33 A-1 [NH4+.(CH3)2PO4-, M(r) = 143.08, monoclinic, P2(1)/c, a = 10.007(1), b = 6.926(1), c = 9.599(2) A, beta = 105.40(1) degrees, V = 641.4(3) A3, Z = 4, F000 = 304, Dx = 1.4815 g.cm-3, mu = 3.726 cm-1]. Least-squares structure refinement assuming Stewart's rigid pseudoatom model (variables including Slater-type radial exponents and electron populations for multipole terms extending to octapoles for C, N, O, and P, and dipoles for H) gave R(F2) = 0.039 for all reflections. The dimethylphosphate anion is in the gauche-gauche conformation and has approximate twofold symmetry. One phosphoryl O atom forms three hydrogen bonds and the other forms one. Neither of the ester O atoms is hydrogen bonded. For the dimethylphosphate anion isolated from the crystal structure, a map of the electrostatic potential obtained using the pseudoatom charge parameters shows that the phosphoryl O atoms are considerably more electronegative than the ester O atoms. The electrostatic potential distribution obtained in this way has been fitted by least squares to a system of atom-centered point charges. The potential calculated from these point charges agrees with the experimental result. It also agrees reasonably well with potentials obtained from three other systems of point charges that are widely used as part of the semiempirical force field for molecular mechanics and molecular dynamics calculations involving nucleic acids.

Electric Conductivity↗

Modeling of long-range electrostatic interactions in DNA.

We used a Monte Carlo approach to compute the statistical properties of closed DNA chains with different descriptions of the electrostatic interactions. We compared these computed results with experimentally measured knotting probabilities, which are very sensitive to intersegment interactions. The calculated results based on the Debye-Hückel approximation of the solution of the Poisson-Boltzmann equation agreed very well with the published experimental data, while potential based on counterion condensation theory was clearly less satisfactory. We then compared the simpler hard-core approximation of electrostatic interactions to the Debye-Hückel potential. The hard-core approximation, specified in terms of a DNA effective diameter, gives the same conformational properties of random coils as the Debye-Hückel approximation. We found clear but relatively small differences between the two potentials for the conformational properties of supercoiled DNA.

DNA↗

Biomass/adsorbent electrostatic interactions in expanded bed adsorption: a zeta potential study.

Expanded bed adsorption is an integrated technology that allows the introduction of particle-containing feedstock without the risk of blocking the bed. The biomass particles contained in the feedstock have to be treated as an integral part of the process and potential interactions between suspended biomass and the adsorbent must be excluded during process design. Because the electrostatic forces dominate the interactions between the biomass and adsorbent, the zeta potential has been studied as a tool to characterize biomass/adsorbent electrostatic interactions. The zeta potentials of four types of biomass (yeast intact cells, yeast homogenate, Escherichia coli intact cells, and E. coli homogenate) and two types of ion exchanger were measured systematically at varying process conditions. Using the cell transmission index from biomass pulse-response experiments as a parameter, the relations between zeta potential and the biomass/adsorbent interaction were evaluated. Combining the influences from zeta potential of adsorbent (zeta(a)), zeta potential of biomass (zeta(b)), and biomass size (d), parameter (-zeta(a)zeta(b)d) was found to be an appropriate indicator of the biomass/adsorbent interactions in expanded beds under various liquid-phase conditions for different types of biomass. The threshold value of parameter (-zeta(a)zeta(b)d) can be defined as 120 mV(2) microm for cell transmission of >90%, which means that systems with (-zeta(a)zeta(b)d) < 120 may have a considerable probability of forming stable expanded beds in a biomass suspension under the particular experimental conditions.

Adsorption↗

Electrostatic effects on the thermodynamics of protonation of reduced plastocyanin.

The L12E, L12K, Q88E, and Q88K variants of spinach plastocyanin have been electrochemically investigated. The effects of insertion of net charges near the metal site on the thermodynamics of protonation and detachment from the copper(I) ion of the His87 ligand have been evaluated. The mutation-induced changes in transition enthalpy cannot be explained by electrostatic considerations. The existence of enthalpy/entropy (H/S) compensation within the protein series indicates that solvent-reorganization effects control the differences in transition thermodynamics. Once these compensating contributions are factorized out, the resulting modest differences in transition enthalpies turn out to be those that can be expected on purely electrostatic grounds. Therefore, this work shows that the acid transition in cupredoxins involves a reorganization of the H-bonding network within the hydration sphere of the molecule in the proximity of the metal center that dominates the observed transition thermodynamics and masks the differences that are due to protein-based effects.

Amino Acid Substitution↗

Free-energy simulations and experiments reveal long-range electrostatic interactions and substrate-assisted specificity in an aminoacyl-tRNA synthetase.

Specific recognition of their cognate amino acid substrates by the aminoacyl-tRNA synthetase enzymes is essential for the correct translation of the genetic code. For aspartyl-tRNA synthetase (AspRS), electrostatic interactions are expected to play an important role, since its three substrates (aspartate, ATP, tRNA) are all electrically charged. We used molecular-dynamics free-energy simulations and experiments to compare the binding of the substrate Asp and its electrically neutral analogue Asn to AspRS. The preference for Asp is found to be very strong, with good agreement between simulations and experiment. The simulations reveal long-range interactions that electrostatically couple the amino acid ligand, ATP, and its associated Mg2+ cations, a histidine side chain (His448) next to the amino acid ligand and a flexible loop that closes over the active site in response to amino acid binding. Closing this loop brings a negatively charged glutamate into the active site; this causes His448 to recruit a labile proton, which interacts favorably with Asp and accounts for most of the Asp/Asn discrimination. Cobinding of the second substrate, ATP, increases specificity for Asp further and makes the system robust towards removal of His448, which is mutated to a neutral amino acid in many organisms. Thus, AspRS specificity is assisted by a labile proton and a cosubstrate, and ATP acts as a mobile discriminator for specific Asp binding to AspRS. In asparaginyl-tRNA synthetase, a close homologue of AspRS, a few binding-pocket differences modify the charge balance so that asparagine binding predominates.

Adenosine Triphosphate↗

Iron twin-coronet porphyrins as models of myoglobin and hemoglobin: amphibious electrostatic effects of overhanging hydroxyl groups for successful CO/O2 discrimination.

Inspired by the observation of polar interactions between CO and O(2) ligands and the peptide residues at the active site of hemoglobin and myoglobin, we synthesized two kinds of superstructured porphyrins: TCP-IM, which contains a linked imidazole ligand, and TCP-PY, which contains a linked pyridine ligand, and examined the thermodynamic, kinetic, and spectroscopic (UV/Vis, IR, NMR, and resonance Raman) properties of their CO and O(2) complexes. On both sides of each porphyrin plane, bulky binaphthyl bridges form hydrophobic cavities that are suitable for the binding of small molecules. In the proximal site, an imidazole or pyridine residue is covalently fixed and coordinates axially to the central iron atom. In the distal site, two naphtholic hydroxyl groups overhang toward the center above the heme. The CO affinities of TCPs are significantly lower than those of other heme models. In contrast, TCPs have moderate O(2) binding ability. Compared with reported model hemes, the binding selectivity of O(2) over CO in TCP-IM and TCP-PY complexes is greatly improved. The high O(2) selectivity of the TCPs is mainly attributable to a low CO affinity. The comparison of k(on)(CO) values of TCPs with those of unhindered hemes indicates the absence of steric hindrance to the intrinsically linear CO coordination to Fe(II) in TCP-IM and TCP-PY. The abnormally large k(off)(CO) values are responsible for the low CO affinities. In contrast, k(off)(O(2)) of TCP-PY is smaller than those of other pyridine-coordinated model hemes. For the CO adducts of TCPs, unusually low nu(Fe-CO) and unusually high nu(C-O) frequencies are observed. These results can be ascribed to decreased back-bonding from the iron atom to the bound CO. The lone pairs of the oxygen atoms of the hydroxyl groups prevent back-bonding by exertion of a strong negative electrostatic interaction. On the other hand, high nu(Fe-O(2)) frequencies are observed for the O(2) adducts of TCPs. In the resonance Raman (RR) spectrum of oxy-TCP-IM, we observed simultaneous enhancement of the Fe-O(2) and O-O stretching modes. Furthermore, direct evidence for hydrogen bonding between the hydroxyl groups and bound dioxygen was obtained by RR and IR spectroscopy. These spectroscopic data strongly suggest that O(2) and CO binding to TCPs is controlled mainly by the two different electrostatic effects exerted by the overhanging OH groups: destabilization of CO binding by decreasing back-bonding and stabilization of O(2) binding by hydrogen bonding.

Algorithms↗

Silicon addition to hydroxyapatite increases nanoscale electrostatic, van der Waals, and adhesive interactions.

The normal intersurface forces between nanosized probe tips functionalized with COO(-)-terminated alkanethiol self-assembling monolayers and dense, polycrystalline silicon-substituted synthetic hydroxyapatite (SiHA) and phase pure hydroxyapatite (HA) were measured via a nanomechanical technique called chemically specific high-resolution force spectroscopy. A significantly larger van der Waals interaction was observed for the SiHA compared to HA; Hamaker constants (A) were found to be A(SiHA) = 35 +/- 27 zJ and A(HA) = 13 +/- 12 zJ. Using the Derjaguin-Landau-Verwey-Overbeek approximation, which assumes linear additivity of the electrostatic double layer and van der Waals components, and the nonlinear Poisson-Boltzmann surface charge model for electrostatic double-layer forces, the surface charge per unit area, sigma (C/m(2)), was calculated as a function of position for specific nanosized areas within individual grains. SiHA was observed to be more negatively charged than HA with sigma(SiHA) = -0.024 +/- 0.013 C/m(2), two times greater than sigma(HA) = -0.011 +/- 0.006 C/m(2). Additionally, SiHA was found to have increased surface adhesion (0.7 +/- 0.3 nN) compared to HA (0.5 +/- 0.3 nN). The characterization of the nanoscale variations in surface forces of SiHA and HA will enable an improved understanding of the initial stages of bone-biomaterial bonding.

Biocompatible Materials↗

Evaluation of the role of electrostatic residues in human epidermal growth factor by site-directed mutagenesis and chemical modification.

Four residues in the carboxy-terminal domain of human epidermal growth factor (hEGF), glutamate 40, glutamine 43, arginine 45, and aspartate 46 were targeted for site-directed mutagenesis to evaluate their potential role in epidermal growth factor (EGF) receptor-ligand interaction. One or more mutations were generated at each of these sites and the altered recombinant hEGF gene products were purified and evaluated by radioreceptor competition binding assay. Charge-conservative replacement of glutamate 40 with aspartate resulted in a decrease in receptor binding affinity to 30% relative to wild-type hEGF. On the other hand, removal of the electrostatic charge by substitution of glutamate 40 with glutamine or alanine resulted in only a slightly greater decrease in receptor binding to 25% relative receptor affinity. The introduction of a positive charge upon substitution of glutamine 43 with lysine had no effect on receptor binding. The substitution of arginine 45 with lysine also showed no effect on receptor binding, unlike the absolute requirement observed for the arginine side-chain at position 41 [Engler DA, Campion SR, Hauser MR, Cook JS, Niyogi, SK: J Biol Chem 267:2274-2281, 1992]. Subsequent elimination of the positive charge of lysine 45 by reaction with potassium cyanate showed that the electrostatic property of the residue at this site, as well as that at lysine 28 and lysine 48, was not required for receptor-ligand association.(ABSTRACT TRUNCATED AT 250 WORDS)

Amino Acid Sequence↗

Can the calculation of ligand binding free energies be improved with continuum solvent electrostatics and an ideal-gas entropy correction?

The prediction of a ligand binding constant requires generating three-dimensional structures of the complex concerned and reliably scoring these structures. Here, the scoring problem is investigated by examining benzamidine-like inhibitors of trypsin, a system for which errors in the structures are small. Precise and consistent binding free energies for the inhibitors are determined experimentally for this test system. To examine possible improvement of scoring methods, we test the suitability of continuum electrostatics to account for solvation effects and use an ideal-gas entropy correction to account for the changes in the degrees of freedom of the ligand. The small observed root-mean-square deviation of 0.55 kcal/mol of the calculated relative to the experimental values indicates that the essentials of the binding process have been captured. Even though all six ligands make the same salt bridge and H-bonds to the protein, the electrostatic contribution varies among the ligands by as much as 2 kcal/mol. Moreover, although the ligands are rigid and similar in size, the entropic terms also significantly affect the relative binding affinities (by up to 2.7 kcal/mol). The present approach to solvation and entropy may allow the ranking of the ligands to be considerably improved at a cost that makes the method applicable to the optimization of lead compounds or to the screening of small collections of ligands.

Journal Article↗

A computationally inexpensive modification of the point dipole electrostatic polarization model for molecular simulations.

We present an approximation, which allows reduction of computational resources needed to explicitly incorporate electrostatic polarization into molecular simulations utilizing empirical force fields. The proposed method is employed to compute three-body energies of molecular complexes with dipolar electrostatic probes, gas-phase dimerization energies, and pure liquid properties for five systems that are important in biophysical and organic simulations-water, methanol, methylamine, methanethiol, and acetamide. In all the cases, the three-body energies agreed with high level ab initio data within 0.07 kcal/mol, dimerization energies-within 0.43 kcal/mol (except for the special case of the CH(3)SH), and computed heats of vaporization and densities differed from the experimental results by less than 2%. Moreover, because the presented method allows a significant reduction in computational cost, we were able to carry out the liquid-state calculations with Monte Carlo technique. Comparison with the full-scale point dipole method showed that the computational time was reduced by 3.5 to more than 20 times, depending on the system in hand and on the desired level of the full-scale model accuracy, while the difference in energetic results between the full-scale and the presented approximate model was not great in the most cases. Comparison with the nonpolarizable OPLS-AA force field for all the substances involved and with the polarizable POL3 and q90 models for water and methanol, respectively, demonstrates that the presented technique allows reduction of computational cost with no sacrifice of accuracy. We hope that the proposed method will be of benefit to research employing molecular modeling technique in the biophysical and physical organic chemistry areas.

Algorithms↗

Electrostatic potentials and covalent radii.

We begin with a brief overview of the electrostatic potential V(r) as a fundamental determinant of the properties of systems of electrons and nuclei. The minimum of V(r) along the internuclear axis between two bonded atoms is a natural and physically meaningful boundary point, at which the electrostatic forces of the two nuclei upon an element of charge exactly cancel. We propose that the distances from nuclei to V(r) bond minima provide the basis for a well-defined set of covalent radii. Density functional calculations at the B3PW91/6-311+G** level were carried out for 59 molecules to locate the V(r) minima in 95 bonds and use these as the basis for determining single- and multiple-bond covalent radii for eight first- and second-row atoms plus hydrogen. It was found to be unrealistic to assign a single covalent radius to each atom; different values are needed for bonds to first- and second-row atoms, as well as to hydrogen. Using these results, we are able to predict the bond lengths of 33 single and multiple bonds with average errors of less than 0.04 A relative to experimental data.

Journal Article↗

On the transferability of hydration-parametrized continuum electrostatics models to solvated binding calculations.

Using molecular mechanics force field partial atomic charges, we show the nonuniqueness of the parametrization of continuum electrostatics models with respect to solute atomic radii and interior dielectric constant based on hydration (vacuum-to-water transfer) free energy data available for small molecules. Moreover, parameter sets that are optimal and equivalent for hydration free energy calculations lead to large variations of calculated absolute and relative electrostatic binding free energies. Hence, parametrization of solvation effects based on hydration data, although a necessary condition, is not sufficient to guarantee its transferability to the calculation of binding free energies in solution.

Journal Article↗

A charge-scaling implementation of the variational electrostatic projection method.

Two new charge-scaling methods for efficient modeling of the solvated macromolecular environment in hybrid QM/MM calculations of biological reactions are presented. The methods are extensions of the variational electrostatic projection (VEP) method, and allows a subset of atomic charges in the external environment to be adjusted to mimic, in the active dynamical region, the electrostatic potential and field due to the large surrounding macromolecule and solvent. The method has the advantages that it offers improved accuracy, does not require the use of a three-dimensional grid or auxiliary set of fitting points, and requires only minor molecular simulation code modifications. The VEP-cs and VEP-RVM+cs methods are able to attain very high accuracy (relative force errors of 10(-7) or better with appropriate choice of control parameters), and take advantage of a recently introduced set of high-order discretization schemes and Gaussian exponents for boundary element solvation and VEP methods. The methods developed here serve as potentially powerful tools in the arsenal of computational techniques used in multiscale computational modeling problems.

Biopolymers↗

The origin of the two-electron/four-centers C--C bond in pi-TCNE(2)2- dimers: electrostatic or dispersion?

The structure and stability of the pi-TCNE(2)2- dimers in K2TCNE2 aggregates is revisited trying to find if the origin of their two-electron/four-centers C--C bond are the electrostatic K+-TCNE- interactions or the dispersion interactions between the anions. The study is done at the HF, B3LYP, CASSCF (2,2), and MCQDPT/CASSCF (2,2) levels using the 6-31+G(d) basis set. Our results show that the only minima of this aggregate that preserves the pi-TCNE(2)2- structure has the two K+ atoms placed in equatorial positions in between the two TCNE- planes. When the K+ atoms are placed along the D2h axis of the anions the structure is not a minimum. The main energetic component responsible for the stability of these aggregates comes from the cation-anion interactions. However, a proper accounting of the dispersion component (as done in the MCQDPT/CASSCF (2,2) calculations) is needed to make the closed-shell singlet more stable than the open-shell singlet. Thus, the bond results from the combination of the electrostatic and dispersion components, being the first the dominant one. The optimum geometry of the closed-shell singlet is very similar to the experimental one found in crystals.

Journal Article↗

Combination of transfection agents and magnetic resonance contrast agents for cellular imaging: relationship between relaxivities, electrostatic forces, and chemical composition.

The purpose of this study was to investigate the changes in electrostatic and magnetic resonance (MR) properties observed when MR contrast agents (CAs) (Feridex, MION-46L, or G5-dendrimer-DOTA-Gd) are combined with transfection agents (TAs) under various conditions for use as a CA-TA complex basis for cellular labeling and MRI. CAs were incubated with various classes of TAs for 0-48 hr in solutions of varying concentrations and pH values. NMR relaxation rates (1/T(1), 1/T(2)), MRI and zeta potential (ZP) of CA-TA solutions were measured. TAs decreased the 1/T(1) and 1/T(2) of G5-DOTA-Gd, Feridex, and MION-46L by 0-95%. Altering the pH of G5-DOTA-Gd-TA decreased the T(1)-weighted signal intensity (SI) on MRI from 0 to 78%. Measured ZP values for G5-DOTA-Gd, Feridex, and MION-46L were -51, -41, and -2.0 mV, respectively. The TA LV had a negative ZP, while the other TAs had ZPs ranging from +20 to +65 mV. The alteration of the ZP and NMR relaxivities of the MR CAs, Feridex, MION-46L, and G5-DOTA-Gd by TAs has been demonstrated. These results enhance our understanding of the relationship between electrostatic and MR properties.

Cells, Cultured↗

Pulsed electrostatic fields (ETG) to reduce hair loss in women undergoing chemotherapy for breast carcinoma: a pilot study.

AIMS: To determine whether specific pulsed electrostatic fields, or electrotrichogenesis (ETG), could potentially prevent or reduce hair loss in patients undergoing adjuvant cyclophosphamide, methotrexate and 5-fluorouracil (CMF) chemotherapy for breast cancer. METHODS: Thirteen women were followed during their adjuvant ETG and chemotherapy treatment to determine the efficacy of ETG. All patients were treated for 12 min, twice weekly with a pulsed electrostatic field. Quantitative hair loss was measured by photographic assessment, and manual hair count. Quality of life assessment was conducted at the end of the study. RESULTS: Twelve out of 13 participants had good hair retention throughout the chemotherapy period and afterwards. There were no reported side effects attributable to ETG. CONCLUSIONS: This study shows encouraging results in an area where no other appropriate treatment is available Reducing alopecia, secondary to chemotherapy has the potential to increase CMF treatment compliance, enhance patient self-esteem, and improve overall quality of life during this stressful period.

Adult↗