In this work, the force-constant disorder in nickel-niobium metallic glass, Ni44Nb56, was studied using the deep inelastic neutron scattering (DINS) technique augmented by isotopic substitution. The distributions of DINS observables (the nuclear kinetic energies, the width of the nuclear momentum distributions, and the effective force constants) were measured in Ni44Nb56 and compared with their counterparts obtained from ab initio harmonic lattice (HLD) simulations for the crystalline forms of nickel, niobium, and the NiNb crystal and from the reverse Monte Carlo (RMC) simulations augmented by effective force fields performed for Ni44Nb56. The force-constant distribution of nickel, obtained from the analysis of the results of the DINS experiments, was found to be two times broader than its counterparts estimated based on the HLD and RMC simulations. In the case of niobium, the force-constant distribution inferred from the DINS experiments is estimated to be an order of magnitude broader than the ab initio HLD prediction in the NiNb crystal. Moreover, no disorder-induced softening (with respect to its crystalline counterparts) of the effective force constants of Ni and Nb in Ni44Nb56 was observed. The lack of disorder-induced softening in Ni44Nb56 is consistent with the correlation between the short-range order, defined by the average coordination number and the interatomic distances, and the magnitudes of the effective force constants. The obtained results are consistent with a picture, whereby disorder induces symmetrical broadening of phonon dispersion curves, and phonon softening is limited to low-energy modes carrying negligible amounts of nuclear kinetic energy. The obtained results have important ramifications for engineering the properties of bulk metallic glasses.