Using citrate synthase from the hyperthermophile Pyrococcus furiosus (PfCS) as our test molecule, we show through guanidine hydrochloride-induced unfolding that the dimer separates into folded, but inactive, monomers before individual subunit unfolding takes place. Given that forces across the dimer interface are vital for thermostability, a robust computational method was derived that uses the University of Houston Brownian Dynamics (UHBD) program to calculate both the hydrophobic and electrostatic contribution to the dimerisation energy at 100°C. The results from computational and experimental determination of the lowered stability of interface mutants were correlated, being both of the same order of magnitude and placing the mutant proteins in the same order of stability. This computational method, optimised for hyperthermophilic molecules and tested in the laboratory, after further testing on other examples, could be of widespread use in the prediction of thermostabilising mutations in other oligomeric proteins for which dissociation is the first step in unfolding.