1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
|
# This example mutates every residue of a chain to alanine.
# load a structure
cmd.load("$TUT/1hpv.pdb")
# activate the mutagensis wizard
cmd.wizard("mutagenesis")
# set target residue type
cmd.get_wizard().set_mode("ALA")
from pymol import stored
stored.list = []
# generate our list of residues (CA-atoms) to mutate
cmd.iterate("1hpv//A//CA",
"stored.list.append(string.join([model,segi,chain,resi,name],'/'))")
# now iterate through each residue
for sele in stored.list:
# pick a residue
cmd.edit(sele)
# notify the wizard about the picked residue
cmd.get_wizard().do_pick(0)
# apply the mutation
cmd.get_wizard().apply()
# update the screen, showing the mutation
cmd.indicate("byres "+sele)
cmd.refresh()
# now close the wizard
cmd.set_wizard()
cmd.deselect()
|
