Skip to content

Commit 17da0c2

Browse files
committed
Use RMG smiles parser for fragments
1 parent c529dae commit 17da0c2

File tree

1 file changed

+10
-1
lines changed

1 file changed

+10
-1
lines changed

src/Parse.jl

Lines changed: 10 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -161,7 +161,16 @@ function getatomdictfromrmg(mol)
161161
molecularweight = mol.get_molecular_weight()
162162
return atmD,nbonds,molecularweight
163163
end
164-
getatomdictsmiles(smiles) = getatomdictfromrdkit(Chem.AddHs(Chem.MolFromSmiles(smiles)))
164+
function getatomdictsmiles(smiles)
165+
if occursin(r"R", smiles) || occursin(r"L", smiles)
166+
mol = fragment.Fragment().from_smiles_like_string(smiles)
167+
mol.assign_representative_molecule()
168+
getatomdictfromrmg(mol.mol_repr)
169+
else
170+
getatomdictfromrdkit(Chem.AddHs(Chem.MolFromSmiles(smiles)))
171+
end
172+
end
173+
165174
export getatomdictsmiles
166175
getatomdictinchi(inchi) = getatomdictfromrdkit(Chem.AddHs(Chem.MolFromInchi(inchi)))
167176
export getatomdictinchi

0 commit comments

Comments
 (0)