Untitled

from rdkit import Chem

from rdkit.Chem import Draw


# Define the reactants and products as SMILES strings
# Reaction (i)
reactant1_i = Chem.MolFromSmiles("CC(=O)CC(=O)C(=O)OCC")
product_i = Chem.MolFromSmiles("CC1CCC(O)C1=O")

# Reaction (ii)
reactant1_ii = Chem.MolFromSmiles("CC(C)C(=O)O")
product_ii = Chem.MolFromSmiles("CC(C(=O)O)C(=O)O")

# Generate reaction depictions
images = Draw.MolsToGridImage(
    [reactant1_i, product_i, reactant1_ii, product_ii],
    legends=["Reaction (i) Reactant", "Reaction (i) Product",
             "Reaction (ii) Reactant", "Reaction (ii) Product"],
    molsPerRow=2
)

imag
es