You are here

Random perturbations to peptide backbone

1 post / 0 new
Random perturbations to peptide backbone
#1

Dear all, 

I am trying to explore peptide local conformations bound to its receptor. I was planning to start by using the peptide docking impletementation in pyrosetta (D100_Docking.py), which works well, except for two things (that are specific for my problem). First, my pose comes from a minimized structure (using relax) so the high resolution protocol doesn't move the peptide at all, It must be because it is already in a deep local minimum and the protocol is designed for guiding the conformations to exactly these kind of enegetic spots. Second, the ab intio version (which allows to escape from minima) makes too drastic changes to the pose, and I only want to explore the local conformational space. So, my plan is to use a small initial random perturber for the peptide and then let the high resolution protocol solve the energy minimization to find new local minimas around my intial configuration. I am looking to use the mover:

pyrosetta.rosetta.protocols.simple_moves.SetTorsion

it has the ability of perform small random perturbations while setting up a temporary fold_tree at an specified residue (which is just great!). I tested on RosettaScripts and works well enough to give it a try. My problem is that is has been difficult to apply in PyRosetta, because I cannot find were to set up which residue(s) to perturb. It doesn't accept any movemap or list of residues, that I am aware of, and I know the Pose() object doesn't have one either. Could anyone be so kind to indicate me how can I execute this mover on my pose? My code now looks something like this:

...
setTorsion = protocols.simple_moves.SetTorsion()
setTorsion.add_perturbation_magnitude(2.0)
setTorsion.add_perturbation_type('gaussian')
setTorsion.set_fold_tree_root(residues[0])
...
setTorsion.apply(test_pose)
...

But nothing happens to the pose. 

Category: 
Post Situation: 
Mon, 2018-12-10 17:16
Martin Floor