|
From: Martin H. <ma...@bl...> - 2011-08-22 10:06:20
|
Hi Troels
It's indeed possible.
The wizards are available through the cmd module.
# Initialize
load yourProtein
cmd.wizard("mutagenesis")
cmd.do("refresh_wizard")
# To get an overview over the wizard API:
for i in dir(cmd.get_wizard(): print i
# lets mutate residue 104 to GLN
cmd.get_wizard().set_mode("GLN")
cmd.get_wizard().do_select("104/")
# Select the rotamer
cmd.frame(11)
# Apply the mutation
cmd.get_wizard().apply()
There is alot more you could do inbetween these lines, but this is what
I use it like.
Martin
On 21.08.11 15:56, Troels Emtekær Linnet wrote:
> Hi.
>
> I wonder if it is possible to interact with the mutagenesis wizard
> from the command line?
>
> I am trying to determine best mutants for FRET labelling, by
> predicting the pka value/reactivity for a possible cysteine.
>
> I would like to make a script that loops over the residues in my
> protein, change the residue to a cysteine, save the molecule and then
> initiate propka <http://pymolwiki.org/index.php/Propka> to get the pka
> for the cysteine.
>
> I have googled myself to:
>
> fetch 4ins, async=0
> wizard mutagenesis
> cmd.get_wizard().do_select("/4ins//A/ASN`18/CB")
> cmd.get_wizard().do_pick("/4ins//A/ASN`18/CB")
>
> But I get:
> Error: please select an atom, not a bond.
>
> Can someone help here?
>
> Best
> Troels
>
>
> ------------------------------------------------------------------------------
> Get a FREE DOWNLOAD! and learn more about uberSVN rich system,
> user administration capabilities and model configuration. Take
> the hassle out of deploying and managing Subversion and the
> tools developers use with it. http://p.sf.net/sfu/wandisco-d2d-2
>
>
> _______________________________________________
> PyMOL-users mailing list (PyM...@li...)
> Info Page: https://lists.sourceforge.net/lists/listinfo/pymol-users
> Archives: http://www.mail-archive.com/pym...@li...
|