RelativeFEPSetup cannot handle list of smiles as input ligands #1155
Description
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)