We present variational results for the ground state of the antiferromagnetic quantum Heisenberg model with frustrating next-nearest-neighbour interactions. The trial wave functions employed are of resonating-valence-bond type, elaborated to account for various geometric motifs of adjacent bond pairs. The calculation is specialized to a square-lattice cluster consisting of just sixteen sites, large enough that the system can accommodate nontrivial singlet dimer correlations but small enough that exhaustive enumeration of states in the total spin zero sector is still feasible. A symbolic computation approach allows us to generate an algebraic expression for the expectation value of any observable and hence to carry out the energy optimization exactly. While we have no measurements that could unambiguously identify a spin liquid state in the controversial region at intermediate frustration, we can say that the bond-bond correlation factors that emerge do not appear to be consistent with the existence of a columnar valence bond crystal. Furthermore, our results suggest that the magnetically disordered region may accommodate two distinct phases.