Avaliable on regioqm.org.

Predict Regioselectivity of electrophilic aromatic substitution reactions in heteroaromatic systems. RegioSQM finds the aromatic CH group with the highest proton affinity estimated by the PM3/COSMO method using the MOPAC program. More information is available at the RegioSQM paper.

- MOPAC (http://openmopac.net/)
- RDKit (http://www.rdkit.org/docs/Install.html)
- obabel (https://openbabel.org/docs/dev/Installation/install.html)

RegioSQM depends on MOPAC for quantum calculations, OpenBabel for some formation convertion and RDkit in the python enviroment for everything else.

The workflow is as following

cd example

# generate conformations from SMILES
python ../regiosqm/regiosqm.py -g example.csv > example.csv

# Run all .mop files with mopac
# or submit them to a cluster
ls *mop | parallel -j4 "mopac {}"

# use the generated csv file to analyse all the 
python ../regiosqm/regiosqm.py -a example.csv example.csv > example_results.csv

The results are now parseable from the results file, or viewable by 2D structures with regioselective indicators (in svg format).