Molecule Alignment

The following two code examples demonstrate how to align molecules either based on their maximum common substructure (Listing 1) or based on substructure matches (Listing 2).

Molecule Alignment Based on MCS

The Figure: Example of depiction without alignment is generated by highlighting the maximum common substructure of two molecules (MCS), but keeping the original orientation of the compounds.

../_images/MCSNoAlignment.png

Example of depiction without alignment

The Listing 1 code example shows how to align these two molecules by their maximum common substructures (i.e. by the highlighted atoms and bonds in Figure: Example of depiction without alignment).

  • First, two molecules are initialized from SMILES string and prepared for depiction.
  • A OEMCSSearch object is initialized with the reference 2D molecule onto which the fit molecule will be aligned.
  • An OEImageGrid object is constructed that allows to render the two molecules in a grid next to each other.
  • The fit molecule is then aligned to the reference by calling the OEPrepareAlignedDepiction function that performs MCS search and then aligns the fit molecule to the reference based on the detected common substructure(s).
  • An OE2DMolDisplayOptions object is initialized that stores the properties determine how the molecules are rendered. In order to render the molecules in equal size the smallest depiction scaling factor is calculated.
  • If the molecule alignment was successful (i.e OEAlignmentResult::IsValid is true), then the correspondence between the two molecules stored in the OEAlignmentResult returned by the OEPrepareAlignedDepiction function. This OEAlignmentResult object can be used to highlight the matched common substructure by invoking the OEAddHighlighting function.
  • After both molecules are rendered into the separate cells of the grid layout, the image is written into a png file.

After initializing the two molecules and the OEMCSSearch object, the OEMCSSearch::Match method is called to perform the search and return an iterator over the identified maximum common substructures. The first match is then utilized to call the OEPrepareAlignedDepiction function that aligns one molecule to the other based on the given match. After the alignment, the common atoms and bonds are highlighted by invoking the OEAddHighlighting function. The image created by Listing 1 is shown in Figure: Example of depiction with alignment based on MCS.

Listing 1: Example of molecule alignment based on MCS

#include <openeye.h>

#include <oesystem.h>
#include <oechem.h>
#include <oedepict.h>

using namespace OESystem;
using namespace OEChem;
using namespace OEDepict;

int main()
{
  OEGraphMol refmol;
  OESmilesToMol(refmol, "c1cc(c2cc(cnc2c1)CCCO)C(=O)CCO");
  OEPrepareDepiction(refmol);

  OEGraphMol fitmol;
  OESmilesToMol(fitmol, "c1cc2ccc(cc2c(c1)C(=O)O)CCO");
  OEPrepareDepiction(fitmol);

  OEMCSSearch mcss(OEMCSType::Approximate);
  unsigned int atomexpr = OEExprOpts::DefaultAtoms;
  unsigned int bondexpr = OEExprOpts::DefaultBonds;
  mcss.Init(refmol, atomexpr, bondexpr);
  mcss.SetMCSFunc(OEMCSMaxBondsCompleteCycles());

  OEAlignmentResult alignres = OEPrepareAlignedDepiction(fitmol, mcss);

  OEImage image(400, 200);

  unsigned int rows = 1;
  unsigned int cols = 2;
  OEImageGrid grid(image, rows, cols);

  OE2DMolDisplayOptions opts(grid.GetCellWidth(), grid.GetCellHeight(), OEScale::AutoScale);
  opts.SetTitleLocation(OETitleLocation::Hidden);

  const double refscale = OEGetMoleculeScale(refmol, opts);
  const double fitscale = OEGetMoleculeScale(fitmol, opts);
  opts.SetScale(std::min(refscale, fitscale));

  OE2DMolDisplay refdisp(mcss.GetPattern(), opts);
  OE2DMolDisplay fitdisp(fitmol, opts);

  if (alignres.IsValid())
  {
    OEAtomBondSet refabset(alignres.GetPatternAtoms(), alignres.GetPatternBonds());
    OEAddHighlighting(refdisp, OEBlueTint, OEHighlightStyle::BallAndStick, refabset);

    OEAtomBondSet fitabset(alignres.GetTargetAtoms(), alignres.GetTargetBonds());
    OEAddHighlighting(fitdisp, OEBlueTint, OEHighlightStyle::BallAndStick, fitabset);
  }

  OEImageBase* refcell = grid.GetCell(1, 1);
  OERenderMolecule(*refcell, refdisp);

  OEImageBase* fitcell = grid.GetCell(1, 2);
  OERenderMolecule(*fitcell, fitdisp);

  OEWriteImage("MCSAlign.png", image);
  return 0;
}
../_images/MCSAlign.png

Example of depiction with alignment based on maximum common substructure

See also