File: MMPA.h

package info (click to toggle)
rdkit 201603.5-2
  • links: PTS, VCS
  • area: main
  • in suites: stretch
  • size: 72,364 kB
  • ctags: 18,217
  • sloc: cpp: 167,966; python: 58,855; java: 5,318; ansic: 5,239; sql: 1,908; yacc: 1,553; lex: 1,131; makefile: 418; xml: 229; sh: 192; fortran: 183; cs: 93
file content (40 lines) | stat: -rw-r--r-- 1,418 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
//
//  Copyright (C) 2015 Novartis Institutes for BioMedical Research
//
//   @@ 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.
//
#pragma once
#include <vector>
#include <string>
#include <stdexcept>
#include "../RDKitBase.h"

namespace RDKit {

namespace MMPA {
//! fragments a Molecule for processing with the Matched Molecular Pairs
//!  MMPA algorithm (Hussain et al)
/*!
  \param mol           Molecule to fragment
  \param result        Vector of Core and Sidechain results from the various
 cuts
  \param maxCuts        Maximum number of times to cut the molecule to generate
                        fragments.  A max cut of 3 will fragment with 1,2 and 3
                        cuts.
  \param maxCutBonds  Set the bond limit for determining which molecules
                        to analyze.  If a molecule has more than
                        this number of cutabble bonds, ignore.

 \return true if the molecule was fragmented, false otherwise.
*/

bool fragmentMol(const ROMol& mol,
                 std::vector<std::pair<ROMOL_SPTR, ROMOL_SPTR> >& result,
                 unsigned int maxCuts = 3, unsigned int maxCutBonds = 20,
                 const std::string& pattern = "[#6+0;!$(*=,#[!#6])]!@!=!#[*]");
}
}  // namespace RDKit