We use a combination of first principles density functional theory (DFT) calculations and the recently developed Prototype Electrostatic Ground State (PEGS) method to predict low energy crystal structures and study phase stability of Li-Zn and Na-Zn mixed-metal borohydride compounds, i.e., NaZn(BH 4 ) 3 , NaZn 2 (BH 4 ) 5 , LiZn(BH 4 ) 3 and LiZn 2 (BH 4 ) 5 .We find the following: (i) DFT+PEGS successfully predicts low-energy structures in these mixed-metal borohydride systems. (ii) DFT calculations show negative mixing energies in both the Li-Zn and Na-Zn borohydride systems, consistent with the observation of mixed-metal ordering in these systems. (iii) Our DFT calculations of the recently reported experimental crystal structures of NaZn 2 (BH 4 ) 5 and NaZn(BH 4 ) 3 show that the former has a negative mixing energy, while the latter has a positive mixing energy, (iv) Using the PEGS approach, we predict a new crystal structure of NaZn(BH 4 ) 3 with negative mixing energy and find that the experimental structure of NaZn 2 (BH 4 ) 5 and the PEGS obtained structure of NaZn(BH 4 ) 3 lie on the ground state convex hull. (v) In the Li-Zn borohydride system, we have used the PEGS+DFT approach to predict a stable crystal structure of new, previously unobserved stoichiometry, LiZn(BH 4 ) 3 . As a consequence of this predicted low-energy compound,