Before you turn this problem in, make sure everything runs as expected. First, restart the kernel (in the menubar, select Kernel$\rightarrow$Restart) and then run all cells (in the menubar, select Cell$\rightarrow$Run All).
Make sure you fill in any place that says YOUR CODE HERE
or "YOUR ANSWER HERE", as well as your name and collaborators below:
NAME = ""
COLLABORATORS = ""
Keywords: pose_from_rcsb(), pdb2pose(), EMapVector()
# Notebook setup
import sys
if 'google.colab' in sys.modules:
!pip install pyrosettacolabsetup
import pyrosettacolabsetup
pyrosettacolabsetup.setup()
print ("Notebook is set for PyRosetta use in Colab. Have fun!")
from pyrosetta import *
from pyrosetta.teaching import *
init()
# From previous section:
sfxn = get_score_function(True)
Analyze the energy between residues Y102 and Q408 in cetuximab (PDB code 1YY9, use the pyrosetta.toolbox.pose_from_rcsb
function to download it and load it into a new Pose
object) by following the steps below.
A. Internally, a Pose object has a list of residues, numbered starting from 1. To find the residue numbers of Y102 of chain D and Q408 of chain A, use the residue chain identifier and the PDB residue number to convert to the pose numbering using the pose2pdb()
method:
pose = pyrosetta.toolbox.pose_from_rcsb("1YY9")
res102 = pose.pdb_info().pdb2pose("D", 102)
res408 = pose.pdb_info().pdb2pose("A", 408)
# get the pose numbers for Y102 (chain D) and Q408 (chain A)
# YOUR CODE HERE
raise NotImplementedError()
B. Score the pose and determine the van der Waals energies and solvation energy between these two residues. Use the following commands to isolate contributions from particular pairs of residues, where rsd102
and rsd408
are the two residue objects of interest from above (not the residue number -- use pose.residue(res_num)
to access the objects):
emap = EMapVector()
sfxn.eval_ci_2b(pose.residue(res102), pose.residue(res408), pose, emap)
print(emap[fa_atr])
print(emap[fa_rep])
print(emap[fa_sol])
# YOUR CODE HERE
raise NotImplementedError()