pdb2pqr

A pdb2pqr task can be used to determine the protonation states of residues based on their calculated pKa values and a specified pH. The only subdirective that is required is pH, which specifies the pH at which the pKa values are calculated. The task will use the PROPKA algorithm to calculate the pKa values and assign ionization states to residues in the task’s incoming structure. It is best to use this task after a psfgen task has processed the external structure file but before any solvation or membrane tasks.

Example

tasks:
  - ... (prior tasks, excluding solvation)
  - pdb2pqr:
      pH: 7.0
  - solvate:
  - ... (subsequent tasks follow)

How it works

This task first takes the incoming structure and prepares an input for pdb2pqr by stripping all hydrogens, naming all histidine residues as HIS and all water residues as HOH. Since this happens after a psfgen task, the input structure’s names are already in CHARMM format. The task then runs pdb2pqr with the PROPKA algorithm to calculate pKa values and assign ionization states to residues based on the specified pH. The output of the command is a PDB file that indicates which residues had their protonation states updated. This task also captures the summary table showing all residue pKa values.

The task then prepares for a new input script for psfgen that will apply all the protonation states as patches. This will use the incoming PSF/PDB combination, NOT the input or output of pdb2pqr. This task can handle the following patches:

  • NNEU: Neutral N-terminus

  • CNEU: Neutral C-terminus

  • TYRO: Deprotonated tyrosine

  • SERD: Deprotonated serine

  • ASPP: Protonate aspartic acid (ASP) on OD1 (ASPP)

  • GLUP: Protonate glutamic acid (GLU) on OE1 (GLUP)

  • LSN : Deprotonated lysine (atom HZ3 is removed)

  • RN2 : Deprotonated arginine (atom HH12 is removed)

Furthermore, the protonation states of the histidine residues are specified via in-segment mutations to the appropriate side-chain. By default, all histidines are interpreted as residue HSD (protonated on atom ND1). pdb2pqr assigns one of HSD, HSE, or HSP as residue names to each histidine residue, depending on the pKa value of the residue and the specified pH. The task will then apply the appropriate mutation to each histidine residue based on its assigned name. For example, if a histidine residue is assigned HSE, then the task will apply a mutation to change the residue from HSD to HSE. If a histidine residue is assigned HSP, then the task will apply a mutation to change the residue from HSD to HSP.

The task then runs VMD with the psfgen input script to apply the mutations and patches to the incoming PSF/PDB combination, and it outputs a new PSF/PDB combination that has the protonation states applied.

References

  • Very Fast Empirical Prediction and Rationalization of Protein pKa Values. Hui Li, Andrew D. Robertson and Jan H. Jensen. PROTEINS: Structure, Function, and Bioinformatics. 61:704-721 DOI 10.1002/prot.20660 (2005)

  • Very Fast Prediction and Rationalization of pKa Values for Protein-Ligand Complexes. Delphine C. Bas, David M. Rogers and Jan H. Jensen. PROTEINS: Structure, Function, and Bioinformatics 73:765-783 DOI 10.1002/prot.22102 (2008)

  • PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical pKa predictions. Mats H.M. Olsson, Chresten R. Sondergard, Michal Rostkowski, and Jan H. Jensen. Journal of Chemical Theory and Computation, 7(2):525-537 DOI 10.1021/ct100578z (2011)

  • Improved Treatment of Ligands and Coupling Effects in Empirical Calculation and Rationalization of pKa Values. Chresten R. Sondergaard, Mats H.M. Olsson, Michal Rostkowski, and Jan H. Jensen. Journal of Chemical Theory and Computation, 7(7):2284-2295 DOI 10.1021/ct200133y (2011)