Reaction stereo

    Chirality

    Reaction stereo centers can be defined in two ways:

    1. as an INVRET (inversion/retention) atom property in the MOL (or RXN, RDF) form

    2. by setting atom parity using up/down-wedge bonds

    In both cases the stereo center should be mapped with the same atom map on the reactant and the product sides.

    The difference is that while reactions defined by the INVRET atom property match both stereo configurations as well as the non-stereo form, the up/down wedge stereo definition will only match the stereo configuration specified in the reactant side of the reaction.

    While the first possibility is more effective, it is not available in SMILES.

    If reaction stereo information is given by the second option and there are at least two unmapped non-hydrogen ligands connected to the stereo center, then the INVRET (inversion/retention) atom property is first determined by preprocessing the reaction and then the reaction processing is done as in the first case. If the ligands of the stereo center are sufficiently mapped (that is, there is at most one unmapped non-hydrogen ligands) then there is no need of this preprocession and the reaction handling is more effective.

    The advantage of the second option is that it is also available in SMILES.

    If the RETENTION atom property is specified then the stereo configuration in the reactants is kept in the products, while in case of INVERSION the stereo configuration is inverted.

    Limitations:

    • by setting atom parity using up/down-wedge bonds (second option) with at least two unmapped non-hydrogen ligands only one stereo center can be specified per reaction

    • the stereo configuration can be processed only if at most one atom changes around the stereo center

      Examples:

    The following reaction defines a stereoselective reaction:

    images/download/attachments/5311879/image2019-8-13_9-55-10.png

    Reactor reacts the first compound according to the above equation, but not the second one, since the reactant is not matching:

    images/download/attachments/5311879/image2019-8-13_9-55-47.png
    images/download/attachments/5311879/image2019-8-13_9-56-12.png

    The inversion/retention flag can be used to define stereoselective reactions producing the appropriate chiral compounds from both enantiomers. The definition:

    images/download/attachments/5311879/image2019-8-13_9-56-55.png

    and the actual stereo reactions inverting the reaction center (please note, that the third reactant is racemic, thus the product is racemic as well; toluene- p -sulphonic acid is not shown in the example):

    images/download/attachments/5311879/image2019-8-13_9-57-54.png
    images/download/attachments/5311879/image2019-8-13_9-58-14.png
    images/download/attachments/5311879/image2019-8-13_9-58-46.png

    Double bond stereo

    To propagate the cis/trans stereo configuration in the reaction definition to the reaction product, the corresponding four atoms should be mapped. For example, the following reaction changes the double bond configuration from cis to trans . Reactor can process double bond isomerisation reactions when the corresponding double bond ligands (at least four atoms) are mapped (the small square on the double bond marks stereo specific matching):

    images/download/attachments/5311879/image2019-8-13_10-0-9.png

    The actual reaction equations. The trans-coumarin is transformed while the cis isomer is not:

    images/download/attachments/5311879/image2019-8-13_12-12-37.png
    images/download/thumbnails/5311879/image2019-8-13_12-13-9.png

    Instead of CIS or TRANS configuration on the product side, you can also specify an undetermined double bond stereo configuration, denoted by a wiggly bond:

    images/download/attachments/5311879/image2019-8-13_12-14-5.png

    In this case you will get the same cis-or-trans configuration in the actual product, denoted by a wiggly bond in the same way as in the reaction equation:

    images/download/attachments/5311879/image2019-8-13_12-14-47.png

    If you also set unspecified double bond stereo on the reactant side of the reaction equation:

    images/download/attachments/5311879/image2019-8-13_12-15-34.png

    then the cis-coumarin will also be transformed:

    images/download/attachments/5311879/image2019-8-13_12-16-14.png

    Double bond stereo configuration on atoms not taking part in the reaction is left unchanged; however, there may be cases when the double bond stereo type cannot be determined: if the double bond has stereo configuration (that is, different ligands attached to at least one of the end-atoms) but there is neither an unmapped, nor a fully mapped A1-A2=A3-A4 configuration along which the configuration could be determined. In such cases an error is reported in command-line interface if the --check parameter is specified, otherwise the double bond stereo type is unpredictable. To prevent this, map a full A1-A2=A3-A4 atom configuration for each double bond on the product side of the reaction that has at least one mapped ligand.

    Examples:

    changing bond type images/download/attachments/5311879/image2019-8-13_12-18-1.png
    orphan atoms images/download/attachments/5311879/image2019-8-13_12-18-29.png
    chirality images/download/attachments/5311879/image2019-8-13_12-19-1.png
    double bond stereo images/download/attachments/5311879/image2019-8-13_12-19-32.png