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 739 records · Page 41Linked to original sources

The refined 1.9-A X-ray crystal structure of D-Phe-Pro-Arg chloromethylketone-inhibited human alpha-thrombin: structure analysis, overall structure, electrostatic properties, detailed active-site geometry, and structure-function relationships.

Thrombin is a multifunctional serine proteinase that plays a key role in coagulation while exhibiting several other key cellular bioregulatory functions. The X-ray crystal structure of human alpha-thrombin was determined in its complex with the specific thrombin inhibitor D-Phe-Pro-Arg chloromethylketone (PPACK) using Patterson search methods and a search model derived from trypsinlike proteinases of known spatial structure (Bode, W., Mayr, I., Baumann, U., Huber, R., Stone, S.R., & Hofsteenge, J., 1989, EMBO J. 8, 3467-3475). The crystallographic refinement of the PPACK-thrombin model has now been completed at an R value of 0.156 (8 to 1.92 A); in particular, the amino- and the carboxy-termini of the thrombin A-chain are now defined and all side-chain atoms localized; only proline 37 was found to be in a cis-peptidyl conformation. The thrombin B-chain exhibits the characteristic polypeptide fold of trypsinlike serine proteinases; 195 residues occupy topologically equivalent positions with residues in bovine trypsin and 190 with those in bovine chymotrypsin with a root-mean-square (r.m.s.) deviation of 0.8 A for their alpha-carbon atoms. Most of the inserted residues constitute novel surface loops. A chymotrypsinogen numbering is suggested for thrombin based on the topological equivalences. The thrombin A-chain is arranged in a boomeranglike shape against the B-chain globule opposite to the active site; it resembles somewhat the propeptide of chymotrypsin(ogen) and is similarly not involved in substrate and inhibitor binding. Thrombin possesses an exceptionally large proportion of charged residues. The negatively and positively charged residues are not distributed uniformly over the whole molecule, but are clustered to form a sandwichlike electrostatic potential; in particular, two extended patches of mainly positively charged residues occur close to the carboxy-terminal B-chain helix (forming the presumed heparin-binding site) and on the surface of loop segment 70-80 (the fibrin[ogen] secondary binding exosite), respectively; the negatively charged residues are more clustered in the ringlike region between both poles, particularly around the active site. Several of the charged residues are involved in salt bridges; most are on the surface, but 10 charged protein groups form completely buried salt bridges and clusters. These electrostatic interactions play a particularly important role in the intrachain stabilization of the A-chain, in the coherence between the A- and the B-chain, and in the surface structure of the fibrin(ogen) secondary binding exosite (loop segment 67-80).(ABSTRACT TRUNCATED AT 400 WORDS)

Amino Acid Chloromethyl Ketones↗

Electrostatic effects in the control of glycogen phosphorylase by phosphorylation.

Electrostatic effects are important in the initial activation mechanism of glycogen phosphorylase by phosphorylation. Analysis of the electrostatic surface potential of glycogen phosphorylase with the program GRASP shows that in the unphosphorylated state, the N-terminal 20 residues, which include a number of basic amino acids, are located close to a position on the surface of the molecule that is highly acidic. Upon phosphorylation by phosphorylase kinase at Ser 14, the N-terminal residues change their position and conformation so that the Ser-P is directed away from the acidic patch and to an intersubunit site where 2 arginines bind the phosphate. This recognition site is created through tertiary and quaternary structural changes that accompany the activation mechanism.

Computer Simulation↗

Identification of electrostatic interaction sites between the regulatory and catalytic subunits of cyclic AMP-dependent protein kinase.

Two classes of molecules inhibit the catalytic subunit (C) of the cyclic AMP-dependent protein kinase (cAPK), the heat-stable protein kinase inhibitors (PKIs) and the regulatory (R) subunits. Basic sites on C, previously identified as important for R/C interaction in yeast TPK1 and corresponding to Lys213, Lys217, and Lys189 in murine C alpha, were replaced with either Ala or Thr and characterized for their kinetic properties and ability to interact with RI and PKI. rC(K213A) and rC(K217A) were both defective in forming holoenzyme with RI but were inhibited readily with PKI. This contrasts with rC(R133A), which is defective in binding PKI but not RI (Wen & Taylor, 1994). Thus, the C-subunit employs two distinct electrostatic surfaces to achieve high-affinity binding with these two types of inhibitory molecules even though all inhibitors share a common consensus site that occupies the active site cleft. Unlike TPK1, mutation of Lys189 had no effect. The mutant C subunits that were defective in binding RI, rC(K213A) and rC(K217A), were then paired with three RI mutants, rRI(D140A), rRI(E143A), and rRI(D258A), shown previously to be defective in recognition of C. Although the mutations at Asp140 and Asp258 in RI were additive with respect to the C mutations. rC(K213A) and rRI(E143A) were compensatory, thus identifying a specific electrostatic interaction site between RI and C. The results are discussed in terms of the RI and C crystal structures and the sequence homology between the yeast and mammalian enzymes.

Alanine↗

Computation of electrostatic complements to proteins: a case of charge stabilized binding.

Recent evidence suggests that the net effect of electrostatics is generally to destabilize protein binding due to large desolvation penalties. A novel method for computing ligand-charge distributions that optimize the tradeoff between ligand desolvation penalty and favorable interactions with a binding site has been applied to a model for barnase. The result is a ligand-charge distribution with a favorable electrostatic contribution to binding due, in part, to ligand point charges whose direct interaction with the binding site is unfavorable, but which make strong intra-molecular interactions that are uncloaked on binding and thus act to lessen the ligand desolvation penalty.

Bacterial Proteins↗

Electrostatic interactions in the acid denaturation of alpha-lactalbumin determined by NMR.

alpha-Lactalbumin (alpha-LA) undergoes a pH-dependent unfolding from the native state to a partially unfolded state (the molten globule state). To understand the role of electrostatic interactions in protein denaturation, NMR and CD pH titration experiments are performed on guinea pig alpha-LA. Variation of pH over the range of 7.0 to 2.0 simultaneously leads to the acid denaturation of the protein and the titration of individual ionizable groups. The pH titrations are interpreted in the context of these coupled events, and indicate that acid denaturation in alpha-LA is a cooperative event that is triggered by the protonation of two ionizable residues. Our NMR results suggest that the critical electrostatic interactions that contribute to the denaturation of alpha-LA are concentrated in the calcium binding region of the protein.

Animals↗

Protein molecular dynamics with electrostatic force entirely determined by a single Poisson-Boltzmann calculation.

The electrostatic force including the intramolecular Coulombic interactions and the electrostatic contribution of solvation effect were entirely calculated by using the finite difference Poisson-Boltzmann method (FDPB), which was incorporated into the GROMOS96 force field to complete a new finite difference stochastic dynamics procedure (FDSD). Simulations were performed on an insulin dimer. Different relative dielectric constants were successively assigned to the protein interior; a value of 17 was selected as optimal for our system. The simulation data were analyzed and compared with those obtained from 500-ps molecular dynamics (MD) simulation with explicit water and a 500-ps conventional stochastic dynamics (SD) simulation without the mean solvent force. The results indicate that the FDSD method with GROMOS96 force field is suitable to study the dynamics and structure of proteins in solution if used with the optimal protein dielectric constant.

Computational Biology↗

Role of the protein side-chain fluctuations on the strength of pair-wise electrostatic interactions: comparing experimental with computed pK(a)s.

The effect of the protein side-chain fluctuations on the strength of electrostatic interactions was studied. The effect was modeled on 7 different crystal structures on the same enzyme as well as on 20 molecular dynamics snapshot structures. It was shown that the side-chain flexibility affects predominantly the magnitude of the strong pair-wise interactions, that is, the pair-wise interaction among ion pairs, and practically does not affect the interactions with the rest of the protein. This was used to suggest a correction function that should be applied to the original pair-wise electrostatic interaction to mimic the effects of the fluctuations. The procedure is applied on three ion pairs identified in lysozyme. It was shown that sampling different side-chain rotamers and modifying the strength of the pair-wise interaction energies makes calculated pK(a)s less sensitive to the fluctuations of the structure and improves the prediction accuracy.

Amino Acids↗

The role of electrostatic interaction in the molecular recognition of selective agonists to metabotropic glutamate receptors.

The influence of electrostatic interactions in determining selectivity for individual subtypes of metabotropic glutamate receptors (mGluRs) is evaluated for a small set of agonists by using the program Delphi and the information thus obtained is compared with docking experiments carried out with AutoDock. The evaluation of the electrostatic component of the free energy of binding for L-Glu, L-AP4, or S-PPG to mGluR1, mGluR2, and mGluR4 subtypes allowed for the detection of subtle differences in the electronic properties of the three subtypes, differences that can account for the observed agonist selectivity.

Amino Acid Sequence↗

What are the dielectric "constants" of proteins and how to validate electrostatic models?

Implicit models for evaluation of electrostatic energies in proteins include dielectric constants that represent effect of the protein environment. Unfortunately, the results obtained by such models are very sensitive to the value used for the dielectric constant. Furthermore, the factors that determine the optimal value of these constants are far from being obvious. This review considers the meaning of the protein dielectric constants and the ways to determine their optimal values. It is pointed out that typical benchmarks for validation of electrostatic models cannot discriminate between consistent and inconsistent models. In particular, the observed pK(a) values of surface groups can be reproduced correctly by models with entirely incorrect physical features. Thus, we introduce a discriminative benchmark that only includes residues whose pK(a) values are shifted significantly from their values in water. We also use the semimacroscopic version of the protein dipole Langevin dipole (PDLD/S) formulation to generate a series of models that move gradually from microscopic to fully macroscopic models. These include the linear response version of the PDLD/S models, Poisson Boltzmann (PB)-type models, and Tanford Kirkwwod (TK)-type models. Using our different models and the discriminative benchmark, we show that the protein dielectric constant, epsilon(p), is not a universal constant but simply a parameter that depends on the model used. It is also shown in agreement with our previous works that epsilon(p) represents the factors that are not considered explicitly. The use of a discriminative benchmark appears to help not only in identifying nonphysical models but also in analyzing effects that are not reproduced in an accurate way by consistent models. These include the effect of water penetration and the effect of the protein reorganization. Finally, we show that the optimal dielectric constant for self-energies is not the optimal constant for charge-charge interactions.

Ions↗

Structure-based prediction of DNA-binding sites on proteins using the empirical preference of electrostatic potential and the shape of molecular surfaces.

Protein-DNA interactions play an essential role in the genetic activities of life. Many structures of protein-DNA complexes are already known, but the common rules on how and where proteins bind to DNA have not emerged. Many attempts have been made to predict protein-DNA interactions using structural information, but the success rate is still about 80%. We analyzed 63 protein-DNA complexes by focusing our attention on the shape of the molecular surface of the protein and DNA, along with the electrostatic potential on the surface, and constructed a new statistical evaluation function to make predictions of DNA interaction sites on protein molecular surfaces. The shape of the molecular surface was described by a combination of local and global average curvature, which are intended to describe the small convex and concave and the large-scale concave curvatures of the protein surface preferentially appearing at DNA-binding sites. Using these structural features, along with the electrostatic potential obtained by solving the Poisson-Boltzmann equation numerically, we have developed prediction schemes with 86% and 96% accuracy for DNA-binding and non-DNA-binding proteins, respectively.

Binding Sites↗

Scoring a diverse set of high-quality docked conformations: a metascore based on electrostatic and desolvation interactions.

Predicting protein-protein interactions involves sampling and scoring docked conformations. Barring some large structural rearrangement, rapidly sampling the space of docked conformations is now a real possibility, and the limiting step for the successful prediction of protein interactions is the scoring function used to reduce the space of conformations from billions to a few, and eventually one high affinity complex. An atomic level free-energy scoring function that estimates in units of kcal/mol both electrostatic and desolvation interactions (plus van der Waals if appropriate) of protein-protein docked conformations is used to rerank the blind predictions (860 in total) submitted for six targets to the community-wide Critical Assessment of PRediction of Interactions (CAPRI; http://capri.ebi.ac.uk). We found that native-like models often have varying intermolecular contacts and atom clashes, making unlikely that one can construct a universal function that would rank all these models as native-like. Nevertheless, our scoring function is able to consistently identify the native-like complexes as those with the lowest free energy for the individual models of 16 (out of 17) human predictors for five of the targets, while at the same time the modelers failed to do so in more than half of the cases. The scoring of high-quality models developed by a wide variety of methods and force fields confirms that electrostatic and desolvation forces are the dominant interactions determining the bound structure. The CAPRI experiment has shown that modelers can predict valuable models of protein-protein complexes, and improvements in scoring functions should soon solve the docking problem for complexes whose backbones do not change much upon binding. A scoring server and programs are available at http://structure.pitt.edu.

Models, Biological↗

Electrostatic calculations and model-building suggest that DNA bound to CAP is sharply bent.

Two observations suggest that DNA, upon binding to E. coli catabolite gene activator protein (CAP), is sharply bent by a total angle of at least 100-150 degrees: (1) The electrostatic potential field of CAP shows regions of positive potential that form a ramp on 3 sides of the protein. (2) The DNA binding site size as determined by DNA ethylation interference with binding, (Majors: "Control of the E. coli Lac Operon at the Molecular Level." Ph.D. Thesis, Harvard University, Cambridge, 1977) and by relative affinities of DNA fragments of various lengths (Liu-Johnson et al.: Cell 47:995-1005, 1986) requires severe bending of the DNA to maintain its favorable electrostatic contact with the protein.

Bacterial Proteins↗

Protein stability and electrostatic interactions between solvent exposed charged side chains.

To investigate the contribution to protein stability of electrostatic interactions between charged surface residues, we have studied the effect of substituting three negatively charged solvent exposed residues with their side-chain amide analogs in bovine calbindin D9k--a small (Mr 8,500) globular protein of the calmodulin superfamily. The free energy of urea-induced unfolding for the wild-type and seven mutant proteins has been measured. The mutant proteins have increased stability towards unfolding relative to the wild-type. The experimental results correlate reasonably well with theoretically calculated relative free energies of unfolding and show that electrostatic interactions between charges on the surface of a protein can have significant effects on protein stability.

Apoproteins↗

Ab initio study of the electrostatic multipole nature of torsional potentials in CH3SSCH3, CH3SSH, and HOOH.

The origin of torsional potentials in H3CSSCH3, H3CSSH, and HOOH and the anisotropy of the local charge distribution has been analyzed in terms of atomic multipoles calculated from the ab initio LCAO-MO-SCF wave function in the 6-31G* basis set. The results indicate that for longer -S-S-bonds the major contribution to these torsional barriers are electrostatic interactions of the atomic multipoles located on two atoms forming the rotated bond. This finding demonstrates the important role of electrostatic 1-2 interatomic interactions, usually neglected in conformational studies. It also opens the possibility to derive directly from accurate ab initio wave functions a simple nonempirical torsional potential involving atomic multipoles of two bonded atoms defining the torsional angle. For shorter -O-O- bonds, use of more precise models and inclusion of 1-3 interactions seems to be necessary.

Chemistry, Organic↗

Topological polymer chemistry by dynamic selection from electrostatic polymer self-assembly.

A collection of recent developments in topological polymer chemistry is presented. First, topological isomerism occurring on randomly coiled, flexible polymer molecules having cyclic and linear structures is discussed. Second, an electrostatic self-assembly and covalent fixation strategy has been developed for the synthesis of polymeric topological isomers. These isomers have double cyclic, manacle-, and theta-shaped constructions, and are prepared by using either linear or star telechelic polymer precursors having moderately strained cyclic ammonium salt groups, which carry multifunctional carboxylate counteranions. A technique of reversed-phase chromatography (RPC) is demonstrated as an effective means to separate polymers with different topologies, especially polymeric topological isomers. A further extension of topological polymer chemistry has been observed by dynamic selection from electrostatic polymer self-assembly to enable the effective formation of tadpole-shaped, cyclic-linear hybrid topologies.

Chromatography, Liquid↗

Electrostatic recognition between enzyme and inhibitor: interaction between papain and leupeptin.

Electrostatic forces are involved in a wide variety of molecular interactions that are of biological interest, including, among others, DNA-Protein interactions, protein folding, and the interactions between enzymes and their substrates and inhibitors. In this work, the interaction between papain and an inhibitor, leupeptin, is analyzed from the point of view of their electrostatic interaction. The computations enable one to suggest that negatively charged amino acids located in the region of the active site are responsible for creating an environment that enables efficient binding of the inhibitor. This binding occurs despite the fact that the net global charge of both molecules is positive; an explanation for this apparent contradiction is proposed.

Binding Sites↗

Redox control in heme proteins: electrostatic substitution in the active site of leghemoglobin.

The effects of electrostatic substitutions on the spectroscopic, ligand binding, and redox properties of the heme in leghemoglobin have been examined by replacement of the proximal leucine 88 residue with an aspartic acid residue (Leu88Asp). Electronic and resonance Raman spectra of the ferric derivative of Leu88Asp indicate a mixture of 6-coordinate, high-spin and 6-coordinate, low-spin hemes, analogous to that observed in the recombinant wild-type protein (rLb). At alkaline pH, formation of hydroxide-bound heme is indicated for Leu88Asp; the pK(a) for this transition (8.7 +/- 0.2, micro = 0.10 M, 25.0 degrees C) is 0.4 pH units higher than for rLb. Equilibrium dissociation constants (sodium phosphate, pH 7.0, micro = 0.10 M, 25.0 +/- 0.1 degrees C) for binding of anionic ligands (N(-)(3), nicotinate) to Leu88Asp are higher (K(d,nicotinate) = 6.8 +/- 0.2 microM; K(d,azide) = 33 +/- 0.6 microM) than the corresponding values for rLb (K(d,nicotinate) = 1.4 +/- 0.3 microM (pH 5.5, micro = 0.10 M, 25.0 +/- 0.1 degrees C); K(d,azide) = 4.8 +/- 0.2 microM). Resonance Raman spectra (sodium phosphate, pH 7.0, micro = 0.10 M) for the ferrous derivatives of Leu88Asp and rLb exhibit a strong nu(Fe-His) stretching frequency at 223 cm(-1) in both cases, indicating that the hydrogen bonding structure on the proximal side is not substantially altered in the variant. The reduction potential of Leu88Asp is -14 +/- 2 mV vs standard hydrogen electrode (SHE) (25.0 degrees C, micro = 0.10 M, pH 7.0), a decrease of 35 mV over the corresponding value for the wild-type protein under the same conditions (21 +/- 3 mV vs SHE). An assessment of these data in terms of electrostatic and hydrogen bonding considerations is presented.

Aspartic Acid↗

A test of the role of electrostatic interactions in determining the CO stretch frequency in carbonmonoxymyoglobin.

The vibrational frequency of CO bound to myoglobin can be varied by up to 60 cm-1 by making site-specific mutations in the distal pocket. These changes may result from specific chemical interactions between distal amino acids and the CO or from changes in the electrostatic field of the distal pocket. In this paper, we separate the relative contributions of these two effects by comparing the IR spectra of the carbonmonoxy complexes of human myoglobin mutants V68N, V68D, and V68E. The effect of replacing valine with these polar amino acids on the electrostatic environment of the distal heme pocket has been independently determined earlier by measurements of the heme reduction potential and electronic absorption spectral band shifts. While all three mutations result in a negative dipole pointing towards the CO ligand, the CO stretch frequency shifts differently in each case. These differences are attributed to specific chemical interactions between the amino acids and the CO ligand.

Electrochemistry↗