def METHOD_NAME():
    from openff.toolkit.topology import Molecule
    # Test something a bit more tricky than water
    molecules = [Molecule.from_smiles("CC(=O)NC1=CC=C(C=C1)O")]
    trajectory, _ = packmol.pack_box(
        molecules, [1], box_size=([20] * 3) * unit.angstrom
    )
    assert trajectory is not None
    assert trajectory.n_chains == 1
    assert trajectory.n_residues == 1
    assert trajectory.n_atoms == 20
    assert trajectory.topology.n_bonds == 20