choderalab / perses

Experiments with expanded ensembles to explore chemical space

Home Page:http://perses.readthedocs.io

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

`RelativeFEPSetup` cannot handle list of smiles as input ligands

ijpulidos opened this issue · comments

According to the documentation, RelativeFEPSetup should be able to handle a list of SMILES strings as input ligands but it makes the interpreter to freeze giving the following warnings:

Warning: Unable to set unknown read format (UNDEFINED) in input stream, using existing format (Canonical stereo SMILES)
Warning: Unable to open "Oc1ccccc1" for reading molecules.

Warning: Error reading molecule "" in Canonical stereo SMILES format.

And example code to reproduce the issue is as follows (Warning, this kinda kills the python interpreter, use at your own risk):

from perses.app.relative_setup import RelativeFEPSetup
from openmm import unit
phenol_smiles = "Oc1ccccc1"
paracetamol_smiles = "CC(=O)Nc1ccc(O)cc1"
input_smiles = [phenol_smiles, paracetamol_smiles]
fe_setup = RelativeFEPSetup(
                ligand_input=input_smiles,
                old_ligand_index=0,
                new_ligand_index=1,
                forcefield_files=["amber14/tip3p.xml"],
                small_molecule_forcefield="gaff-2.11",
                phases=["solvent"],
                solvent_padding=1.1 * unit.nanometers)