annotate mayachemtools/bin/PathLengthFingerprints.pl @ 9:ab29fa5c8c1f draft default tip

Uploaded
author deepakjadmin
date Thu, 15 Dec 2016 14:18:03 -0500
parents 73ae111cf86f
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1 #!/usr/bin/perl -w
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
2 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
3 # $RCSfile: PathLengthFingerprints.pl,v $
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
4 # $Date: 2015/02/28 20:46:20 $
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
5 # $Revision: 1.50 $
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
6 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
7 # Author: Manish Sud <msud@san.rr.com>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
8 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
9 # Copyright (C) 2015 Manish Sud. All rights reserved.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
10 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
11 # This file is part of MayaChemTools.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
12 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
13 # MayaChemTools is free software; you can redistribute it and/or modify it under
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
14 # the terms of the GNU Lesser General Public License as published by the Free
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
15 # Software Foundation; either version 3 of the License, or (at your option) any
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
16 # later version.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
17 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
18 # MayaChemTools is distributed in the hope that it will be useful, but without
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
19 # any warranty; without even the implied warranty of merchantability of fitness
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
20 # for a particular purpose. See the GNU Lesser General Public License for more
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
21 # details.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
22 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
23 # You should have received a copy of the GNU Lesser General Public License
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
24 # along with MayaChemTools; if not, see <http://www.gnu.org/licenses/> or
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
25 # write to the Free Software Foundation Inc., 59 Temple Place, Suite 330,
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
26 # Boston, MA, 02111-1307, USA.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
27 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
28
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
29 use strict;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
30 use FindBin; use lib "$FindBin::Bin/../lib";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
31 use Getopt::Long;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
32 use File::Basename;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
33 use Text::ParseWords;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
34 use Benchmark;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
35 use FileUtil;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
36 use TextUtil;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
37 use SDFileUtil;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
38 use MoleculeFileIO;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
39 use FileIO::FingerprintsSDFileIO;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
40 use FileIO::FingerprintsTextFileIO;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
41 use FileIO::FingerprintsFPFileIO;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
42 use AtomTypes::AtomicInvariantsAtomTypes;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
43 use AtomTypes::FunctionalClassAtomTypes;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
44 use Fingerprints::PathLengthFingerprints;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
45
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
46 my($ScriptName, %Options, $StartTime, $EndTime, $TotalTime);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
47
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
48 # Autoflush STDOUT
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
49 $| = 1;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
50
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
51 # Starting message...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
52 $ScriptName = basename($0);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
53 print "\n$ScriptName: Starting...\n\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
54 $StartTime = new Benchmark;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
55
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
56 # Get the options and setup script...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
57 SetupScriptUsage();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
58 if ($Options{help} || @ARGV < 1) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
59 die GetUsageFromPod("$FindBin::Bin/$ScriptName");
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
60 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
61
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
62 my(@SDFilesList);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
63 @SDFilesList = ExpandFileNames(\@ARGV, "sdf sd");
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
64
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
65 # Process options...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
66 print "Processing options...\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
67 my(%OptionsInfo);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
68 ProcessOptions();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
69
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
70 # Setup information about input files...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
71 print "Checking input SD file(s)...\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
72 my(%SDFilesInfo);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
73 RetrieveSDFilesInfo();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
74
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
75 # Process input files..
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
76 my($FileIndex);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
77 if (@SDFilesList > 1) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
78 print "\nProcessing SD files...\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
79 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
80 for $FileIndex (0 .. $#SDFilesList) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
81 if ($SDFilesInfo{FileOkay}[$FileIndex]) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
82 print "\nProcessing file $SDFilesList[$FileIndex]...\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
83 GeneratePathLengthFingerprints($FileIndex);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
84 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
85 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
86 print "\n$ScriptName:Done...\n\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
87
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
88 $EndTime = new Benchmark;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
89 $TotalTime = timediff ($EndTime, $StartTime);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
90 print "Total time: ", timestr($TotalTime), "\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
91
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
92 ###############################################################################
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
93
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
94 # Generate fingerprints for a SD file...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
95 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
96 sub GeneratePathLengthFingerprints {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
97 my($FileIndex) = @_;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
98 my($CmpdCount, $IgnoredCmpdCount, $SDFile, $MoleculeFileIO, $Molecule, $PathLengthFingerprints, $NewFPSDFileIO, $NewFPTextFileIO, $NewFPFileIO);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
99
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
100 $SDFile = $SDFilesList[$FileIndex];
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
101
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
102 # Setup output files...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
103 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
104 ($NewFPSDFileIO, $NewFPTextFileIO, $NewFPFileIO) = SetupAndOpenOutputFiles($FileIndex);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
105
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
106 $MoleculeFileIO = new MoleculeFileIO('Name' => $SDFile);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
107 $MoleculeFileIO->Open();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
108
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
109 $CmpdCount = 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
110 $IgnoredCmpdCount = 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
111
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
112 COMPOUND: while ($Molecule = $MoleculeFileIO->ReadMolecule()) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
113 $CmpdCount++;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
114
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
115 # Filter compound data before calculating fingerprints...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
116 if ($OptionsInfo{Filter}) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
117 if (CheckAndFilterCompound($CmpdCount, $Molecule)) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
118 $IgnoredCmpdCount++;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
119 next COMPOUND;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
120 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
121 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
122
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
123 $PathLengthFingerprints = GenerateMoleculeFingerprints($Molecule);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
124 if (!$PathLengthFingerprints) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
125 $IgnoredCmpdCount++;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
126 ProcessIgnoredCompound('FingerprintsGenerationFailed', $CmpdCount, $Molecule);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
127 next COMPOUND;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
128 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
129
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
130 WriteDataToOutputFiles($FileIndex, $CmpdCount, $Molecule, $PathLengthFingerprints, $NewFPSDFileIO, $NewFPTextFileIO, $NewFPFileIO);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
131 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
132 $MoleculeFileIO->Close();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
133
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
134 if ($NewFPSDFileIO) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
135 $NewFPSDFileIO->Close();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
136 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
137 if ($NewFPTextFileIO) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
138 $NewFPTextFileIO->Close();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
139 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
140 if ($NewFPFileIO) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
141 $NewFPFileIO->Close();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
142 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
143
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
144 WriteFingerprintsGenerationSummaryStatistics($CmpdCount, $IgnoredCmpdCount);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
145 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
146
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
147 # Process compound being ignored due to problems in fingerprints geneation...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
148 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
149 sub ProcessIgnoredCompound {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
150 my($Mode, $CmpdCount, $Molecule) = @_;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
151 my($CmpdID, $DataFieldLabelAndValuesRef);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
152
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
153 $DataFieldLabelAndValuesRef = $Molecule->GetDataFieldLabelAndValues();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
154 $CmpdID = SetupCmpdIDForOutputFiles($CmpdCount, $Molecule, $DataFieldLabelAndValuesRef);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
155
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
156 MODE: {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
157 if ($Mode =~ /^ContainsNonElementalData$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
158 warn "\nWarning: Ignoring compound record number $CmpdCount with ID $CmpdID: Compound contains atom data corresponding to non-elemental atom symbol(s)...\n\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
159 next MODE;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
160 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
161
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
162 if ($Mode =~ /^ContainsNoElementalData$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
163 warn "\nWarning: Ignoring compound record number $CmpdCount with ID $CmpdID: Compound contains no atom data...\n\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
164 next MODE;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
165 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
166
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
167 if ($Mode =~ /^FingerprintsGenerationFailed$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
168 warn "\nWarning: Ignoring compound record number $CmpdCount with ID $CmpdID: Fingerprints generation didn't succeed...\n\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
169 next MODE;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
170 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
171 warn "\nWarning: Ignoring compound record number $CmpdCount with ID $CmpdID: Fingerprints generation didn't succeed...\n\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
172 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
173 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
174
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
175 # Check and filter compounds....
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
176 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
177 sub CheckAndFilterCompound {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
178 my($CmpdCount, $Molecule) = @_;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
179 my($ElementCount, $NonElementCount);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
180
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
181 ($ElementCount, $NonElementCount) = $Molecule->GetNumOfElementsAndNonElements();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
182
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
183 if ($NonElementCount) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
184 ProcessIgnoredCompound('ContainsNonElementalData', $CmpdCount, $Molecule);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
185 return 1;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
186 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
187
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
188 if (!$ElementCount) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
189 ProcessIgnoredCompound('ContainsNoElementalData', $CmpdCount, $Molecule);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
190 return 1;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
191 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
192
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
193 return 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
194 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
195
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
196 # Write out compounds fingerprints generation summary statistics...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
197 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
198 sub WriteFingerprintsGenerationSummaryStatistics {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
199 my($CmpdCount, $IgnoredCmpdCount) = @_;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
200 my($ProcessedCmpdCount);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
201
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
202 $ProcessedCmpdCount = $CmpdCount - $IgnoredCmpdCount;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
203
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
204 print "\nNumber of compounds: $CmpdCount\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
205 print "Number of compounds processed successfully during fingerprints generation: $ProcessedCmpdCount\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
206 print "Number of compounds ignored during fingerprints generation: $IgnoredCmpdCount\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
207 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
208
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
209 # Open output files...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
210 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
211 sub SetupAndOpenOutputFiles {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
212 my($FileIndex) = @_;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
213 my($NewFPSDFile, $NewFPFile, $NewFPTextFile, $NewFPSDFileIO, $NewFPTextFileIO, $NewFPFileIO, %FingerprintsFileIOParams);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
214
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
215 ($NewFPSDFileIO, $NewFPTextFileIO, $NewFPFileIO) = (undef) x 3;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
216
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
217 # Setup common parameters for fingerprints file IO objects...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
218 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
219 %FingerprintsFileIOParams = ();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
220 if ($OptionsInfo{Mode} =~ /^PathLengthBits$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
221 %FingerprintsFileIOParams = ('Mode' => 'Write', 'Overwrite' => $OptionsInfo{OverwriteFiles}, 'FingerprintsStringMode' => 'FingerprintsBitVectorString', 'BitStringFormat' => $OptionsInfo{BitStringFormat}, 'BitsOrder' => $OptionsInfo{BitsOrder});
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
222 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
223 elsif ($OptionsInfo{Mode} =~ /^PathLengthCount$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
224 %FingerprintsFileIOParams = ('Mode' => 'Write', 'Overwrite' => $OptionsInfo{OverwriteFiles}, 'FingerprintsStringMode' => 'FingerprintsVectorString', 'VectorStringFormat' => $OptionsInfo{VectorStringFormat});
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
225 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
226
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
227 if ($OptionsInfo{SDOutput}) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
228 $NewFPSDFile = $SDFilesInfo{SDOutFileNames}[$FileIndex];
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
229 print "Generating SD file $NewFPSDFile...\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
230 $NewFPSDFileIO = new FileIO::FingerprintsSDFileIO('Name' => $NewFPSDFile, %FingerprintsFileIOParams, 'FingerprintsFieldLabel' => $OptionsInfo{FingerprintsLabel});
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
231 $NewFPSDFileIO->Open();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
232 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
233
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
234 if ($OptionsInfo{FPOutput}) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
235 $NewFPFile = $SDFilesInfo{FPOutFileNames}[$FileIndex];
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
236 print "Generating FP file $NewFPFile...\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
237 $NewFPFileIO = new FileIO::FingerprintsFPFileIO('Name' => $NewFPFile, %FingerprintsFileIOParams);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
238 $NewFPFileIO->Open();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
239 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
240
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
241 if ($OptionsInfo{TextOutput}) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
242 my($ColLabelsRef);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
243
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
244 $NewFPTextFile = $SDFilesInfo{TextOutFileNames}[$FileIndex];
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
245 $ColLabelsRef = SetupFPTextFileCoulmnLabels($FileIndex);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
246
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
247 print "Generating text file $NewFPTextFile...\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
248 $NewFPTextFileIO = new FileIO::FingerprintsTextFileIO('Name' => $NewFPTextFile, %FingerprintsFileIOParams, 'DataColLabels' => $ColLabelsRef, 'OutDelim' => $OptionsInfo{OutDelim}, 'OutQuote' => $OptionsInfo{OutQuote});
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
249 $NewFPTextFileIO->Open();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
250 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
251
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
252 return ($NewFPSDFileIO, $NewFPTextFileIO, $NewFPFileIO);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
253 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
254
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
255 # Write fingerpritns and other data to appropriate output files...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
256 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
257 sub WriteDataToOutputFiles {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
258 my($FileIndex, $CmpdCount, $Molecule, $PathLengthFingerprints, $NewFPSDFileIO, $NewFPTextFileIO, $NewFPFileIO) = @_;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
259 my($DataFieldLabelAndValuesRef);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
260
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
261 $DataFieldLabelAndValuesRef = undef;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
262 if ($NewFPTextFileIO || $NewFPFileIO) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
263 $DataFieldLabelAndValuesRef = $Molecule->GetDataFieldLabelAndValues();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
264 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
265
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
266 if ($NewFPSDFileIO) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
267 my($CmpdString);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
268
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
269 $CmpdString = $Molecule->GetInputMoleculeString();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
270 $NewFPSDFileIO->WriteFingerprints($PathLengthFingerprints, $CmpdString);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
271 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
272
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
273 if ($NewFPTextFileIO) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
274 my($ColValuesRef);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
275
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
276 $ColValuesRef = SetupFPTextFileCoulmnValues($FileIndex, $CmpdCount, $Molecule, $DataFieldLabelAndValuesRef);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
277 $NewFPTextFileIO->WriteFingerprints($PathLengthFingerprints, $ColValuesRef);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
278 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
279
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
280 if ($NewFPFileIO) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
281 my($CompoundID);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
282
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
283 $CompoundID = SetupCmpdIDForOutputFiles($CmpdCount, $Molecule, $DataFieldLabelAndValuesRef);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
284 $NewFPFileIO->WriteFingerprints($PathLengthFingerprints, $CompoundID);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
285 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
286 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
287
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
288 # Generate approriate column labels for FPText output file...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
289 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
290 sub SetupFPTextFileCoulmnLabels {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
291 my($FileIndex) = @_;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
292 my($Line, @ColLabels);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
293
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
294 @ColLabels = ();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
295 if ($OptionsInfo{DataFieldsMode} =~ /^All$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
296 push @ColLabels, @{$SDFilesInfo{AllDataFieldsRef}[$FileIndex]};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
297 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
298 elsif ($OptionsInfo{DataFieldsMode} =~ /^Common$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
299 push @ColLabels, @{$SDFilesInfo{CommonDataFieldsRef}[$FileIndex]};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
300 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
301 elsif ($OptionsInfo{DataFieldsMode} =~ /^Specify$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
302 push @ColLabels, @{$OptionsInfo{SpecifiedDataFields}};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
303 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
304 elsif ($OptionsInfo{DataFieldsMode} =~ /^CompoundID$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
305 push @ColLabels, $OptionsInfo{CompoundIDLabel};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
306 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
307 # Add fingerprints label...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
308 push @ColLabels, $OptionsInfo{FingerprintsLabel};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
309
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
310 return \@ColLabels;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
311 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
312
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
313 # Generate column values FPText output file..
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
314 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
315 sub SetupFPTextFileCoulmnValues {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
316 my($FileIndex, $CmpdCount, $Molecule, $DataFieldLabelAndValuesRef) = @_;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
317 my(@ColValues);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
318
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
319 @ColValues = ();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
320 if ($OptionsInfo{DataFieldsMode} =~ /^CompoundID$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
321 push @ColValues, SetupCmpdIDForOutputFiles($CmpdCount, $Molecule, $DataFieldLabelAndValuesRef);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
322 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
323 elsif ($OptionsInfo{DataFieldsMode} =~ /^All$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
324 @ColValues = map { exists $DataFieldLabelAndValuesRef->{$_} ? $DataFieldLabelAndValuesRef->{$_} : ''} @{$SDFilesInfo{AllDataFieldsRef}[$FileIndex]};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
325 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
326 elsif ($OptionsInfo{DataFieldsMode} =~ /^Common$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
327 @ColValues = map { exists $DataFieldLabelAndValuesRef->{$_} ? $DataFieldLabelAndValuesRef->{$_} : ''} @{$SDFilesInfo{CommonDataFieldsRef}[$FileIndex]};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
328 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
329 elsif ($OptionsInfo{DataFieldsMode} =~ /^Specify$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
330 @ColValues = map { exists $DataFieldLabelAndValuesRef->{$_} ? $DataFieldLabelAndValuesRef->{$_} : ''} @{$OptionsInfo{SpecifiedDataFields}};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
331 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
332
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
333 return \@ColValues;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
334 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
335
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
336 # Generate compound ID for FP and FPText output files..
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
337 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
338 sub SetupCmpdIDForOutputFiles {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
339 my($CmpdCount, $Molecule, $DataFieldLabelAndValuesRef) = @_;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
340 my($CmpdID);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
341
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
342 $CmpdID = '';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
343 if ($OptionsInfo{CompoundIDMode} =~ /^MolNameOrLabelPrefix$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
344 my($MolName);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
345 $MolName = $Molecule->GetName();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
346 $CmpdID = $MolName ? $MolName : "$OptionsInfo{CompoundID}${CmpdCount}";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
347 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
348 elsif ($OptionsInfo{CompoundIDMode} =~ /^LabelPrefix$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
349 $CmpdID = "$OptionsInfo{CompoundID}${CmpdCount}";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
350 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
351 elsif ($OptionsInfo{CompoundIDMode} =~ /^DataField$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
352 my($SpecifiedDataField);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
353 $SpecifiedDataField = $OptionsInfo{CompoundID};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
354 $CmpdID = exists $DataFieldLabelAndValuesRef->{$SpecifiedDataField} ? $DataFieldLabelAndValuesRef->{$SpecifiedDataField} : '';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
355 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
356 elsif ($OptionsInfo{CompoundIDMode} =~ /^MolName$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
357 $CmpdID = $Molecule->GetName();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
358 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
359 return $CmpdID;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
360 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
361
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
362 # Generate fingerprints for molecule...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
363 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
364 sub GenerateMoleculeFingerprints {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
365 my($Molecule) = @_;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
366 my($PathLengthFingerprints);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
367
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
368 if ($OptionsInfo{KeepLargestComponent}) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
369 $Molecule->KeepLargestComponent();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
370 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
371 if ($OptionsInfo{IgnoreHydrogens}) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
372 $Molecule->DeleteHydrogens();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
373 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
374
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
375 if ($OptionsInfo{DetectAromaticity}) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
376 if (!$Molecule->DetectRings()) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
377 return undef;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
378 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
379 $Molecule->SetAromaticityModel($OptionsInfo{AromaticityModel});
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
380 $Molecule->DetectAromaticity();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
381 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
382
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
383 $PathLengthFingerprints = undef;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
384 if ($OptionsInfo{Mode} =~ /^PathLengthBits$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
385 $PathLengthFingerprints = GeneratePathLengthBitsFingerprints($Molecule);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
386 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
387 elsif ($OptionsInfo{Mode} =~ /^PathLengthCount$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
388 $PathLengthFingerprints = GeneratePathLengthCountFingerprints($Molecule);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
389 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
390 else {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
391 die "Error: The value specified, $Options{mode}, for option \"-m, --mode\" is not valid. Allowed values: PathLengthBits or PathLengthCount\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
392 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
393
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
394 return $PathLengthFingerprints;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
395 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
396
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
397 # Generate pathlength bits finerprints for molecule...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
398 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
399 sub GeneratePathLengthBitsFingerprints {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
400 my($Molecule) = @_;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
401 my($PathLengthFingerprints);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
402
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
403 $PathLengthFingerprints = new Fingerprints::PathLengthFingerprints('Molecule' => $Molecule, 'Type' => 'PathLengthBits', 'AtomIdentifierType' => $OptionsInfo{AtomIdentifierType}, 'NumOfBitsToSetPerPath' => $OptionsInfo{NumOfBitsToSetPerPath}, 'Size' => $OptionsInfo{Size}, 'MinLength' => $OptionsInfo{MinPathLength}, 'MaxLength' => $OptionsInfo{MaxPathLength}, 'AllowRings' => $OptionsInfo{AllowRings}, 'AllowSharedBonds' => $OptionsInfo{AllowSharedBonds}, 'UseBondSymbols' => $OptionsInfo{UseBondSymbols}, 'UseUniquePaths' => $OptionsInfo{UseUniquePaths}, 'UsePerlCoreRandom' => $OptionsInfo{UsePerlCoreRandom});
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
404
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
405 # Set atom identifier type...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
406 SetAtomIdentifierTypeValuesToUse($PathLengthFingerprints);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
407
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
408 # Generate fingerprints...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
409 $PathLengthFingerprints->GenerateFingerprints();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
410
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
411 # Make sure fingerprints generation is successful...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
412 if (!$PathLengthFingerprints->IsFingerprintsGenerationSuccessful()) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
413 return undef;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
414 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
415
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
416 if ($OptionsInfo{Fold}) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
417 my($CheckSizeValue) = 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
418 $PathLengthFingerprints->FoldFingerprintsBySize($OptionsInfo{FoldedSize}, $CheckSizeValue);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
419 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
420
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
421 return $PathLengthFingerprints;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
422 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
423
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
424 # Generate pathlength count finerprints for molecule...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
425 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
426 sub GeneratePathLengthCountFingerprints {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
427 my($Molecule) = @_;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
428 my($PathLengthFingerprints);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
429
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
430 $PathLengthFingerprints = new Fingerprints::PathLengthFingerprints('Molecule' => $Molecule, 'Type' => 'PathLengthCount', 'AtomIdentifierType' => $OptionsInfo{AtomIdentifierType}, 'MinLength' => $OptionsInfo{MinPathLength}, 'MaxLength' => $OptionsInfo{MaxPathLength}, 'AllowRings' => $OptionsInfo{AllowRings}, 'AllowSharedBonds' => $OptionsInfo{AllowSharedBonds}, 'UseBondSymbols' => $OptionsInfo{UseBondSymbols}, 'UseUniquePaths' => $OptionsInfo{UseUniquePaths});
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
431
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
432 # Set atom identifier type...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
433 SetAtomIdentifierTypeValuesToUse($PathLengthFingerprints);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
434
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
435 # Generate fingerprints...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
436 $PathLengthFingerprints->GenerateFingerprints();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
437
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
438 # Make sure fingerprints generation is successful...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
439 if (!$PathLengthFingerprints->IsFingerprintsGenerationSuccessful()) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
440 return undef;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
441 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
442 return $PathLengthFingerprints;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
443 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
444
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
445 # Set atom identifier type to use for generating path strings...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
446 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
447 sub SetAtomIdentifierTypeValuesToUse {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
448 my($PathLengthFingerprints) = @_;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
449
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
450 if ($OptionsInfo{AtomIdentifierType} =~ /^AtomicInvariantsAtomTypes$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
451 $PathLengthFingerprints->SetAtomicInvariantsToUse(\@{$OptionsInfo{AtomicInvariantsToUse}});
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
452 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
453 elsif ($OptionsInfo{AtomIdentifierType} =~ /^FunctionalClassAtomTypes$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
454 $PathLengthFingerprints->SetFunctionalClassesToUse(\@{$OptionsInfo{FunctionalClassesToUse}});
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
455 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
456 elsif ($OptionsInfo{AtomIdentifierType} =~ /^(DREIDINGAtomTypes|EStateAtomTypes|MMFF94AtomTypes|SLogPAtomTypes|SYBYLAtomTypes|TPSAAtomTypes|UFFAtomTypes)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
457 # Nothing to do for now...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
458 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
459 else {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
460 die "Error: The value specified, $Options{atomidentifiertype}, for option \"-a, --AtomIdentifierType\" is not valid. Supported atom identifier types in current release of MayaChemTools: AtomicInvariantsAtomTypes, DREIDINGAtomTypes, EStateAtomTypes, FunctionalClassAtomTypes, MMFF94AtomTypes, SLogPAtomTypes, SYBYLAtomTypes, TPSAAtomTypes, UFFAtomTypes\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
461 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
462 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
463
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
464 # Retrieve information about SD files...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
465 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
466 sub RetrieveSDFilesInfo {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
467 my($SDFile, $Index, $FileDir, $FileExt, $FileName, $OutFileRoot, $TextOutFileExt, $SDOutFileExt, $FPOutFileExt, $NewSDFileName, $NewFPFileName, $NewTextFileName, $CheckDataField, $CollectDataFields, $AllDataFieldsRef, $CommonDataFieldsRef);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
468
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
469 %SDFilesInfo = ();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
470 @{$SDFilesInfo{FileOkay}} = ();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
471 @{$SDFilesInfo{OutFileRoot}} = ();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
472 @{$SDFilesInfo{SDOutFileNames}} = ();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
473 @{$SDFilesInfo{FPOutFileNames}} = ();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
474 @{$SDFilesInfo{TextOutFileNames}} = ();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
475 @{$SDFilesInfo{AllDataFieldsRef}} = ();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
476 @{$SDFilesInfo{CommonDataFieldsRef}} = ();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
477
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
478 $CheckDataField = ($OptionsInfo{TextOutput} && ($OptionsInfo{DataFieldsMode} =~ /^CompoundID$/i) && ($OptionsInfo{CompoundIDMode} =~ /^DataField$/i)) ? 1 : 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
479 $CollectDataFields = ($OptionsInfo{TextOutput} && ($OptionsInfo{DataFieldsMode} =~ /^(All|Common)$/i)) ? 1 : 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
480
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
481 FILELIST: for $Index (0 .. $#SDFilesList) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
482 $SDFile = $SDFilesList[$Index];
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
483
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
484 $SDFilesInfo{FileOkay}[$Index] = 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
485 $SDFilesInfo{OutFileRoot}[$Index] = '';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
486 $SDFilesInfo{SDOutFileNames}[$Index] = '';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
487 $SDFilesInfo{FPOutFileNames}[$Index] = '';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
488 $SDFilesInfo{TextOutFileNames}[$Index] = '';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
489
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
490 $SDFile = $SDFilesList[$Index];
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
491 if (!(-e $SDFile)) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
492 warn "Warning: Ignoring file $SDFile: It doesn't exist\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
493 next FILELIST;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
494 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
495 if (!CheckFileType($SDFile, "sd sdf")) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
496 warn "Warning: Ignoring file $SDFile: It's not a SD file\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
497 next FILELIST;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
498 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
499
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
500 if ($CheckDataField) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
501 # Make sure data field exists in SD file..
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
502 my($CmpdString, $SpecifiedDataField, @CmpdLines, %DataFieldValues);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
503
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
504 @CmpdLines = ();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
505 open SDFILE, "$SDFile" or die "Error: Couldn't open $SDFile: $! \n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
506 $CmpdString = ReadCmpdString(\*SDFILE);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
507 close SDFILE;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
508 @CmpdLines = split "\n", $CmpdString;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
509 %DataFieldValues = GetCmpdDataHeaderLabelsAndValues(\@CmpdLines);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
510 $SpecifiedDataField = $OptionsInfo{CompoundID};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
511 if (!exists $DataFieldValues{$SpecifiedDataField}) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
512 warn "Warning: Ignoring file $SDFile: Data field value, $SpecifiedDataField, using \"--CompoundID\" option in \"DataField\" \"--CompoundIDMode\" doesn't exist\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
513 next FILELIST;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
514 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
515 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
516
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
517 $AllDataFieldsRef = '';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
518 $CommonDataFieldsRef = '';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
519 if ($CollectDataFields) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
520 my($CmpdCount);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
521 open SDFILE, "$SDFile" or die "Error: Couldn't open $SDFile: $! \n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
522 ($CmpdCount, $AllDataFieldsRef, $CommonDataFieldsRef) = GetAllAndCommonCmpdDataHeaderLabels(\*SDFILE);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
523 close SDFILE;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
524 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
525
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
526 # Setup output file names...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
527 $FileDir = ""; $FileName = ""; $FileExt = "";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
528 ($FileDir, $FileName, $FileExt) = ParseFileName($SDFile);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
529
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
530 $TextOutFileExt = "csv";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
531 if ($Options{outdelim} =~ /^tab$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
532 $TextOutFileExt = "tsv";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
533 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
534 $SDOutFileExt = $FileExt;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
535 $FPOutFileExt = "fpf";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
536
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
537 if ($OptionsInfo{OutFileRoot} && (@SDFilesList == 1)) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
538 my ($RootFileDir, $RootFileName, $RootFileExt) = ParseFileName($OptionsInfo{OutFileRoot});
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
539 if ($RootFileName && $RootFileExt) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
540 $FileName = $RootFileName;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
541 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
542 else {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
543 $FileName = $OptionsInfo{OutFileRoot};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
544 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
545 $OutFileRoot = $FileName;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
546 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
547 else {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
548 $OutFileRoot = "${FileName}PathLengthFP";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
549 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
550
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
551 $NewSDFileName = "${OutFileRoot}.${SDOutFileExt}";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
552 $NewFPFileName = "${OutFileRoot}.${FPOutFileExt}";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
553 $NewTextFileName = "${OutFileRoot}.${TextOutFileExt}";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
554
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
555 if ($OptionsInfo{SDOutput}) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
556 if ($SDFile =~ /$NewSDFileName/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
557 warn "Warning: Ignoring input file $SDFile: Same output, $NewSDFileName, and input file names.\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
558 print "Specify a different name using \"-r --root\" option or use default name.\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
559 next FILELIST;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
560 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
561 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
562
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
563 if (!$OptionsInfo{OverwriteFiles}) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
564 # Check SD, FP and text outout files...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
565 if ($OptionsInfo{SDOutput}) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
566 if (-e $NewSDFileName) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
567 warn "Warning: Ignoring file $SDFile: The file $NewSDFileName already exists\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
568 next FILELIST;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
569 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
570 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
571 if ($OptionsInfo{FPOutput}) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
572 if (-e $NewFPFileName) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
573 warn "Warning: Ignoring file $SDFile: The file $NewFPFileName already exists\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
574 next FILELIST;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
575 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
576 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
577 if ($OptionsInfo{TextOutput}) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
578 if (-e $NewTextFileName) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
579 warn "Warning: Ignoring file $SDFile: The file $NewTextFileName already exists\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
580 next FILELIST;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
581 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
582 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
583 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
584
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
585 $SDFilesInfo{FileOkay}[$Index] = 1;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
586
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
587 $SDFilesInfo{OutFileRoot}[$Index] = $OutFileRoot;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
588 $SDFilesInfo{SDOutFileNames}[$Index] = $NewSDFileName;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
589 $SDFilesInfo{FPOutFileNames}[$Index] = $NewFPFileName;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
590 $SDFilesInfo{TextOutFileNames}[$Index] = $NewTextFileName;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
591
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
592 $SDFilesInfo{AllDataFieldsRef}[$Index] = $AllDataFieldsRef;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
593 $SDFilesInfo{CommonDataFieldsRef}[$Index] = $CommonDataFieldsRef;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
594 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
595 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
596
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
597 # Process option values...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
598 sub ProcessOptions {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
599 %OptionsInfo = ();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
600
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
601 $OptionsInfo{Mode} = $Options{mode};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
602 $OptionsInfo{AromaticityModel} = $Options{aromaticitymodel};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
603 $OptionsInfo{PathMode} = $Options{pathmode};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
604
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
605 ProcessAtomIdentifierTypeOptions();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
606
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
607 $OptionsInfo{BitsOrder} = $Options{bitsorder};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
608 $OptionsInfo{BitStringFormat} = $Options{bitstringformat};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
609
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
610 $OptionsInfo{CompoundIDMode} = $Options{compoundidmode};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
611 $OptionsInfo{CompoundIDLabel} = $Options{compoundidlabel};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
612 $OptionsInfo{DataFieldsMode} = $Options{datafieldsmode};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
613
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
614 my(@SpecifiedDataFields);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
615 @SpecifiedDataFields = ();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
616
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
617 @{$OptionsInfo{SpecifiedDataFields}} = ();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
618 $OptionsInfo{CompoundID} = '';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
619
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
620 if ($Options{datafieldsmode} =~ /^CompoundID$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
621 if ($Options{compoundidmode} =~ /^DataField$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
622 if (!$Options{compoundid}) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
623 die "Error: You must specify a value for \"--CompoundID\" option in \"DataField\" \"--CompoundIDMode\". \n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
624 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
625 $OptionsInfo{CompoundID} = $Options{compoundid};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
626 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
627 elsif ($Options{compoundidmode} =~ /^(LabelPrefix|MolNameOrLabelPrefix)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
628 $OptionsInfo{CompoundID} = $Options{compoundid} ? $Options{compoundid} : 'Cmpd';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
629 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
630 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
631 elsif ($Options{datafieldsmode} =~ /^Specify$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
632 if (!$Options{datafields}) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
633 die "Error: You must specify a value for \"--DataFields\" option in \"Specify\" \"-d, --DataFieldsMode\". \n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
634 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
635 @SpecifiedDataFields = split /\,/, $Options{datafields};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
636 push @{$OptionsInfo{SpecifiedDataFields}}, @SpecifiedDataFields;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
637 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
638
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
639 if ($Options{atomidentifiertype} !~ /^AtomicInvariantsAtomTypes$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
640 if ($Options{detectaromaticity} =~ /^No$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
641 die "Error: The value specified, $Options{detectaromaticity}, for option \"--DetectAromaticity\" is not valid. No value is only allowed during AtomicInvariantsAtomTypes value for \"-a, --AtomIdentifierType\" \n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
642 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
643 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
644 $OptionsInfo{DetectAromaticity} = ($Options{detectaromaticity} =~ /^Yes$/i) ? 1 : 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
645
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
646 $OptionsInfo{Filter} = ($Options{filter} =~ /^Yes$/i) ? 1 : 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
647
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
648 $OptionsInfo{FingerprintsLabel} = $Options{fingerprintslabel} ? $Options{fingerprintslabel} : 'PathLengthFingerprints';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
649
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
650 my($Size, $MinSize, $MaxSize);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
651 $MinSize = 32;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
652 $MaxSize = 2**32;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
653 $Size = $Options{size};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
654 if (!(IsPositiveInteger($Size) && $Size >= $MinSize && $Size <= $MaxSize && IsNumberPowerOfNumber($Size, 2))) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
655 die "Error: Invalid size value, $Size, for \"-s, --size\" option. Allowed values: power of 2, >= minimum size of $MinSize, and <= maximum size of $MaxSize.\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
656 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
657 $OptionsInfo{Size} = $Size;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
658
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
659 $OptionsInfo{Fold} = ($Options{fold} =~ /^Yes$/i) ? 1 : 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
660 my($FoldedSize);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
661 $FoldedSize = $Options{foldedsize};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
662 if ($Options{fold} =~ /^Yes$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
663 if (!(IsPositiveInteger($FoldedSize) && $FoldedSize < $Size && IsNumberPowerOfNumber($FoldedSize, 2))) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
664 die "Error: Invalid folded size value, $FoldedSize, for \"--FoldedSize\" option. Allowed values: power of 2, >= minimum size of $MinSize, and < size value of $Size.\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
665 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
666 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
667 $OptionsInfo{FoldedSize} = $FoldedSize;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
668
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
669 $OptionsInfo{IgnoreHydrogens} = ($Options{ignorehydrogens} =~ /^Yes$/i) ? 1 : 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
670 $OptionsInfo{KeepLargestComponent} = ($Options{keeplargestcomponent} =~ /^Yes$/i) ? 1 : 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
671
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
672 my($MinPathLength, $MaxPathLength);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
673 $MinPathLength = $Options{minpathlength};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
674 $MaxPathLength = $Options{maxpathlength};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
675 if (!IsPositiveInteger($MinPathLength)) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
676 die "Error: Invalid path length value, $MinPathLength, for \"--MinPathLength\" option. Allowed values: > 0\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
677 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
678 if (!IsPositiveInteger($MaxPathLength)) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
679 die "Error: Invalid path length value, $MaxPathLength, for \"--MinPathLength\" option. Allowed values: > 0\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
680 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
681 if ($MinPathLength >= $MaxPathLength) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
682 die "Error: Invalid minimum and maximum path length values, $MinPathLength and $MaxPathLength, for \"--MinPathLength\" and \"--MaxPathLength\"options. Allowed values: minimum path length value must be smaller than maximum path length value.\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
683 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
684 $OptionsInfo{MinPathLength} = $MinPathLength;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
685 $OptionsInfo{MaxPathLength} = $MaxPathLength;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
686
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
687 my($NumOfBitsToSetPerPath);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
688 $NumOfBitsToSetPerPath = $Options{numofbitstosetperpath};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
689 if (!IsPositiveInteger($MaxPathLength)) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
690 die "Error: Invalid value, $NumOfBitsToSetPerPath, for \"-n, --NumOfBitsToSetPerPath\" option. Allowed values: > 0\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
691 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
692 if ($NumOfBitsToSetPerPath >= $Size) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
693 die "Error: Invalid value, $NumOfBitsToSetPerPath, for \"-n, --NumOfBitsToSetPerPath\" option. Allowed values: It must be less than the size, $Size, of the fingerprint bit-string.\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
694 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
695 $OptionsInfo{NumOfBitsToSetPerPath} = $NumOfBitsToSetPerPath;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
696
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
697 $OptionsInfo{Output} = $Options{output};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
698 $OptionsInfo{SDOutput} = ($Options{output} =~ /^(SD|All)$/i) ? 1 : 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
699 $OptionsInfo{FPOutput} = ($Options{output} =~ /^(FP|All)$/i) ? 1 : 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
700 $OptionsInfo{TextOutput} = ($Options{output} =~ /^(Text|All)$/i) ? 1 : 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
701
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
702 $OptionsInfo{OutDelim} = $Options{outdelim};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
703 $OptionsInfo{OutQuote} = ($Options{quote} =~ /^Yes$/i) ? 1 : 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
704
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
705 $OptionsInfo{OverwriteFiles} = $Options{overwrite} ? 1 : 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
706 $OptionsInfo{OutFileRoot} = $Options{root} ? $Options{root} : 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
707
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
708 $OptionsInfo{UseBondSymbols} = ($Options{usebondsymbols} =~ /^Yes$/i) ? 1 : 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
709
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
710 $OptionsInfo{UsePerlCoreRandom} = ($Options{useperlcorerandom} =~ /^Yes$/i) ? 1 : 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
711
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
712 $OptionsInfo{UseUniquePaths} = ($Options{useuniquepaths} =~ /^Yes$/i) ? 1 : 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
713
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
714 $OptionsInfo{VectorStringFormat} = $Options{vectorstringformat};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
715
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
716 # Setup parameters used during generation of fingerprints by PathLengthFingerprints class...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
717 my($AllowRings, $AllowSharedBonds);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
718 $AllowRings = 1;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
719 $AllowSharedBonds = 1;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
720 MODE: {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
721 if ($Options{pathmode} =~ /^AtomPathsWithoutRings$/i) { $AllowSharedBonds = 0; $AllowRings = 0; last MODE;}
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
722 if ($Options{pathmode} =~ /^AtomPathsWithRings$/i) { $AllowSharedBonds = 0; $AllowRings = 1; last MODE;}
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
723 if ($Options{pathmode} =~ /^AllAtomPathsWithoutRings$/i) { $AllowSharedBonds = 1; $AllowRings = 0; last MODE;}
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
724 if ($Options{pathmode} =~ /^AllAtomPathsWithRings$/i) { $AllowSharedBonds = 1; $AllowRings = 1; last MODE;}
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
725 die "Error: ProcessOptions: mode value, $Options{pathmode}, is not supported.\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
726 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
727 $OptionsInfo{AllowRings} = $AllowRings;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
728 $OptionsInfo{AllowSharedBonds} = $AllowSharedBonds;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
729 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
730
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
731 # Process atom identifier type and related options...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
732 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
733 sub ProcessAtomIdentifierTypeOptions {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
734
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
735 $OptionsInfo{AtomIdentifierType} = $Options{atomidentifiertype};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
736
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
737 if ($Options{atomidentifiertype} =~ /^AtomicInvariantsAtomTypes$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
738 ProcessAtomicInvariantsToUseOption();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
739 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
740 elsif ($Options{atomidentifiertype} =~ /^FunctionalClassAtomTypes$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
741 ProcessFunctionalClassesToUse();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
742 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
743 elsif ($OptionsInfo{AtomIdentifierType} =~ /^(DREIDINGAtomTypes|EStateAtomTypes|MMFF94AtomTypes|SLogPAtomTypes|SYBYLAtomTypes|TPSAAtomTypes|UFFAtomTypes)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
744 # Nothing to do for now...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
745 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
746 else {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
747 die "Error: The value specified, $Options{atomidentifiertype}, for option \"-a, --AtomIdentifierType\" is not valid. Supported atom identifier types in current release of MayaChemTools: AtomicInvariantsAtomTypes, DREIDINGAtomTypes, EStateAtomTypes, FunctionalClassAtomTypes, MMFF94AtomTypes, SLogPAtomTypes, SYBYLAtomTypes, TPSAAtomTypes, UFFAtomTypes\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
748 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
749 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
750
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
751 # Process specified atomic invariants to use...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
752 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
753 sub ProcessAtomicInvariantsToUseOption {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
754 my($AtomicInvariant, $AtomSymbolSpecified, @AtomicInvariantsWords);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
755
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
756 @{$OptionsInfo{AtomicInvariantsToUse}} = ();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
757 if (IsEmpty($Options{atomicinvariantstouse})) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
758 die "Error: Atomic invariants value specified using \"--AtomicInvariantsToUse\" option is empty\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
759 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
760 $AtomSymbolSpecified = 0;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
761 @AtomicInvariantsWords = split /\,/, $Options{atomicinvariantstouse};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
762 for $AtomicInvariant (@AtomicInvariantsWords) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
763 if (!AtomTypes::AtomicInvariantsAtomTypes::IsAtomicInvariantAvailable($AtomicInvariant)) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
764 die "Error: Atomic invariant specified, $AtomicInvariant, using \"--AtomicInvariantsToUse\" option is not valid...\n ";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
765 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
766 if ($AtomicInvariant =~ /^(AS|AtomSymbol)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
767 $AtomSymbolSpecified = 1;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
768 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
769 push @{$OptionsInfo{AtomicInvariantsToUse}}, $AtomicInvariant;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
770 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
771 if (!$AtomSymbolSpecified) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
772 die "Error: Atomic invariant, AS or AtomSymbol, must be specified as using \"--AtomicInvariantsToUse\" option...\n ";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
773 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
774 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
775
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
776 # Process specified functional classes invariants to use...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
777 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
778 sub ProcessFunctionalClassesToUse {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
779 my($FunctionalClass, @FunctionalClassesToUseWords);
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
780
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
781 @{$OptionsInfo{FunctionalClassesToUse}} = ();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
782 if (IsEmpty($Options{functionalclassestouse})) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
783 die "Error: Functional classes value specified using \"--FunctionalClassesToUse\" option is empty\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
784 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
785 @FunctionalClassesToUseWords = split /\,/, $Options{functionalclassestouse};
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
786 for $FunctionalClass (@FunctionalClassesToUseWords) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
787 if (!AtomTypes::FunctionalClassAtomTypes::IsFunctionalClassAvailable($FunctionalClass)) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
788 die "Error: Functional class specified, $FunctionalClass, using \"--FunctionalClassesToUse\" option is not valid...\n ";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
789 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
790 push @{$OptionsInfo{FunctionalClassesToUse}}, $FunctionalClass;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
791 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
792 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
793
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
794 # Setup script usage and retrieve command line arguments specified using various options...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
795 sub SetupScriptUsage {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
796
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
797 # Retrieve all the options...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
798 %Options = ();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
799
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
800 $Options{aromaticitymodel} = 'MayaChemToolsAromaticityModel';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
801
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
802 $Options{atomidentifiertype} = 'AtomicInvariantsAtomTypes';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
803 $Options{atomicinvariantstouse} = 'AS';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
804
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
805 $Options{functionalclassestouse} = 'HBD,HBA,PI,NI,Ar,Hal';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
806
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
807 $Options{bitsorder} = 'Ascending';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
808 $Options{bitstringformat} = 'HexadecimalString';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
809
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
810 $Options{compoundidmode} = 'LabelPrefix';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
811 $Options{compoundidlabel} = 'CompoundID';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
812 $Options{datafieldsmode} = 'CompoundID';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
813 $Options{detectaromaticity} = 'Yes';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
814
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
815 $Options{filter} = 'Yes';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
816
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
817 $Options{fold} = 'No';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
818 $Options{foldedsize} = 256;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
819
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
820 $Options{ignorehydrogens} = 'Yes';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
821 $Options{keeplargestcomponent} = 'Yes';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
822
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
823 $Options{mode} = 'PathLengthBits';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
824 $Options{pathmode} = 'AllAtomPathsWithRings';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
825
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
826 $Options{minpathlength} = 1;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
827 $Options{maxpathlength} = 8;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
828
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
829 $Options{numofbitstosetperpath} = 1;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
830
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
831 $Options{output} = 'text';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
832 $Options{outdelim} = 'comma';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
833 $Options{quote} = 'yes';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
834
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
835 $Options{size} = 1024;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
836
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
837 $Options{usebondsymbols} = 'yes';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
838 $Options{useperlcorerandom} = 'yes';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
839 $Options{useuniquepaths} = 'yes';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
840
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
841 $Options{vectorstringformat} = 'IDsAndValuesString';
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
842
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
843 if (!GetOptions(\%Options, "aromaticitymodel=s", "atomidentifiertype|a=s", "atomicinvariantstouse=s", "functionalclassestouse=s", "bitsorder=s", "bitstringformat|b=s", "compoundid=s", "compoundidlabel=s", "compoundidmode=s", "datafields=s", "datafieldsmode|d=s", "detectaromaticity=s", "filter|f=s", "fingerprintslabel=s", "fold=s", "foldedsize=i", "help|h", "ignorehydrogens|i=s", "keeplargestcomponent|k=s", "mode|m=s", "minpathlength=i", "maxpathlength=i", "numofbitstosetperpath|n=i", "outdelim=s", "output=s", "overwrite|o", "pathmode|p=s", "quote|q=s", "root|r=s", "size|s=i", "usebondsymbols|u=s", "useperlcorerandom=s", "useuniquepaths=s", "vectorstringformat|v=s", "workingdir|w=s")) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
844 die "\nTo get a list of valid options and their values, use \"$ScriptName -h\" or\n\"perl -S $ScriptName -h\" command and try again...\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
845 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
846 if ($Options{workingdir}) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
847 if (! -d $Options{workingdir}) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
848 die "Error: The value specified, $Options{workingdir}, for option \"-w --workingdir\" is not a directory name.\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
849 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
850 chdir $Options{workingdir} or die "Error: Couldn't chdir $Options{workingdir}: $! \n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
851 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
852 if (!Molecule::IsSupportedAromaticityModel($Options{aromaticitymodel})) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
853 my(@SupportedModels) = Molecule::GetSupportedAromaticityModels();
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
854 die "Error: The value specified, $Options{aromaticitymodel}, for option \"--AromaticityModel\" is not valid. Supported aromaticity models in current release of MayaChemTools: @SupportedModels\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
855 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
856 if ($Options{atomidentifiertype} !~ /^(AtomicInvariantsAtomTypes|DREIDINGAtomTypes|EStateAtomTypes|FunctionalClassAtomTypes|MMFF94AtomTypes|SLogPAtomTypes|SYBYLAtomTypes|TPSAAtomTypes|UFFAtomTypes)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
857 die "Error: The value specified, $Options{atomidentifiertype}, for option \"-a, --AtomIdentifierType\" is not valid. Supported atom identifier types in current release of MayaChemTools: AtomicInvariantsAtomTypes, DREIDINGAtomTypes, EStateAtomTypes, FunctionalClassAtomTypes, MMFF94AtomTypes, SLogPAtomTypes, SYBYLAtomTypes, TPSAAtomTypes, UFFAtomTypes\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
858 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
859 if ($Options{bitsorder} !~ /^(Ascending|Descending)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
860 die "Error: The value specified, $Options{bitsorder}, for option \"--BitsOrder\" is not valid. Allowed values: Ascending or Descending\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
861 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
862 if ($Options{bitstringformat} !~ /^(BinaryString|HexadecimalString)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
863 die "Error: The value specified, $Options{bitstringformat}, for option \"-b, --bitstringformat\" is not valid. Allowed values: BinaryString or HexadecimalString\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
864 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
865 if ($Options{compoundidmode} !~ /^(DataField|MolName|LabelPrefix|MolNameOrLabelPrefix)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
866 die "Error: The value specified, $Options{compoundidmode}, for option \"--CompoundIDMode\" is not valid. Allowed values: DataField, MolName, LabelPrefix or MolNameOrLabelPrefix\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
867 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
868 if ($Options{datafieldsmode} !~ /^(All|Common|Specify|CompoundID)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
869 die "Error: The value specified, $Options{datafieldsmode}, for option \"-d, --DataFieldsMode\" is not valid. Allowed values: All, Common, Specify or CompoundID\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
870 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
871 if ($Options{detectaromaticity} !~ /^(Yes|No)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
872 die "Error: The value specified, $Options{detectaromaticity}, for option \"--DetectAromaticity\" is not valid. Allowed values: Yes or No\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
873 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
874 if ($Options{filter} !~ /^(Yes|No)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
875 die "Error: The value specified, $Options{filter}, for option \"-f, --Filter\" is not valid. Allowed values: Yes or No\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
876 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
877 if ($Options{fold} !~ /^(Yes|No)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
878 die "Error: The value specified, $Options{fold}, for option \"--fold\" is not valid. Allowed values: Yes or No\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
879 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
880 if (!IsPositiveInteger($Options{foldedsize})) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
881 die "Error: The value specified, $Options{foldedsize}, for option \"--FoldedSize\" is not valid. Allowed values: > 0 \n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
882 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
883 if ($Options{ignorehydrogens} !~ /^(Yes|No)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
884 die "Error: The value specified, $Options{ignorehydrogens}, for option \"-i, --IgnoreHydrogens\" is not valid. Allowed values: Yes or No\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
885 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
886 if ($Options{keeplargestcomponent} !~ /^(Yes|No)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
887 die "Error: The value specified, $Options{keeplargestcomponent}, for option \"-k, --KeepLargestComponent\" is not valid. Allowed values: Yes or No\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
888 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
889 if ($Options{mode} !~ /^(PathLengthBits|PathLengthCount)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
890 die "Error: The value specified, $Options{mode}, for option \"-m, --mode\" is not valid. Allowed values: PathLengthBits or PathLengthCount\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
891 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
892 if (!IsPositiveInteger($Options{minpathlength})) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
893 die "Error: The value specified, $Options{minpathlength}, for option \"--MinPathLength\" is not valid. Allowed values: > 0 \n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
894 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
895 if (!IsPositiveInteger($Options{numofbitstosetperpath})) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
896 die "Error: The value specified, $Options{NumOfBitsToSetPerPath}, for option \"--NumOfBitsToSetPerPath\" is not valid. Allowed values: > 0 \n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
897 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
898 if (!IsPositiveInteger($Options{maxpathlength})) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
899 die "Error: The value specified, $Options{maxpathlength}, for option \"--MaxPathLength\" is not valid. Allowed values: > 0 \n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
900 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
901 if ($Options{output} !~ /^(SD|FP|text|all)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
902 die "Error: The value specified, $Options{output}, for option \"--output\" is not valid. Allowed values: SD, FP, text, or all\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
903 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
904 if ($Options{outdelim} !~ /^(comma|semicolon|tab)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
905 die "Error: The value specified, $Options{outdelim}, for option \"--outdelim\" is not valid. Allowed values: comma, tab, or semicolon\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
906 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
907 if ($Options{pathmode} !~ /^(AtomPathsWithoutRings|AtomPathsWithRings|AllAtomPathsWithoutRings|AllAtomPathsWithRings)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
908 die "Error: The value specified, $Options{pathmode}, for option \"-m, --PathMode\" is not valid. Allowed values: AtomPathsWithoutRings, AtomPathsWithRings, AllAtomPathsWithoutRings or AllAtomPathsWithRings\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
909 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
910 if ($Options{quote} !~ /^(Yes|No)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
911 die "Error: The value specified, $Options{quote}, for option \"-q --quote\" is not valid. Allowed values: Yes or No\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
912 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
913 if ($Options{outdelim} =~ /semicolon/i && $Options{quote} =~ /^No$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
914 die "Error: The value specified, $Options{quote}, for option \"-q --quote\" is not allowed with, semicolon value of \"--outdelim\" option: Fingerprints string use semicolon as delimiter for various data fields and must be quoted.\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
915 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
916
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
917 if (!IsPositiveInteger($Options{size})) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
918 die "Error: The value specified, $Options{size}, for option \"-s, --size\" is not valid. Allowed values: > 0 \n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
919 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
920 if ($Options{usebondsymbols} !~ /^(Yes|No)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
921 die "Error: The value specified, $Options{usebondsymbols}, for option \"-u, --UseBondSymbols\" is not valid. Allowed values: Yes or No\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
922 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
923 if ($Options{useperlcorerandom} !~ /^(Yes|No)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
924 die "Error: The value specified, $Options{useperlcorerandom}, for option \"--UsePerlCoreRandom\" is not valid. Allowed values: Yes or No\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
925 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
926 if ($Options{useuniquepaths} !~ /^(Yes|No)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
927 die "Error: The value specified, $Options{useuniquepaths}, for option \"--UseUniquePaths\" is not valid. Allowed values: Yes or No\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
928 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
929 if ($Options{vectorstringformat} !~ /^(IDsAndValuesString|IDsAndValuesPairsString|ValuesAndIDsString|ValuesAndIDsPairsString)$/i) {
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
930 die "Error: The value specified, $Options{vectorstringformat}, for option \"-v, --VectorStringFormat\" is not valid. Allowed values: IDsAndValuesString, IDsAndValuesPairsString, ValuesAndIDsString or ValuesAndIDsPairsString\n";
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
931 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
932 }
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
933
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
934 __END__
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
935
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
936 =head1 NAME
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
937
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
938 PathLengthFingerprints.pl - Generate atom path length based fingerprints for SD files
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
939
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
940 =head1 SYNOPSIS
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
941
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
942 PathLengthFingerprints.pl SDFile(s)...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
943
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
944 PathLengthFingerprints.pl [B<--AromaticityModel> I<AromaticityModelType>]
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
945 [B<-a, --AtomIdentifierType> I<AtomicInvariantsAtomTypes>]
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
946 [B<--AtomicInvariantsToUse> I<"AtomicInvariant1,AtomicInvariant2...">]
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
947 [B<--FunctionalClassesToUse> I<"FunctionalClass1,FunctionalClass2...">]
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
948 [B<--BitsOrder> I<Ascending | Descending>] [B<-b, --BitStringFormat> I<BinaryString | HexadecimalString>]
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
949 [B<--CompoundID> I<DataFieldName or LabelPrefixString>] [B<--CompoundIDLabel> I<text>]
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
950 [B<--CompoundIDMode> I<DataField | MolName | LabelPrefix | MolNameOrLabelPrefix>]
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
951 [B<--DataFields> I<"FieldLabel1,FieldLabel2,... ">] [B<-d, --DataFieldsMode> I<All | Common | Specify | CompoundID>]
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
952 [B<--DetectAromaticity> I<Yes | No>] [B<-f, --Filter> I<Yes | No>] [B<--FingerprintsLabel> I<text>]
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
953 [B<--fold> I<Yes | No>] [B<--FoldedSize> I<number>] [B<-h, --help>]
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
954 [B<-i, --IgnoreHydrogens> I<Yes | No>] [B<-k, --KeepLargestComponent> I<Yes | No>]
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
955 [B<-m, --mode> I<PathLengthBits | PathLengthCount>]
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
956 [B<--MinPathLength> I<number>] [B<--MaxPathLength> I<number>] [B<-n, --NumOfBitsToSetPerPath> I<number>]
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
957 [B<--OutDelim> I<comma | tab | semicolon>]
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
958 [B<--output> I<SD | FP | text | all>] [B<-q, --quote> I<Yes | No>] [B<-r, --root> I<RootName>]
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
959 [B<-p, --PathMode> I<AtomPathsWithoutRings | AtomPathsWithRings | AllAtomPathsWithoutRings | AllAtomPathsWithRings>]
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
960 [B<-s, --size> I<number>] [B<-u, --UseBondSymbols> I<Yes | No>] [B<--UsePerlCoreRandom> I<Yes | No>]
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
961 [B<--UseUniquePaths> I<Yes | No>] [B<-q, --quote> I<Yes | No>] [B<-r, --root> I<RootName>]
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
962 [B<-v, --VectorStringFormat> I<IDsAndValuesString | IDsAndValuesPairsString | ValuesAndIDsString | ValuesAndIDsPairsString>]
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
963 [B<-w, --WorkingDir> dirname] SDFile(s)...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
964
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
965 =head1 DESCRIPTION
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
966
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
967 Generate atom path length fingerprints for I<SDFile(s)> and create appropriate SD, FP or
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
968 CSV/TSV text file(s) containing fingerprints bit-vector or vector strings corresponding to
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
969 molecular fingerprints.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
970
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
971 Multiple SDFile names are separated by spaces. The valid file extensions are I<.sdf>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
972 and I<.sd>. All other file names are ignored. All the SD files in a current directory
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
973 can be specified either by I<*.sdf> or the current directory name.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
974
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
975 The current release of MayaChemTools supports generation of path length fingerprints
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
976 corresponding to following B<-a, --AtomIdentifierTypes>:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
977
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
978 AtomicInvariantsAtomTypes, DREIDINGAtomTypes, EStateAtomTypes,
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
979 FunctionalClassAtomTypes, MMFF94AtomTypes, SLogPAtomTypes,
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
980 SYBYLAtomTypes, TPSAAtomTypes, UFFAtomTypes
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
981
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
982 Based on the values specified for B<-p, --PathMode>, B<--MinPathLength> and B<--MaxPathLength>,
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
983 all appropriate atom paths are generated for each atom in the molecule and collected in a list and
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
984 the list is filtered to remove any structurally duplicate paths as indicated by the value of
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
985 B<--UseUniquePaths> option.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
986
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
987 For each atom path in the filtered atom paths list, an atom path string is created using value of
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
988 B<-a, --AtomIdentifierType> and specified values to use for a particular atom identifier type.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
989 Value of B<-u, --UseBondSymbols> controls whether bond order symbols are used during generation
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
990 of atom path string. For each atom path, only lexicographically smaller atom path strings are kept.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
991
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
992 For I<PathLengthBits> value of B<-m, --mode> option, each atom path is hashed to a 32 bit unsigned
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
993 integer key using B<TextUtil::HashCode> function. Using the hash key as a seed for a random number
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
994 generator, a random integer value between 0 and B<--Size> is used to set corresponding bits
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
995 in the fingerprint bit-vector string. Value of B<--NumOfBitsToSetPerPath> option controls the number
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
996 of time a random number is generated to set corresponding bits.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
997
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
998 For I< PathLengthCount> value of B<-m, --mode> option, the number of times an atom path appears
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
999 is tracked and a fingerprints count-string corresponding to count of atom paths is generated.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1000
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1001 Example of I<SD> file containing path length fingerprints string data:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1002
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1003 ... ...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1004 ... ...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1005 $$$$
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1006 ... ...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1007 ... ...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1008 ... ...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1009 41 44 0 0 0 0 0 0 0 0999 V2000
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1010 -3.3652 1.4499 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1011 ... ...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1012 2 3 1 0 0 0 0
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1013 ... ...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1014 M END
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1015 > <CmpdID>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1016 Cmpd1
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1017
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1018 > <PathLengthFingerprints>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1019 FingerprintsBitVector;PathLengthBits:AtomicInvariantsAtomTypes:MinLengt
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1020 h1:MaxLength8;1024;HexadecimalString;Ascending;9c8460989ec8a49913991a66
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1021 03130b0a19e8051c89184414953800cc2151082844a201042800130860308e8204d4028
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1022 00831048940e44281c00060449a5000ac80c894114e006321264401600846c050164462
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1023 08190410805000304a10205b0100e04c0038ba0fad0209c0ca8b1200012268b61c0026a
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1024 aa0660a11014a011d46
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1025
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1026 $$$$
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1027 ... ...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1028 ... ...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1029
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1030 Example of I<FP> file containing path length fingerprints string data:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1031
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1032 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1033 # Package = MayaChemTools 7.4
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1034 # ReleaseDate = Oct 21, 2010
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1035 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1036 # TimeStamp = Mon Mar 7 15:14:01 2011
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1037 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1038 # FingerprintsStringType = FingerprintsBitVector
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1039 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1040 # Description = PathLengthBits:AtomicInvariantsAtomTypes:MinLength1:...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1041 # Size = 1024
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1042 # BitStringFormat = HexadecimalString
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1043 # BitsOrder = Ascending
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1044 #
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1045 Cmpd1 9c8460989ec8a49913991a6603130b0a19e8051c89184414953800cc21510...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1046 Cmpd2 000000249400840040100042011001001980410c000000001010088001120...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1047 ... ...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1048 ... ..
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1049
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1050 Example of CSV I<Text> file containing pathlength fingerprints string data:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1051
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1052 "CompoundID","PathLengthFingerprints"
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1053 "Cmpd1","FingerprintsBitVector;PathLengthBits:AtomicInvariantsAtomTypes
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1054 :MinLength1:MaxLength8;1024;HexadecimalString;Ascending;9c8460989ec8a4
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1055 9913991a6603130b0a19e8051c89184414953800cc2151082844a20104280013086030
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1056 8e8204d402800831048940e44281c00060449a5000ac80c894114e006321264401..."
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1057 ... ...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1058 ... ...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1059
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1060 The current release of MayaChemTools generates the following types of path length
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1061 fingerprints bit-vector and vector strings:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1062
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1063 FingerprintsBitVector;PathLengthBits:AtomicInvariantsAtomTypes:MinLeng
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1064 th1:MaxLength8;1024;BinaryString;Ascending;001000010011010101011000110
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1065 0100010101011000101001011100110001000010001001101000001001001001001000
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1066 0010110100000111001001000001001010100100100000000011000000101001011100
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1067 0010000001000101010100000100111100110111011011011000000010110111001101
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1068 0101100011000000010001000011000010100011101100001000001000100000000...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1069
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1070 FingerprintsBitVector;PathLengthBits:AtomicInvariantsAtomTypes:MinLeng
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1071 th1:MaxLength8;1024;HexadecimalString;Ascending;48caa1315d82d91122b029
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1072 42861c9409a4208182d12015509767bd0867653604481a8b1288000056090583603078
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1073 9cedae54e26596889ab121309800900490515224208421502120a0dd9200509723ae89
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1074 00024181b86c0122821d4e4880c38620dab280824b455404009f082003d52c212b4e6d
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1075 6ea05280140069c780290c43
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1076
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1077 FingerprintsVector;PathLengthCount:AtomicInvariantsAtomTypes:MinLength
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1078 1:MaxLength8;432;NumericalValues;IDsAndValuesPairsString;C.X1.BO1.H3 2
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1079 C.X2.BO2.H2 4 C.X2.BO3.H1 14 C.X3.BO3.H1 3 C.X3.BO4 10 F.X1.BO1 1 N.X
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1080 2.BO2.H1 1 N.X3.BO3 1 O.X1.BO1.H1 3 O.X1.BO2 2 C.X1.BO1.H3C.X3.BO3.H1
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1081 2 C.X2.BO2.H2C.X2.BO2.H2 1 C.X2.BO2.H2C.X3.BO3.H1 4 C.X2.BO2.H2C.X3.BO
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1082 4 1 C.X2.BO2.H2N.X3.BO3 1 C.X2.BO3.H1:C.X2.BO3.H1 10 C.X2.BO3.H1:C....
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1083
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1084 FingerprintsVector;PathLengthCount:DREIDINGAtomTypes:MinLength1:MaxLen
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1085 gth8;410;NumericalValues;IDsAndValuesPairsString;C_2 2 C_3 9 C_R 22 F_
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1086 1 N_3 1 N_R 1 O_2 2 O_3 3 C_2=O_2 2 C_2C_3 1 C_2C_R 1 C_2N_3 1 C_2O_3
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1087 1 C_3C_3 7 C_3C_R 1 C_3N_R 1 C_3O_3 2 C_R:C_R 21 C_R:N_R 2 C_RC_R 2 C
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1088 _RF_ 1 C_RN_3 1 C_2C_3C_3 1 C_2C_R:C_R 2 C_2N_3C_R 1 C_3C_2=O_2 1 C_3C
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1089 _2O_3 1 C_3C_3C_3 5 C_3C_3C_R 2 C_3C_3N_R 1 C_3C_3O_3 4 C_3C_R:C_R ...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1090
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1091 FingerprintsVector;PathLengthCount:EStateAtomTypes:MinLength1:MaxLengt
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1092 h8;454;NumericalValues;IDsAndValuesPairsString;aaCH 14 aasC 8 aasN 1 d
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1093 O 2 dssC 2 sCH3 2 sF 1 sOH 3 ssCH2 4 ssNH 1 sssCH 3 aaCH:aaCH 10 aaCH:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1094 aasC 8 aasC:aasC 3 aasC:aasN 2 aasCaasC 2 aasCdssC 1 aasCsF 1 aasCssNH
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1095 1 aasCsssCH 1 aasNssCH2 1 dO=dssC 2 dssCsOH 1 dssCssCH2 1 dssCssNH 1
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1096 sCH3sssCH 2 sOHsssCH 2 ssCH2ssCH2 1 ssCH2sssCH 4 aaCH:aaCH:aaCH 6 a...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1097
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1098 FingerprintsVector;PathLengthCount:FunctionalClassAtomTypes:MinLength1
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1099 :MaxLength8;404;NumericalValues;IDsAndValuesPairsString;Ar 22 Ar.HBA 1
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1100 HBA 2 HBA.HBD 3 HBD 1 Hal 1 NI 1 None 10 Ar.HBA:Ar 2 Ar.HBANone 1 Ar:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1101 Ar 21 ArAr 2 ArHBD 1 ArHal 1 ArNone 2 HBA.HBDNI 1 HBA.HBDNone 2 HBA=NI
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1102 1 HBA=None 1 HBDNone 1 NINone 1 NoneNone 7 Ar.HBA:Ar:Ar 2 Ar.HBA:ArAr
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1103 1 Ar.HBA:ArNone 1 Ar.HBANoneNone 1 Ar:Ar.HBA:Ar 1 Ar:Ar.HBANone 2 ...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1104
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1105 FingerprintsVector;PathLengthCount:MMFF94AtomTypes:MinLength1:MaxLengt
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1106 h8;463;NumericalValues;IDsAndValuesPairsString;C5A 2 C5B 2 C=ON 1 CB 1
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1107 8 COO 1 CR 9 F 1 N5 1 NC=O 1 O=CN 1 O=CO 1 OC=O 1 OR 2 C5A:C5B 2 C5A:N
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1108 5 2 C5ACB 1 C5ACR 1 C5B:C5B 1 C5BC=ON 1 C5BCB 1 C=ON=O=CN 1 C=ONNC=O 1
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1109 CB:CB 18 CBF 1 CBNC=O 1 COO=O=CO 1 COOCR 1 COOOC=O 1 CRCR 7 CRN5 1 CR
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1110 OR 2 C5A:C5B:C5B 2 C5A:C5BC=ON 1 C5A:C5BCB 1 C5A:N5:C5A 1 C5A:N5CR ...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1111
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1112 FingerprintsVector;PathLengthCount:SLogPAtomTypes:MinLength1:MaxLength
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1113 8;518;NumericalValues;IDsAndValuesPairsString;C1 5 C10 1 C11 1 C14 1 C
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1114 18 14 C20 4 C21 2 C22 1 C5 2 CS 2 F 1 N11 1 N4 1 O10 1 O2 3 O9 1 C10C1
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1115 1 C10N11 1 C11C1 2 C11C21 1 C14:C18 2 C14F 1 C18:C18 10 C18:C20 4 C18
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1116 :C22 2 C1C5 1 C1CS 4 C20:C20 1 C20:C21 1 C20:N11 1 C20C20 2 C21:C21 1
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1117 C21:N11 1 C21C5 1 C22N4 1 C5=O10 1 C5=O9 1 C5N4 1 C5O2 1 CSO2 2 C10...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1118
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1119 FingerprintsVector;PathLengthCount:SYBYLAtomTypes:MinLength1:MaxLength
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1120 8;412;NumericalValues;IDsAndValuesPairsString;C.2 2 C.3 9 C.ar 22 F 1
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1121 N.am 1 N.ar 1 O.2 1 O.3 2 O.co2 2 C.2=O.2 1 C.2=O.co2 1 C.2C.3 1 C.2C.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1122 ar 1 C.2N.am 1 C.2O.co2 1 C.3C.3 7 C.3C.ar 1 C.3N.ar 1 C.3O.3 2 C.ar:C
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1123 .ar 21 C.ar:N.ar 2 C.arC.ar 2 C.arF 1 C.arN.am 1 C.2C.3C.3 1 C.2C.ar:C
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1124 .ar 2 C.2N.amC.ar 1 C.3C.2=O.co2 1 C.3C.2O.co2 1 C.3C.3C.3 5 C.3C.3...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1125
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1126 FingerprintsVector;PathLengthCount:TPSAAtomTypes:MinLength1:MaxLength8
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1127 ;331;NumericalValues;IDsAndValuesPairsString;N21 1 N7 1 None 34 O3 2 O
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1128 4 3 N21:None 2 N21None 1 N7None 2 None:None 21 None=O3 2 NoneNone 13 N
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1129 oneO4 3 N21:None:None 2 N21:NoneNone 2 N21NoneNone 1 N7None:None 2 N7N
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1130 one=O3 1 N7NoneNone 1 None:N21:None 1 None:N21None 2 None:None:None 20
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1131 None:NoneNone 12 NoneN7None 1 NoneNone=O3 2 NoneNoneNone 8 NoneNon...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1132
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1133 FingerprintsVector;PathLengthCount:UFFAtomTypes:MinLength1:MaxLength8;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1134 410;NumericalValues;IDsAndValuesPairsString;C_2 2 C_3 9 C_R 22 F_ 1 N_
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1135 3 1 N_R 1 O_2 2 O_3 3 C_2=O_2 2 C_2C_3 1 C_2C_R 1 C_2N_3 1 C_2O_3 1 C_
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1136 3C_3 7 C_3C_R 1 C_3N_R 1 C_3O_3 2 C_R:C_R 21 C_R:N_R 2 C_RC_R 2 C_RF_
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1137 1 C_RN_3 1 C_2C_3C_3 1 C_2C_R:C_R 2 C_2N_3C_R 1 C_3C_2=O_2 1 C_3C_2O_3
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1138 1 C_3C_3C_3 5 C_3C_3C_R 2 C_3C_3N_R 1 C_3C_3O_3 4 C_3C_R:C_R 1 C_3...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1139
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1140 =head1 OPTIONS
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1141
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1142 =over 4
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1143
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1144 =item B<--AromaticityModel> I<MDLAromaticityModel | TriposAromaticityModel | MMFFAromaticityModel | ChemAxonBasicAromaticityModel | ChemAxonGeneralAromaticityModel | DaylightAromaticityModel | MayaChemToolsAromaticityModel>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1145
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1146 Specify aromaticity model to use during detection of aromaticity. Possible values in the current
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1147 release are: I<MDLAromaticityModel, TriposAromaticityModel, MMFFAromaticityModel,
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1148 ChemAxonBasicAromaticityModel, ChemAxonGeneralAromaticityModel, DaylightAromaticityModel
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1149 or MayaChemToolsAromaticityModel>. Default value: I<MayaChemToolsAromaticityModel>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1150
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1151 The supported aromaticity model names along with model specific control parameters
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1152 are defined in B<AromaticityModelsData.csv>, which is distributed with the current release
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1153 and is available under B<lib/data> directory. B<Molecule.pm> module retrieves data from
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1154 this file during class instantiation and makes it available to method B<DetectAromaticity>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1155 for detecting aromaticity corresponding to a specific model.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1156
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1157 This option is ignored during I<No> value of B<--DetectAromaticity> option.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1158
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1159 =item B<-a, --AtomIdentifierType> I<AtomicInvariantsAtomTypes | DREIDINGAtomTypes | EStateAtomTypes | FunctionalClassAtomTypes | MMFF94AtomTypes | SLogPAtomTypes | SYBYLAtomTypes | TPSAAtomTypes | UFFAtomTypes>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1160
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1161 Specify atom identifier type to use for assignment of atom types to hydrogen and/or
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1162 non-hydrogen atoms during calculation of atom types fingerprints. Possible values in the
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1163 current release are: I<AtomicInvariantsAtomTypes, DREIDINGAtomTypes, EStateAtomTypes,
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1164 FunctionalClassAtomTypes, MMFF94AtomTypes, SLogPAtomTypes, SYBYLAtomTypes,
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1165 TPSAAtomTypes, UFFAtomTypes>. Default value: I<AtomicInvariantsAtomTypes>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1166
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1167 =item B<-a, --AtomIdentifierType> I<AtomicInvariantsAtomTypes | DREIDINGAtomTypes | EStateAtomTypes | FunctionalClassAtomTypes | MMFF94AtomTypes | SLogPAtomTypes | SYBYLAtomTypes | TPSAAtomTypes | UFFAtomTypes>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1168
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1169 Specify atom identifier type to use during generation of atom path strings
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1170 corresponding to path length fingerprints. Possible values in the current release are:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1171 I<AtomicInvariantsAtomTypes, DREIDINGAtomTypes, EStateAtomTypes,
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1172 FunctionalClassAtomTypes, MMFF94AtomTypes, SLogPAtomTypes, SYBYLAtomTypes,
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1173 TPSAAtomTypes, UFFAtomTypes>. Default value: I<AtomicInvariantsAtomTypes>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1174
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1175 =item B<--AtomicInvariantsToUse> I<"AtomicInvariant1,AtomicInvariant2...">
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1176
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1177 This value is used during I<AtomicInvariantsAtomTypes> value of B<a, --AtomIdentifierType>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1178 option. It's a list of comma separated valid atomic invariant atom types.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1179
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1180 Possible values for atomic invariants are: I<AS, X, BO, LBO, SB, DB, TB,
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1181 H, Ar, RA, FC, MN, SM>. Default value: I<AS>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1182
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1183 The atomic invariants abbreviations correspond to:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1184
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1185 AS = Atom symbol corresponding to element symbol
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1186
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1187 X<n> = Number of non-hydrogen atom neighbors or heavy atoms
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1188 BO<n> = Sum of bond orders to non-hydrogen atom neighbors or heavy atoms
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1189 LBO<n> = Largest bond order of non-hydrogen atom neighbors or heavy atoms
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1190 SB<n> = Number of single bonds to non-hydrogen atom neighbors or heavy atoms
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1191 DB<n> = Number of double bonds to non-hydrogen atom neighbors or heavy atoms
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1192 TB<n> = Number of triple bonds to non-hydrogen atom neighbors or heavy atoms
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1193 H<n> = Number of implicit and explicit hydrogens for atom
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1194 Ar = Aromatic annotation indicating whether atom is aromatic
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1195 RA = Ring atom annotation indicating whether atom is a ring
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1196 FC<+n/-n> = Formal charge assigned to atom
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1197 MN<n> = Mass number indicating isotope other than most abundant isotope
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1198 SM<n> = Spin multiplicity of atom. Possible values: 1 (singlet), 2 (doublet) or
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1199 3 (triplet)
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1200
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1201 Atom type generated by AtomTypes::AtomicInvariantsAtomTypes class corresponds to:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1202
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1203 AS.X<n>.BO<n>.LBO<n>.<SB><n>.<DB><n>.<TB><n>.H<n>.Ar.RA.FC<+n/-n>.MN<n>.SM<n>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1204
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1205 Except for AS which is a required atomic invariant in atom types, all other atomic invariants are
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1206 optional. Atom type specification doesn't include atomic invariants with zero or undefined values.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1207
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1208 In addition to usage of abbreviations for specifying atomic invariants, the following descriptive words
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1209 are also allowed:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1210
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1211 X : NumOfNonHydrogenAtomNeighbors or NumOfHeavyAtomNeighbors
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1212 BO : SumOfBondOrdersToNonHydrogenAtoms or SumOfBondOrdersToHeavyAtoms
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1213 LBO : LargestBondOrderToNonHydrogenAtoms or LargestBondOrderToHeavyAtoms
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1214 SB : NumOfSingleBondsToNonHydrogenAtoms or NumOfSingleBondsToHeavyAtoms
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1215 DB : NumOfDoubleBondsToNonHydrogenAtoms or NumOfDoubleBondsToHeavyAtoms
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1216 TB : NumOfTripleBondsToNonHydrogenAtoms or NumOfTripleBondsToHeavyAtoms
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1217 H : NumOfImplicitAndExplicitHydrogens
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1218 Ar : Aromatic
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1219 RA : RingAtom
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1220 FC : FormalCharge
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1221 MN : MassNumber
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1222 SM : SpinMultiplicity
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1223
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1224 Examples:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1225
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1226 B<Benzene>: Using value of I<AS> for B<--AtomicInvariantsToUse>, I<Yes> for B<UseBondSymbols>,
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1227 and I< AllAtomPathsWithRings> for B<-p, --PathMode>, atom path strings generated are:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1228
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1229 C C:C C:C:C C:C:C:C C:C:C:C:C C:C:C:C:C:C C:C:C:C:C:C:C
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1230
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1231 And using I<AS,X,BO> for B<--AtomicInvariantsToUse> generates following atom path
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1232 strings:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1233
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1234 C.X2.BO3 C.X2.BO3:C.X2.BO3 C.X2.BO3:C.X2.BO3:C.X2.BO3
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1235 C.X2.BO3:C.X2.BO3:C.X2.BO3:C.X2.BO3
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1236 C.X2.BO3:C.X2.BO3:C.X2.BO3:C.X2.BO3:C.X2.BO3
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1237 C.X2.BO3:C.X2.BO3:C.X2.BO3:C.X2.BO3:C.X2.BO3:C.X2.BO3
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1238 C.X2.BO3:C.X2.BO3:C.X2.BO3:C.X2.BO3:C.X2.BO3:C.X2.BO3:C.X2.BO3
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1239
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1240 B<Urea>: Using value of I<AS> for B<--AtomicInvariantsToUse>, I<Yes> for B<UseBondSymbols>,
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1241 and I< AllAtomPathsWithRings> for B<-p, --PathMode>, atom path strings are:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1242
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1243 C N O C=O CN NC=O NCN
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1244
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1245 And using I<AS,X,BO> for B<--AtomicInvariantsToUse> generates following atom path
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1246 strings:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1247
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1248 C.X3.BO4 N.X1.BO1 O.X1.BO2 C.X3.BO4=O.X1.BO2
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1249 C.X3.BO4N.X1.BO1 N.X1.BO1C.X3.BO4=O.X1.BO2
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1250 N.X1.BO1C.X3.BO4N.X1.BO1
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1251
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1252 =item B<--FunctionalClassesToUse> I<"FunctionalClass1,FunctionalClass2...">
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1253
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1254 This value is used during I<FunctionalClassAtomTypes> value of B<a, --AtomIdentifierType>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1255 option. It's a list of comma separated valid functional classes.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1256
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1257 Possible values for atom functional classes are: I<Ar, CA, H, HBA, HBD, Hal, NI, PI, RA>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1258 Default value [ Ref 24 ]: I<HBD,HBA,PI,NI,Ar,Hal>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1259
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1260 The functional class abbreviations correspond to:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1261
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1262 HBD: HydrogenBondDonor
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1263 HBA: HydrogenBondAcceptor
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1264 PI : PositivelyIonizable
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1265 NI : NegativelyIonizable
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1266 Ar : Aromatic
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1267 Hal : Halogen
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1268 H : Hydrophobic
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1269 RA : RingAtom
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1270 CA : ChainAtom
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1271
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1272 Functional class atom type specification for an atom corresponds to:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1273
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1274 Ar.CA.H.HBA.HBD.Hal.NI.PI.RA
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1275
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1276 I<AtomTypes::FunctionalClassAtomTypes> module is used to assign functional class atom
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1277 types. It uses following definitions [ Ref 60-61, Ref 65-66 ]:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1278
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1279 HydrogenBondDonor: NH, NH2, OH
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1280 HydrogenBondAcceptor: N[!H], O
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1281 PositivelyIonizable: +, NH2
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1282 NegativelyIonizable: -, C(=O)OH, S(=O)OH, P(=O)OH
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1283
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1284 =item B<--BitsOrder> I<Ascending | Descending>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1285
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1286 Bits order to use during generation of fingerprints bit-vector string for I<PathLengthBits> value of
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1287 B<-m, --mode> option. Possible values: I<Ascending, Descending>. Default: I<Ascending>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1288
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1289 I<Ascending> bit order which corresponds to first bit in each byte as the lowest bit as
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1290 opposed to the highest bit.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1291
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1292 Internally, bits are stored in I<Ascending> order using Perl vec function. Regardless
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1293 of machine order, big-endian or little-endian, vec function always considers first
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1294 string byte as the lowest byte and first bit within each byte as the lowest bit.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1295
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1296 =item B<-b, --BitStringFormat> I<BinaryString | HexadecimalString>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1297
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1298 Format of fingerprints bit-vector string data in output SD, FP or CSV/TSV text file(s) specified by
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1299 B<--output> used during I<PathLengthBits> value of B<-m, --mode> option. Possible
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1300 values: I<BinaryString, HexadecimalString>. Default value: I<HexadecimalString>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1301
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1302 I<BinaryString> corresponds to an ASCII string containing 1s and 0s. I<HexadecimalString>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1303 contains bit values in ASCII hexadecimal format.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1304
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1305 Examples:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1306
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1307 FingerprintsBitVector;PathLengthBits:AtomicInvariantsAtomTypes:MinLeng
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1308 th1:MaxLength8;1024;BinaryString;Ascending;001000010011010101011000110
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1309 0100010101011000101001011100110001000010001001101000001001001001001000
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1310 0010110100000111001001000001001010100100100000000011000000101001011100
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1311 0010000001000101010100000100111100110111011011011000000010110111001101
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1312 0101100011000000010001000011000010100011101100001000001000100000000...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1313
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1314 FingerprintsBitVector;PathLengthBits:AtomicInvariantsAtomTypes:MinLeng
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1315 th1:MaxLength8;1024;HexadecimalString;Ascending;48caa1315d82d91122b029
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1316 42861c9409a4208182d12015509767bd0867653604481a8b1288000056090583603078
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1317 9cedae54e26596889ab121309800900490515224208421502120a0dd9200509723ae89
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1318 00024181b86c0122821d4e4880c38620dab280824b455404009f082003d52c212b4e6d
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1319 6ea05280140069c780290c43
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1320
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1321 =item B<--CompoundID> I<DataFieldName or LabelPrefixString>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1322
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1323 This value is B<--CompoundIDMode> specific and indicates how compound ID is generated.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1324
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1325 For I<DataField> value of B<--CompoundIDMode> option, it corresponds to datafield label name
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1326 whose value is used as compound ID; otherwise, it's a prefix string used for generating compound
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1327 IDs like LabelPrefixString<Number>. Default value, I<Cmpd>, generates compound IDs which
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1328 look like Cmpd<Number>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1329
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1330 Examples for I<DataField> value of B<--CompoundIDMode>:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1331
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1332 MolID
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1333 ExtReg
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1334
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1335 Examples for I<LabelPrefix> or I<MolNameOrLabelPrefix> value of B<--CompoundIDMode>:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1336
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1337 Compound
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1338
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1339 The value specified above generates compound IDs which correspond to Compound<Number>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1340 instead of default value of Cmpd<Number>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1341
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1342 =item B<--CompoundIDLabel> I<text>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1343
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1344 Specify compound ID column label for FP or CSV/TSV text file(s) used during I<CompoundID> value
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1345 of B<--DataFieldsMode> option. Default: I<CompoundID>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1346
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1347 =item B<--CompoundIDMode> I<DataField | MolName | LabelPrefix | MolNameOrLabelPrefix>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1348
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1349 Specify how to generate compound IDs and write to FP or CSV/TSV text file(s) along with generated
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1350 fingerprints for I<FP | text | all> values of B<--output> option: use a I<SDFile(s)> datafield value;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1351 use molname line from I<SDFile(s)>; generate a sequential ID with specific prefix; use combination
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1352 of both MolName and LabelPrefix with usage of LabelPrefix values for empty molname lines.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1353
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1354 Possible values: I<DataField | MolName | LabelPrefix | MolNameOrLabelPrefix>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1355 Default: I<LabelPrefix>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1356
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1357 For I<MolNameAndLabelPrefix> value of B<--CompoundIDMode>, molname line in I<SDFile(s)> takes
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1358 precedence over sequential compound IDs generated using I<LabelPrefix> and only empty molname
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1359 values are replaced with sequential compound IDs.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1360
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1361 This is only used for I<CompoundID> value of B<--DataFieldsMode> option.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1362
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1363 =item B<--DataFields> I<"FieldLabel1,FieldLabel2,... ">
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1364
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1365 Comma delimited list of I<SDFiles(s)> data fields to extract and write to CSV/TSV text file(s) along
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1366 with generated fingerprints for I<text | all> values of B<--output> option.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1367
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1368 This is only used for I<Specify> value of B<--DataFieldsMode> option.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1369
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1370 Examples:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1371
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1372 Extreg
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1373 MolID,CompoundName
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1374
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1375 =item B<-d, --DataFieldsMode> I<All | Common | Specify | CompoundID>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1376
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1377 Specify how data fields in I<SDFile(s)> are transferred to output CSV/TSV text file(s) along
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1378 with generated fingerprints for I<text | all> values of B<--output> option: transfer all SD
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1379 data field; transfer SD data files common to all compounds; extract specified data fields;
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1380 generate a compound ID using molname line, a compound prefix, or a combination of both.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1381 Possible values: I<All | Common | specify | CompoundID>. Default value: I<CompoundID>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1382
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1383 =item B<--DetectAromaticity> I<Yes | No>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1384
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1385 Detect aromaticity before generating fingerprints. Possible values: I<Yes or No>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1386 Default value: I<Yes>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1387
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1388 I<No> B<--DetectAromaticity> forces usage of atom and bond aromaticity values
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1389 from I<SDFile(s)> and skips the step which detects and assigns aromaticity.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1390
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1391 I<No> B<--DetectAromaticity> value is only allowed uring I<AtomicInvariantsAtomTypes>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1392 value of B<-a, --AtomIdentifierType> options; for all possible values B<-a, --AtomIdentifierType>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1393 values, it must be I<Yes>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1394
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1395 =item B<-f, --Filter> I<Yes | No>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1396
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1397 Specify whether to check and filter compound data in SDFile(s). Possible values: I<Yes or No>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1398 Default value: I<Yes>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1399
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1400 By default, compound data is checked before calculating fingerprints and compounds containing
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1401 atom data corresponding to non-element symbols or no atom data are ignored.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1402
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1403 =item B<--FingerprintsLabel> I<text>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1404
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1405 SD data label or text file column label to use for fingerprints string in output SD or
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1406 CSV/TSV text file(s) specified by B<--output>. Default value: I<PathLenghFingerprints>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1407
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1408 =item B<--fold> I<Yes | No>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1409
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1410 Fold fingerprints to increase bit density during I<PathLengthBits> value of
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1411 B<-m, --mode> option. Possible values: I<Yes or No>. Default value: I<No>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1412
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1413 =item B<--FoldedSize> I<number>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1414
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1415 Size of folded fingerprint during I<PathLengthBits> value of B<-m, --mode> option. Default
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1416 value: I<256>. Valid values correspond to any positive integer which is less than
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1417 B<-s, --size> and meets the criteria for its value.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1418
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1419 Examples:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1420
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1421 128
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1422 512
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1423
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1424 =item B<-h, --help>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1425
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1426 Print this help message
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1427
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1428 =item B<-i, --IgnoreHydrogens> I<Yes | No>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1429
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1430 Ignore hydrogens during fingerprints generation. Possible values: I<Yes or No>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1431 Default value: I<Yes>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1432
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1433 For I<yes> value of B<-i, --IgnoreHydrogens>, any explicit hydrogens are also used for
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1434 generation of atoms path lengths and fingerprints; implicit hydrogens are still ignored.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1435
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1436 =item B<-k, --KeepLargestComponent> I<Yes | No>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1437
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1438 Generate fingerprints for only the largest component in molecule. Possible values:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1439 I<Yes or No>. Default value: I<Yes>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1440
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1441 For molecules containing multiple connected components, fingerprints can be generated
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1442 in two different ways: use all connected components or just the largest connected
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1443 component. By default, all atoms except for the largest connected component are
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1444 deleted before generation of fingerprints.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1445
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1446 =item B<-m, --mode> I<PathLengthBits | PathLengthCount>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1447
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1448 Specify type of path length fingerprints to generate for molecules in I<SDFile(s)>. Possible
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1449 values: I<PathLengthBits, PathLengthCount>. Default value: I<PathLengthBits>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1450
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1451 For I<PathLengthBits> value of B<-m, --mode> option, a fingerprint bit-vector string containing
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1452 zeros and ones is generated and for I<PathLengthCount> value, a fingerprint vector string
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1453 corresponding to number of atom paths is generated.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1454
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1455 =item B<--MinPathLength> I<number>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1456
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1457 Minimum atom path length to include in fingerprints. Default value: I<1>. Valid values:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1458 positive integers and less than B<--MaxPathLength>. Path length of 1 correspond to
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1459 a path containing only one atom.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1460
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1461 =item B<--MaxPathLength> I<number>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1462
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1463 Maximum atom path length to include in fingerprints. Default value: I<8>. Valid values:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1464 positive integers and greater than B<--MinPathLength>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1465
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1466 =item B<-n, --NumOfBitsToSetPerPath> I<number>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1467
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1468 Number of bits to set per path during generation of fingerprints bit-vector string for I<PathLengthBits>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1469 value of B<-m, --mode> option. Default value: I<1>. Valid values: positive integers.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1470
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1471 =item B<--OutDelim> I<comma | tab | semicolon>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1472
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1473 Delimiter for output CSV/TSV text file(s). Possible values: I<comma, tab, or semicolon>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1474 Default value: I<comma>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1475
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1476 =item B<--output> I<SD | FP | text | all>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1477
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1478 Type of output files to generate. Possible values: I<SD, FP, text, or all>. Default value: I<text>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1479
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1480 =item B<-o, --overwrite>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1481
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1482 Overwrite existing files.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1483
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1484 =item B<-p, --PathMode> I<AtomPathsWithoutRings | AtomPathsWithRings | AllAtomPathsWithoutRings | AllAtomPathsWithRings>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1485
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1486 Specify type of atom paths to use for generating pathlength fingerprints for molecules in
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1487 I<SDFile(s)>. Possible values:I<AtomPathsWithoutRings, AtomPathsWithRings,
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1488 AllAtomPathsWithoutRings, AllAtomPathsWithRings>. Default value: I<AllAtomPathsWithRings>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1489
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1490 For molecules with no rings, first two and last two options are equivalent and generate
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1491 same set of atom paths starting from each atom with length between B<--MinPathLength>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1492 and B<--MaxPathLength>. However, all these four options can result in the same set of
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1493 final atom paths for molecules containing fused, bridged or spiro rings.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1494
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1495 For molecules containing rings, atom paths starting from each atom can be traversed in
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1496 four different ways:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1497
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1498 I<AtomPathsWithoutRings> - Atom paths containing no rings and without sharing of bonds
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1499 in traversed paths.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1500
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1501 I<AtomPathsWithRings> - Atom paths containing rings and without any sharing of bonds in
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1502 traversed paths.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1503
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1504 I<AllAtomPathsWithoutRings> - All possible atom paths containing no rings and without any
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1505 sharing of bonds in traversed paths.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1506
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1507 I<AllAtomPathsWithRings> - All possible atom paths containing rings and with sharing of
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1508 bonds in traversed paths.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1509
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1510 Atom path traversal is terminated at the ring atom.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1511
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1512 Based on values specified for for B<-p, --PathMode>, B<--MinPathLength> and
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1513 B<--MaxPathLength>, all appropriate atom paths are generated for each atom in the molecule
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1514 and collected in a list.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1515
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1516 For each atom path in the filtered atom paths list, an atom path string is created using value of
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1517 B<-a, --AtomIdentifierType> and specified values to use for a particular atom identifier type.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1518 Value of B<-u, --UseBondSymbols> controls whether bond order symbols are used during generation
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1519 of atom path string. Atom symbol corresponds to element symbol and characters used to represent
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1520 bond order are: I<1 - None; 2 - '='; 3 - '#'; 1.5 or aromatic - ':'; others: bond order value>. By default,
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1521 bond symbols are included in atom path strings. Exclusion of bond symbols in atom path strings
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1522 results in fingerprints which correspond purely to atom paths without considering bonds.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1523
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1524 B<UseUniquePaths> controls the removal of structurally duplicate atom path strings are removed
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1525 from the list.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1526
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1527 For I<PathLengthBits> value of B<-m, --mode> option, each atom path is hashed to a 32 bit unsigned
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1528 integer key using B<TextUtil::HashCode> function. Using the hash key as a seed for a random number
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1529 generator, a random integer value between 0 and B<--Size> is used to set corresponding bits
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1530 in the fingerprint bit-vector string. Value of B<--NumOfBitsToSetPerPaths> option controls the number
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1531 of time a random number is generated to set corresponding bits.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1532
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1533 For I< PathLengthCount> value of B<-m, --mode> option, the number of times an atom path appears
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1534 is tracked and a fingerprints count-string corresponding to count of atom paths is generated.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1535
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1536 For molecule containing rings, combination of B<-p, --PathMode> and B<--UseBondSymbols> allows
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1537 generation of up to 8 different types of atom path length strings:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1538
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1539 AllowSharedBonds AllowRings UseBondSymbols
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1540
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1541 0 0 1 - AtomPathsNoCyclesWithBondSymbols
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1542 0 1 1 - AtomPathsWithCyclesWithBondSymbols
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1543
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1544 1 0 1 - AllAtomPathsNoCyclesWithBondSymbols
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1545 1 1 1 - AllAtomPathsWithCyclesWithBondSymbols
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1546 [ DEFAULT ]
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1547
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1548 0 0 0 - AtomPathsNoCyclesNoBondSymbols
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1549 0 1 0 - AtomPathsWithCyclesNoBondSymbols
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1550
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1551 1 0 0 - AllAtomPathsNoCyclesNoBondSymbols
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1552 1 1 0 - AllAtomPathsWithCyclesNoWithBondSymbols
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1553
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1554 Default atom path length fingerprints generation for molecules containing rings with
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1555 I<AllAtomPathsWithRings> value for B<-p, --PathMode>, I<Yes> value for B<--UseBondSymbols>,
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1556 I<2> value for B<--MinPathLength> and I<8> value for B<--MaxPathLength> is the most time
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1557 consuming. Combinations of other options can substantially speed up fingerprint generation
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1558 for molecules containing complex ring systems.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1559
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1560 Additionally, value for option B<-a, --AtomIdentifierType> in conjunction with corresponding specified
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1561 values for atom types changes the nature of atom path length strings and the fingerprints.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1562
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1563 =item B<-q, --quote> I<Yes | No>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1564
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1565 Put quote around column values in output CSV/TSV text file(s). Possible values:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1566 I<Yes or No>. Default value: I<Yes>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1567
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1568 =item B<-r, --root> I<RootName>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1569
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1570 New file name is generated using the root: <Root>.<Ext>. Default for new file
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1571 names: <SDFileName><PathLengthFP>.<Ext>. The file type determines <Ext> value.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1572 The sdf, fpf, csv, and tsv <Ext> values are used for SD, FP, comma/semicolon, and tab
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1573 delimited text files, respectively.This option is ignored for multiple input files.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1574
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1575 =item B<-s, --size> I<number>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1576
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1577 Size of fingerprints. Default value: I<1024>. Valid values correspond to any positive
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1578 integer which satisfies the following criteria: power of 2, >= 32 and <= 2 ** 32.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1579
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1580 Examples:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1581
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1582 256
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1583 512
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1584 2048
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1585
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1586 =item B<-u, --UseBondSymbols> I<Yes | No>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1587
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1588 Specify whether to use bond symbols for atom paths during generation of atom path strings.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1589 Possible values: I<Yes or No>. Default value: I<Yes>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1590
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1591 I<No> value option for B<-u, --UseBondSymbols> allows the generation of fingerprints corresponding
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1592 purely to atoms disregarding all bonds.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1593
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1594 =item B<--UsePerlCoreRandom> I<Yes | No>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1595
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1596 Specify whether to use Perl CORE::rand or MayaChemTools MathUtil::random function
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1597 during random number generation for setting bits in fingerprints bit-vector strings. Possible
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1598 values: I<Yes or No>. Default value: I<Yes>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1599
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1600 I<No> value option for B<--UsePerlCoreRandom> allows the generation of fingerprints
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1601 bit-vector strings which are same across different platforms.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1602
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1603 The random number generator implemented in MayaChemTools is a variant of
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1604 linear congruential generator (LCG) as described by Miller et al. [ Ref 120 ].
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1605 It is also referred to as Lehmer random number generator or Park-Miller
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1606 random number generator.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1607
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1608 Unlike Perl's core random number generator function rand, the random number
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1609 generator implemented in MayaChemTools, MathUtil::random, generates consistent
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1610 random values across different platforms for a specific random seed and leads
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1611 to generation of portable fingerprints bit-vector strings.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1612
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1613 =item B<--UseUniquePaths> I<Yes | No>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1614
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1615 Specify whether to use structurally unique atom paths during generation of atom path strings.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1616 Possible values: I<Yes or No>. Default value: I<Yes>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1617
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1618 I<No> value option for B<--UseUniquePaths> allows usage of all atom paths generated by
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1619 B<-p, --PathMode> option value for generation of atom path strings leading to duplicate
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1620 path count during I<PathLengthCount> value of B<-m, --mode> option. It doesn't change fingerprint
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1621 string generated during I<PathLengthBits> value of B<-m, --mode>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1622
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1623 For example, during I<AllAtomPathsWithRings> value of B<-p, --PathMode> option, benzene has
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1624 12 linear paths of length 2 and 12 cyclic paths length of 7, but only 6 linear paths of length 2 and
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1625 1 cyclic path of length 7 are structurally unique.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1626
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1627 =item B<-v, --VectorStringFormat> I<IDsAndValuesString | IDsAndValuesPairsString | ValuesAndIDsString | ValuesAndIDsPairsString>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1628
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1629 Format of fingerprints vector string data in output SD, FP or CSV/TSV text file(s) specified by
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1630 B<--output> used during I<PathLengthCount> value of B<-m, --mode> option. Possible
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1631 values: I<IDsAndValuesString | IDsAndValuesPairsString | ValuesAndIDsString |
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1632 ValuesAndIDsPairsString>. Defaultvalue: I<IDsAndValuesString>.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1633
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1634 Examples:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1635
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1636 FingerprintsVector;PathLengthCount:AtomicInvariantsAtomTypes:MinLength
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1637 1:MaxLength8;432;NumericalValues;IDsAndValuesPairsString;C.X1.BO1.H3 2
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1638 C.X2.BO2.H2 4 C.X2.BO3.H1 14 C.X3.BO3.H1 3 C.X3.BO4 10 F.X1.BO1 1 N.X
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1639 2.BO2.H1 1 N.X3.BO3 1 O.X1.BO1.H1 3 O.X1.BO2 2 C.X1.BO1.H3C.X3.BO3.H1
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1640 2 C.X2.BO2.H2C.X2.BO2.H2 1 C.X2.BO2.H2C.X3.BO3.H1 4 C.X2.BO2.H2C.X3.BO
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1641 4 1 C.X2.BO2.H2N.X3.BO3 1 C.X2.BO3.H1:C.X2.BO3.H1 10 C.X2.BO3.H1:C....
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1642
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1643 FingerprintsVector;PathLengthCount:EStateAtomTypes:MinLength1:MaxLengt
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1644 h8;454;NumericalValues;IDsAndValuesPairsString;aaCH 14 aasC 8 aasN 1 d
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1645 O 2 dssC 2 sCH3 2 sF 1 sOH 3 ssCH2 4 ssNH 1 sssCH 3 aaCH:aaCH 10 aaCH:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1646 aasC 8 aasC:aasC 3 aasC:aasN 2 aasCaasC 2 aasCdssC 1 aasCsF 1 aasCssNH
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1647 1 aasCsssCH 1 aasNssCH2 1 dO=dssC 2 dssCsOH 1 dssCssCH2 1 dssCssNH 1
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1648 sCH3sssCH 2 sOHsssCH 2 ssCH2ssCH2 1 ssCH2sssCH 4 aaCH:aaCH:aaCH 6 a...
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1649
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1650 =item B<-w, --WorkingDir> I<DirName>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1651
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1652 Location of working directory. Default: current directory.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1653
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1654 =back
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1655
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1656 =head1 EXAMPLES
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1657
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1658 To generate path length fingerprints corresponding to all unique paths from length 1
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1659 through 8 in hexadecimal bit-vector string format of size 1024 and create a
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1660 SamplePLFPHex.csv file containing sequential compound IDs along with fingerprints
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1661 bit-vector strings data, type:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1662
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1663 % PathLengthFingerprints.pl -o -r SamplePLFPHex Sample.sdf
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1664
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1665 To generate path length fingerprints corresponding to all unique paths from length 1
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1666 through 8 in hexadecimal bit-vector string format of size 1024 and create SamplePLFPHex.sdf,
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1667 SamplePLFPHex.fpf, and SamplePLFPHex.csv files containing sequential compound IDs
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1668 in CSV file along with fingerprints bit-vector strings data, type:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1669
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1670 % PathLengthFingerprints.pl --output all -o -r SamplePLFPHex Sample.sdf
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1671
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1672 To generate path length fingerprints corresponding to all unique paths from length 1
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1673 through 8 in binary bit-vector string format of size 1024 and create a
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1674 SamplePLFPBin.csv file containing sequential compound IDs along with fingerprints
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1675 bit-vector strings data, type:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1676
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1677 % PathLengthFingerprints.pl --BitStringFormat BinaryString --size 2048
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1678 -o -r SamplePLFPBin Sample.sdf
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1679
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1680 To generate path length fingerprints corresponding to count of all unique paths from
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1681 length 1 through 8 in IDsAndValuesString format and create a SamplePLFPCount.csv file
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1682 containing sequential compound IDs along with fingerprints vector strings data, type:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1683
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1684 % PathLengthFingerprints.pl -m PathLengthCount -o -r SamplePLFPCount
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1685 Sample.sdf
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1686
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1687 To generate path length fingerprints corresponding to count of all unique paths from
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1688 length 1 through 8 in IDsAndValuesString format using E-state atom types and
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1689 create a SamplePLFPCount.csv file containing sequential compound IDs along with fingerprints
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1690 vector strings data, type:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1691
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1692 % PathLengthFingerprints.pl -m PathLengthCount --AtomIdentifierType
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1693 EStateAtomTypes -o -r SamplePLFPCount Sample.sdf
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1694
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1695 To generate path length fingerprints corresponding to count of all unique paths from
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1696 length 1 through 8 in IDsAndValuesString format using SLogP atom types and
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1697 create a SamplePLFPCount.csv file containing sequential compound IDs along with fingerprints
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1698 vector strings data, type:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1699
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1700 % PathLengthFingerprints.pl -m PathLengthCount --AtomIdentifierType
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1701 SLogPAtomTypes -o -r SamplePLFPCount Sample.sdf
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1702
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1703 To generate path length fingerprints corresponding to count of all unique paths from
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1704 length 1 through 8 in IDsAndValuesString format and create a SamplePLFPCount.csv file
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1705 containing sequential compound IDs along with fingerprints vector strings data, type:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1706
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1707 % PathLengthFingerprints.pl -m PathLengthCount --VectorStringFormat
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1708 ValuesAndIDsPairsString -o -r SamplePLFPCount Sample.sdf
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1709
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1710 To generate path length fingerprints corresponding to count of all unique paths from
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1711 length 1 through 8 in IDsAndValuesString format using AS,X,BO as atomic invariants and
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1712 create a SamplePLFPCount.csv file containing sequential compound IDs along with fingerprints
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1713 vector strings data, type:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1714
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1715 % PathLengthFingerprints.pl -m PathLengthCount --AtomIdentifierType
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1716 AtomicInvariantsAtomTypes --AtomicInvariantsToUse "AS,X,BO" -o
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1717 -r SamplePLFPCount Sample.sdf
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1718
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1719 To generate path length fingerprints corresponding to count of all paths from
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1720 length 1 through 8 in IDsAndValuesString format and create a SamplePLFPCount.csv file
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1721 containing compound IDs from MolName line along with fingerprints vector strings data, type:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1722
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1723 % PathLengthFingerprints.pl -m PathLengthCount --UseUniquePaths No
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1724 -o --CompoundIDMode MolName -r SamplePLFPCount --UseUniquePaths No
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1725 Sample.sdf
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1726
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1727 To generate path length fingerprints corresponding to all unique paths from length 1
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1728 through 8 in hexadecimal bit-vector string format of size 512 after folding and create
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1729 SamplePLFPHex.sdf, SamplePLFPHex.fpf, and SamplePLFPHex.sdf files containing sequential
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1730 compound IDs along with fingerprints bit-vector strings data, type:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1731
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1732 % PathLengthFingerprints.pl --output all --Fold Yes --FoldedSize 512
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1733 -o -r SamplePLFPHex Sample.sdf
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1734
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1735 To generate path length fingerprints corresponding to all unique paths from length 1
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1736 through 8 containing no rings and without sharing of bonds in hexadecimal bit-vector
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1737 string format of size 1024 and create a SamplePLFPHex.csv file containing sequential
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1738 compound IDs along with fingerprints bit-vector strings data and all data fields, type:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1739
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1740 % PathLengthFingerprints.pl -p AtomPathsWithoutRings --DataFieldsMode All
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1741 -o -r SamplePLFPHex Sample.sdf
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1742
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1743 To generate path length fingerprints corresponding to all unique paths from length 1
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1744 through 8 containing rings and without sharing of bonds in hexadecimal bit-vector
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1745 string format of size 1024 and create a SamplePLFPHex.tsv file containing compound IDs
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1746 derived from combination of molecule name line and an explicit compound prefix
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1747 along with fingerprints bit-vector strings data and all data fields, type:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1748
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1749 % PathLengthFingerprints.pl -p AtomPathsWithRings --DataFieldsMode
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1750 CompoundID --CompoundIDMode MolnameOrLabelPrefix --CompoundID Cmpd
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1751 --CompoundIDLabel MolID --FingerprintsLabel PathLengthFP --OutDelim Tab
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1752 -r SamplePLFPHex -o Sample.sdf
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1753
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1754 To generate path length fingerprints corresponding to count of all unique paths from
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1755 length 1 through 8 in IDsAndValuesString format and create a SamplePLFPCount.csv file
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1756 containing sequential compound IDs along with fingerprints vector strings data using
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1757 aromaticity specified in SD file, type:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1758
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1759 % PathLengthFingerprints.pl -m PathLengthCount --DetectAromaticity No
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1760 -o -r SamplePLFPCount Sample.sdf
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1761
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1762 To generate path length fingerprints corresponding to all unique paths from length 2
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1763 through 6 in hexadecimal bit-vector string format of size 1024 and create a
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1764 SamplePLFPHex.csv file containing sequential compound IDs along with fingerprints
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1765 bit-vector strings data, type:
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1766
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1767 % PathLengthFingerprints.pl --MinPathLength 2 --MaxPathLength 6
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1768 -o -r SamplePLFPHex Sample.sdf
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1769
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1770 =head1 AUTHOR
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1771
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1772 Manish Sud <msud@san.rr.com>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1773
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1774 =head1 SEE ALSO
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1775
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1776 InfoFingerprintsFiles.pl, SimilarityMatricesFingerprints.pl, AtomNeighborhoodsFingerprints.pl,
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1777 ExtendedConnectivityFingerprints.pl, MACCSKeysFingerprints.pl,
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1778 TopologicalAtomPairsFingerprints.pl, TopologicalAtomTorsionsFingerprints.pl,
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1779 TopologicalPharmacophoreAtomPairsFingerprints.pl, TopologicalPharmacophoreAtomTripletsFingerprints.pl
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1780
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1781 =head1 COPYRIGHT
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1782
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1783 Copyright (C) 2015 Manish Sud. All rights reserved.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1784
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1785 This file is part of MayaChemTools.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1786
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1787 MayaChemTools is free software; you can redistribute it and/or modify it under
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1788 the terms of the GNU Lesser General Public License as published by the Free
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1789 Software Foundation; either version 3 of the License, or (at your option)
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1790 any later version.
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1791
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1792 =cut