In the present thermodynamic model, the optimization of bulk metallic glass forming compositions in Zr-Ti-Cu-Ni-Be system using enthalpy of chemical mixing ('H chem ) as thermodynamic, mismatch entropy ('S V /k B ) as topological and configurational entropy ('S config /R) as statistical parameters were studied. The product of 'H chem and 'S V /k B which is termed as 'P HS in the 'S config /R range of 0.9 to 1.0 can be strongly correlated to glass forming ability. Using the contributions of 'P HS values of all ternary compositions possible in Zr-Ti-Cu-Ni-Be system, the final quinary composition is designed. Composition obtained from present model was found to closely match with compositions reported in literature.