Symplectic quaternion scheme for biophysical molecular dynamics