Calculator Plugins

Contents

Marvin contains a framework for integrating chemical computations into the drawing/viewing application environment. These tools - called plugins - are loaded dynamically upon request. ChemAxon provides various tools for calculating charge, pKa, logP, etc. The available calculator plugins are located in the Tools menu. Batch processing is available using Calculator (see the list of calculations accessible from Calculator).

Some plugins (Charge, Polarizability, Polar Surface Area, Molecular Surface Area and Hydrogen Bond Donor-Acceptor) optionally perform their computation on the major microspecies of the input molecule taken at a specified pH. This option together with the corresponding pH can also be set in the parameter panel.

Elemental Analysis Plugin

Basic molecular values related to the elemental composition of the molecule are displayed by the Elemental Analysis plugin.

By default, molecules are handled separately if more than one molecules are drawn in the sketcher. However, sometimes a single molecule consists of more fragments (e.g. salt molecules), in which case the fragments should be treated as one molecule. This behavior can be reached by switching off the "Single fragment mode" option in the Elemental Analysis Options panel.

Isotopes of Hdrogen are displayed in formulas as 2H and 3H by default, but using D and T symbols instead is possible by switching on the "Use D/T symbols for Deuterium/Tritium" option in the Elemental Analysis Options panel.

IUPAC Naming Plugin

IUPAC naming plugin generates preferred IUPAC names, that conform (when possible) to the IUPAC Provisional Recommendations for the Nomenclature of Organic Chemistry published in 2004. Traditional names of the molecules can also be generated optionally.

You can generate either the "Traditional Name" or the "Preferred IUPAC Name" of the molecules; you can change between these otion in the IUPAC Name Options panel. By default, the "Traditional Name" option is set, but even in this case, if the traditional name cannot be generated, the preferred IUPAC Name will be generated automatically.

In case of the IUPAC name generator plugin, by default, the molecules are handled together (as a mixture), and the generated names are separated by a semicolon. To handle the fragments and get their names separately, the "Single fragment mode" option in the IUPAC Name Options panel has to be switched on.

Learn more about IUPAC name generation in the IUPAC name generator documentation.

Protonation

pKa Plugin

Most molecules contain some specific functional groups likely to loose or gain proton under specific circumstances. Each ionization equilibrium between the protonated or deprotonated forms of the molecule can be described with a constant value called pKa. The pKa plugin calculates the pKa values of all proton gaining or loosing atoms on the basis of the partial charge distribution.
Learn more about how the plugin calculates pKa.

The chart shows the microspecies distribution curves by pH. The microspecies images are shown in the legend. When clicking on an image, the corresponding microspecies molecule is displayed in the upper-left viewer. The viewer can be detached from the chart panel by double clicking in it, or else by selecting Open Viewer from the View menu. The original molecule with the pKa values is shown when clicking on the chart outside of the legend image areas, or else when selecting pKa Values from the View menu.

Major Microspecies Plugin

Determines the major protonation form at a specified pH.

The pH can be set in the Major Microspecies Options panel, the default pH is 7.4.

Isoelectric Point Plugin

Gross charge of an ionizable molecule is zero at certain pH. This pH is called as isoelectric point. Isoelectric point plugin calculates gross charge distribution of a molecule as function of pH.

Partitioning

logP Plugin

The logP plugin calculates the octanol/water partition coefficient, which is used in QSAR analysis and rational drug design as a measure of molecular hydrophobicity. The calculation method is based on the publication of Viswanadhan et al. (see Ref.1.) The logP value of a molecule is composed of the increment values of its atoms. The algorithm described in the paper was modified at several points. Many atomic types were redefined to accommodate electron delocalization. Contributions of ionic forms were added. The logP value of zwitterions are calculated from the logD value at the isoelectric point. The effect of hydrogen bonds on logP is considered if there is a chance to form a six membered ring between suitable donor and acceptor atoms. New atom types were introduced especially for sulfur, carbon, nitrogen, and metal atoms.
Learn more about how the plugin calculates logP.

The result window shows the logP increments for each atom. The numbers in brackets refer to the logP increment sums of implicit H atoms, and displayed only if the "Increment of Hs" option is switched on in the logP Options panel.

logD Plugin

Compounds having ionizable groups exist in solution as a mixture of different ionic forms. The ionization of those groups, thus the ratio of the ionic forms depends on the pH. Since logP describes the hydrophobicity of one form only, the apparent logP value can be different. The octanol-water distribution coefficient, logD represents the compounds at any pH value (see Ref.2.).
Learn more about how the plugin calculates logD.

The chart shows the logD(pH) curves for each molecule drawn in the sketcher. The molecule images are shown in the legend. When clicking on an image, the corresponding molecule is displayed in the upper-left viewer. The viewer can be detached from the chart panel by double clicking in it, or else by selecting Open Viewer from the View menu. The reference logD values originally shown can be restored by either clicking on the chart outside of the legend image areas, or else by selecting logD at reference pH-s from the View menu.

Charge

Charge Plugin

The partial charge distribution determines many physico-chemical properties of a molecule, such as ionization constants, reactivity and pharmacophore pattern. Use Charge plugin to compute the partial charge value of each atom. Total charge is calculated from sigma and pi charge components, and any of these three charge values can be displayed. You can change between these by setting the "Type" option in the Charge Options panel; by default, the total charge is displayed. You can also
Learn more about how the plugin calculates the partial charge.

The numbers in brackets refer to the charge sums of implicit H atoms, and displayed only if the "Increment of Hs" option is switched on in the Charge Options panel.

Polarizability Plugin

The electric field generated by partial charges of a molecule spread through intermolecular cavities and the solvent that the molecule is solved within. The induced partial charge (induced dipole) has a tendency to diminish the external electric field. This phenomenon is called as polarizability. The more stable each ionized site is the more its vicinity is polarizable. This is why atomic polarizability is an important factor in the determination of pKa and why it is considered in our pKa calculation plugin. Atomic polarizability is altered by partial charges of atoms. Our calculation is based on Ref.3., and takes into account the effect of partial charge upon atomic polarizability.

Orbital Electronegativity Plugin

Partial charge distribution of the molecule is governed by the orbital electronegativity of the atoms contained in the molecule.
Learn more about how the plugin calculates orbital electronegativity.

Isomers

Tautomerization Plugin

Tautomers are structural isomers which are in dynamic equilibrium due to the migration of protons. All tautomers of a compound can be determined with the help of the Tautomerization plugin. For example the next structures (on the left) are the calculated tautomers of pyrimidine-2,4-diol. The tautomers shown on the result window below were generated taking into consideration the pH effect, and the dominant tautomer distibution is shown for each tautomer. The plugin is also able to generate the canonical tautomeric form (the generated canonical tautomer of pyrimidine-2,4-diol can be seen on the right). The choices between the generation of 'basic' tautomers, pH-effected tautomers and only the canonical tautomer can be set by the different options of the Tautomers Options panel.

Resonance Plugin

The Resonance plugin generates all resonance forms of a molecule. The major contributors of the resonance structures can be calculated separately. For example the next two structures on the left are the major resonance contributors of diazomethane, while the structure on the right is the canonical form.

Stereoisomer Plugin

The Stereoisomer plugin produces all possible stereoisomers of a given compound. The plugin handles both tetrahedral and double bond stereo centers.

In case the "Filter invalid 3D structures" option is switched on in the Stereoisomers Options panel, the stereoisomers can also be displayed in 3D.

Conformation

Conformer Plugin

Conformer plugin can generate selected number of conformers or the lowest energy conformer of a molecule. For conformer calculation Dreiding force field is used. The 3D structure (conformation) strongly affects the properties and the reactions of the molecules.

Molecular Dynamics Plugin

The molecular dynamics plugin calculates the configurations of the system by integrating Newton's laws of motion.

The calculation and the display options can be set on the Molecular Dynamics Options panel:

Geometry

Topology Analysis Plugin

The Topology Analysis plugin provides characteristic values related to the topological structure of a molecule.

Simple

Ring

Path based

Distance based

Other

Geometry Plugin

The Geometry plugin provides characteristic values related to the geometrical structure of a molecule. It can calculate interatomic distances, bond angles, dihedral angles, steric hindrance and Dreiding energy. The calculation can predict and use the lowest energy conformer of the input structure.

Calculations

Polar Surface Area Plugin (2D)

Polar surface area (PSA) is formed by polar atoms of a molecule. It is a descriptor that shows good correlation with passive molecular transport through membranes, and so allows estimation of transport properties of drugs. Estimation of topoligical polar surface area (TPSA) is based on the method given in Ref.4.. The method provides results which are practically identical with the 3D PSA, while calculation time of TPSA is approximately 100-times faster. This method is more suitable for fast bioavailability screening of large virtual libraries. The TPSA value can be calculated both for the neutral form and the major microspecies.

Molecular Surface Area Plugin (3D)

There are two types of available molecular surface area calculations: Van der Waals and solvent accessible. Calculation method is based on the publication of Ferrara et al. (Ref.6.).

Other

Hydrogen Bond Donor-Acceptor Plugin

Hydrogen Bond Donor-Acceptor calculates atomic hydrogen bond donor and acceptor inclination. Atomic data and overall hydrogen bond donor and acceptor multiplicity are displayed for the input molecule (or its physiological microspecies at a given pH). The weighted average hydrogen bond donor and acceptor multiplicities taken over the microspecies th proportions of their occurrences are computed for different pH-s and displayed in a chart.

Huckel Analysis Plugin

Localization energies L(+) and L(-) for electrophilic and nucleophilic attack at an aromatic center are calculated by the Hückel method. The smaller L(+) or L(-) means more reactive atomic location. Order of atoms in E(+) or in Nu(-) attack are adjusted according to their localization energies. The total pi energy, the pi electron density and the total electron density are also calculated by the Hückel method.

Theoretical background is given in Neil S. Isaacs5.

Refractivity Plugin

Our calculation is based on the atomic method proposed by Viswanadhan et al.1. Molar refractivity is strongly related to the volume of the molecules and to London dispersive forces that has important effect in drug-receptor interaction.

The numbers in brackets refer to the refractivity sums of the implicit H atoms.

Markush Enumeration Plugin

A Markush structure is a description of a compound class by generic notations, primarily used in patent claims and the description of combinatorial libraries. The library of a Markush structure is the total set of specific molecules that are described by the Markush structure.

The Markush enumeration plugin can be used to generate a whole or a subset of the library of a generic Markush structure. It is also capable of calculating the total number of specific structures present in a Markush library. The plugin is accessible from the Marvin GUI (Tools->Markush Enumeration), through the cxcalc command-line program (See this link for the detailed usage of the plugin in command line.), via API and in the Chemical Terms functions in JChem.

Markush features

Currently, the Enumeration plugin supports the following features that describe mostly Markush structures in combinatorial libraries:

Name Description Example Example Markush library member
R-groups R-groups (also referred to as "substituent variation") are the most widely known Markush generic features. The variable part of the structure is denoted by an R-atom (eg. R1), and the definitions are given separately. In each definition the connection points must be defined to show where the bonds of the R-atom are linked. R-atoms can appear in both rings and chains and can have up to two attachments points. The same R-atom can appear multiple times, and the different occurrences are handled as different cases. (So they can be substituted with different definitions.) R-group nesting in R-group definitions is allowed to any depth, but without recursion. (An R-group definition cannot use the R-atom it is defining, not even through the use of other embedding R-atom(s).) In Marvin and the plugin, R-groups up to number R32767 can be used.
Atom lists Atom lists are another example of substituent variation. They define lists of atom types at a given position. There is no restriction for the length of the list and for bond count of atom lists.
Bond lists The following bond lists (generic bond types) are supported by the plugin: single or double, any(single, double or triple), single or aromatic, double or aromatic.
Link nodes Link nodes are atoms that may repeat between two of their designated bonds (called outer bonds, denoted by brackets). All other substituents (if exist) repeat together with the atom. In the results, the new bonds between the repeating atoms will have the bond type of the lower order outer bond.

Functionality of the plugin

The plugin allows the following functionality. Examples are given using Marvin GUI.
Sequential enumeration
Enumerates members of the Markush library in a sequential manner (by substituting the first definition of the first variable, etc). The results are specific structures. The plugin user interface allows the enumeration of all library members, or a specified number.

Random enumeration
This mode generates a random subset of the Markush library to give a quick sampling. It is especially helpful for huge libraries, where full enumeration is impossible. In random mode variable parts are chosen randomly, and the substitution probability of each definition is proportionate with the fragment library size that the given definition generates. This ensures the generation of a representative random subset over the Markush library space.

Calculate library size
The size of the Markush library can be calculated by arbitrary precision. On the user interface, the exact value is displayed until 20 digits, above that only the magnitude is shown (for example, 10^28). The calculated number is the size of the whole library, and does not consider the valence check filter. (See below.)

Selected part enumeration
If part of the Markush structure is selected, only the generic features in the selected part are considered for enumeration/calculation. This allows focusing on a particular area of the Markush structure. Enumeration of selected parts only may result in generating (more specific) Markush structures.

Valence filter
If the Markush structure is not properly (or too generally) formulated, it is possible that it describes structures with valence errors. In this case, the valence filter setting is useful to filter out the offending result structures. The default value is on.

References

  1. Viswanadhan, V. N.; Ghose, A. K.; Revankar, G. R.; Robins, R. K., J. Chem. Inf. Comput. Sci., 1989, 29, 163-172
  2. Csizmadia, F; Tsantili-Kakoulidou, A.; Pander, I.; Darvas, F., J. Pharm. Sci., 1997, 86, 865-871
  3. Miller, K. J.; Savchik, J. A., J. Am. Chem. Soc., 1979, 101, 7206-7213
  4. Ertl, P.; Rohde, B.; Selzer, P., J. Med. Chem., 2000, 43, 3714-3717
  5. Neil S. Isaacs, Physical Organic Chemistry, John Wiley & Sons, Inc., New York, 1987.
  6. Ferrara, P; Apostolakis J; Caflisch A., Proteins 2002, 46, 24-33

Return to Top

Back to Marvin User's Guide

Copyright © 1998-2008 ChemAxon Ltd.

http://www.chemaxon.com/marvin