Skip Navigation

Nucleic Acids Research 2006 34(Web Server issue):W43-W47; doi:10.1093/nar/gkl165
This Article
Right arrow Abstract Freely available
Right arrow Print PDF (92K) Freely available
Right arrow Screen PDF (94K) Freely available
Right arrow Supplementary Material
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Commercial Re-use Guidelines
for Open Access NAR Content
Google Scholar
Right arrow Articles by Kantardjiev, A. A.
Right arrow Articles by Atanasov, B. P.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Kantardjiev, A. A.
Right arrow Articles by Atanasov, B. P.
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2006. Published by Oxford University Press. All rights reserved
The online version of this article has been published under an open access model. Users are entitled to use, reproduce, disseminate, or display the open access version of this article for non-commercial purposes provided that: the original authorship is properly and fully attributed; the Journal and Oxford University Press are attributed as the original place of publication with the correct citation details given; if an article is subsequently reproduced or disseminated not in its entirety but only in part or as a derivative work this must be clearly indicated. For commercial re-use, please contact journals.permissions@oxfordjournals.org

Article

PHEPS: web-based pH-dependent Protein Electrostatics Server

Alexander A. Kantardjiev and Boris P. Atanasov*

Biophysical Chemistry Group, Institute of Organic Chemistry, Bulgarian Academy of Sciences Sofia-1113, Bulgaria

*To whom correspondence should be addressed. Tel: +359-2 960 6123; Fax: +359-2 870 0225; Email: boris{at}orgchm.bas.bg

Received February 10, 2006. Revised March 1, 2006. Accepted March 20, 2006.


   ABSTRACT
TOP
 ABSTRACT
INTRODUCTION
 METHODS
 IMPLEMENTATION
 RESULTS
 CONCLUSION AND FUTURE...
 SUPPLEMENTARY DATA
 REFERENCES
 
PHEPS (pH-dependent Protein Electrostatics Server) is a web service for fast prediction and experiment planning support, as well as for correlation and analysis of experimentally obtained results, reflecting charge-dependent phenomena in globular proteins. Its implementation is based on long-term experience (PHEI package) and the need to explain measured physicochemical characteristics at the level of protein atomic structure. The approach is semi-empirical and based on a mean field scheme for description and evaluation of global and local pH-dependent electrostatic properties: protein proton binding; ionic sites proton population; free energy electrostatic term; ionic groups proton affinities (pKa,i) and their Coulomb interaction with whole charge multipole; electrostatic potential of whole molecule at fixed pH and pH-dependent local electrostatic potentials at user-defined set of points. The speed of calculation is based on fast determination of distance-dependent pair charge-charge interactions as empirical three exponential function that covers charge–charge, charge–dipole and dipole–dipole contributions. After atomic coordinates input, all standard parameters are used as defaults to facilitate non-experienced users. Special attention was given to interactive addition of non-polypeptide charges, extra ionizable groups with intrinsic pKas or fixed ions. The output information is given as plain-text, readable by ‘RasMol’, ‘Origin’ and the like. The PHEPS server is accessible at http://pheps.orgchm.bas.bg/home.html.


   INTRODUCTION
TOP
 ABSTRACT
 INTRODUCTION
METHODS
 IMPLEMENTATION
 RESULTS
 CONCLUSION AND FUTURE...
 SUPPLEMENTARY DATA
 REFERENCES
 
Electrostatic phenomena are widely manifested as a fundamental feature of protein structure–function relationships (1–5). Protein molecules are very complex dielectric systems but can be treated as ‘solid state’ nanometer particles, immersed in buffered water solutions. Many approaches [‘macroscopic’—continuum dielectrics (6,7) and ‘microscopic’—polarizability (8,10)] were developed with different degree of validity: from simple TK-‘dielectric cavity models’ and analytical solution of Poisson–Boltzmann equation to their non-linear numerical and sophisticated empirical generalized Born solutions (11–13). The application of detailed and complex model description leads to increased difficulties for experimentalists to understand and use such sophisticated models. At present there are number of popular program packages [(14,15) and others] and web servers (16,17) but they are of limited significance for everyday problems of experimentalists. To the best of our knowledge, there is not available web server for fast pH-dependent calculation and analysis of protein electrostatic properties. Such software is needed because proteins are polyelectrolytes and their system of ionizable groups is pH-dependent. Many programs and servers compute pKas at ‘neutral pH’ yielding pH-independent pKas, which leads to erroneous results and distorts our view of principal properties of protein molecules, important for their functions. It is known that ‘pKa’ is directly related to free energy change of the corresponding protolytic reaction ({Delta}Ga = RTpKa) and that this {Delta}Ga is pH-dependent. It is well known that protein pKas are also pH-dependent, because ionic groups are closely arranged in the molecule. There are excellent theoretical works describing pH-dependent protein electrostatics [(7,18,19) and so on] but they are not straightforward for application by experimentalists. For many years a method addressing aforementioned requirements was developed and applied successfully in Biophysical Chemistry Laboratory at IOCh. The theoretical results are unequivocally validated by comparison with experimental studies as shown in a number of peer-reviewed publications over the years. Typical examples for this are pKas prediction of lysozyme, BPTI and cytochrome c (21–23); spectrophometric titration prediction and infrared carboxylic groups titration (24), enthalpy of protein ionization prediction (25), pH-dependent protein ultrasonic compressibility analysis (26); local electrostatic potentials (27); electrostatic contribution to a protein crystal lattice energy (28) and so on. The method was also applied to clarify enzyme mechanisms (29) and proved to be an invaluable tool for fast evaluation of electrostatic interactions and their analysis in large biomolecular immunochemical complexes (30). We hope that our server http://pheps.orgchm.bas.bg/home.html will be useful for experimentalist (protein scientist in need for fast evaluation of pH-dependent properties, enzymologists in need of pK values, spectroscopists and the like) as well as for in silico analysis by structural biologists and bioinformaticians. Being fast and easy to use this server is suitable for first acquaintance and training in the field.


   METHODS
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
IMPLEMENTATION
 RESULTS
 CONCLUSION AND FUTURE...
 SUPPLEMENTARY DATA
 REFERENCES
 
Protein self-consistent electrostatics
It is generally accepted that a model for protein electrostatics can be build on the assumption of continuum medium description, fixed atom approximation, protein–solvent boundary numerically described by atomic static accessibilities, SAi [variants of Lee-Richards algorithm (31)] and two type of charges: (i) permanent (pH-independent) partial charges (par) and (ii) proton-binding sites with pH-dependent titratable charges (tit). The model accepts experimentally measured pKa of model compounds (e.g. N-acetyl amides of each i-th ionogenic amino acids) (pKmod,i) and evaluates work for charge transfer from highly polar water solvent ({varepsilon}w = 80) to protein macromolecule (4 < {varepsilon}p,i < 40). Exposure to the solvent is evaluated by SAi at absence of other ionic groups. Born term, which is proportional to [1 – ({varepsilon}p,i/{varepsilon}w)], is roughly estimated as (1 – SAi). Partial charges assume values from AMBER and PARSE parameterization sets. Since the ratio of the number of ionic AAR (Nion) to total number of AAR (Ntot) Rel = Nion/Ntot is relative high for protein particles with small radii (Rp), the pairwise interaction between any i-th and j-th ionic groups counts contributions from charge–charge, charge–dipole and dipole–dipole interactions which can be simulated by an empirical three exponential curve:

Formula
where k = 1 for long-range (Columbic) interactions; k = 2 for mid-range, charge–dipole interactions; and k = 3 for short-range dipole–dipole interactions. The ak were estimated by a non-linear procedure by minimizing the functional F(a1,a2,a3) (21):

Formula
where the values of Zexp are taken from experimental data and Zth are the calculated values of the protein net charge as a function of pH. The initial values of the coefficients ak are obtained by numerical approximation of W(rij). Through extensive testing, using large dataset of structures, it was found that a1, a2 and a3 values are practically constants for a great number of proteins.

The pH-dependence of the electrostatic potential {Phi}el,i (pH) at the i-th proton binding site in PHEI was evaluated according to the following equation:

Formula
where Qj(pH) is defined by degree of dissociation or statistical mechanical proton population of given H+-binding site; Qj (pH) = (1 – <sj>) and Qj (pH) = –<sj> for basic and acidic groups respectively, where

Formula
Thus using partial titration of each j-th group we can find the pH-dependent net-charge of the whole molecule, Z(pH), i.e. potentiometric titration curve:

Formula
By definition if Z = 0 than pH = pI, i.e. the isoelectric point (the only pH at which the dipole moment of a protein molecule can be evaluated).

Thus starting with pKint,i = pKmod,i + {Delta}pKBorn,i + {Delta}pKpar,i, where pKmod,i is the pKa of the i-th site according to model compounds—see set of pKmod,i in (21,22,29); {Delta}pKBorn,i is the Born self-energy of the i-th site buried within the ‘uncharged’ protein, and {Delta}pKpar,i is the contribution of the i-th site interacting with the set of partial (permanent, fixed) atomic charges (see above).

Formula
where C is the Debye–Hückel term for ionic strength (Is). The term pKtit,i is the pKa shift of the i-th site caused by interactions with all other proton-binding groups and is evaluated according to efficient self-consistent iterative procedure (32). Coming to self consisted pH-dependent ionization the free energy term Gel(pH) is calculated as follows.

Formula
as well as pH-dependent Coulomb energy of each i-th ionic group with whole charge multipole:

Formula
After applying this iterative algorithm the electrostatic system is converged and all basic pH-dependent properties are reported.


   IMPLEMENTATION
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 IMPLEMENTATION
RESULTS
 CONCLUSION AND FUTURE...
 SUPPLEMENTARY DATA
 REFERENCES
 
The web sever is a front end of our program package PHEI, developed over many years in our Biophysical Chemistry Lab. Its current version is written in PERL and C/C++ by one of us (A.K.). Our package is capable of much more functionality and only basic electrostatic properties are presented online, the rest being under consideration for the next release (Conclusions and Future). The web implementation is driven by CGI/PERL routines. The only input file is a coordinate file in Protein Data Bank (PDB) format (33)—either user supplied or just as a PDB ID, following retrieval from our local PDB database. Following submission, the user is given some basic information about the protein molecule (chains; number of residues; ratio of ionogenic to all groups, Rel) and warned about certain inconsistencies in structure, related to subsequent calculation (interruption in residue numbering which might influence appearance of terminal charges). The user is given the possibility to edit initial setup of ionogenic groups (attention to CYS in SSBONDs and excluding covalently modified groups). This is accomplished by convenient interactive selection of used set of groups. This gives opportunity for simulation of ‘electrostatic mutagenesis’. Full ‘charge mutant analysis’ is supposed for next versions. The same screen visualizes the PDB file in a text field which allows for direct editing: adding missing terminal charges, fixed (non-titratable) whole or partial charges and titratable groups with user defined pKa intrinsic. All other parameters used as input are predefined or automatically calculated. After initial setup completion the calculation proceeds through aforementioned steps—evaluation of accessibilities and Born term {Delta}pKBorn,i, perturbation of pKa by partial charges {Delta}pKpar,i and finally the iterative procedure for self-consistent evaluation of titratable {Delta}pKtit,i.

To calculate {Delta}Gel(pH), the following energy conversion units were used: 1 kcal = 4.186 kJ = 1.68RT units (at 298 K) = 0.735 pKa units. The units of {phi}i(pH) are kcal/mol·e = 43.176 mV or 30.24 mC/m2. All calculations are provided at ionic strength (Is) 0.1.

The obtained results are organized in two groups: (i) GLOBAL [Z(pH), {Delta}Gel(pH) and {Phi}j at fixed pH] and (ii) LOCAL [si(pH), pKa,i(pH), Eel,i(pH) and {phi}i(pH)]. For each of them there is a link to own page. The contents of each page is comprised of the result itself, related derivatives (e.g. pI, {partial}Z/{partial}pH, pK1/2 and so on) as well as a short description and examples for visualization of this type of data. All output data files are in standard plain text format. Visualization is straightforward with any 2D plotting software and molecular graphics programs (RasMol, JMol, PyMol and so on).


   RESULTS
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 IMPLEMENTATION
 RESULTS
CONCLUSION AND FUTURE...
 SUPPLEMENTARY DATA
 REFERENCES
 
Global pH-characteristics
(1) pH-dependent protein net charge Z(pH) and its derivatives: Isoelectric point pI/Z = 0 and protein buffer capacity ß = {partial}pI/{partial}pH at three pH: (pI – 1.5), pI and (pI + 1.5).

It is equivalent to experimental potentiometric titration curve (34) and reflects basic global electrostatic characteristic of protein proton binding (35). The definition of pI is pH at which Z = 0. Protein buffer capacity (ß) is an important parameter for design of precise ion-exchange (36) and electrophoresis (37) experiments. The difference between two Z(pH) of analogous but perturbed states [e.g. apo-holo (30), oxidized-reduced, free-liganded and the like] can be useful in analysis the nature of such perturbation and identify pH region where it has maximal effect on proton binding. Other relevant issues are: the net charge of protein under condition of electro-spray mass-spectrometry (38); the critical Z-values at extreme pH in water (39,40) and in vacuum (41) at which protein ‘denature’ and many others (Supplementary Figure S1).

(2) pH-dependent electrostatic free energy term [{Delta}Gel(pH)] and its derivatives: {Delta}Gel,min, {Delta}Gel,pI; pHa and pHb at {Delta}Gel = 0 for acid and alkaline/basic denaturation, respectively.

Quantitative estimate for charge dependent stability {Delta}Gel(pH) is basic electrostatic characteristic of protein molecules (2). By evaluating {Delta}{Delta}Gel,ion(pH) = {Delta}Gel,holo(pH) – {Delta}Gel,apo(pH) it is possible to determine pH-dependent specific ion and/or cofactor binding (30). Similarly ‘electron affinity’ can be evaluated from difference {Delta}{Delta}Gel,e(pH) = {Delta}Gel,red(pH) – {Delta}Gel,oxid(pH) (42) and the like. It is easy to obtain experimental values for pHd,a and pHd,b and compare with calculated by our server (43). Another option is estimation of stability of pH-induced conformational states and evaluation of energetic barrier between them (44). Presence or absence of stricture ruled charge asymmetry is reflected in {Delta}Gel,min – {Delta}Gel,pI difference (also from their pHmin – pI) (Supplementary Figure S2).

(3) Electrostatic potential, EP({Phi}el) at user selected pH for all j-th protein non-hydrogen atoms in a PDB-formatted file and can be visualized in color scale by RasMol.

The electrostatic potential at each point within (45), on the molecular surface (46) and at near vicinity in solvent (47) for a protein molecule is its fundamental electrostatic characteristic (8,48). In fact all above quantities are derivative of {Phi}el = f(pH, ligands). Using present PHEPS version output file, it is straightforward to visualize {Phi}el (or EP) at each protein non-hydrogen atoms by switching on ‘color by temperature’ using color scale (dark blue: positive EP; green: zeroed EP; and red: negative EP) applicable to entire variety of RasMol model representations (Supplementary Figure S3).

Local pH-characteristics
(4) pH-dependent proton population or degree of ionization of each i-th ionic group (Si).

The results Si(pH) for ionic groups in order of increasing sequence numbers are presented in the form of column formatted file (all in one table).

Si(pH) can be related to NMR pH-dependent chemical shifts, {delta}i(pH) (49,50) or other individual titration characteristics—FTIR carboxylate titration (51); differential Tyr UV-titration (52); calorimetric/enthalpy titration [{Delta}Hi(pH)] (53) and so on (Supplementary Figure S4).

(5) pH-dependent proton affinity pKa,i(pH) at each individual i-th ionic site: The results pKa,i(pH) for ionic groups in order of increasing sequence numbers are presented in the form of column formatted file (all in one table). The set of pK1/2,i for each i-th ionizable group [their pKas at si(pH) = 0.5] is available in another table.

Predicted pK1/2,i can be compared directly to experimentally obtained. Plotting pKa,i(pH) is a fast way to differentiate ‘normal’ and ‘abnormal’ titration groups (19,53); functionally important ionogenic sites (54); {Delta}pKda = (pKdonor – pKacceptor) as function of pH to necessary for description pf H+-transfer processes (Hydrogen Bonded Networks, Brønsted's relations and the like) (55,56) (Supplementary Figure S5).

(6) pH-dependent electrostatic energy Eel,i(pH) of interaction of each i-th ionic group with whole multipole of partial and protonic/ionic charges—individual sites and their sum. The results Eel,i(pH) for ionic groups in order of increasing sequence numbers are presented in the form of column formatted file (all in one table).

The pH-dependent Columbic energy of interaction of given ionic group with whole charge multipole is evident characteristic reflecting specific electrostatic site property: influence on charge stability (57), participation in charge-driven processes (58), through space interactions with introduced charged systems in protein complexes (30) and the like (Supplementary Figure S6).

(7) pH-dependent local electrostatic potential, {phi}i(pH) at each i-th point within protein molecule or its close surrounding. The user is supposed to define points in PDB format. It is recommended that number of points do not exceed 20. The results are presented as pH-dependence of electrostatic potential at each point.

Knowledge of local EP at user-defined points is a great tool for elucidating electrostatic response of these sites to intra/inter-molecular interactions with charged groups (ions and dipoles of different kind of ligands in static and dynamic manner) (59,60). This characteristic is of indispensable use for evaluation of the effect of whole protein electrostatic field on crucially important sites (e.g. for understanding its role in intermediate species of enzyme catalytic cycles (59), protein stability (60,62) and many others (Supplementary Figure S7).

Tested proteins. All these features of our program package PHEI were developed and have been tested for many years. The method was applied to numerous proteins [Supplementary Table 1, (20–30)]. Calculated pK1/2-s was compared with experimental estimates of pKas (21,22) and correlation was made of calculated Z(pH) to published experimental curves.


   CONCLUSION AND FUTURE DEVELOPMENT
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 IMPLEMENTATION
 RESULTS
 CONCLUSION AND FUTURE...
SUPPLEMENTARY DATA
 REFERENCES
 
We hope this server will be useful to anyone who needs fast and detailed analysis of pH-dependent properties of a protein with known atomic structure and a tool for protein electrostatic design (61,62). We are ready to share our experience in the field with other protein scientists and are open for discussion.

Features in preparation for next PHEPS version are as follows:

  • For each AA in sequence order n (backbone, side chain and residue) with respective (B-factor)n, static accessibility <SA>n and {varphi}el,n.—now implemented.
  • 3D-contour EP map generation (in static and dynamic regime)—search of saddle and other critical points on multidimensional maps.
  • Correct determination of dipole (at pI) and electric (at any other pH) moments (µd and µe, respectively) using 3D-EP grid—their scalar and vector values.
  • Thorough ‘electrostatic mutation’ analysis with ‘{Delta} = mut – wild’ as function of pH—data for all mentioned above characteristics: {Delta}Z, {Delta}Si, {Delta}pKi, {Delta}{Delta}Gel, {Delta}Eel and {Delta}{Phi}el.
  • EP gradients (electrostatic forces, EF) at pH control, located at defined atoms and sites (user selected fragments, domains, subunits).

Many of these features are implemented in our program package PHEI, but their online access will be realized after extensive testing.


   SUPPLEMENTARY DATA
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 IMPLEMENTATION
 RESULTS
 CONCLUSION AND FUTURE...
 SUPPLEMENTARY DATA
REFERENCES
 
Supplementary Data are available at NAR Online.


   ACKNOWLEDGEMENTS
 
We thank Profs B. Honig and E. Alexov for kind donation of computers, one of which hosts our server. We thank L. Roumenina for her kind help in text edition and correction. This work is partially supported by grant X-1310 of National Fund ‘Scientific Research’, Sofia, Bulgaria. The Open Access publication charges for this article were waived by Oxford University Press

Conflict of interest statement. None declared.


   REFERENCES
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 IMPLEMENTATION
 RESULTS
 CONCLUSION AND FUTURE...
 SUPPLEMENTARY DATA
 REFERENCES
 

  1. Perutz, M.F. (1978) Electrostatic effects in proteins Science, 201, 1178–1191 .

  2. Honig, B. and Nicholls, A. (1995) Classical electrostatics in biology and chemistry Science, 268, 1144–1149[Abstract/Free Full Text] .

  3. Warshel, A. (1981) Electrostatic basis of structure–function correlation in proteins Acc. Chem. Res, . 14, 284–290[CrossRef][Web of Science] .

  4. Bashford, D. and Karplus, M. (1990) pKas of ionization groups in proteins: atomic detail from a continuum electrostatic model Biochemistry, 29, 10219–10225[CrossRef][Medline] .

  5. Antosiewicz, J., McCammon, J.A., and Gilson, M.K. (1994) Prediction of pH-dependent properties of proteins J. Mol. Biol, . 238, 415–436[CrossRef][Web of Science][Medline] .

  6. Simonson, T. (2001) Macromolecular electrostatics: continuum models and their grooving pains Curr. Opin. Struct. Biol, . 11, 243–253[CrossRef][Web of Science][Medline] .

  7. Bashford, D. (2004) Macroscopic electrostatic models for protonation state in proteins Front. Biosci, . 9, 1082–1099[Web of Science][Medline] .

  8. Warshel, A. and Papazyan, A. (1998) Electrostatic effects in macromolecules: fundamental concepts and practical modeling Curr. Opin. Struct. Biol, . 8, 211–217[CrossRef][Web of Science][Medline] .

  9. Braun-Sand, S. and Warshel, A. (2005) ‘Electrostatics of proteins: models and principles’ In Buchner, J. and Kiefhaber, T. (Eds.). Protein Folding Handbook, Part I, New York Wiley-VCH GmbH & Co .

  10. Aquist, J. (1996) Calculation of absolute binding free energies for charged ligands and effects of long-range electrostatic interactions J. Comput. Chem, . 17, 1587–1597[CrossRef][Web of Science] .

  11. Matthews, J.B. and Gurd, F.R.N. (1986) Calculation of electrostatic interactions in proteins Methods Enzymol, . 130, 413–453[Web of Science][Medline] .

  12. Petersen, M.T.N. and Petersen, S.B. (2001) How to lipases and esterases work: the electrostatic contribution J. Biotechnol, . 85, 115–147[CrossRef][Web of Science][Medline] .

  13. Feig, M., Onufriev, A., Lee, M.S., Im, W., Case, E.A., Brooks, C.L., III. (2004) Performance comparison of generalized Born and Poisson methods in the calculation of electrostatic solvation energies for protein structures J. Comput. Chem, . 25, 265–284[CrossRef][Web of Science][Medline] .

  14. Nayal, M., Hitz, B.C., Honig, B. (1999) GRASS: a server for the graphical representation and analysis of structures Protein Sci, . 8, 676–679[Web of Science][Medline] .

  15. Dolinsky, D.J., Nielson, J.E., McCammon, A.J., Backer, N.A. (2004) PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatic calculations Nucleic Acid Res, . 32, W665–W667[Abstract/Free Full Text] .

  16. Gordon, J.C., Meyers, J.B., Folta, T., Shoja, V., Heath, L.S., Onufriev, A. (2005) H++: a server for estimating pKas and adding missing hydrogens to macromolecules Nucleic Acid Res, . 33, W368–W371[Abstract/Free Full Text] .

  17. Miteva, M.A., Tuffery, P., Villoutreix, B.O. (2005) PCE: web tools to compute protein continuum electrostatics Nucleic Acid Res, . 33, W372–W375[Abstract/Free Full Text] .

  18. Georgesku, R.E., Alexov, E., Gunner, M. (2002) Combining conformational flexibility and continuum electrostatics for calculating pK(a)s in proteins Biophys. J, . 83, 1731–1748[Web of Science][Medline] .

  19. Spassov, V.Z. and Bashford, D. (1998) Electrostatic coupling to pH-titration sites are a source of cooperativity in protein-ligand binding Protein Sci, . 7, 2012–2025[Web of Science][Medline] .

  20. Atanasov, B. and Karshikov, A. (1985) Semi-empirical method for calculation of electrostatic interactions in proteins Studia Biophysica, 105, 11–22 .

  21. Spassov, V.Z., Karshikov, A.D., Atanasov, B.P. (1989) Electrostatic interactions in proteins: a theoretical analysis of lysozyme ionization Biochim. Biophys. Acta, 999, 1–6 .

  22. Karshikov, A.D., Engh, R., Bode, W., Atanasov, B.P. (1989) Electrostatic interactions in proteins: Calculations of the electrostatic term of free energy and the electrostatic potential field Eur. Biophys. J, . 17, 287–297 .

  23. Kossekova, G.P., Miteva, M.A., Atanasov, B.P. (1997) Characterization of pyridoxal phosphate as an optical label for measuring electrostatic potentials in proteins J. Photochem. Photobiol. B, 37, 74–83[CrossRef][Medline] .

  24. Miteva, M., Alexov, E., Atanasov, B. (1998) Numerical simulation of aldolase tetramer stability Eur. Biophys J, . 28, 67–73 .

  25. Atanasov, B. and Miteva, M. (1997) Prediction and structural analysis of the enthalpy of ionisation of proteins Thermochim. Acta, 291, 141–153[CrossRef] .

  26. Miteva, M.A., Mishonova, E.I., Atanasov, B.P. (1996) A theoretical model for evaluation of the ultrasonic velocimetric titration of proteins Comput. Rend. Acad. Bulg. Sci, . 49, 101–104 .

  27. Miteva, M.A., Kossekova, G.P., Villoutreix, B.O., Atanasov, B.P. (1997) Local electrostatic potentials in pyridoxal phosphate labeled horse heart cytochrome c J. Photochem. Photobiol. B, 37, 74–83 .

  28. Alexov, E. and Atanasov, B. (1994) Analysis of electrostatic interactions in ribonuclease A monoclinic crystal Biochim. Biophys. Acta, 1206, 55–62[CrossRef][Medline] .

  29. Atanasov, B., Mustafi, D., Makinen, M.W. (2000) Protonation of the ß-lactam nitrogen is the trigger event in the catalytic action of class A ß-lactamases Proc. Natl Acad. Sci. USA, 97, 3160–3165[Abstract/Free Full Text] .

  30. Roumenina, L.T., Kantardjiev, A.A., Atanasov, B.P., Waters, P., Gadjeva, M., Reid, K.B.M., Mantovani, A., Kishore, U., Kojouharova, M.S. (2005) Role of Ca2+ in the electrostatic stability and the functional activity of the globular domain of the human C1q Biochemistry, 44, 14097–14109[CrossRef][Medline] .

  31. Lee, B. and Richards, F.M. (1971) The interpretation of protein structures: estimation of static accessibility J. Mol. Biol, . 55, 379–400[CrossRef][Web of Science][Medline] .

  32. Karshikoff, A. (1995) A simple algorithm for calculation of multiple site titration curves Protein Eng, . 8, 243–248[Abstract/Free Full Text] .

  33. Westbrook, J., Feng, Z., Chen, L., Yang, H., Berman, H.M. (2003) The Protein Data Bank and structural genomics Nucleic Acid Res, . 31, 489–491[Abstract/Free Full Text] .

  34. Tanford, C. (1962) The interpretation of hydrogen ion titration curves of proteins Adv. Protein Chem, . 17, 69–165[Web of Science] .

  35. Nitta, K. and Sugai, Sh. (1972) Potentiometric titration studies on globular proteins Biopolymers, 11, 1893–1901[CrossRef][Web of Science][Medline] .

  36. Egmond, M.R., Antheunisse, W.P., van Bemmel, C.J., Ravestein, P., de Vlieg, J., Peters, H., Branner, S. (1994) Engineering surface charges in a subtilisin: the effects on electrophoretic and ion-exchange behavior Protein Eng, . 7, 793–800[Abstract/Free Full Text] .

  37. Palusinski, O.A., Su, Y., Fife, P.C. (1990) Numerical technique and computational procedure for isotachophoresis Electrophoresis, 11, 903–907[CrossRef][Web of Science][Medline] .

  38. Kaltashov, I. and Eyles, S.J. (2002) Studies of biomolecular conformations and conformational dynamics by Mass Spectrometry Mass Spectrom. Rev, . 21, 37–71[CrossRef][Web of Science][Medline] .

  39. Oliveberg, M., Vuilleumier, S., Fersht, A.R. (1994) Thermodynamic study of the acid denaturation of barnase and its dependence on ionic strength: evidence for residual electrostatic interactions in the acid/thermally denatured state Biochemistry, 33, 8826–8832[CrossRef][Medline] .

  40. Anderson, D., Becktel, W.J., Dahlquist, F.W. (1990) pH-induced denaturation of proteins: A single salt bridge contributes 3–5 kcal/mol to the free energy of folding of T4 lysozyme Biochemistry, 29, 2403–2408[CrossRef][Medline] .

  41. Valentine, S.J., Counterman, A.E., Clemmer, D.E. (1997) Conformer dependent proton transfer reactions of Ubiquitin ions J. Am. Soc. Mass Spectrom, . 8, 954–961[CrossRef][Web of Science] .

  42. Salgueiro, C.A., daCosta, P.N., Turner, D.L., Messias, A.C., van Dongen, W.M., Saraiva, L.M., Xavier, A.V. (2001) Effect of hydrogen-bond networks in controlling reduction potentials in Desulfovibrio vulgaris (Hildenborough) cytochrome C3 probed by site-specific mutagenesis Biochemistry, 40, 9709–9716[CrossRef][Medline] .

  43. Meiering, E.M., Serrano, L., Fersht, A.R. (1992) Effect of active site residues in barnase on activity and stability J. Mol. Biol, . 225, 585–589[CrossRef][Web of Science][Medline] .

  44. Yang, A.S. and Honig, B. (1994) Structural origin of pH and ionic-strength effects of protein stability: acid denaturation of Sperm Whale Myoglobin J. Mol. Biol, . 237, 602–612[CrossRef][Web of Science][Medline] .

  45. Franklin, J.C. and Cafiso, D.S. (1993) Internal electrostatic potentials in bilayers: measuring and controlling dipole potentials in lipid vesicles Biophys. J, . 65, 289–299[Web of Science][Medline] .

  46. Hirono, S., Liu, Q., Moriguchi, I. (1991) High correlation between hydrophobic free energy and molecular surface area characterized by electrostatic potential Chem. Pharm. Bull. (Tokyo), 39, 3106–3109 .

  47. Rashin, A.A., Iofin, M., Honig, B. (1991) Internal cavities and buried waters in globular proteins Biochemistry, 25, 3619–3625 .

  48. Welinder, K.G., Bjornholm, B., Dunford, H.B. (1995) Functions of electrostatic potentials and conserved distal and proximal His- Asp H-bonding networks in haem peroxidases Biochem. Soc. Trans, . 23, 257–262[Web of Science][Medline] .

  49. Craescu, C.T., Schaeffer, C., Mispelter, J., Garin, J., Rosa, J. (1986) High-resolution NMR studies of histidine-substituted and histidine-perturbed hemoglobin variants. Histidine assignments, electrostatic interactions at the protein surface, and implications for hemoglobin S polymerization J. Biol. Chem, . 261, 7894–7901[Abstract/Free Full Text] .

  50. Cocco, M.J., Kao, Y.H., Phillips, A.T., Lecomte, J.T. (1992) Structural comparison of apomyoglobin and metaquomyoglobin: pH titration of histidines by NMR spectroscopy Biochemistry, 31, 6481–6491[CrossRef][Medline] .

  51. Nadolny, C., Kempf, I., Zundel, G. (1993) Specific interactions of the allosteric effector 2,3-bisphosphoglycerate with human hemoglobin—a difference FTIR study Biol. Chem. Hoppe Seyler, 374, 403–407[Web of Science][Medline] .

  52. Honore, B. and Brodersen, R. (1992) Ionization of Tyrosine residues in human serum albumin and in its complexes with bilirubin and laurate Int. J. Pept. Protein Res, . 39, 24–28[Web of Science][Medline] .

  53. Sham, Y.Y., Chu, Z.T., Warshel, A. (1997) Consistent calculations of pKa of ionizable residues in proteins—semimicroscopic and microscopic approaches J. Phys. Chem. B, 101, 4458–4472[CrossRef] .

  54. Davoodi, J., Wakarchuk, W.W., Campbell, R.L., Carey, P.R., Surewicz, W.K. (1995) Abnormally high pKa of an active-site glutamic acid residue in Bacillus circulans xylanase. The role of electrostatic interactions Eur. J. Biochem, . 323, 839–843 .

  55. Wade, R.C. and Goodford, P.J. (1989) The role of hydrogen bonds in drug binding Prog. Clin. Biol. Res, . 289, 433–444[Medline] .

  56. Zhang, T. and Koshland, D.E., Jr. (1997) Computational method for relative binding energies of enzyme–substrate complexes Protein Sci, . 5, 348–356[Web of Science] .

  57. Akke, M. and Forsen, S. (1990) Protein stability and electrostatic interactions between solvent exposed charged side chains Proteins, 8, 23–29[CrossRef][Web of Science][Medline] .

  58. Swartz, P.D., Beck, B.W., Ichiye, T. (1996) Structural origin of redox potentials in Fe-S proteins—electrostatic potentials of crystal structures Biophys. J, . 71, 2958–2969[Web of Science][Medline] .

  59. Lokhart, D.J. and Kim, P.S. (1992) Internal Stark effect measurements of the electric field at the amino terminus of a {alpha}-helix Science, 257, 947–951[Abstract/Free Full Text] .

  60. Friedrich, K. and Wooley, P. (1988) Electrostatic potential of macromolecules measured by pKa shift of a fluorophore. I. The 3' terminus of 16 S RNA Eur. J. Biochem, . 173, 227–231[Web of Science][Medline] .

  61. Makhatadze, G.I., Loladze, V.V., Ermolenko, D.N., Chen, XF., Thomas, S.T. (2003) Contribution of surface salt bridges to protein stability: guidelines for protein engineering J. Mol. Biol, . 327, 1135–1148[CrossRef][Web of Science][Medline] .

  62. Makhatadze, G.I., Loladze, V.V., Gribenko, A.V., Lopez, M.M. (2004) Mechanism of thermostabilization in a designed cold shock protein with optimized surface electrostatic interactions J. Mol. Biol, . 336, 929–942[CrossRef][Web of Science][Medline] .


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
 Nucleic Acids ResHome page
A. A. Kantardjiev and B. P. Atanasov
PHEMTO: protein pH-dependent electric moment tools
Nucleic Acids Res., July 1, 2009; 37(suppl_2): W422 - W427.
[Abstract] [Full Text] [PDF]

This Article
Right arrow Abstract Freely available
Right arrow Print PDF (92K) Freely available
Right arrow Screen PDF (94K) Freely available
Right arrow Supplementary Material
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Commercial Re-use Guidelines
for Open Access NAR Content
Google Scholar
Right arrow Articles by Kantardjiev, A. A.
Right arrow Articles by Atanasov, B. P.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Kantardjiev, A. A.
Right arrow Articles by Atanasov, B. P.
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?