We use classical molecular dynamics (MD) to study electron-ion temperature equilibration in two-component plasmas in regimes for which the presence of coupled collective modes has been predicted to substantively reduce the equilibration rate. Guided by previous kinetic theory work, we examine hydrogen plasmas at a density of n = 10 26 cm −3 , T i = 10 5 K, and 10 7 K < T e < 10 9 K. The non-equilibrium classical MD simulations are performed with inter-particle interactions modeled by quantum statistical potentials (QSPs). Our MD results indicate: 1. a large effect from time-varying potential energy, which we quantify by appealing to an adiabatic 2-temperature equation of state, and 2. a notable deviation in the energy equilibration rate when compared to calculations from classical Lenard-Balescu theory including the QSPs. In particular, it is shown that the energy equilibration rates from MD are more similar to those of the theory when coupled modes are neglected. We suggest possible reasons for this surprising result, and propose directions of further research along these lines.