................ SHORT DOC .............................................
RADI: Computation of the Radius and Diameter of a spatial set
Reference:
M. Petitjean, Applications of the Radius-Diameter Diagram to
the Classification of Topological and Geometrical Shapes of
Chemical Compounds, J.Chem.Inf.Comput.Sci. 1992,32[4],331-337
Author email: petitjean@itodys.jussieu.fr
RADI reads the cartesian coordinates of the molecule. It extracts
the convex hull of the atomic positions, then computes the diameter
and the radius of the hull. Three shape coefficients are computed,
all related to the convex hull of the set.
Input data and parameters:
-------------------------
INPUT FORMAT:
CAS : Reserved for internal purposes
HIN : Hyperchem-type files
MDL : Cambridge Crystallographic Model files
ML2 : SYBYL Mol2 files
PDB : Protein Data Bank or Nucleic Acid Data Bank files
(only HEADER, ATOM, ENDMDL and END records are recognized)
BIO : Biosym (MSI) files
ISU : Reserved for internal purposes
INPUT MOLEC FILE NAME: name of the input file containing the molecule
Output results:
--------------
3D CONVEX HULL: the K atoms of the hull, i.e. the K vertices of the
smallest convex polyhedron containing the set of the
N atomic positions
3D MOL. MEAN POINT: arithmetic mean of the N atomic positions
3D HULL MEAN POINT: arithmetic mean of the K atoms of the hull
The SURFACE and VOLUME of the full 3D convex hull (no cavity is considered)
DIAMETER: the distance D of the longest atom-pair of the molecule
RADIUS: the radius R of the smallest sphere containing the N atoms
MIN SPHERE CENTER: the center of this minimal sphere
GEOMETRICAL CONVEX-SHAPE COEFFICIENT : the ratio (D-R)/R
SURFACIC CONVEX-SHAPE COEFFICIENT :
The ratio: surface of the hull, divided by the surface of the minimal sphere,
followed by the square root of this ratio
VOLUMIC CONVEX-SHAPE COEFFICIENT :
The ratio: volume of the hull, divided by the volume of the minimal sphere,
followed by the cubic root of this ratio
Remarks:
-------
The number of atoms is currently limited to 15000 for each molecule.
The source has to be recompiled to read larger molecules.
Perfectly planar molecules are not handled. When there are at least
four atoms in the molecule, it is suggested to slightly modify one
atomic coordinate to get a spatial set.
The computing time is worst-case bounded by O(K**5), K being the number
of vertices of the hull, and K is at worst equal to N.
The worst situation is not expected in chemistry, even for proteins.
................END SHORT DOC .....................................
|