annotate mayachemtool/mayachemtools/lib/Fingerprints/EStateIndiciesFingerprints.pm @ 0:a4a2ad5a214e draft default tip

Uploaded
author deepakjadmin
date Thu, 05 Nov 2015 02:37:56 -0500
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
1 package Fingerprints::EStateIndiciesFingerprints;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
2 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
3 # $RCSfile: EStateIndiciesFingerprints.pm,v $
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
4 # $Date: 2015/02/28 20:48:54 $
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
5 # $Revision: 1.19 $
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
6 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
7 # Author: Manish Sud <msud@san.rr.com>
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
8 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
9 # Copyright (C) 2015 Manish Sud. All rights reserved.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
10 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
11 # This file is part of MayaChemTools.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
12 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
13 # MayaChemTools is free software; you can redistribute it and/or modify it under
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
14 # the terms of the GNU Lesser General Public License as published by the Free
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
15 # Software Foundation; either version 3 of the License, or (at your option) any
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
16 # later version.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
17 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
18 # MayaChemTools is distributed in the hope that it will be useful, but without
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
19 # any warranty; without even the implied warranty of merchantability of fitness
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
20 # for a particular purpose. See the GNU Lesser General Public License for more
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
21 # details.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
22 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
23 # You should have received a copy of the GNU Lesser General Public License
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
24 # along with MayaChemTools; if not, see <http://www.gnu.org/licenses/> or
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
25 # write to the Free Software Foundation Inc., 59 Temple Place, Suite 330,
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
26 # Boston, MA, 02111-1307, USA.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
27 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
28
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
29 use strict;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
30 use Carp;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
31 use Exporter;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
32 use Text::ParseWords;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
33 use TextUtil ();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
34 use FileUtil ();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
35 use MathUtil ();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
36 use Fingerprints::Fingerprints;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
37 use Molecule;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
38 use AtomTypes::EStateAtomTypes;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
39 use AtomicDescriptors::EStateValuesDescriptors;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
40
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
41 use vars qw(@ISA @EXPORT @EXPORT_OK %EXPORT_TAGS);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
42
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
43 @ISA = qw(Fingerprints::Fingerprints Exporter);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
44 @EXPORT = qw();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
45 @EXPORT_OK = qw();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
46
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
47 %EXPORT_TAGS = (all => [@EXPORT, @EXPORT_OK]);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
48
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
49 # Setup class variables...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
50 my($ClassName);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
51 _InitializeClass();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
52
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
53 # Overload Perl functions...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
54 use overload '""' => 'StringifyEStateIndiciesFingerprints';
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
55
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
56 # Class constructor...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
57 sub new {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
58 my($Class, %NamesAndValues) = @_;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
59
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
60 # Initialize object...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
61 my $This = $Class->SUPER::new();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
62 bless $This, ref($Class) || $Class;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
63 $This->_InitializeEStateIndiciesFingerprints();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
64
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
65 $This->_InitializeEStateIndiciesFingerprintsProperties(%NamesAndValues);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
66
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
67 return $This;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
68 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
69
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
70 # Initialize object data...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
71 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
72 sub _InitializeEStateIndiciesFingerprints {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
73 my($This) = @_;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
74
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
75 # EStateIndicies is a vector containing sum of E-state values for E-state atom types
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
76 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
77 $This->{Type} = 'EStateIndicies';
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
78
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
79 # EStateAtomTypesSetToUse for EStateIndicies:
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
80 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
81 # ArbitrarySize - Corrresponds to only E-state atom types detected in molecule
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
82 # FixedSize - Corresponds to fixed number of E-state atom types previously defined [ Ref 77 ]
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
83 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
84 # The default EStateAtomTypesSetToUse value for EStateIndicies fingerprints type: ArbitrarySize.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
85 # Possible values: ArbitrarySize or FixedSize.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
86 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
87 $This->{EStateAtomTypesSetToUse} = '';
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
88
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
89 # Assigned E-state atom types...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
90 %{$This->{EStateAtomTypes}} = ();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
91
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
92 # Vector values precision for real values during E-state indicies...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
93 $This->{ValuesPrecision} = 3;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
94
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
95 # Calculated E-state values and indicies for generating E-state indicies fingerprints...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
96 %{$This->{EStateValues}} = ();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
97 %{$This->{EStateIndicies}} = ();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
98 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
99
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
100 # Initialize class ...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
101 sub _InitializeClass {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
102 #Class name...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
103 $ClassName = __PACKAGE__;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
104
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
105 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
106
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
107 # Initialize object properties....
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
108 sub _InitializeEStateIndiciesFingerprintsProperties {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
109 my($This, %NamesAndValues) = @_;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
110
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
111 my($Name, $Value, $MethodName);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
112 while (($Name, $Value) = each %NamesAndValues) {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
113 $MethodName = "Set${Name}";
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
114 $This->$MethodName($Value);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
115 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
116
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
117 # Make sure molecule object was specified...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
118 if (!exists $NamesAndValues{Molecule}) {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
119 croak "Error: ${ClassName}->New: Object can't be instantiated without specifying molecule...";
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
120 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
121
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
122 $This->_InitializeEstateIndicies();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
123
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
124 return $This;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
125 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
126
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
127 # Initialize E-state indicies...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
128 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
129 sub _InitializeEstateIndicies {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
130 my($This) = @_;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
131
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
132 # Set default EStateAtomTypesSetToUse...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
133 if (!$This->{EStateAtomTypesSetToUse}) {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
134 $This->{EStateAtomTypesSetToUse} = 'ArbitrarySize';
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
135 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
136
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
137 # Vector type...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
138 $This->{VectorType} = 'FingerprintsVector';
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
139
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
140 if ($This->{EStateAtomTypesSetToUse} =~ /^FixedSize$/i) {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
141 $This->{FingerprintsVectorType} = 'OrderedNumericalValues';
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
142 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
143 else {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
144 $This->{FingerprintsVectorType} = 'NumericalValues';
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
145 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
146
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
147 $This->_InitializeFingerprintsVector();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
148
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
149 return $This;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
150 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
151
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
152 # Disable set size method...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
153 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
154 sub SetSize {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
155 my($This, $Type) = @_;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
156
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
157 croak "Error: ${ClassName}->SetSize: Can't change size: It's not allowed...";
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
158 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
159
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
160 # Set E-state atom types set to use...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
161 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
162 sub SetEStateAtomTypesSetToUse {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
163 my($This, $Value) = @_;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
164
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
165 if ($This->{EStateAtomTypesSetToUse}) {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
166 croak "Error: ${ClassName}->SetEStateAtomTypesSetToUse: Can't change size: It's already set...";
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
167 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
168
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
169 if ($Value !~ /^(ArbitrarySize|FixedSize)/i) {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
170 croak "Error: ${ClassName}->SetEStateAtomTypesSetToUse: Unknown EStateAtomTypesSetToUse value: $Value; Supported values: ArbitrarySize or FixedSize";
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
171 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
172
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
173 $This->{EStateAtomTypesSetToUse} = $Value;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
174
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
175 return $This;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
176 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
177
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
178 # Set vector values precision for real values for E-state indicies...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
179 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
180 sub SetValuesPrecision {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
181 my($This, $Value) = @_;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
182
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
183 if (!TextUtil::IsPositiveInteger($Value)) {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
184 croak "Error: ${ClassName}->SetValuesPrecision: ValuesPrecision value, $Value, is not valid: It must be a positive integer...";
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
185 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
186 $This->{ValuesPrecision} = $Value;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
187
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
188 return $This;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
189 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
190
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
191 # Generate fingerprints description...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
192 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
193 sub GetDescription {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
194 my($This) = @_;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
195
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
196 # Is description explicity set?
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
197 if (exists $This->{Description}) {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
198 return $This->{Description};
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
199 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
200
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
201 # Generate fingerprints description...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
202
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
203 return "$This->{Type}:$This->{EStateAtomTypesSetToUse}";
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
204 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
205
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
206 # Generate electrotopological state indicies (E-state) [ Ref 75-78 ] fingerprints for
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
207 # non-hydrogen atoms in a molecule...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
208 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
209 # EStateIndicies fingerprints constitute a vector containing sum of E-state values
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
210 # for E-state atom types. Two types of E-state atom types set size are allowed:
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
211 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
212 # ArbitrarySize - Corrresponds to only E-state atom types detected in molecule
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
213 # FixedSize - Corresponds to fixed number of E-state atom types previously defined
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
214 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
215 # Module AtomTypes::EStateAtomTypes.pm is used to assign E-state atom types to
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
216 # non-hydrogen atoms in the molecule which is able to assign atom types to any valid
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
217 # atom group. However, for FixedSize value of EStateAtomTypesSetToUse, only a fixed
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
218 # set of E-state atom types corresponding to specific atom groups [ Appendix III in
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
219 # Ref 77 ] are used for fingerprints.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
220 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
221 # The fixed size E-state atom type set size used during generation of fingerprints corresponding
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
222 # FixedSize value of EStateAtomTypesSetToUse contains 87 E-state non-hydrogen atom types
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
223 # in EStateAtomTypes.csv data file distributed with MayaChemTools.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
224 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
225 # Combination of Type and EStateAtomTypesSetToUse allow generation of 2 different types of
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
226 # E-state indicies fingerprints:
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
227 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
228 # Type EStateAtomTypesSetToUse
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
229 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
230 # EStateIndicies ArbitrarySize [ default fingerprints ]
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
231 # EStateIndicies FixedSize
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
232 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
233 # The default is generate EStateIndicies type fingeprints corresponding to ArbitrarySize as
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
234 # EStateAtomTypesSetToUse value.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
235 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
236 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
237 sub GenerateFingerprints {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
238 my($This) = @_;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
239
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
240 # Cache appropriate molecule data...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
241 $This->_SetupMoleculeDataCache();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
242
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
243 # Assign E-state atom types...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
244 if (!$This->_AssignEStateAtomTypes()) {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
245 carp "Warning: ${ClassName}->GenerateFingerprints: $This->{Type} fingerprints generation didn't succeed: Couldn't assign valid E-state atom types to all atoms...";
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
246 return $This;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
247 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
248
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
249 # Calculate E-state indicies...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
250 if (!$This->_CalculateEStateIndicies()) {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
251 carp "Warning: ${ClassName}->GenerateFingerprints: $This->{Type} fingerprints generation didn't succeed: Couldn't calculate E-state values for all atoms...";
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
252 return $This;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
253 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
254
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
255 # Set final fingerprints...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
256 $This->_SetFinalFingerprints();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
257
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
258 # Clear cached molecule data...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
259 $This->_ClearMoleculeDataCache();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
260
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
261 return $This;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
262 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
263
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
264 # Assign E-state atom types...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
265 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
266 sub _AssignEStateAtomTypes {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
267 my($This) = @_;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
268 my($EStateAtomTypes, $Atom, $AtomID, $AtomType);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
269
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
270 %{$This->{EStateAtomTypes}} = ();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
271
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
272 # Assign E-state atom types...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
273 $EStateAtomTypes = new AtomTypes::EStateAtomTypes('Molecule' => $This->{Molecule}, 'IgnoreHydrogens' => 1);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
274 $EStateAtomTypes->AssignAtomTypes();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
275
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
276 # Make sure atom types assignment is successful...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
277 if (!$EStateAtomTypes->IsAtomTypesAssignmentSuccessful()) {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
278 return undef;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
279 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
280
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
281 # Collect assigned atom types...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
282 for $Atom (@{$This->{Atoms}}) {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
283 $AtomID = $Atom->GetID();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
284
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
285 $AtomType = $EStateAtomTypes->GetAtomType($Atom);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
286 $This->{EStateAtomTypes}{$AtomID} = $AtomType;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
287 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
288 return $This;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
289 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
290
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
291 # Calculate E-state indicies by summing up E-state values for specific
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
292 # E-state atom types...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
293 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
294 sub _CalculateEStateIndicies {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
295 my($This) = @_;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
296 my($Atom, $AtomID, $AtomType, $EStateValue);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
297
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
298 # Calculate E-state values to generate E-state indicies...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
299 if (!$This->_CalculateEStateValuesDescriptors()) {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
300 return undef;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
301 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
302
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
303 # Calculate E-state indicies...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
304 for $Atom (@{$This->{Atoms}}) {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
305 $AtomID = $Atom->GetID();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
306
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
307 $AtomType = $This->{EStateAtomTypes}{$AtomID};
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
308 $EStateValue = $This->{EStateValues}{$AtomID};
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
309
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
310 if (!exists $This->{EStateIndicies}{$AtomType}) {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
311 $This->{EStateIndicies}{$AtomType} = 0;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
312 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
313
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
314 $This->{EStateIndicies}{$AtomType} += $EStateValue;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
315 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
316 return $This;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
317 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
318
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
319 # Calculate E-state values for E-state indicies...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
320 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
321 sub _CalculateEStateValuesDescriptors {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
322 my($This) = @_;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
323 my($EStateValuesDescriptors, $Atom, $AtomID, $EStateValue);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
324
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
325 %{$This->{EStateValues}} = ();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
326
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
327 # Calculate and assign E-state values...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
328 $EStateValuesDescriptors = new AtomicDescriptors::EStateValuesDescriptors('Molecule' => $This->{Molecule});
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
329 $EStateValuesDescriptors->GenerateDescriptors();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
330
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
331 # Make sure E-state values calculation is successful...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
332 if (!$EStateValuesDescriptors->IsDescriptorsGenerationSuccessful()) {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
333 return undef;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
334 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
335
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
336 # Collect assigned E-state values...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
337 for $Atom (@{$This->{Atoms}}) {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
338 $AtomID = $Atom->GetID();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
339 $EStateValue = $EStateValuesDescriptors->GetDescriptorValue($Atom);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
340 $This->{EStateValues}{$AtomID} = $EStateValue;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
341 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
342 return $This;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
343 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
344
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
345 # Set final final fingerpritns for E-state indicies...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
346 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
347 sub _SetFinalFingerprints {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
348 my($This) = @_;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
349 my($AtomType, $ValuesPrecision, $EStateAtomTypesDataRef, @Values, @IDs);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
350
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
351 # Mark successful generation of fingerprints...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
352 $This->{FingerprintsGenerated} = 1;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
353
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
354 @Values = ();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
355 @IDs = ();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
356
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
357 $ValuesPrecision = $This->{ValuesPrecision};
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
358
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
359 if ($This->{EStateAtomTypesSetToUse} =~ /^FixedSize$/i) {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
360 # Use fixed size E-state atom types set for non-hydrogen atoms...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
361 for $AtomType (@{AtomTypes::EStateAtomTypes::GetAllPossibleEStateNonHydrogenAtomTypes()}) {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
362 push @IDs, "S${AtomType}";
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
363 push @Values, exists($This->{EStateIndicies}{$AtomType}) ? MathUtil::round($This->{EStateIndicies}{$AtomType}, $ValuesPrecision) : 0;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
364 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
365 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
366 else {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
367 for $AtomType (sort keys %{$This->{EStateIndicies}}) {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
368 push @IDs, "S${AtomType}";
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
369 push @Values, MathUtil::round($This->{EStateIndicies}{$AtomType}, $ValuesPrecision);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
370 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
371 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
372
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
373 # Add IDs and values to fingerprint vector...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
374 if (@IDs) {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
375 $This->{FingerprintsVector}->AddValueIDs(\@IDs);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
376 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
377 $This->{FingerprintsVector}->AddValues(\@Values);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
378
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
379 return $This;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
380 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
381
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
382 # Cache appropriate molecule data...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
383 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
384 sub _SetupMoleculeDataCache {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
385 my($This) = @_;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
386
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
387 # Get all non-hydrogen atoms...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
388 my($NegateAtomCheckMethod);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
389 $NegateAtomCheckMethod = 1;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
390 @{$This->{Atoms}} = $This->GetMolecule()->GetAtoms("IsHydrogen", $NegateAtomCheckMethod);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
391
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
392 return $This;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
393 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
394
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
395 # Clear cached molecule data...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
396 #
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
397 sub _ClearMoleculeDataCache {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
398 my($This) = @_;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
399
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
400 @{$This->{Atoms}} = ();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
401
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
402 return $This;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
403 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
404
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
405 # Return a string containg data for EStateIndiciesFingerprints object...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
406 sub StringifyEStateIndiciesFingerprints {
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
407 my($This) = @_;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
408 my($EStateIndiciesFingerprintsString);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
409
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
410 # Type of Keys...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
411 $EStateIndiciesFingerprintsString = "Type: $This->{Type}; EStateAtomTypesSetToUse: $This->{EStateAtomTypesSetToUse}";
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
412
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
413 # Fingerprint vector...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
414 $EStateIndiciesFingerprintsString .= "; FingerprintsVector: < $This->{FingerprintsVector} >";
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
415
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
416 return $EStateIndiciesFingerprintsString;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
417 }
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
418
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
419 1;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
420
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
421 __END__
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
422
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
423 =head1 NAME
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
424
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
425 EStateIndiciesFingerprints
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
426
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
427 =head1 SYNOPSIS
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
428
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
429 use Fingerprints::EStateIndiciesFingerprints;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
430
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
431 use Fingerprints::EStateIndiciesFingerprints qw(:all);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
432
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
433 =head1 DESCRIPTION
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
434
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
435 B<EStateIndiciesFingerprints> [ Ref 75-78 ] class provides the following methods:
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
436
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
437 new, GenerateFingerprints, GetDescription, SetEStateAtomTypesSetToUse,
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
438 SetValuesPrecision, StringifyEStateIndiciesFingerprints
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
439
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
440 B<EStateIndiciesFingerprints> is derived from B<Fingerprints> class which in turn
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
441 is derived from B<ObjectProperty> base class that provides methods not explicitly defined
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
442 in B<AtomNeighborhoodsFingerprints>, B<Fingerprints> or B<ObjectProperty> classes using Perl's
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
443 AUTOLOAD functionality. These methods are generated on-the-fly for a specified object property:
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
444
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
445 Set<PropertyName>(<PropertyValue>);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
446 $PropertyValue = Get<PropertyName>();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
447 Delete<PropertyName>();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
448
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
449 E-state atom types are assigned to all non-hydrogen atoms in a molecule using module
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
450 AtomTypes::EStateAtomTypes.pm and E-state values are calculated using module
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
451 AtomicDescriptors::EStateValues.pm. Using E-state atom types and E-state values,
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
452 B<EStateIndiciesFingerprints> constituting sum of E-state values for E-sate atom types
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
453 are generated.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
454
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
455 Two types of E-state atom types set size are allowed:
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
456
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
457 ArbitrarySize - Corresponds to only E-state atom types detected
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
458 in molecule
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
459 FixedSize - Corresponds to fixed number of E-state atom types previously
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
460 defined
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
461
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
462 Module AtomTypes::EStateAtomTypes.pm, used to assign E-state atom types to
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
463 non-hydrogen atoms in the molecule, is able to assign atom types to any valid
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
464 atom group. However, for I<FixedSize> value of B<EStateAtomTypesSetToUse>, only a
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
465 fixed set of E-state atom types corresponding to specific atom groups [ Appendix III in
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
466 Ref 77 ] are used for fingerprints.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
467
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
468 The fixed size E-state atom type set size used during generation of fingerprints contains
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
469 87 E-state non-hydrogen atom types in EStateAtomTypes.csv data file distributed with
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
470 MayaChemTools.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
471
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
472 Combination of Type and EStateAtomTypesSetToUse allow generation of 2 different types of
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
473 E-state indicies fingerprints:
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
474
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
475 Type EStateAtomTypesSetToUse
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
476
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
477 EStateIndicies ArbitrarySize [ default fingerprints ]
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
478 EStateIndicies FixedSize
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
479
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
480 The current release of MayaChemTools generates the following types of E-state
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
481 fingerprints vector strings:
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
482
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
483 FingerprintsVector;EStateIndicies:ArbitrarySize;11;NumericalValues;IDs
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
484 AndValuesString;SaaCH SaasC SaasN SdO SdssC SsCH3 SsF SsOH SssCH2 SssN
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
485 H SsssCH;24.778 4.387 1.993 25.023 -1.435 3.975 14.006 29.759 -0.073 3
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
486 .024 -2.270
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
487
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
488 FingerprintsVector;EStateIndicies:FixedSize;87;OrderedNumericalValues;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
489 ValuesString;0 0 0 0 0 0 0 3.975 0 -0.073 0 0 24.778 -2.270 0 0 -1.435
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
490 4.387 0 0 0 0 0 0 3.024 0 0 0 0 0 0 0 1.993 0 29.759 25.023 0 0 0 0 1
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
491 4.006 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
492 0 0 0 0 0 0 0 0 0 0 0 0 0 0
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
493
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
494 FingerprintsVector;EStateIndicies:FixedSize;87;OrderedNumericalValues;
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
495 IDsAndValuesString;SsLi SssBe SssssBem SsBH2 SssBH SsssB SssssBm SsCH3
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
496 SdCH2 SssCH2 StCH SdsCH SaaCH SsssCH SddC StsC SdssC SaasC SaaaC Sssss
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
497 C SsNH3p SsNH2 SssNH2p SdNH SssNH SaaNH StN SsssNHp SdsN SaaN SsssN Sd
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
498 0 0 0 0 0 0 0 3.975 0 -0.073 0 0 24.778 -2.270 0 0 -1.435 4.387 0 0 0
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
499 0 0 0 3.024 0 0 0 0 0 0 0 1.993 0 29.759 25.023 0 0 0 0 14.006 0 0 0 0
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
500 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0...
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
501
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
502 =head2 METHODS
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
503
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
504 =over 4
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
505
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
506 =item B<new>
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
507
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
508 $EStateIndiciesFingerprints = new EStateIndiciesFingerprints(%NamesAndValues);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
509
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
510 Using specified I<EStateIndiciesFingerprints> property names and values hash, B<new> method creates a new object
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
511 and returns a reference to newly created B<PathLengthFingerprints> object. By default, the
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
512 following properties are initialized:
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
513
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
514 Molecule = '';
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
515 Type = 'EStateIndicies'
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
516 EStateAtomTypesSetToUse = 'ArbitrarySize'
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
517 ValuesPrecision = 3
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
518
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
519 Examples:
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
520
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
521 $EStateIndiciesFingerprints = new AtomTypesFingerprints(
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
522 'Molecule' => $Molecule,
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
523 'EStateAtomTypesSetToUse' =>
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
524 'ArbitrarySize');
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
525
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
526 $EStateIndiciesFingerprints = new AtomTypesFingerprints(
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
527 'Molecule' => $Molecule,
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
528 'EStateAtomTypesSetToUse' =>
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
529 'FixedSize');
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
530
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
531 $EStateIndiciesFingerprints->GenerateFingerprints();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
532 print "$EStateIndiciesFingerprints\n";
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
533
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
534 =item B<GenerateFingerprints>
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
535
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
536 $EStateIndiciesFingerprints = $EStateIndiciesFingerprints->
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
537 GenerateEStateIndiciesFingerprints();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
538
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
539 Generates EState keys fingerprints and returns I<EStateIndiciesFingerprints>.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
540
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
541 =item B<GetDescription>
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
542
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
543 $Description = $EStateIndiciesFingerprints->GetDescription();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
544
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
545 Returns a string containing description of EState keys fingerprints.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
546
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
547 =item B<SetEStateAtomTypesSetToUse>
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
548
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
549 $EStateIndiciesFingerprints->SetEStateAtomTypesSetToUse($Value);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
550
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
551 Sets I<Value> of I<EStateAtomTypesSetToUse> and returns I<EStateIndiciesFingerprints>.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
552 Possible values: I<ArbitrarySize or FixedSize>. Default value: I<ArbitrarySize>.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
553
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
554 =item B<SetValuesPrecision>
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
555
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
556 $EStateIndiciesFingerprints->SetValuesPrecision($Precision);
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
557
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
558 Sets precesion of E-state values to use during generation of E-state indices fingerprints
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
559 and returns I<EStateIndiciesFingerprints>. Possible values: I<Positive integers>.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
560 Default value: I<3>.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
561
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
562 =item B<StringifyEStateIndiciesFingerprints>
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
563
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
564 $String = $EStateIndiciesFingerprints->StringifyEStateIndiciesFingerprints();
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
565
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
566 Returns a string containing information about I<EStateIndiciesFingerprints> object.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
567
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
568 =back
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
569
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
570 =head1 AUTHOR
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
571
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
572 Manish Sud <msud@san.rr.com>
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
573
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
574 =head1 SEE ALSO
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
575
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
576 Fingerprints.pm, FingerprintsStringUtil.pm, AtomNeighborhoodsFingerprints.pm, AtomTypesFingerprints.pm,
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
577 ExtendedConnectivityFingerprints.pm, MACCSKeys.pm, PathLengthFingerprints.pm,
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
578 TopologicalAtomPairsFingerprints.pm, TopologicalAtomTripletsFingerprints.pm,
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
579 TopologicalAtomTorsionsFingerprints.pm, TopologicalPharmacophoreAtomPairsFingerprints.pm,
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
580 TopologicalPharmacophoreAtomTripletsFingerprints.pm
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
581
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
582 =head1 COPYRIGHT
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
583
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
584 Copyright (C) 2015 Manish Sud. All rights reserved.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
585
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
586 This file is part of MayaChemTools.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
587
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
588 MayaChemTools is free software; you can redistribute it and/or modify it under
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
589 the terms of the GNU Lesser General Public License as published by the Free
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
590 Software Foundation; either version 3 of the License, or (at your option)
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
591 any later version.
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
592
a4a2ad5a214e Uploaded
deepakjadmin
parents:
diff changeset
593 =cut