GPU Accelerated Side-Chain Optimization
Amino acid side-chain optimization algorithms that utilize discrete rotamer libraries of optimal side-chain conformations are used to find the global minimum energy conformation of a protein. Under a many-body energy expansion, Force Field X can optimize side-chain conformations up to four-body side-chain interactions. Using this aggressive many-body expansion is computationally expensive; through the use of Graphical Processing Units (GPUs) and the OpenMM software package, Force Field X is capable of computing the energy of side-chain conformations 25x faster than the first implementation of side-chain optimization.
Best Practices Overview
The following example demonstrates global optimization of Chignolin side-chains (PDB 1UAO). The best practices protocol is as follows:
- Minimize the model to eliminate obvious steric clashes.
- Optimize the model using all-with-rotamer-elimination global optimization. If the protein is large, use the sliding box method approximate global optimization.
- Minimize the model to allow rigid rotamer conformations to establish ideal intermolecular contacts.
Minimization of a Short Peptide
To perform a minimization on the peptide, use the command below:
ffxc Minimize 1uao.pdb
Output from the minimization will show each iteration of the algorithm and the resulting energy of the protein model. Pictured below are iterations 76-80 and the energy output. This simulation should take less than one minute of computing time on a standard laptop.
76 -242.8328 1.1926 0.3876 0.0032 75.14 77 0.025 77 -243.3192 1.3723 0.4864 0.0041 67.24 78 0.027 78 -243.7189 1.1017 0.3997 0.0068 75.87 79 0.034 79 -243.8713 1.4691 0.1524 0.0028 66.85 80 0.035 80 -244.1078 0.8571 0.2365 0.0010 60.08 81 0.031 Optimization achieved convergence criteria: 0.85712 Computed Potential Energy Bond Stretching 9.04900448 141 0.000 ( 0.01322) Angle Bending 50.63504369 249 0.000 ( 3.54833) Stretch-Bend -0.39314453 229 0.000 Out-of-Plane Bend 4.38249086 108 0.000 Torsional Angle 8.76885188 369 0.000 Pi-Orbital Torsion 2.77304417 25 0.000 Torsion-Torsion -4.58642973 7 0.000 Van der Waals 63.64719498 9063 0.003 Atomic Multipoles -316.37910863 9453 Polarization -62.00472529 9453 0.037 Total Potential -244.10777813 (Kcal/mole) 0.042 (sec) Saving 1uao.pdb_2
The resulting PDB file is automatically saved as 1uao.pdb_2.
Side-Chain Optimization of a Short Peptide Without a GPU
To perform an all-with-rotamer-elimination side-chain optimization on the minimized peptide without use of a GPU, use the command below:
ffxc ManyBody -a 2 -l 2 -x 1 1uao.pdb_2
Output from the ManyBody algorithm will show energy evaluations of self and pair side-chain energies along with a printout overall system energy. The final rotamers and overall energy are shown below. This simulation should take less than one minute of computing time on a standard laptop.
Final rotamers: -------------------------------------------------------------------------------------------- Residue | Chi 1 | Chi 2 | Chi 3 | Chi 4 | Energy | -------------------------------------------------------------------------------------------- 1 A ( 2-TYR, 0) | -172.5 0.0 | 90.2 0.0 | 114.2 0.0 | 0.0 0.0 | -19.4350 | 2 A ( 3-ASP, 0) | 176.3 0.0 | 0.0 0.0 | 0.0 0.0 | 0.0 0.0 | -28.3183 | 3 A ( 5-GLU, 0) | -165.5 0.0 | -84.1 0.0 | 22.3 0.0 | 0.0 0.0 | 9.5233 | 4 A ( 6-THR, 5) | 62.0 0.0 | -60.0 0.0 | 0.0 0.0 | 0.0 0.0 | -8.2411 | 5 A ( 8-THR,12) | -175.0 0.0 | -120.0 0.0 | 0.0 0.0 | 0.0 0.0 | -3.3993 | 6 A ( 9-TRP, 0) | -41.5 0.0 | 87.5 0.0 | 0.0 0.0 | 0.0 0.0 | 5.9391 | -------------------------------------------------------------------------------------------- Final Minimum Energy Computed Potential Energy Bond Stretching 8.77023477 141 0.000 ( 0.01299) Angle Bending 50.20607206 249 0.000 ( 3.53036) Stretch-Bend -0.50772738 229 0.000 Out-of-Plane Bend 4.40102720 108 0.000 Torsional Angle 8.77155306 369 0.000 Pi-Orbital Torsion 2.76777372 25 0.000 Torsion-Torsion -4.58901318 7 0.000 Van der Waals 71.10186683 9063 0.000 Atomic Multipoles -331.31181960 9453 Polarization -63.51240025 9453 0.012 Total Potential -253.90243277 (Kcal/mole) 0.014 (sec) Saving 1uao.pdb_3
The resulting PDB file is automatically saved as 1uao.pdb_3.
Alternatively, the previous command can be excuted using OpenMM. Please be sure OpenMM is configured by following the instructions here: OpenMM Installation.
ffxc ManyBody -a 2 -l 2 -x 1 1uao.pdb_2 -Dplatform=omm
Optimization of a Deafness Related Protein With a GPU
This next example is based on the deafness related COCH protein. Start by installing OpenMM by following the instructions here: OpenMM Installation. If OpenMM is not installed, the program will run without a GPU. Next, download the coordinate file and rename it to COCH.pdb (PDB COCH). Perform a minimization and a GPU accelerated side-chain optimization by using the commands below:
ffxc Minimize COCH.pdb
ffxc ManyBody -Dplatform=OMM -a 2 -l 2 -tC 1000 -x 1 -p 1 COCH.pdb_2
When using a GPU with the OpenMM platform, energy evaluations on side-chains will occur approximately 25x faster than without a GPU. This simulation should finish within one hour on a GPU. The resulting PDB file is automatically saved as COCH.pdb_3.
A description of all of the flags used in the ManyBody command is below:
Command | Description |
---|---|
-a 2 | Algorithm using Goldstein Elimination. |
-l 2 | Richardson rotamer library. |
-tC 1000 | Cutoff distance for two-body side-chain interactions (1000 Angstroms). |
-x 1 | Start simulation at first amino acid in protein. |
-p 1 | Prune side-chain conformations that clash with protein backbone. |
-Dplatform=OMM | Use OpenMM to run the algorithm on a GPU. |