This manual gives you a walk-through on how to use the Stereo Analysis module:
Stereochemistry involves the study of the relative spatial arrangement of atoms and their transformations. Stereochemistry defines different isomerisms that describe the molecule's behaviour in the space. These isomer groups cover the entire space of organic molecules. We do not give further details on stereochemistry here, instead we refer to Marvin's stereochemistry summary.
It can be important to determine a molecule's stereochemical behaviour. This can be achieved by calculating stereochemical descriptors for the molecule.
Chemaxon's Stereo Analysis module is able to calculate stereochemical descriptors for a molecule, giving an analysis of the stereochemical properties present in the input.
The stereo analysis functionality is built in different Chemaxon products.
Stereo Analysis module is integrated into the cxcalc command line tool. The general command syntax is
cxcalc [general options] [input files/strings] stereoanalysis [stereoanalysis options] [input files/strings]
where the stereoanalysis options are the following:
stereoanalysis options: -h --help this help message -T --type stereo descriptor type [tetrahedral | cistrans | axial | atrop] (default: not set)
Here are some examples of how to use the stereo analysis module via cxcalc :
Calculating tetrahedral stereo descriptors for the following molecule:
cxcalc stereoanalysis "[H][C@](cl)(Br)I"
The output shows the recognized tetrahedral stereo descriptors:
TETRAHEDRAL [1] - [4, 3, 2, 0] : S
Calculating cis-trans stereo descriptors:
cxcalc stereoanalysis "Cl/C=C\C=C\Cl"
The output shows the recognized cis-trans stereo descriptors:
CISTRANS [1, 2] - [0, 3] : Z CISTRANS [3, 4] - [2, 5] : E
Calculating axial stereo descriptors for a molecule:
cxcalc stereoanalysis CC=C=CC
The output shows the recognized axial stereo descriptors:
AXIAL [1, 3] - [0, 4] : WIGGLY
Calculating atrop stereo descriptor for a molecule:
cxcalc stereoanalysis "CCC1=C(C(C)=CC=C1)C1=CC=CC=C1C"
The output shows the recognized atrop stereo descriptors:
ATROP [3, 9] - [4, 2, 14, 10] : UNKNOWN
Calculating all stereo descriptors for a molecule with different stereo features:
cxcalc stereoanalysis "ClC(Br)C1CCCC(C1)/C=C/Cl"
The output shows the recognized stereo descriptors:
TETRAHEDRAL [7] - [9, 8, 6] : UNKNOWN TETRAHEDRAL [3] - [1, 8, 4] : UNKNOWN TETRAHEDRAL [1] - [2, 0, 3] : UNKNOWN CISTRANS [9, 10] - [7, 11] : E
Listing just the cis-trans stereo descriptors for the previous molecule:
cxcalc stereoanalysis -T cistrans "ClC(Br)C1CCCC(C1)/C=C/Cl"
The output is:
CISTRANS [9, 10] - [7, 11] : E
The stereo analysis module is also integrated into Chemaxon's Chemical Terms language. The stereo descriptors can be calculated by two different functions in Chemical Terms:
stereoanalysis() : calculates and returns stereo descriptors of all types for the input molecule.
stereoanalysis(T) : calculates and returns stereo centers of type T . T can be any stereo function defined in the API.
The following examples show how the functions above can be used with the Chemical Terms evaluator command line tool:
Calculating all stereo descriptors for geldanamycin:
evaluate -e "stereoanalysis()" geldanamycin
The output is:
TETRAHEDRAL [26] - [27, 23, 31] : S CISTRANS [18, 20] - [16, 19, 21] : E TETRAHEDRAL [36] - [37, 2, 34] : R TETRAHEDRAL [4] - [3, 6, 5] : R CISTRANS [21, 22] - [20, 23] : Z TETRAHEDRAL [34] - [36, 33, 35] : S TETRAHEDRAL [23] - [24, 26, 22] : S TETRAHEDRAL [2] - [1, 36, 3] : S CISTRANS [31, 33] - [26, 32, 34] : E
Calculating only tetrahedral stereo centers for the same molecule:
evaluate -e "stereoanalysis('tetrahedralStereoCenters')" geldanamycin
The output is:
2;4;23;26;34;36
Calculating the number of unknown, resolved, tetrahedral stereocenters for the following molecule with the count() function:
evaluate -e "count(stereoanalysis('unknownResolvedStereoCenters'))" "C[C@@H]1OP(O)(=O)O[C@@H](c)C1C |o1:7,r|"
Finding only the carbon tetrahedral stereocenters can be done using the filter() function for the same molecule:
evaluate -e "filter(stereoanalysis('tetrahedralStereoCenters'), \"match('C')\")" "C[C@@H]1OP(O)(=O)O[C@@H](c)C1C |o1:7,r|"