High-fidelity kinetic modeling of instabilities and gyromotion physics in nonuniform low-beta plasmas