Disorder-induced zero-bias anomaly in the Anderson-Hubbard model: Numerical and analytical calculations