1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
|
//
// Copyright (C) 2003-2006 Rational Discovery LLC
//
// @@ All Rights Reserved @@
// This file is part of the RDKit.
// The contents are covered by the terms of the BSD license
// which is included in the file license.txt, found at the root
// of the RDKit source tree.
//
#include <RDGeneral/export.h>
#ifndef _RD_FRAG_CATALOG_UTILS_H_
#define _RD_FRAG_CATALOG_UTILS_H_
#include <GraphMol/Subgraphs/Subgraphs.h>
#include <GraphMol/RDKitBase.h>
#include <GraphMol/Substruct/SubstructMatch.h>
#include "FragCatParams.h"
#include <iostream>
namespace RDKit {
// get the functional groups from file or stream
// each functional groups is read in as a molecule with queryatoms and
// querybonds
RDKIT_FRAGCATALOG_EXPORT MOL_SPTR_VECT readFuncGroups(std::string fileName);
RDKIT_FRAGCATALOG_EXPORT MOL_SPTR_VECT readFuncGroups(std::istream &inStream, int nToRead = -1);
// REVIEW: should this return a vector of pairs or a map?
// mark the functional groups of interest on the molecule
// and return a vector os std::pair <aid, fid>
// aid - is the atom id in mol that connect to a functional (fid)
// fid - the functional groups in the list maintained in params
// ARGUMENTS:
// mol - molecule of interest
// params - fragment catalog paramter object (contains a list of functional
// groups of interest
// fgBonds - container for bondIds in mol that are part of the functional
// groups
// the connection bond is included. these need to be chopped from
// the molecule later
RDKIT_FRAGCATALOG_EXPORT MatchVectType findFuncGroupsOnMol(const ROMol &mol, const FragCatParams *params,
INT_VECT &fgBonds);
// This functions is called before either adding the fragments from a molecule
// to a fragment catalog or generating the fincgerprint for this molecule
// using a fragment catalog. These are the things this function does
// - recognize the function groups (and their location) on the molecule
// - chop these functional groups of the molecule to create a core molecule
// "coreMol"
// - map the function group locations onto this "coreMol" (bacause the atom ids
// on coreMol are different from the original molecule
// - return coreMol to the caller of this function and the enter the atom ids to
// func
// group ids mapping into aToFmap argument
RDKIT_FRAGCATALOG_EXPORT ROMol *prepareMol(const ROMol &mol, const FragCatParams *fparams,
MatchVectType &aToFmap);
}
#endif
|