Deciphering Chemomechanical Couplings In Proteins Using Microsecond-Level Molecular Dynamics Simulations