File: FragCatalogUtils.h

package info (click to toggle)
rdkit 201809.1%2Bdfsg-6
  • links: PTS, VCS
  • area: main
  • in suites: buster
  • size: 123,688 kB
  • sloc: cpp: 230,509; python: 70,501; java: 6,329; ansic: 5,427; sql: 1,899; yacc: 1,739; lex: 1,243; makefile: 445; xml: 229; fortran: 183; sh: 123; cs: 93
file content (60 lines) | stat: -rw-r--r-- 2,492 bytes parent folder | download
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