cdkr icon indicating copy to clipboard operation
cdkr copied to clipboard

Only "0" returned for AromaticAtomsCountDescriptor and AromaticBondsCountDescriptor

Open bsedio opened this issue 1 year ago • 4 comments

Hello rcdk,

I am querying the CDK chemical properties using smiles generated using the CSI:FingerID module of Sirius. For some reason, I get "0" for number of aromatic atoms and bonds for every metabolite in my dataset, which includes many aromatic compounds (flavonoids, alkaloids from plant extracts).

Here is a small example:

dn.arom = c("org.openscience.cdk.qsar.descriptors.molecular.AromaticBondsCountDescriptor", "org.openscience.cdk.qsar.descriptors.molecular.AromaticAtomsCountDescriptor")
mol.aromtest = parse.smiles("C=CC1C2CC3C4=C(CCN3C(=O)C2=COC1OC5C(C(C(C(O5)CO)O)O)O)C6=CC=CC=C6N4")
mol.aromtest = parse.smiles(c("C=CC1C2CC3C4=C(CCN3C(=O)C2=COC1OC5C(C(C(C(O5)CO)O)O)O)C6=CC=CC=C6N4","CC(=NO)C1=CC(=C(C=C1CC2=NC=CC3=CC(=C(C=C32)OC)OC)OC)OC"))
desc.arom.test = eval.desc(mol.aromtest, dn.arom)
head(desc.arom.test)
                                                                    nAromBond naAromAtom
C=CC1C2CC3C4=C(CCN3C(=O)C2=COC1OC5C(C(C(C(O5)CO)O)O)O)C6=CC=CC=C6N4         0          0
CC(=NO)C1=CC(=C(C=C1CC2=NC=CC3=CC(=C(C=C32)OC)OC)OC)OC                      0          0

I am using R version 4.3.1 "Beagle Scouts" on Mac OS 12.6 Monterey and java version "1.8.0_391" Java(TM) SE Runtime Environment (build 1.8.0_391-b13) Java HotSpot(TM) 64-Bit Server VM (build 25.391-b13, mixed mode)

sessionInfo() R version 4.3.1 (2023-06-16) Platform: x86_64-apple-darwin20 (64-bit) Running under: macOS Monterey 12.6

Matrix products: default BLAS: /Library/Frameworks/R.framework/Versions/4.3-x86_64/Resources/lib/libRblas.0.dylib LAPACK: /Library/Frameworks/R.framework/Versions/4.3-x86_64/Resources/lib/libRlapack.dylib; LAPACK version 3.11.0

locale: [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

time zone: America/Chicago tzcode source: internal

attached base packages: [1] stats graphics grDevices utils datasets methods base

other attached packages: [1] rcdk_3.8.1 rcdklibs_2.8 rJava_1.0-6

loaded via a namespace (and not attached): [1] compiler_4.3.1 parallel_4.3.1 fingerprint_3.5.7 iterators_1.0.14 itertools_0.1-3 png_0.1-8

Thank you very much for any assistance, Brian

bsedio avatar Nov 08 '23 23:11 bsedio