OESubSearchScreenType

Attention

This API is currently available in C++ and Python.

The OESubSearchScreen objects are typed bit-vectors. The type of a screen determines what kind of queries in mind the screen was designed.

In OEChem TK, a query molecule (OEQMolBase) can be constructed in three different ways.

Screen Type

Description

Initialize

See Also

OESubSearchScreenType::MDL

query read from MDL query file

OEReadMDLQueryFile

Substructure Search with MDL Queries

OESubSearchScreenType::Molecule

query built from molecule

BuildExpressions

OESubSearchScreenType::SMARTS

query built from SMARTS pattern

OEParseSmarts

SMARTS Pattern Matching

The expression trees of the queries that are initialized from different sources can be significantly different. For example, the “language” of SMARTS is the most expressive one allowing complex logical operations and recursive patterns. A SMARTS such as [NX3;H2,H1;!$(NC=O)] (primary or secondary amine that is not in amide) defines a substructure that can not be expressed using MDL query features. This also makes queries built from SMARTS the most complex ones. The other major difference between SMARTS and MDL queries is the way how they are expected to match aromaticity. For these reasons, OEChem TK provides three different screen types especially designed for the three different ways that are available to built query molecules using OEChem TK.

Note

For performance reason, OEChem TK does not support screening-based fast substructure search with explicit hydrogens. The OEMakeSubSearchTargetScreen function automatically suppresses explicit hydrogens of the molecule prior to the screen generation. The OEMakeSubSearchQueryScreen function will return false if the query molecule has any explicit hydrogen expression.

Default

The default screen type is OESubSearchScreenType::MDL.

MDL

The following code snippet shows how to create OESubSearchScreenType::MDL screens.

OEQMol qmol;
const auto queryopts = OEMDLQueryOpts::Default | OEMDLQueryOpts::SuppressExplicitH;
OEReadMDLQueryFile(ifs, qmol.QMol(), queryopts);

OESubSearchScreen qscreen;
if (!OEMakeSubSearchQueryScreen(qscreen, qmol, OESubSearchScreenType::MDL))
  std::cout << "query MDL screen cannot be generated" << std::endl;

OEGraphMol tmol;
OESmilesToMol(tmol, "c1c[nH]cc1");

OESubSearchScreen tscreen;
if (!OEMakeSubSearchTargetScreen(tscreen, tmol, OESubSearchScreenType::MDL))
  std::cout << "target MDL screen cannot be generated" << std::endl;

Molecule

The following code snippet shows how to create OESubSearchScreenType::Molecule screens. The query molecule is created from a “real” molecule using the BuildExpressions method.

OEQMol qmol;
OESmilesToMol(qmol, "c1cnccc1");
qmol.BuildExpressions(OEExprOpts::DefaultAtoms, OEExprOpts::DefaultBonds);

OESubSearchScreen qscreen;
if (!OEMakeSubSearchQueryScreen(qscreen, qmol, OESubSearchScreenType::Molecule))
  std::cout << "query Molecule screen cannot be generated" << std::endl;

OEGraphMol tmol;
OESmilesToMol(tmol, "c1ccc2c(c1)cccn2");

OESubSearchScreen tscreen;
if (!OEMakeSubSearchTargetScreen(tscreen, tmol, OESubSearchScreenType::Molecule))
  std::cout << "target Molecule screen cannot be generated" << std::endl;

Note

SMARTS

The following code snippet shows how to create OESubSearchScreenType::SMARTS screens.

OEQMol qmol;
OEParseSmarts(qmol, "C1CC[O,N]CC1");

// generate screens
OESubSearchScreen qscreen;
if (!OEMakeSubSearchQueryScreen(qscreen, qmol, OESubSearchScreenType::SMARTS))
  std::cout << "query SMARTS screen cannot be generated" << std::endl;

OEGraphMol tmol;
OESmilesToMol(tmol, "C1CCOCC1");

OESubSearchScreen tscreen;
if (!OEMakeSubSearchTargetScreen(tscreen, tmol, OESubSearchScreenType::SMARTS))
  std::cout << "target SMARTS screen cannot be generated" << std::endl;

Hint

In case of a query molecule of unknown source, the OEIsComplementaryScreenType function can be used to check whether the query molecule is suited for a specific screen type.