Nonadiabatic vibronic effects in single-molecule junctions: A theoretical study using the hierarchical equations of motion approach