  We develop a mathematical model for the energy landscape of polyhedral
supramolecular cages recently synthesized by self-assembly [Sun et al., Science
2010]. Our model includes two essential features of the experiment: (i)
geometry of the organic ligands and metallic ions; and (ii) combinatorics. The
molecular geometry is used to introduce an energy that favors square-planar
vertices (modeling $\mathrm{Pd}^{2+}$ ions) and bent edges with one of two
preferred opening angles (modeling boomerang-shaped ligands of two types). The
combinatorics of the model involve $2$-colorings of edges of polyhedra with
$4$-valent vertices. The set of such $2$-colorings, quotiented by the
octahedral symmetry group, has a natural graph structure, and is called the
combinatorial configuration space. The energy landscape of our model is the
energy of each state in the combinatorial configuration space.
  The challenge in the computation of the energy landscape is a combinatorial
explosion in the number of $2$-colorings of edges. We describe sampling methods
based on the symmetries of the configurations and connectivity of the
configuration graph. When the two preferred opening angles encompass the
geometrically ideal angle, the energy landscape exhibits a very low-energy
minimum for the most symmetric configuration at equal mixing of the two angles,
even when the average opening angle does not match the ideal angle.
