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 901 records · Page 50Linked to original sources
PubMed1993

Brownian dynamics study of the influences of electrostatic interaction and diffusion on protein-protein association kinetics.

A unified model is presented for protein-protein association processes that are under the influences of electrostatic interaction and diffusion (e.g., protein oligomerization, enzyme catalysis, electron and energy transfer). The proteins are modeled as spheres that bear point charges and undergo translational and rotational Brownian motion. Before association can occur the two spheres have to be aligned properly to form a reaction complex via diffusion. The reaction complex can either go on to form the product or it can dissociate into the separate reactants through diffusion. The electrostatic interaction, like diffusion, influences every step except the one that brings the reaction complex into the product. The interaction potential is obtained by extending the Kirkwood-Tanford protein model (Tanford, C., and J. G. Kirkwood. 1957. J. Am. Chem. Soc. 79:5333-5339) to two charge-embedded spheres and solving the consequent equations under a particular basis set. The time-dependent association rate coefficient is then obtained through Brownian dynamics simulations according an algorithm developed earlier (Zhou, H.-X. 1990. J. Phys. Chem. 94:8794-8800). This method is applied to a model system of the cytochrome c and cytochrome c peroxidase association process and the results confirm the experimental dependence of the association rate constant on the solution ionic strength. An important conclusion drawn from this study is that when the product is formed by very specific alignment of the reactants, as is often the case, the effect of the interaction potential is simply to scale the association rate constant by a Boltzmann factor. This explains why mutations in the interface of the reaction complex have strong influences on the association rate constant whereas those away from the interface have minimal effects. It comes about because the former mutations change the interaction potential of the reaction complex significantly and the latter ones do not.

Algorithms↗
PubMed1994

Calculation of electrostatic effects at the amino terminus of an alpha helix.

It is generally believed that the electrostatic field arising from the dipolar charge distribution in alpha helices is important for protein structure and function. We report a calculation of the electrostatic potential and field at the amino terminus of an alpha helix in water, obtained from a finite difference solution to the Poisson-Boltzmann equation. This method takes into account the detailed helix shape and charge distribution, as well as solvent, and generalized ionic strength effects. The calculated potential and field are found to be in good agreement with the experimentally observed helix-induced Stark effect and pKa shifts of a probe at the N-terminus of a stable, monomeric alpha-helical peptide (Lockhart and Kim, 1992, 1993). Ionic screening effects are reproduced at low salt concentrations. Deviations at higher salt concentrations may result from specific ion effects (specific ion-solute and/or ion-solvent interactions). The FDPB method was used to analyze the contributions from each residue, charged side chains, and solvent to the helix potential and field. Backbone contributions come primarily from the first one to two helical turns. Charged side chains contribute to helix-induced pKa shifts for certain probe-peptide combinations, even at relatively large distances from the probe (> 14 A).

Amino Acid Sequence↗
PubMed1995

Localized contact formation by erythrocyte membranes: electrostatic effects.

The topology of the contact seam of human erythrocytes adhered by dextran, an uncharged polymer, has been examined. Particular attention has been paid to the influence of electrostatic intermembrane interactions since their magnitude and range can be accurately estimated. Normal cells formed a continuous seam, whereas erythrocytes with pronase-modified glycocalices formed localized contact points on adhesion in 72 kDa dextran in buffered 145 mM NaCl. The dependence of the inter-contact distance lambda on dextran concentration [D] over the range 2-6% w/v, was given by lambda = C[D]-0.62, where C was a constant. The index of [D] was independent of dextran molecular mass over the range 20 to 450 kDa. The inter-contact distance for pronase-pretreated cells in 6% w/v 72 kDa dextran increased from 0.78 to 1.4 microns as [NaCl] was reduced through the range 145 to 90 mM and the suspending phase was maintained at isotonicity by using sorbitol to replace NaCl. The formation and lateral separation of the contact points are discussed from the perspective of linear interfacial instability theory. The theory allows a quantitative explanation for the experimentally observed dependence of inter-contact distance and of disturbance growth rate on change in electrostatic interaction. The results suggest that the dominant wavelength, determining the inter-contact distance, is established on approaching membranes when the layers of cell surface charge are separated by a perpendicular distance of < 14 nm (bilayer separation of 24 nm).

Cell Adhesion↗
PubMed1995

Conformation and hydrogen ion titration of proteins: a continuum electrostatic model with conformational flexibility.

A new method for including local conformational flexibility in calculations of the hydrogen ion titration of proteins using macroscopic electrostatic models is presented. Intrinsic pKa values and electrostatic interactions between titrating sites are calculated from an ensemble of conformers in which the positions of titrating side chains are systematically varied. The method is applied to the Asp, Glu, and Tyr residues of hen lysozyme. The effects of different minimization and/or sampling protocols for both single-conformer and multi-conformer calculations are studied. For single-conformer calculations it is found that the results are sensitive to the choice of all-hydrogen versus polar-hydrogen-only atomic models and to the minimization protocol chosen. The best overall agreement of single-conformer calculations with experiment is obtained with an all-hydrogen model and either a two-step minimization process or minimization using a high dielectric constant. Multi-conformational calculations give significantly improved agreement with experiment, slightly smaller shifts between model compound pKa values and calculated intrinsic pKa values, and reduced sensitivity of the intrinsic pKa calculations to the initial details of the structure compared to single-conformer calculations. The extent of these improvements depends on the type of minimization used during the generation of conformers, with more extensive minimization giving greater improvements. The ordering of the titrations of the active-site residues, Glu-35 and Asp-52, is particularly sensitive to the minimization and sampling protocols used. The balance of strong site-site interactions in the active site suggests a need for including site-site conformational correlations.

Animals↗
PubMed1995

Computer modeling of electrostatic steering and orientational effects in antibody-antigen association.

Brownian dynamics simulations are performed to investigate the role of long-range electrostatic forces in the association of the monoclonal antibody HyHEL-5 with hen egg lysozyme. The electrostatic field of the antibody is obtained from a solution of the nonlinear Poisson-Boltzmann using the x-ray crystal coordinates of this protein. The lysozyme is represented as an asymmetric dumbell consisting of two spheres of unequal size, an arrangement that allows for the modeling of the orientational requirements for docking. Calculations are done with the wild-type antibody and several point mutants at different ionic strengths. Changes in the charge distribution of the lysozyme are also considered. Results are compared with experiment and a simpler model in which the lysozyme is approximately by a single charged sphere.

Animals↗
PubMed1996

The pore domain of the nicotinic acetylcholine receptor: molecular modeling, pore dimensions, and electrostatics.

The pore domain of the nicotinic acetylcholine receptor has been modeled as a bundle of five kinked M2 helices. Models were generated via molecular dynamics simulations incorporating restraints derived from 9-A resolution cryoelectron microscopy data (Unwin, 1993; 1995), and from mutagenesis data that identify channel-lining side chains. Thus, these models conform to current experimental data but will require revision as higher resolution data become available. Models of the open and closed states of a homopentameric alpha 7 pore are compared. The minimum radius of the closed-state model is less than 2 A; the minimum radius of the open-state models is approximately 6 A. It is suggested that the presence of "bound" water molecules within the pore may reduce the effective minimum radii below these values by up to approximately 3 A. Poisson-Boltzmann calculations are used to obtain a first approximation to the potential energy of a monovalent cation as it moves along the pore axis. The differences in electrostatic potential energy profiles between the open-state models of alpha 7 and of a mutant of alpha 7 are consistent with the experimentally observed change in ion selectivity from cationic to anionic. Models of the open state of the heteropentameric Torpedo nicotinic acetylcholine receptor pore domain are also described. Relatively small differences in pore radius and electrostatic potential energy profiles are seen when the Torpedo and alpha 7 models are compared.

Amino Acid Sequence↗
PubMed1996

Structural origins of redox potentials in Fe-S proteins: electrostatic potentials of crystal structures.

Redox potentials often differ dramatically for homologous proteins that have identical redox centers. For two types of iron-sulfur proteins, the rubredoxins and the high-potential iron-sulfur proteins (HiPIPs), no structural explanations for these differences have been found. We calculated the classical electrostatic potential at the redox site using static crystal structures of four rubredoxins and four HiPIPs to identify important structural determinants of their redox potentials. The contributions from just the backbone and polar side chains are shown to explain major features of the experimental redox potentials. For instance, in the rubredoxins, the presence of Val 44 versus Ala 44 causes a backbone shift that explains a approximately 50 mV lower redox potential in one of the four rubredoxins. This result is consistent with experimental redox potentials of five additional rubredoxins with known sequence. Also, we attribute the unusually lower redox potentials of two of the HiPIPs studied to a less positive electrostatic environment around their redox sites. Finally, molecular dynamics simulations of solvent around static rubredoxin crystal structures indicate that water alone is a major factor in dampening the contribution of charged side chains, in accord with experiments showing that mutations of surface charges produce relatively little effect on redox potentials.

Alanine↗
PubMed1996

Theory of electrostatic effects in soft biological interfaces using atomic force microscopy.

We calculated the electrostatic force between a planar interface, such as a planar-supported lipid bilayer membrane, and the tip of a stylus on which another lipid bilayer or some other biomacromolecular system might be deposited. We considered styli with rounded tips as well as conical tips. To take into account the effect of dynamical hydrogen-bonded structures in the aqueous phase, we used a theory of nonlocal electrostatics. We used the Derjaguin approximation and identified the systems for which its use is valid. We pointed out where our approach differs from previous calculations and to what extent the latter are inadequate. We found that 1) the nonlocal interactions have significant effects over distances of 10-15 A from the polar zone and that, at the surface of this zone, the effect on the calculated force can be some orders of magnitude; 2) the lipid dipoles and charges are located a distance L from the hydrophobic layer in the aqueous medium and this can have consequences that may not be appreciated if it is ignored; 3) dipoles, located in the aqueous region, can give rise to forces even though the polar layer is unchanged, and if this is ignored the interpretation of force data can be erroneous if an attempt is made to rationalize an observed force with a knowledge of an uncharged surface; 4) the shape of the stylus tip can be very important, and a failure to take this into account can result in incorrect conclusions, a point made by other workers; and 5) when L is nonzero, the presence of charges and dipoles can yield a force that can be nonmonotonic as a function of ionic concentration.

Biophysical Phenomena↗
PubMed1997

The height of biomolecules measured with the atomic force microscope depends on electrostatic interactions.

In biological applications of atomic force microscopy, the different surface properties of the biological sample and its support become apparent. Observed height differences between the biomolecule and its supporting surface are thus not only of structural origin, but also depend on the different sample-tip and support-tip interactions. This can result in negative or positive contributions to the measured height, effects that are described by the DLVO (Derjaguin, Landau, Verwey, Overbeek) theory. Experimental verification shows that the electrostatic interactions between tip and sample can strongly influence the result obtained. To overcome this problem, pH and electrolyte concentration of the buffer solution have to be adjusted to screen out electrostatic forces. Under these conditions, the tip comes into direct contact with the surface of support and biological system, even when low forces required to prevent sample deformation are applied. In this case, the measured height can be related to the thickness of the native biological structure. The observed height dependence of the macromolecules on electrolyte concentration makes it possible to estimate surface charge densities.

Aquaporin 1↗
PubMed1997

The contribution of electrostatic and van der Waals interactions to the stereospecificity of the reaction catalyzed by lactate dehydrogenase.

Continuum electrostatic calculations in conjunction with molecular dynamics simulations have been used to investigate the source of the stereospecificity in the hydride transfer reaction catalyzed by lactate dehydrogenase (LDH). These studies show that favorable electrostatic interactions between the carboxamide group of the reduced nicotinamide adenine dinucleotide coenzyme and protein residues of the active site of LDH can account for much if not all of the stereospecificity of the LDH-catalyzed reaction, with A-side hydride transfer more than 10(7) times greater than B-side transfer. Unfavorable steric interactions within the binding complex for B-side transfer are not found.

Animals↗
PubMed1997

Electrostatic interactions between transmembrane segments mediate folding of Shaker K+ channel subunits.

In voltage-dependent Shaker K+ channels, charged residues E293 in transmembrane segment S2 and R365, R368, and R371 in S4 contribute significantly to the gating charge movement that accompanies activation. Using an intragenic suppression strategy, we have now probed for structural interaction between transmembrane segments S2, S3, and S4 in Shaker channels. Charge reversal mutations of E283 in S2 and K374 in S4 disrupt maturation of the protein. Maturation was specifically and efficiently rescued by second-site charge reversal mutations, indicating that electrostatic interactions exist between E283 in S2 and R368 and R371 in S4, and between K374 in S4 and E293 in S2 and D316 in S3. Rescued subunits were incorporated into functional channels, demonstrating that a native structure was restored. Our data indicate that K374 interacts with E293 and D316 within the same subunit. These electrostatic interactions mediate the proper folding of the protein and are likely to persist in the native structure. Our results raise the possibility that the S4 segment is tilted relative to S2 and S3 in the voltage-sensing domain of Shaker channels. Such an arrangement might provide solvent access to voltage-sensing residues, which we find to be highly tolerant of mutations.

Animals↗
PubMed1998

The coupling of electron transfer and proton translocation: electrostatic calculations on Paracoccus denitrificans cytochrome c oxidase.

We have calculated the electrostatic potential and interaction energies of ionizable groups and analyzed the response of the protein environment to redox changes in Paracoccus denitrificans cytochrome c oxidase by using a continuum dielectric model and finite difference technique. Subsequent Monte Carlo sampling of protonation states enabled us to calculate the titration curves of all protonatable groups in the enzyme complex. Inclusion of a model membrane allowed us to restrict the calculations to the functionally essential subunits I and II. Some residues were calculated to have complex titration curves, as a result of strong electrostatic coupling, desolvation, and dipolar interactions. Around the heme a3-CuB binuclear center, we have identified a cluster of 18 strongly interacting residues that account for most of the proton uptake linked to electron transfer. This was calculated to be between 0.7 and 1.1 H+ per electron, depending on the redox transition considered. A hydroxide ion bound to CuB was determined to become protonated to form water upon transfer of the first electron to the binuclear site. The bulk of the protonation changes linked to further reduction of the heme a3-CuB center was calculated to be due to proton uptake by the interacting cluster and Glu(II-78). Upon formation of the three-electron reduced state (P1), His325, modeled in an alternative orientation away from CuB, was determined to become protonated. The agreement of these results with experiment and their relevance in the light of possible mechanisms of redox-coupled proton transfer are discussed.

Catalysis↗
PubMed1998

Binding of basic peptides to membranes produces lateral domains enriched in the acidic lipids phosphatidylserine and phosphatidylinositol 4,5-bisphosphate: an electrostatic model and experimental results.

Direct fluorescence digital imaging microscopy observations demonstrate that a basic peptide corresponding to the effector region of the myristoylated alanine-rich C kinase substrate (MARCKS) self-assembles into membrane domains enriched in the acidic phospholipids phosphatidylserine (PS) and phosphatidylinositol 4,5-bisphosphate (PIP2). We show here that pentalysine, which corresponds to the first five residues of the MARCKS effector region peptide and binds to membranes through electrostatic interactions, also forms domains enriched in PS and PIP2. We present a simple model of domain formation that represents the decrease in the free energy of the system as the sum of two contributions: the free energy of mixing of neutral and acidic lipids and the electrostatic free energy. The first contribution is always positive and opposes domain formation, whereas the second contribution may become negative and, at low ionic strength, overcome the first contribution. Our model, based on Gouy-Chapman-Stern theory, makes four predictions: 1) multivalent basic ligands, for which the membrane binding is a steep function of the mole fraction of acidic lipid, form domains enriched in acidic lipids; domains break up at high concentrations of either 2) basic ligand or 3) monovalent salt; and 4) if multivalent anionic lipids (e.g., PIP2) are present in trace concentrations in the membrane, they partition strongly into the domains. These predictions agree qualitatively with experimental data obtained with pentalysine and spermine, another basic ligand.

4-Chloro-7-nitrobenzofurazan↗
PubMed1998

Electrostatic contributions to heat capacity changes of DNA-ligand binding.

Significant heat capacity changes (DeltaCp) often accompany protein unfolding, protein binding, and specific DNA-ligand binding reactions. Such changes are widely used to analyze contributions arising from hydrophobic and polar hydration. Current models relate the magnitude of DeltaCp to the solvent accessible surface area (ASA) of the molecule. However, for many binding systems-particularly those involving non-peptide ligands-these models predict a DeltaCp that is significantly different from the experimentally measured value. Electrostatic interactions provide a potential source of heat capacity changes and do not scale with ASA. Using finite-difference Poisson-Boltzmann methods (FDPB), we have determined the contribution of electrostatics to the DeltaCp associated with binding for DNA binding reactions involving the ligands DAPI, netropsin, lexitropsin, and the lambda repressor binding domain.

Antineoplastic Agents↗
PubMed1998

The pH-induced release of iron from transferrin investigated with a continuum electrostatic model.

A reduction in pH induces the release of iron from transferrin in a process that involves a conformational change in the protein from a closed to an open form. Experimental evidence suggests that there must be changes in the protonation states of certain, as yet not clearly identified, residues in the protein accompanying this conformational change. Such changes in protonation states of residues and the consequent changes in electrostatic interactions are assumed to play a large part in the mechanism of release of iron from transferrin. Using the x-ray crystal structures of human ferri- and apo-lactoferrin, we calculated the pKa values of the titratable residues in both the closed (iron-loaded) and open (iron-free) conformations with a continuum electrostatic model. With the knowledge of a residue's pKa value, its most probable protonation state at any specified pH may be determined. The preliminary results presented here are in good agreement with the experimental observation that the binding of ferric iron and the synergistic anion bicarbonate/carbonate results in the release of approximately three H+ ions. It is suggested that the release of these three H+ ions may be accounted for, in most part, by the deprotonation of the bicarbonate and residues Tyr-92, Lys-243, Lys-282, and Lys-285 together with the protonation of residues Asp-217 and Lys-277.

Apoproteins↗
PubMed1999

Electrostatically balanced subnanometer imaging of biological specimens by atomic force microscope.

To achieve high-resolution topographs of native biological macromolecules in aqueous solution with the atomic force microscope (AFM) interactions between AFM tip and sample need to be considered. Short-range forces produce the submolecular information of high-resolution topographs. In contrast, no significant high-resolution information is provided by the long-range electrostatic double-layer force. However, this force can be adjusted by pH and electrolytes to distribute the force applied to the AFM tip over a large sample area. As demonstrated on fragile biological samples, adjustment of the electrolyte solution results in a local reduction of both vertical and lateral forces between the AFM tip and proteinous substructures. Under such electrostatically balanced conditions, the deformation of the native protein is minimized and the sample surface can be reproducibly contoured at a lateral resolution of 0.6 nm.

Bacterial Proteins↗
PubMed1999

Simulation of electron-proton coupling with a Monte Carlo method: application to cytochrome c3 using continuum electrostatics.

A new method is presented for simulating the simultaneous binding equilibrium of electrons and protons on protein molecules, which makes it possible to study the full equilibrium thermodynamics of redox and protonation processes, including electron-proton coupling. The simulations using this method reflect directly the pH and electrostatic potential of the environment, thus providing a much closer and realistic connection with experimental parameters than do usual methods. By ignoring the full binding equilibrium, calculations usually overlook the twofold effect that binding fluctuations have on the behavior of redox proteins: first, they affect the energy of the system by creating partially occupied sites; second, they affect its entropy by introducing an additional empty/occupied site disorder (here named occupational entropy). The proposed method is applied to cytochrome c3 of Desulfovibrio vulgaris Hildenborough to study its redox properties and electron-proton coupling (redox-Bohr effect), using a continuum electrostatic method based on the linear Poisson-Boltzmann equation. Unlike previous studies using other methods, the full reduction order of the four hemes at physiological pH is successfully predicted. The sites more strongly involved in the redox-Bohr effect are identified by analysis of their titration curves/surfaces and the shifts of their midpoint redox potentials and pKa values. Site-site couplings are analyzed using statistical correlations, a method much more realistic than the usual analysis based on direct interactions. The site found to be more strongly involved in the redox-Bohr effect is propionate D of heme I, in agreement with previous studies; other likely candidates are His67, the N-terminus, and propionate D of heme IV. Even though the present study is limited to equilibrium conditions, the possible role of binding fluctuations in the concerted transfer of protons and electrons under nonequilibrium conditions is also discussed. The occupational entropy contributions to midpoint redox potentials and pKa values are computed and shown to be significant.

Binding Sites↗
PubMed1999

The electrostatics of lipid surfaces.

Charged lipids constitute a substantial fraction of all membrane lipids. Their charges vary in quantity and distribution within their headgroup regions. In long range interactions, their charges' value and electrostatic potential in the vicinity of the membrane surface can be approximated by the Guy-Chapman theory. This theory treats the interface as a charged structureless plain surrounded by uniform environments. However, if one considers intermolecular interactions, such assumptions need to be revised. The interface is in reality a thick region containing the residual charges of lipid headgroups. Their arrangement depends on the type of lipid present in the membrane. The variety of lipids and their biological functions suggests that charge distribution determines the extent and type of interaction with surface associated molecules. Numerous examples show that protein behavior at the lipid bilayer surface is determined by the type of lipid present, indicating protein specificity towards certain surface locations and local properties (determined by lipid composition) of a particular type. Such specificity is achieved by a combination of electrostatic, hydrophobic and enthropic effects. Comparing lipid biological activity, it can be stated that residual charge distribution is one of the factors of intermolecular recognition leading to the specific interaction of lipid molecules and selected proteins in various processes, particularly those involved with signal transduction pathways. Such specificity enables a variety of processes occurring simultaneously on the same membrane surface to function without cross-reaction interference.

Animals↗