Simple 3D Conformer Generation with Smi23D

Rich Apodaca — Wed, 12 Dec 2007 10:32:00 +0000

Three-dimensional conformer generation is a common problem in cheminformatics. The most convenient and generally-useful method for creating chemical structures is the 2D chemical structure editor; applications that require three-dimensional representations need a way to generate reasonable coordinates from 2D user input. Until recently, there were no options for doing so with Open Source software. This article shows how the Open Source package smi23d can be used to convert ordinary SMILES strings into three-dimensional molfile representations.

About smi23d

smi23d uses a two-stage process to generate 3D coordinates.; an initial pass with smi2sdf generates rough coordinates and subsequent refinement by mengine results in the final coordinates. The package was originally written in C by Kevin Gilbert and updated by Rajarshi Guha. As part of what appears to be a growing trend in cheminformatics, smi23d is licensed under the highly-permissive Apache License.

On a related note, the source code for a program called Frog is reportedly on its way into the Open Babel project.

Prerequisites

To build smi23d, you'll need to install Scons, a Make-like build utility written in Python. I was able to install the Scons rpm on my Linux system without a problem. smi23d uses no other dependencies.

Download smi23d

smi23d can be downloaded with Subversion:

$ svn co https://cicc-grid.svn.sourceforge.net/svnroot/cicc-grid/cicc-grid/smi23d/trunk smi23d

Building smi23d

With the source code in place, compilation is just a matter of running Scons:

$ cd smi23d
$ scons
...

Once the sources are compiled, we'll want to configure our system a bit:

$ cd build
$ ls
mmff94.prm  mmxconst.prm

$ cp ../src/smi2sdf/smi2sdf .
$ cp ../src/mengine/mengine .

The two files mmff94.prm and mmxconst.prm are parameter files needed by both smi2sdf and mengine.

With smi2sdf and mengine both in the build directory, we can create a simple test with the SMILES for Ivabradine:

$ vi test.smi
...

$ cat test.smi
CN(CCCN1CCC2=CC(=C(C=C2CC1=O)OC)OC)C[C@H]3CC4=CC(=C(C=C34)OC)OC

With everything ready to go, we can begin Stage one:

$ ./smi2sdf test.smi
Found 1 structures in test.smi
 field : MMX
 Atom Types: 169
 Bonds: 580 Bond3: 0 Bond4: 0 Bond5: 0
 Angle: 434 Angle3: 41 Angle4: 60 Angle5: 0
 Torsion: 697  Torsion4: 58 Torsion5: 0
 Vdw: 172 OOP: 91 Dipole: 474 Charge: 0 Improper: 0
 STBN: 26 ANGANG: 0 STRTOR: 0 VDWPR: 4


Input file  = test.smi
Output file = output.sdf
Param file  = mmxconst.prm
Log file    = error.log
Inorganic file = test_inorg.smi

Structure: 0 CN(CCCN1CCC2=CC(=C(C=C2CC1=O)OC)OC)C[C@H]3CC4=CC(=C(C=C34)OC)OC

You can view the result in an application like Jmol:

It's not much to look at, but we're not quite done yet.

Stage two is accomplished by using the output of Stage one as input to mengine:

$ ./mengine -o optimized.sdf output.sdf
 field : MMX
 Atom Types: 169
 Bonds: 580 Bond3: 0 Bond4: 0 Bond5: 0
 Angle: 434 Angle3: 41 Angle4: 60 Angle5: 0
 Torsion: 697  Torsion4: 58 Torsion5: 0
 Vdw: 172 OOP: 91 Dipole: 474 Charge: 0 Improper: 0
 STBN: 26 ANGANG: 0 STRTOR: 0 VDWPR: 4
 field : MMFF94
 Atom Types: 181
 Bonds: 448 Bond3: 0 Bond4: 0 Bond5: 0
 Angle: 1801 Angle3: 21 Angle4: 61 Angle5: 0
 Torsion: 674  Torsion4: 38 Torsion5: 95
 Vdw: 182 OOP: 112 Dipole: 0 Charge: 0 Improper: 0
 STBN: 286 ANGANG: 0 STRTOR: 0 VDWPR: 0

We now have a file called output.sdf. As you can see, it's a pretty good 3D representation of Ivabradine:

Conclusions

In this tutorial, we've seen how the Open Source program smi23d can be used to assign reasonable 3D coordinates to an arbitrary SMILES string. One very practical use of smi23d would be to process the output of 2D chemical structure editors prior to use in a 3D program. Future articles will discuss some of the possibilities.