Hello, I have some 2D sdf structure and I need to generate 3D conformer, but
when I used OBBuilder class, it generated a 3D conformer with high energy. I
used ForceField to lower the energy but it seem not to get there yet. Here
is my script:
import pybel
import openbabel
conv = openbabel.OBConversion()
mol = openbabel.OBMol()
conv.SetInAndOutFormats("sdf","sdf")
conv.ReadFile(mol,"4.sdf")
builder = openbabel.OBBuilder()
builder.Build(mol)
conv.Convert()
ff = openbabel.OBForceField.FindForceField("mmff94")
ff.Setup(mol)
print ff.Energy()
ff.WeightedRotorSearch(200,25)
print ff.Energy()
ff.UpdateCoordinates(mol)
conv.WriteFile(mol,"4b.sdf")
I tried to use SystematicRotorSearch() but it took long time for molecule
with 5 rotatable bonds.
Please comment how to modify it.
Best wishes,
Phuong Do
--
View this message in context:
http://forums.openbabel.org/Generate-low-energy-3D-conformer-with-python-tp4659929.html
Sent from the General discussion mailing list archive at Nabble.com.
------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
OpenBabel-discuss mailing list
OpenBabel-discuss@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/openbabel-discuss
