Rotamer Libraries in Pyrosetta

I am under the impression that this functionality is not particularly well-exposed for this sort of use, but you can certainly get at the existence of rotamers in PyRosetta. The classes you’ll be looking at are probably those in core/pack/dunbrack/. This code is moderately heavily templated (it’s very speed-optimized) so it may be tricky to use with pyrosetta.

The Dunbrack library isn’t prescriptive in the sense of “these rotamers exist” – it’s set up to have a score for all phi/psi/chi(1-4) combinations. Are you wanting to query the likelihood of rotamers you have, or have Rosetta generate a list of most likely rotamers (the latter duplicating what the packer does, which may be possible to extract into pyrosetta)?

What are you trying to do?

If you can find it in Rosetta src and it’s not in the app directory, 99+% of the time it will work in PyRosetta. You just got to know how to get it to work! What you want is to do:
from rosetta.core.pack.dunbrack import *

This will import all classes from core::pack::dunbrack namespace. You may also want to look in core.pack.rotamer_set. To only import certain classes you want to use (for example):
from rosetta.core.pack.dunbrack import RotamerLibrary

Have you been through the PyRosetta tutorials?

I don’t have much experience working or coding with individual rotamers within Rosetta, so I couldn’t help you there. I was able to construct a FixbbRotamerSets object and assume that most of the objects in core.pack could be constructed in PyRosetta, and hopefully all their functions should work. Make sure to check out the pyrosetta.org faq where it talks about how to construct vector1 objects and use the Rosetta options system within PyRosetta. You might want to use the 2010 Rotamer library, as there are non/semi-rotameric dihedrals not found in the 2002 library – but this also makes some code in Rosetta specific for the 2010 libary, such as the create_semi_rotameric_dunlib function in RotamerLibary just to forwarn you. If you go into a bit more specifics, we may be able to help you furthe.

Note that even in the Rotamer Libraries that you can download from the Dunbrack website, there are no xyz coordinates. They are dihedral angles for each chi. You can manually set an existing side chain’s dihedral angle using the pose.set_chi function. If you can’t get Rosetta to give you or make you a specific set of rotamers to try out or set on a pose, you may be able to do it manually in code through reading the textfile Rotamer Libary and writing a few functions in Python which would use the pose.set_chi function (in rosetta_database/rotamer or downloadable from the Dunbrack website). Just a thought.
Hope at least some of this is useful!

I can smooth your efforts by pointing out that the rotamer library itself is in rosetta_database/rotamer; bbdep02.May.sortlib is the 2002 library, and bbdep02.May.sortlib-correct.12.2010 is the 2010 library (the 2008 library was found to be flawed and deleted). You can surely get them directly from Roland’s lab, as well (http://dunbrack.fccc.edu/bbdep2010/) (and probably should, if you’re going to use it sans-Rosetta – he has a separate license)

These are non-prescriptive as I described before. Rosetta prescribes rotamers by setting a probability cutoff (one of the obvious columns in the file is probability); I guess you can do something similar?

bbdep02.May.sortlib-correct.12.2010 is *NOT* the Dunbrack 2010 library. It’s the “-correct” adjustments (from Dec 2010) to the Dunbrack 2002 library.

The Dunbrack 2010 library is contained in rosetta_database/rotamer/ExtendedOpt1-5 (at least for the most recent versions of PyRosetta )

Regarding generating xyz coordinates for rotamers, I’d recommend looking at core.pack.rotamer_set.RotamerSetFactory and core.pack.rotamer_set.RotamerSet, particularly the build_rotamers(), num_rotamers() and rotamer() member functions of the latter. Note that for the build_rotamers function you’ll need a PackerTask, with which you can set the sampling level (also called the “ex level”) for the rotamers. The default settings should get you just the “on rotamer” conformation. — Technically, you won’t get *all* the rotamers, just the ones that have reasonable probabilities. These can be controlled with the -packing:dunbrack_prob_buried, -packing:dunbrack_prob_nonburied and -packing:dunbrack_prob_nonburied_semirotameric command line flags (I don’t know how to change these from within PyRosetta), with closer to 1.0 being more rare rotamers (defaults are 0.98, 0.95, 0.95, respectively).

Great! Please let us know if you run into any problems. Some bindings of specific variable types (templates) are not directly available in PyRosetta, so you may run into issues on the Rotamer set side. If you do, Sergey Lyskov MAY be able to configure those bindings to work, as he has for the common vector1 objects, C++ std:pair, etc. Then, you would need to wait for the updated binaries, which get released from the Gray lab pretty much every week.

There is a great deal of additional information on the text file at our website, http://dunbrack.fccc.edu/bbdep2010/Tutorial.php , which should help if you have not been there before. If you decide to use the text file separately then the rotamers that are built in Rosetta, the academic licence as Steven suggested on the site can be filled out pretty quickly.