| <<< Previous | Next >>> | |
| Introduction |
APBSCommand Tutorial |
Reference |
In this section, we walk through the determination of the electrostatic
potential of a protein selected from the PDB. There are three steps; creating a PQR file, setting up the calculation in PMV and executing the APBS calculation.
1. Creating a PQR file
Suppose we are interested in Rubredoxin, PDB ID 1CAA. APBS requires the charge and radius information for every atom so currently APBSCommand requires that a pqr file be loaded. The PDB2PQR service is very helpful in generating pqr files. Go to this site and provide the PDB ID (or upload your PDB of interest), select the AMBER, CHARMM, or PARSE parameter set and click Submit. It should be noted that although the PARSE [Sitkoff et al., 1994] parameters have been optimized for continuum calculations (see Further Reading for more details) the 0 radii contained therein are currently incompatible with the cubic spline surfaces (srfm=spl2). The molecular or harmonic smoothed surfaces (srfm=mol/smol) should be used with PARSE. Save the returned file to an appropriate name such as 1caa.pqr, and take note of where you put it. We are now prepared to set up an APBS calculation.
2. Calculation setup
Open PMV and click APBS, Setup, Calculation. There are three windows in the APBSCommand GUI.
|
The Calculation Window
The calculation window sets the calculation type and conditions, selects the molecule, and specifies the output. Click Select Molecule 1..., Load Molecule, select 1caa.pqr from the folder you saved it to, then OK. All of the files created in the APBS run will be stored in the Project folder; unique project folders within the current working directory are automatically used if no others are specified. By default, the electrostatic potential and solvent accessible surface files will be created in OpenDX format in the project folder. This can be changed by clicking on Output files but these will be sufficient for our purposes as will the parameters in the Mathematics section. Note that you can save the Grid, Physics and Mathematics settings for future calculations or load settings from previous calculations. |
![]() |
|
The Grid Window
The grid window allows the user to specify the grid dimensions, resolution, and center while visualizing the grid. Click on Show Coarse/Fine Grid then Autocenter and Autosize. If you have the psize script in your path you can click on Psize Auto Params to generate the recommended grid parameters. Be sure that the coarse grid always extends significantly beyond the molecule in all dimensions and that the coarse grid is at least 1 Å larger than the fine grid. Since Rubredoxin is so small you will need to adjust one of your grids after using Psize Auto Params. Also be sure that the machine you are working on has enough memory to work with the grid you are interested in. For the best performance, you should have a bit more physical memory free than the Memory to be allocated indicated at the bottom of the Grid tab (otherwise your operating system will use your hard disk as extra memory, slowing down the calculation). |
![]() |
|
The Physics Window
The physics window allows the user to set additional physical parameters and to add ions to their calculation. It is important to note that quantitative results (e.g., solvation and binding energies) will depend sensitively on the protein dielectric, and the surface definition. The commonly used solvent probe radius is 1.4 Å for solvent and 0.0 Å for vacuum calculations. The solution is by default ion-free. To add ions, click Add, and then specify the properties of your ion. Make sure your solution is neutral. We have finished setting up the calculation. |
![]() |
| <<< Previous | APBSCommand | Next >>> |
| Introduction | Reference |