OEBestOverlay¶

class OEBestOverlay

Note

This deprecated class will not be developed further. It is being replaced by OEOverlay and OEROCS which provides more control and greater expandability.

This class is used to optimize the overlap between 2 molecules. Both molecules can contain one or more conformers. The reference molecule is held rigid and the fit molecule orientation is calculated to maximize overlap. Each conformer of the fit molecule is optimized against each conformer of the reference molecule, resulting in a large number of results that can be returned. The OEBestOverlayResults and OEBestOverlayScore classes are used to make handling this large amount of results easier.

Additionally, the reference or fit shape can be a grid instead of a molecule.

Since an OEBestOverlay optimization can result in numerous final results, one for each starting position, the OEBestOverlayResults class is a container to hold all the results for a fit conformer against a reference conformer. This class is essentially a container for one or more OEBestOverlayScore instances. The actual number of scores depends on the OEBOOrientation used in the OEBestOverlay.

For OEBOOrientation::Inertial, there will be 4, 8 or 20 scores per ref-fit conformer pair, depending on the symmetry of each conformer. For OEBOOrientation::AsIs, there will be just 1 OEBestOverlayScore inside each OEBestOverlayResults instance. And for OEBOOrientation::Random, then there will be 1 OEBestOverlayScore for each of the N random starts.

When dealing with the results from OEBestOverlay, the user can either use a double loop, i.e. loop over all the OEBestOverlayResults instances for each ref-fit conformer pair then loop over each OEBestOverlayScore inside each OEBestOverlayResults instance. In this case, the results also come out in conformer order.

OEIter<OEBestOverlayResults> resiter = best.Overlay(fitmol);
for (;resiter;++resiter)
{
    OEIter<OEBestOverlayScore> scoreiter = resiter->GetScores();
    for (;scoreiter;++scoreiter)
    {
        // do something with score
    }
}

Or, there is a free function OESortOverlayScores, that takes an iterator of OEBestOverlayResults and returns a single, sorted iterator of OEBestOverlayScores that can be used in a single loop.

OEIter<OEBestOverlayResults> resiter = best.Overlay(fitmol);
OEIter<OEBestOverlayScore> scoreiter;
OESortOverlayScores(scoreiter, resiter, OEHighestTanimotoCombo());
for (;scoreiter;++scoreiter)
{
    // do something with score
}

Constructors¶

OEBestOverlay()
OEBestOverlay(const OEBestOverlay &)
OEBestOverlay(const OEChem::OEMCMolBase &refmol)
OEBestOverlay(const OESystem::OEScalarGrid &refgrid, float interpolate=0.5f)

An empty OEBestOverlay class instance can be created, or a new instance can take the reference molecule or grid as an argument. Note that by default, grids that have a resolution larger that 0.5 Å will be interpolated to that resolution.

operator=¶

OEBestOverlay &operator=(const OEBestOverlay &)

The assignment operator.

ClearColorForceField¶

void ClearColorForceField()

Clear out any color force field. No color scores will be calculated.

ClearUserStarts¶

void ClearUserStarts()

Set the number of User Starts back to zero.

GetCarbonRadius¶

float GetCarbonRadius() const

Return the current value for the carbon radius approximation.

GetInitialOrientation¶

unsigned int GetInitialOrientation() const

Get the current value for initial orientation. Possible values are from the OEBOOrientation namespace.

GetMaxRandomTranslation¶

float GetMaxRandomTranslation() const

Get the value in Å of maximum random translation.

GetMethod¶

unsigned int GetMethod() const

Return the current value of overlap method.

GetMinimizeType¶

unsigned int GetMinimizeType() const

Get the current value for minimize type. Possible values are from the OEBOMinType namespace.

GetNumRandomStarts¶

unsigned int GetNumRandomStarts() const

Get the current number of random starts.

GetNumUserStarts¶

unsigned int GetNumUserStarts() const

Get the number of user starts currently set. Note that these are not actually used unless OEBOOrientation::UserInertialStarts is also set.

GetRadiiApproximation¶

unsigned int GetRadiiApproximation() const

Return the current value of the radii approximation.

GetRandomSeed¶

unsigned int GetRandomSeed() const

Get the current value of the random seed.

GetRefGrid¶

const OESystem::OEScalarGrid *GetRefGrid() const

GetRefMol¶

const OEChem::OEMCMolBase *GetRefMol() const

GetRefSelfColor¶

float GetRefSelfColor()

Return the self color score of the reference molecule.

GetRefSymmetry¶

unsigned int GetRefSymmetry(unsigned int confIdx) const

Return the current representation level.

GetSymmetryThreshold¶

float GetSymmetryThreshold() const

GetUseHydrogens¶

bool GetUseHydrogens() const

Return the status of hydrogen use in OEBestOverlay.

GetUserStarts¶

bool GetUserStarts(float *xyz) const

Extract the coordinates of the user starts. xyz should be sized to 3 * GetNumUserStarts().

Overlay¶

OESystem::OEIterBase<OEBestOverlayResults> *
  Overlay(const OEChem::OEMCMolBase &fitmol)

Perform the calculation and return the an iterator of the results.

SetCarbonRadius¶

void SetCarbonRadius(float cradius)

Set the radius to use when using OEOverlapRadii::Carbon. By default this is set to 1.7 Å.

SetColorForceField¶

bool SetColorForceField(unsigned int type)
bool SetColorForceField(OEPlatform::oeistream &is)
bool SetColorForceField(const std::string &filename)
bool SetColorForceField(const OEColorForceField &cff)

Set the color force field to be used. Once set, color scores will be calculated and included in the results. Color gradients will not be included in the optimization unless specifically set using SetColorOptimize(true).

SetColorOptimize¶

void SetColorOptimize(bool state)

Add color gradients to shape gradients in the optimization. Has no effect unless a color force field is also set via OEBestOverlay::SetColorForceField.

SetInitialOrientation¶

void SetInitialOrientation(unsigned int orient)

Determines the initial orientation (starting position) for each optimization. The default is OEBOOrientation::Inertial. Alternatives are defined in the OEBOOrientation namespace.

SetMaxRandomTranslation¶

void SetMaxRandomTranslation(float trans)

If using random starts, this set the maximum distance (in Å) that the center of mass of the fit molecule will be moved away from the center calculated for inertial frame alignment.

SetMethod¶

void SetMethod(unsigned int method)

Set the method used to calculate overlap. The default for OEBestOverlay is OEOverlapMethod::Grid. Alternatives are defined in the OEOverlapMethod namespace.

SetMinimizeType¶

void SetMinimizeType(unsigned int type)

Set the score to use in the optimization. Options are in the OEBOMinType namespace.

SetNumRandomStarts¶

void SetNumRandomStarts(unsigned int n)

If SetInitialOrientation is set to OEBOOrientation::Random, this method sets the number of random starting positions that will be used.

SetRadiiApproximation¶

void SetRadiiApproximation(unsigned int type)

Set the radius approximation used to calculate overlap. The default for OEBestOverlay is OEOverlapRadii::Carbon. Alternatives are defined in the OEOverlapRadii namespace.

SetRandomSeed¶

void SetRandomSeed(unsigned int seed)

Set a random seed value to allow reproducible random searches.

SetRefGrid¶

bool SetRefGrid(const OESystem::OEScalarGrid &refgrid, float interpolate=0.5f)

Set a reference grid for the calculation. OEBestOverlay makes an internal copy. Pre-existing reference molecules or grids are replaced.

SetRefMol¶

bool SetRefMol(const OEChem::OEMCMolBase &refmol)

Set a reference molecule for the calculation. OEBestOverlay makes an internal copy. Pre-existing reference molecules or grids are replaced.

SetSymmetryThreshold¶

void SetSymmetryThreshold(float threshold)

SetUseHydrogens¶

void SetUseHydrogens(bool state)

Boolean to determine whether hydrogens are included in the shape calculation. By default this is false and hydrogens are ignored.

SetUserStarts¶

void SetUserStarts(const float *xyz, unsigned int nstarts)

Set a set of 3D coordinates for doing separate inertial starts. xyz should be 3 * nstarts. These are not used unless OEBOOrientation::UserInertialStarts is also set.