| Title: | R interface to a subset of OpenBabel functionalities |
|---|---|
| Description: | ChemmineOB provides an R interface to a subset of cheminformatics functionalities implemented by the OpelBabel C++ project. OpenBabel is an open source cheminformatics toolbox that includes utilities for structure format interconversions, descriptor calculations, compound similarity searching and more. ChemineOB aims to make a subset of these utilities available from within R. For non-developers, ChemineOB is primarily intended to be used from ChemmineR as an add-on package rather than used directly. |
| Authors: | Kevin Horan, Thomas Girke |
| Maintainer: | Thomas Girke <[email protected]> |
| License: | Artistic-2.0 |
| Version: | 1.43.0 |
| Built: | 2024-07-12 06:20:11 UTC |
| Source: | https://github.com/bioc/ChemmineOB |
Computes a re-arrangement required to transform the atom numbering of the given compound into the canonical atom numbering. This function uses the OBGraphSym and CanonicalLabels classes of Open Babel to compute the re-arrangement.
canonicalNumbering_OB(obmolRefs)canonicalNumbering_OB(obmolRefs)
obmolRefs |
A list of OBMol references ( of class '_p_OpenBabel__OBMol') representing the compounds. |
A list of vectors of index values. Each item in the list corresponds to one of the given compounds. The values of a list item are the re-arrangement of the atoms. For example, if the value in item 1, column 1 is 25, that means that atom number 1 in the original compound should become atom number 25 in the canonical version of that compound.
Kevin Horan
http://openbabel.org/api/2.3/canonical_code_algorithm.shtml
convertFormat to return a new compound in the canonical format.
## Not run: library(ChemmineR) data(sdfsample) labels = canonicalNumbering_OB(obmol(sdfsample[[1]])) ## End(Not run)## Not run: library(ChemmineR) data(sdfsample) labels = canonicalNumbering_OB(obmol(sdfsample[[1]])) ## End(Not run)
Converts compound data from one format to another. Some formats are not supported on windows. This currently includes "INCHI" and "INCHIKEY".
convertFormat(from, to, source,options=data.frame(names="gen2D",args=""))convertFormat(from, to, source,options=data.frame(names="gen2D",args=""))
from |
The format that |
to |
The format to convert |
source |
The initial compound format, as a string. The format of the string should be identical to the file format of the same name. Tabs and newlines may be represented with \t and \n, repsectivly. |
options |
This is a data frame listing OpenBabel Options that should be applied while converting each compound. The "names" column should be the name of each option and the "args" column sould be the arguments to the corresponding option. For example, the default value shown above will cause 2D coordinates to be generated from the input compound and saved in the output compound, assuming the ouptut format supports it. A full list of available options can be found by executing "obabel -L ops" at the command line. The current list is: 0xout <file.xxx> Additional file output AddInIndex Append input index to title AddPolarH Adds hydrogen to polar atoms only canonical Canonicalize the atom order energy ForceField Energy Evaluation (not displayed in GUI) fillUC <param> Fill the unit cell (strict or keepconnect) gen2D Generate 2D coordinates gen3D Generate 3D coordinates minimize ForceField Energy Minimization (not displayed in GUI) partialcharge <method> Calculate partial charges by specified method readconformer Adjacent conformers combined into a single molecule s Isomorphism filter(-s, -v options replacement)(not displayed in GUI) sort <desc> Sort by descriptor(~desc for reverse) unique [param] remove duplicates by descriptor;default inchi v Isomorphism filter(-s, -v options replacement)(not displayed in GUI) |
Returns the compound given in source in the format specified by to.
Kevin Horan
OpenBabel http://openbabel.org
#create an SDF from a SMILES string and put the atoms in the canonical order ## Not run: sdfStr = convertFormat("SMI","SDF","CC(=O)OC1=CC=CC=C1C(=O)O\ttest_name", options=data.frame(names="canonical",args="")) ## End(Not run)#create an SDF from a SMILES string and put the atoms in the canonical order ## Not run: sdfStr = convertFormat("SMI","SDF","CC(=O)OC1=CC=CC=C1C(=O)O\ttest_name", options=data.frame(names="canonical",args="")) ## End(Not run)
Convert a file from one format to another. Some formats are not supported on windows. This currently includes "INCHI" and "INCHIKEY".
convertFormatFile(from, to, fromFile, toFile,options=data.frame(names="gen2D",args=""))convertFormatFile(from, to, fromFile, toFile,options=data.frame(names="gen2D",args=""))
from |
The format that |
to |
The format to convert |
fromFile |
The name of the file to be converted |
toFile |
The name of the new file to be created or overwritten |
options |
This is a data frame listing OpenBabel Options that should be applied while converting each compound. The "names" column should be the name of each option and the "args" column sould be the arguments to the corresponding option. For example, the default value shown above will cause 2D coordinates to be generated from the input compound and saved in the output compound, assuming the ouptut format supports it. A full list of available options can be found by executing "obabel -L ops" at the command line. The current list is: 0xout <file.xxx> Additional file output AddInIndex Append input index to title AddPolarH Adds hydrogen to polar atoms only canonical Canonicalize the atom order energy ForceField Energy Evaluation (not displayed in GUI) fillUC <param> Fill the unit cell (strict or keepconnect) gen2D Generate 2D coordinates gen3D Generate 3D coordinates minimize ForceField Energy Minimization (not displayed in GUI) partialcharge <method> Calculate partial charges by specified method readconformer Adjacent conformers combined into a single molecule s Isomorphism filter(-s, -v options replacement)(not displayed in GUI) sort <desc> Sort by descriptor(~desc for reverse) unique [param] remove duplicates by descriptor;default inchi v Isomorphism filter(-s, -v options replacement)(not displayed in GUI) |
No value is returned. toFile will be created with the compound in the new format.
Kevin Horan
OpenBabel http://openbabel.org
## Not run: convertFormatFile("SMI","SDF","test.smiles","test.sdf") ## End(Not run)## Not run: convertFormatFile("SMI","SDF","test.smiles","test.sdf") ## End(Not run)
Converts compound data to an image format
convertToImage(from,to,source,toFile,options=data.frame(names="gen2D",args=""), out_options=data.frame(names=c("p"),args=c(300)))convertToImage(from,to,source,toFile,options=data.frame(names="gen2D",args=""), out_options=data.frame(names=c("p"),args=c(300)))
from |
The format that |
to |
The image format to convert |
source |
The initial compound format, as a string. The format of the string should be identical to the file format of the same name. Tabs and newlines may be represented with \t and \n, repsectivly. |
toFile |
The name of the output file. |
options |
This is a data frame listing OpenBabel Options that should be applied while converting each compound. The "names" column should be the name of each option and the "args" column sould be the arguments to the corresponding option. For example, the default value shown above will cause 2D coordinates to be generated from the input compound and saved in the output compound, assuming the ouptut format supports it. A full list of available options can be found by executing "obabel -L ops" at the command line. The current list is: 0xout <file.xxx> Additional file output AddInIndex Append input index to title AddPolarH Adds hydrogen to polar atoms only canonical Canonicalize the atom order energy ForceField Energy Evaluation (not displayed in GUI) fillUC <param> Fill the unit cell (strict or keepconnect) gen2D Generate 2D coordinates gen3D Generate 3D coordinates minimize ForceField Energy Minimization (not displayed in GUI) partialcharge <method> Calculate partial charges by specified method readconformer Adjacent conformers combined into a single molecule s Isomorphism filter(-s, -v options replacement)(not displayed in GUI) sort <desc> Sort by descriptor(~desc for reverse) unique [param] remove duplicates by descriptor;default inchi v Isomorphism filter(-s, -v options replacement)(not displayed in GUI) |
out_options |
Similar to the For the "PNG" format, possible options are: "p" for the number of pixels on the square output file. For the "SVG" format, see the OpenBabel docs for a full list (https://openbabel.org/docs/dev/FileFormats/SVG_2D_depiction.html). |
Generates a square image with the given compounds arranged in a grid format. The width (and height, which must be the same)
can be given using the out_options with argument "p". Writes the resulting image to a file specified by to_file.
Kevin Horan
OpenBabel http://openbabel.org
# Convert the given smiles string to a PNG image with size 500x500. ## Not run: sdfStr = convertToImage("SMI","PNG","CC(=O)OC1=CC=CC=C1C(=O)O\ttest_name","out_image.png", out_options=data.frame(names="p",args="500")) ## End(Not run)# Convert the given smiles string to a PNG image with size 500x500. ## Not run: sdfStr = convertToImage("SMI","PNG","CC(=O)OC1=CC=CC=C1C(=O)O\ttest_name","out_image.png", out_options=data.frame(names="p",args="500")) ## End(Not run)
Computes the exact mass of each compound given.
exactMass_OB(obmolRefs)exactMass_OB(obmolRefs)
obmolRefs |
A list of OBMol references ( of class '_p_OpenBabel__OBMol') representing the compounds. |
A vector of mass values.
Kevin Horan
## Not run: molRefs = forEachMol("SMILES","C1CCCCC1\ttest-compound-name",identity) exactMass_OB(molRefs) ## End(Not run)## Not run: molRefs = forEachMol("SMILES","C1CCCCC1\ttest-compound-name",identity) exactMass_OB(molRefs) ## End(Not run)
Generates fingerprints using OpenBabel. The compound format can be specified as anything supported by OpenBabel. The fingerprint name can also be specified.
fingerprint_OB(obmolRefs, fingerprintName,reverse=FALSE)fingerprint_OB(obmolRefs, fingerprintName,reverse=FALSE)
obmolRefs |
A list of OBMol references ( of class '_p_OpenBabel__OBMol') representing the molecules to
compute properites for. If you have your molecules in string format, you can create a list
of OBMol references using the |
fingerprintName |
The name of the fingerprint to generate. A list of available names can be found with "obabel -L fingerprints". Currently that list is: "FP2", "FP3", "FP4", and "MACCS". Currently "MACCS" is not supported on windows. |
reverse |
Reverse the order of bits in each fingerprint. The fingerprints
returned by this function are in reverse order to those returned by
obabel when using the same fingerprint calculation. As long as all
fingerprings being compared came from the same source (ie, all from
obabel, or all from ChemmineOB) the results should be correct. But if
you need to mix fingerprints from different sources, care muse be
taken to make sure they are all in the same direction. To make
ChemmineOB match obabel output, set We cannot set this to TRUE by default at this point as that could break code written prior to this change. |
A matrix of binary values is returned. There is a row for each compound. The length of a row is determined by the fingerprint specified.
Kevin Horan
## Not run: molRefs = forEachMol("SMILES","C1CCCCC1\ttest-compound-name",identity) fingerprint_OB(molRefs,"FP3") ## End(Not run)## Not run: molRefs = forEachMol("SMILES","C1CCCCC1\ttest-compound-name",identity) fingerprint_OB(molRefs,"FP3") ## End(Not run)
Reads in molecules from the given string in the given
format and calls function f on each molecule.
The results are then combinded using the reduce
function, if given.
forEachMol(inFormat, inString, f, reduce)forEachMol(inFormat, inString, f, reduce)
inFormat |
Format of string in |
inString |
The compounds to generate fingerprints for. The format should be exactly what would be in a file of the same format. Newlines can be represented with "\n". |
f |
A function taking one OBMol reference and possibly returning a result. |
reduce |
This function will be passed to the |
The result will be a List of return values from the
f function if not reduce function was
given. Otherwise it will be the result of the
reduce function applied to the results of the
f function.
Kevin Horan
## Not run: molRefs = forEachMol("SMILES","C1CCCCC1\ttest-compound-name", identity,c) ## End(Not run)## Not run: molRefs = forEachMol("SMILES","C1CCCCC1\ttest-compound-name", identity,c) ## End(Not run)
These are methods generated by SWIG for R. These should not generally be used outside of ChemmineOB, they are listed here to quite some warnings.
Generates the following descriptors: "cansmi", "cansmiNS", "formula", "HBA1", "HBA2", "HBD", "InChI", "InChIKey", "logP", "MR", "MW", "nF","title", "TPSA".
prop_OB(obmolRefs)prop_OB(obmolRefs)
obmolRefs |
A list of OBMol references ( of class '_p_OpenBabel__OBMol') representing the molecules to
compute properites for. If you have your molecules in string format, you can create a list
of OBMol references using the |
Returns a data frame with the following OpenBabel descriptors: "cansmi", "cansmiNS", "formula", "HBA1", "HBA2", "HBD", "InChI", "logP", "MR", "MW", "nF", "title", "TPSA". The "InChI" descriptor is not supported on windows currently, so that column will always be blank.
Kevin Horan
## Not run: # remove when ubuntu 16.04 bug fixed molRefs = forEachMol("SMILES","C1CCCCC1\ttest-compound-name",identity) prop_OB(molRefs) ## End(Not run)## Not run: # remove when ubuntu 16.04 bug fixed molRefs = forEachMol("SMILES","C1CCCCC1\ttest-compound-name",identity) prop_OB(molRefs) ## End(Not run)
Returns the number of matches found for each compound given.
smartsSearch_OB(obmolRefs, smartsPattern, uniqueMatches = TRUE)smartsSearch_OB(obmolRefs, smartsPattern, uniqueMatches = TRUE)
obmolRefs |
A list of OBMol references ( of class '_p_OpenBabel__OBMol') representing the molecules to
compute properites for. If you have your molecules in string format, you can create a list
of OBMol references using the |
smartsPattern |
Any valid SMARTS pattern. |
uniqueMatches |
Should only unique matches be counted? |
A vector of counts.
Kevin Horan
## Not run: molRefs = forEachMol("SMILES","C1CCCCC1\ttest-compound-name",identity) smartsSearch_OB(molRefs,"[CH3X4]") ## End(Not run)## Not run: molRefs = forEachMol("SMILES","C1CCCCC1\ttest-compound-name",identity) smartsSearch_OB(molRefs,"[CH3X4]") ## End(Not run)