Skip Navigation



Nucleic Acids Research Advance Access published online on May 28, 2008
Nucleic Acids Research, doi:10.1093/nar/gkn314
This Article
Right arrow Abstract Freely available
Right arrow Print PDF (705K) Freely available
Right arrow Screen PDF (246K) Freely available
Right arrow Supplementary Data
Right arrowOA All Versions of this Article:
36/suppl_2/W270    most recent
gkn314v1
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 Jo, S.
Right arrow Articles by Im, W.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Jo, S.
Right arrow Articles by Im, W.
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© 2008 The Author(s)
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/2.0/uk/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Web Server Issue

PBEQ-Solver for online visualization of electrostatic potential of biomolecules

Sunhwan Jo1, Miklos Vargyas2, Judit Vasko-Szedlar2, Benoît Roux3 and Wonpil Im4,*

1Department of Chemistry, The University of Kansas, 2030 Becker Drive, Lawrence, KS 66047, USA, 2ChemAxon Kft., Máramaros köz 3/a, Budapest, 1037 Hungary, 3Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Sciences, W323B 929 East 57th Street, Chicago, IL 60637 and 4Department of Molecular Biosciences and Center for Bioinformatics, The University of Kansas, 2030 Becker Drive, Lawrence, KS 66047, USA

*To whom correspondence should be addressed. Tel: (785) 864 1993; Fax: (785) 864 5558; Email: wonpil{at}ku.edu

Received January 6, 2008. Revised April 25, 2008. Accepted May 7, 2008.


   ABSTRACT
TOP
 ABSTRACT
INTRODUCTION
 METHODS AND ILLUSTRATIVE RESULTS
 CONCLUDING DISCUSSION & FUTURE...
 SUPPLEMENTARY DATA
 REFERENCES
 
PBEQ-Solver provides a web-based graphical user interface to read biomolecular structures, solve the Poisson-Boltzmann (PB) equations and interactively visualize the electrostatic potential. PBEQ-Solver calculates (i) electrostatic potential and solvation free energy, (ii) protein–protein (DNA or RNA) electrostatic interaction energy and (iii) pKa of a selected titratable residue. All the calculations can be performed in both aqueous solvent and membrane environments (with a cylindrical pore in the case of membrane). PBEQ-Solver uses the PBEQ module in the biomolecular simulation program CHARMM to solve the finite-difference PB equation of molecules specified by users. Users can interactively inspect the calculated electrostatic potential on the solvent-accessible surface as well as iso-electrostatic potential contours using a novel online visualization tool based on MarvinSpace molecular visualization software, a Java applet integrated within CHARMM-GUI (http://www.charmm-gui.org). To reduce the computational time on the server, and to increase the efficiency in visualization, all the PB calculations are performed with coarse grid spacing (1.5 Å before and 1 Å after focusing). PBEQ-Solver suggests various physical parameters for PB calculations and users can modify them if necessary. PBEQ-Solver is available at http://www.charmm-gui.org/input/pbeqsolver.


   INTRODUCTION
TOP
 ABSTRACT
 INTRODUCTION
METHODS AND ILLUSTRATIVE RESULTS
 CONCLUDING DISCUSSION & FUTURE...
 SUPPLEMENTARY DATA
 REFERENCES
 
Implicit solvent treatments are approximate methods that attempt to incorporate the average influence of the molecular environment on a system of interest without having to explicitly simulate the molecules constituting this environment (1). Implicit solvent methods have emerged as a popular strategy to approximate bulk solvent or membrane environments and have been applied successfully to protein–protein or protein–ligand binding thermodynamics, scoring of protein conformations in structure prediction, peptide and protein folding/unfolding studies and ion channels (1–4). In particular, Poisson-Boltzmann (PB) continuum electrostatics, in which the solvent is represented as a featureless dielectric material, is the most rigorous and popular method to estimate the electrostatic solvation energy of a solute with an arbitrary shape, and particular successes in applications to complex biological problems are evident (2,5–7). The characterization of the electrostatic potential on the macromolecular surface by solving the PB equation is becoming a routine practice in structural biology (5).

Over the last two decades, considerable efforts have been made to generalize and enhance the computational methodologies and techniques to solve the PB equation and visualize the calculated electrostatic potential (7–13). As a result, various user-friendly programs that provide numerical solutions of the PB equation using finite-difference or finite-element methods with a discretized grid are now available as standalone software such as APBS (7), MEAD (14), Qnifft (15), DelPhi (16) and Zap (17), or as modules in biomolecular modeling and simulation programs such as PB solver (18) in AMBER (19), PBEQ (12,20–22) in CHARMM (23), APBS (7) in AMBER (7), CHARMM (23), TINKER (24) and NAMD (25), and PB solver (26) in Jaguar (http://www.schrodinger.com). PyMOL (27), VMD (28), GRASP (29), PMV (30) and DINO (http://www.dino3d.org) (31) provide various visualization tools for calculated electrostatic potentials. There are also several useful web-based interfaces to setup and perform PB calculations such as PDB2PQR (32), PCE (33) and PDB|Hydro (34).

We have developed PBEQ-Solver (http://www.charmm-gui.org/input/pbeqsolver) to provide a web-based graphical user interface (GUI) to read biomolecular protein data bank (PDB) structures (35) through PDB Reader at the CHARMM-GUI website (http://www.charmm-gui.org), solve the PB equations using PBEQ module (12,20–22) in CHARMM (23) and interactively visualize the electrostatic potential on the solvent-accessible surface as well as iso-electrostatic potential contours using the MarvinSpace molecular visualization software (http://www.chemaxon.com/product/mspace.html), a Java applet integrated within CHARMM-GUI. In addition to the calculations of electrostatic potential and solvation free energy, PBEQ-Solver also computes protein–protein (DNA or RNA) electrostatic interaction energy and Formula of a selected titratable residue. All the calculations can be performed in both aqueous solvent and membrane environments (with a cylindrical pore in the case of membrane).


   METHODS AND ILLUSTRATIVE RESULTS
TOP
 ABSTRACT
 INTRODUCTION
 METHODS AND ILLUSTRATIVE RESULTS
CONCLUDING DISCUSSION & FUTURE...
 SUPPLEMENTARY DATA
 REFERENCES
 
PB theory
The electrostatic solvation energy {Delta}Gelec is the work required to assemble the charges {q{alpha}} of the solute in the solvent (1), which can be expressed in terms of the reaction field potential Formula ,

Formula 1

(1)
Based on continuum electrostatics, the reaction field potential, Formula can be computed by solving the PB equation twice for the reference electrostatic potential Formula and the electrostatic potential in the solvent environment Formula ,

Formula 2

(2)
where Formula , Formula and Formula are the dielectric constant, the modified Debye–Hückel screening factor and the fixed charge density of the solute, respectively. For a given dielectric constant in the interior of the solute (Formula ), Formula is calculated by setting the dielectric constant to Formula (dielectric constant of reference environment) at all points in the solvent region, while Formula is calculated by setting the dielectric constant to Formula (solvent dielectric constant) and Formula (from input concentration) in the solvent region. The influence of biological membranes and cylindrical pore can be incorporated into Formula by setting the dielectric constants to Formula (membrane dielectric constant) and Formula (pore dielectric constant) in the membrane and pore regions, respectively.

In the case of a biomolecular complex containing proteins, DNA and RNA, the electrostatic interaction energy in solvent between a selected chain [A] and the rest [B] in the complex can be computed as

Formula 3

(3)
where Formula , Formula and Formula are the electrostatic potentials of the complex, [A] and [B] in solvent or membrane environments, respectively. The Formula of a residue in a protein, Formula , is estimated based on calculating, for both the protonated (p) and unprotonated (u) states of the residue, the difference between its electrostatic free energy when it is in the protein environment (Formula ) and its electrostatic free energy when it is isolated in solution (Formula ). The Formula shift (Formula ) relative to the Formula of the same amino acid isolated in solution (Formula ) is given by

Formula 4

(4)
where

Formula 5

(5)

PDB Reader
Since PBEQ-Solver uses the PBEQ module in CHARMM, its first step is to read a PDB file into CHARMM, which is generally considered not to be straightforward due to complexity of PDB files. In general, this is also the first difficulty that any simulation program users may face. It is typically even harder to introduce different protonation states of titratable residues, disulfide bonds or other posttranslational modifications such as phosphorylation. For a seamless, efficient procedure, PBEQ-Solver first uses PDB Reader in CHARMM-GUI that provides a flexible web-based GUI to convert a PDB file [downloaded from RCSB (35), http://www.rcsb.org, or uploaded from user's computer] into CHARMM readable files with the following options; (i) partial selection of protein chains as well as model selection in the case of NMR structures, (ii) modification of engineered residues, (iii) terminal group selection, (iv) protonation selection, (v) disulfide bond selection, (vi) phosphorylation selection, (vii) generation of a biologically functional unit and (viii) generation of a crystal packing. The terminal patch residues available in PDB Reader are listed in Supplementary Table S1. If indicated in a PDB file, PDB Reader automatically detects the disulfide bonds and displays them. Users can always add or remove them in the list. PDB Reader also automatically detects some of engineered residues listed in Supplementary Table S2 and converts them to corresponding natural residues. For example, an engineered residue ‘Sep’ in a PDB file represents a phosphorylated Ser residue. In such case, PDB Reader automatically converts Sep to Ser and turns on phosphorylation of the Ser residue (Supplementary Figure S3). However, due to the complexity of PDB files with various heteroatoms and other engineered residues, it is users that need to check if all the residues to be read do exist in a CHARMM topology file (currently, top|all27|prot|na.rtf). If there are undetermined coordinates in selected chains, PDB Reader simply builds them using a predetermined internal coordinate table (‘IC BUILD’ command in CHARMM). The usages of PDB Reader is well illustrated in a video demo, ‘PDB Reader Tutorial’, available at the CHARMM-GUI website (http://www.charmm-gui.org/?doc=demo).

It should be stressed that RCSB PDB structures do not contain orientation information of a membrane(-bound) protein relative to lipid bilayers. Therefore, in order to perform reasonable PB calculations in a membrane environment, users need to validate if the protein structure is properly oriented with respect to membranes. It is assumed that the membrane normal is parallel to the Z-axis and its center by default is located at Z = 0 (Table 1). Users can preorient the structure in their local machine and upload it, or use preoriented protein structures from the OPM database (http://opm.phar.umich.edu) (36) by selecting OPM as PDB download source (Supplementary Figure S1).


View this table:
[in this window]
[in a new window]

  		Table 1. Physical parameters in PBEQ-Solver

 
PBEQ-Solver
After PDB Reader, PBEQ-Solver displays the options for three types of PB calculations as well as adjustable physical parameters summarized in Table 1 (see also Supplementary Figure S4). Note that, to reduce the computational time on the server and to increase the efficiency in visualization, all the PB calculations are presently performed with (unadjustable) coarse grid spacing (Dcel|c = 1.5 Å before and Dcel|f = 1 Å after focusing (37)). All the PB calculations uses the molecular surface (with a probe radius of 1.4 Å) to setup the dielectric boundary based on a set of atomic radii for proteins, DNA and RNA that were optimized to closely reproduce the charging free energies calculated by molecular dynamics free energy simulations for standard amino acids (20) and nucleic acids (22) (filename: step2|radii.str, which is downloadable as shown in Supplementary Figure S5).

PBEQ-Solver currently offers three types of PB calculations for (i) electrostatic potential and solvation free energy [based on Equation 1], (ii) protein–protein (DNA or RNA) electrostatic interaction energy [based on Equation 3] and (iii) pKa of a selected titratable residue [based on Equation 4]. As mentioned before, all the calculations can be performed in both aqueous solvent and membrane environments (with a cylindrical pore in the case of membrane). Based on the user's inputs for the physical parameters listed in Table 1, PBEQ-Solver uses the PBEQ module in CHARMM to solve the finite-difference PB equation of molecules specified by users through PDB Reader. Note that a PDB file must contain at least two chains for the electrostatic interaction calculations, such that the interaction energy is calculated between a selected chain and the rest (Supplementary Figure S8). In the case that only one chain is selected in PDB Reader, the interaction energy option is not available.

Tests and illustrations
We have tested PBEQ-Solver with 360 PDB structures whose number of residues range from about 50 to 4000 to examine the computational time to calculate the electrostatic solvation free energy (see Supplementary Table S3 for the full list of PDB IDs). Each test took about 1–14 min on the server depending on the size of the biomolecule. In addition, to validate if the CHARMM input used in PBEQ-Solver is reliable, we have compared the PB solvation energies calculated based on the PBEQ-Solver input with those from other PB programs that are published by Feig et al. (38). As shown in Supplementary Table S4, the PBEQ-Solver input yields similar results to other PB programs when the same parameters are used. This also demonstrates that the visualization of calculated potentials with the input and the coarse grid in PBEQ-Solver is also reliable.

We have illustrated how to use PBEQ-Solver with an example of PDB:1KDX (complex of KIX and phosphorylated KID domains) in Supplementary Material (Figures S1–S11), including all the snapshots with some useful annotations for the three types of PB calculations. Figures 1 and 2 illustrate various molecular graphics that the users can generate using PBEQ-Solver. As shown in Figure 1, the protein complex (PDB:1KDX) of the KIX and phosphorylated KID domains that play an important role in regulation by posttranslational modification (39) can be viewed with its ribbon representation, solvent-accessible surface representation with electrostatic potential and iso-electrostatic potential contours. Note that the unit for the bounds on electrostatic potential visualization is kcal/(mol·e), where e is the unit charge. The iso-contour map helps users to examine how the electrostatic potentials are distributed in a distal place (other than on the surface of the biomolecule). In Figure 2, the protein–DNA complex (PDB:1E3M) of DNA mismatch repair protein MutS and G–T mismatch DNA (40) is shown in its surface electrostatic potentials with and without DNA. It is clearly shown that the DNA binding site has the high positive electrostatic surface for negatively charged DNA.


Figure 1
View larger version (24K):
[in this window]
[in a new window]
[Download PowerPoint slide]
  		Figure 1. Molecular graphics views of the protein complex (PDB:1KDX) of the KIX domain and one of its co-activators, the phosphorylated KID domain, that play an important role in regulation by posttranslational modification (39). PBEQ-Solver provides a tool for online visualization of its (left) ribbon representation as well as (middle) solvent-accessible surface representation with electrostatic potential [+2 kcal/(mol·e) in blue to –2 kcal/(mol·e) in red] and (right) iso-electrostatic potential contours [+1 kcal/(mol·e) in blue and –1 kcal/(mol·e) in red].

 

Figure 2
View larger version (63K):
[in this window]
[in a new window]
[Download PowerPoint slide]
  		Figure 2. Molecular graphics views of the protein–DNA complex (PDB:1E3M) of DNA mismatch repair protein MutS and G–T mismatch DNA (40) with the surface electrostatic potentials (left) with and (right) without DNA [+2 kcal/(mol·e) in blue to –2 kcal/(mol·e) in red]. To generate the surface potential without DNA, one has to unselect DNA during the PDB reading step.

 
The web interfaces for PBEQ-Solver and PDB Reader in CHARMM-GUI have been developed following de facto web standards and extensively tested under Firefox (version 2) and Safari (version 2). Although all the functionalities work just fine with different web browsers, we have so far found some minor inconsistency of user interfaces in Internet Explorer (IE) and the beta version of Safari (version 3). In the case of IE versions 6 and 7, some user interface elements look slightly different than those in Firefox and Safari. In the case of beta version of Safari (version 3), the (job) monitoring panel does not work properly.


   CONCLUDING DISCUSSION & FUTURE DIRECTIONS
TOP
 ABSTRACT
 INTRODUCTION
 METHODS AND ILLUSTRATIVE RESULTS
 CONCLUDING DISCUSSION & FUTURE...
SUPPLEMENTARY DATA
 REFERENCES
 
We have described the functionalities of PBEQ-Solver at the CHARMM-GUI website with illustrations. There are several programs available for visualization of electrostatic potential of biomolecules, calculated by solving the PB equations, so that users can use them in their local machines or generate images through web-based visualization applications. However, we believe that PBEQ-Solver provides unique and interactive web-based GUI for online visualization. We recently found that Honig and coworkers have made similar efforts through the Mark-US server, a functional annotation server (http://luna.bioc.columbia.edu/honiglab/mark-us). In particular, such a GUI development will help non-expert users, especially experimentalist, to perform various PB calculations and visualize calculated electrostatic potentials of their own systems. In a similar way, such an effort will be also useful for educational purposes.

The development of PBEQ-Solver is an ongoing project. Its present drawback is to use relatively coarse grid spacing [1.5 Å before and 1.0 Å after focusing (37)] for efficient PB calculations and visualization. However, it is recommended to use 1–1.5 Å before and at least 0.5 Å after focusing in practical applications for solvation energy calculations. One way to overcome this problem is to download all the CHARMM input files and redo the calculations with finer grid spacing in users’ local machines. We also expect to provide the finer grid calculations as the capability of our server increases. The calculation of transmembrane potential using the modified PB equation (21) and its visualization will be incorporated into PBEQ-Solver.


   SUPPLEMENTARY DATA
TOP
 ABSTRACT
 INTRODUCTION
 METHODS AND ILLUSTRATIVE RESULTS
 CONCLUDING DISCUSSION & FUTURE...
 SUPPLEMENTARY DATA
REFERENCES
 
Supplementary data are available at NAR Online.


   ACKNOWLEDGEMENTS
 
The authors are grateful to Martin Karplus for his support, to Nathan Baker for sharing useful information and to Michael Feig for sharing the results of other PB programs. Sunhwan Jo is the recipient of the Undergraduate Research Assistant Fund from the University of Kansas. Wonpil Im is 2007 Alfred P. Sloan Research Fellow. This work was supported by institutional funding from the University of Kansas (to W.I.) and grant 0415784 from the National Science Foundation (to B.R.). Funding to pay the Open Access publication charges for this article was provided by National Science Foundation grant MCB-0415784.

Conflict of interest statement. None declared.


   REFERENCES
TOP
 ABSTRACT
 INTRODUCTION
 METHODS AND ILLUSTRATIVE RESULTS
 CONCLUDING DISCUSSION & FUTURE...
 SUPPLEMENTARY DATA
 REFERENCES
 

  1. Roux B, Simonson T. Implicit solvent models. Biophys. Chem. (1999) 78:1–20.[CrossRef][Web of Science][Medline]

  2. Roux B, Bernèche S, Im W. Ion channels, permeaion, and electrostatics: insight into the function of KcsA. Biochemistry (2000) 39:13295–13306.[CrossRef][Web of Science][Medline]

  3. Im W, Chen J, Brooks C.L. III. Peptide and protein folding and conformational equilibria: theoretical treatment of electrostatics and hydrogen bonding with implicit solvent models. Adv. Protein Chem. (2005) 72:173–198.[Web of Science][Medline]

  4. Feig M, Brooks C.L. III. Recent advances in the development and application of implicit solvent models in biomolecule simulations. Curr. Opin. Struct. Biol. (2004) 14:217–224.[CrossRef][Web of Science][Medline]

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

  6. Murray D, Honig B. Electrostatic control of the membrane targeting of C2 domains. Mol. Cell. (2002) 9:145–154.[CrossRef][Web of Science][Medline]

  7. Baker NA, Sept D, Joseph S, Holst MJ, McCammon JA. Electrostatics of nanosystems: application to microtubules and the ribosome. Proc. Natl Acad. Sci. USA (2001) 98:10037–10041.[Abstract/Free Full Text]

  8. Warwicker J, Watson HC. Calculation of the electric potential in the active site cleft due to alpha-helix dipoles. J. Mol. Biol. (1982) 157:671–679.[CrossRef][Web of Science][Medline]

  9. Klapper I, Hagstrom R, Fine R, Sharp K, Honig B. Focusing of electric fields in the active site of Cu-Zn superoxide dismutase: effects of ionic strength and amino-acid modification. Proteins (1986) 1:47–59.[CrossRef][Medline]

  10. Bashford D, Karplus M. pKa's of ionizable groups in proteins: atomic detail from a continuum electrostatic model. Biochemistry (1990) 29:10219–10225.[CrossRef][Web of Science][Medline]

  11. Gilson MK, McCammon JA, Madura JD. Molecular dynamics simulation with a continuum electrostatic model of the solvent. J. Comput. Chem. (1995) 16:1081–1095.[CrossRef][Web of Science]

  12. Im W, Beglov D, Roux B. Continuum solvation model: electrostatic forces from numerical solutions to the Poisson-Bolztmann equation. Comput. Phys. Comm. (1998) 111:59–75.[CrossRef][Web of Science]

  13. Zhou YC, Feig M, Wei GW. Highly accurate biomolecular electrostatics in continuum dielectric environments. J. Comput. Chem. (2007) 29:87–97.[Web of Science]

  14. Bashford D. An object-oriented programming suite for electrostatic effects in biological molecules. In: Scientific Computing in Object-Oriented Parallel Environments—Ishikawa Y, Oldehoeft RR, Reynders JVW, eds. (1997) Springer, Berlin, Vol. 1343, pp. 233–240.

  15. Chin K, Sharp KA, Honig B, Pyle AM. Calculating the electrostatic properties of RNA provides new insights into molecular interactions and function. Nat. Struct. Biol. (1999) 6:1055–1061.[CrossRef][Web of Science][Medline]

  16. Rocchia W, Sridharan S, Nicholls A, Alexov E, Chiabrera A, Honig B. Rapid grid-based construction of the molecular surface and the use of induced surface charge to calculate reaction field energies: applications to the molecular systems and geometric objects. J. Comput. Chem. (2002) 23:128–137.[CrossRef][Web of Science][Medline]

  17. Grant JA, Pickup BT, Nicholls A. A smooth permittivity function for Poisson-Boltzmann solvation methods. J. Comput. Chem. (2001) 22:608–640.[CrossRef][Web of Science]

  18. Hsieh MJ, Luo R. Physical scoring function based on AMBER force field and Poisson-Boltzmann implicit solvent for protein structure prediction. Proteins (2004) 56:475–486.[CrossRef][Web of Science][Medline]

  19. Case DA, Cheatham T.E. III, Darden T, Gohlke H, Luo R, Merz K.M. Jr, Onufriev A, Simmerling C, Wang B, Woods RJ. The Amber biomolecular simulation programs. J. Comput. Chem. (2005) 26:1668–1688.[CrossRef][Web of Science][Medline]

  20. Nina M, Beglov D, Roux B. Atomic radii for continuum electrostatics calculations based on molecular dynamics free energy simulations. J. Phys. Chem. B. (1997) 101:5239–5248.

  21. Roux B. The influence of the membrane potential on the free energy of an intrinsic protein. Biophys. J. (1997) 73:2980–2989.[Web of Science][Medline]

  22. Banavali NK, Roux B. Atomic radii for continuum electrostatics calculations on nucleic acids. J. Phys. Chem. B. (2002) 106:11026–11035.

  23. Brooks BR, Bruccoleri RE, Olafson BD, States DJ, Swaminathan S, Karplus M. CHARMM: a program for macromolecular energy minimization and dynamics calculations. J. Comput. Chem. (1983) 4:187–217.[CrossRef][Web of Science]

  24. Ren P, Ponder JW. Consistent treatment of inter- and intramolecular polarization in molecular mechanics calculations. J. Comput. Chem. (2002) 23:1497–1506.[CrossRef][Web of Science][Medline]

  25. Phillips JC, Braun R, Wang W, Gumbart J, Tajkhorshid E, Villa E, Chipot C, Skeel RD, Kale L, Schulten K. Scalable molecular dynamics with NAMD. J. Comput. Chem. (2005) 26:1781–1802.[CrossRef][Web of Science][Medline]

  26. Marten B, Kim K, Cortis C, Friesner RA, Murphy RB, Ringnalda MN, Sitkoff D, Honig B. New model for calculation of solvation free energies: correction of self-consistent reaction field continuum dielectric theory for short-range hydrogen-bonding effects. J. Phys. Chem. (1996) 100:11775–11788.[CrossRef][Web of Science]

  27. DeLano W. The PyMOL Molecular Graphics System. (2002) San Carlos, CA: De Lano Scientific.

  28. Humphrey W, Dalke A, Schulten K. VMD: visual molecular dynamics. J. Mol. Graph. (1996) 14:33–38.[CrossRef][Web of Science][Medline]

  29. Nicholls A, Sharp KA, Honig B. Protein folding and association – insights from the interfacial and thermodynamic properties of hydrocarbons. Prot. Struct. Func. Genet. (1991) 11:281–296.[CrossRef]

  30. Sanner MF. Python: a programming language for software integration and development. J. Mol. Graph. Model. (1999) 17:57–61.[Web of Science][Medline]

  31. Philippsen A. DINO: Visualizing Structural Biology (2001).

  32. Dolinsky TJ, Nielsen JE, McCammon JA, Baker NA. PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatics calculations. Nucleic Acids Res. (2004) 32:W665–W667.[Abstract/Free Full Text]

  33. Miteva MA, Tuffery P, Villoutreix BO. PCE: web tools to compute protein continuum electrostatics. Nucleic Acids Res. (2005) 33:W372–W375.[Abstract/Free Full Text]

  34. Azuara C, Lindahl E, Koehl P, Orland H, Delarue M. PDB|Hydro: incorporating dipolar solvents with variable density in the Poisson-Boltzmann treatment of macromolecule electrostatics. Nucleic Acids Res. (2006) 34:W38–W42.[Abstract/Free Full Text]

  35. Berman HM, Battistuz T, Bhat TN, Bluhm WF, Bourne PE, Burkhardt K, Feng Z, Gilliland GL, Iype L, Jain S, et al. The protein data bank. Acta Crystallogr. D Biol. Crystallogr. (2002) 58:899–907.[CrossRef][Medline]

  36. Lomize MA, Lomize AL, Pogozheva ID, Mosberg HI. OPM: orientations of proteins in membranes database. Bioinformatics (2006) 22:623–625.[Abstract/Free Full Text]

  37. Gilson MK, Sharp K, Honig B. Calculating the electrostatic potential of molecules in solution: method and error assessment. J. Comput. Chem. (1987) 9:327–335.[CrossRef][Web of Science]

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

  39. Radhakrishnan I, Perez-Alvarado GC, Parker D, Dyson HJ, Montminy MR, Wright PE. Solution structure of the KIX domain of CBP bound to the transactivation domain of CREB: a model for activator:coactivator interactions. Cell (1997) 91:741–752.[CrossRef][Web of Science][Medline]

  40. Lamers MH, Perrakis A, Enzlin JH, Winterwerp HH, de Wind N, Sixma TK. The crystal structure of DNA mismatch repair protein MutS binding to a G x T mismatch. Nature (2000) 407:711–717.[CrossRef][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
 Proc. Natl. Acad. Sci. USAHome page
S. Lhuillier, M. Gallopin, B. Gilquin, S. Brasiles, N. Lancelot, G. Letellier, M. Gilles, G. Dethan, E. V. Orlova, J. Couprie, et al.
From the Cover: Structure of bacteriophage SPP1 head-to-tail connection reveals mechanism for viral DNA gating
PNAS, May 26, 2009; 106(21): 8507 - 8512.
[Abstract] [Full Text] [PDF]

This Article
Right arrow Abstract Freely available
Right arrow Print PDF (705K) Freely available
Right arrow Screen PDF (246K) Freely available
Right arrow Supplementary Data
Right arrowOA All Versions of this Article:
36/suppl_2/W270    most recent
gkn314v1
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 Jo, S.
Right arrow Articles by Im, W.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Jo, S.
Right arrow Articles by Im, W.
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?