12 e = modeller.environ()
13 e.edat.dynamic_sphere =
False
14 e.libs.topology.read(
'${LIB}/top_heav.lib')
15 e.libs.parameters.read(
'${LIB}/par.lib')
16 modmodel = modeller.model(e)
17 modmodel.build_sequence(
'GGCC')
20 feat = modeller.features.distance(modmodel.atoms[0], modmodel.atoms[-1])
21 r = modeller.forms.gaussian(feature=feat, mean=10.0, stdev=1.0,
22 group=modeller.physical.xy_distance)
23 modmodel.restraints.add(r)
36 print m.evaluate(
False)