Binary interaction parameters used in the UNIQUAC activity coefficient model are found to be dependent on each other and related by a linear relation termed as the closure equation. For a ternary system, six binary interaction parameters are related by one closure equation. Similarly for quaternary systems, three independent closure equations are obtained for the twelve binary interaction parameters and for quinary systems there are six closure equations for twenty parameters. Each closure equation consists of six parameters. The binary interaction parameters that do not satisfy the closure equations may lead to a less accurate prediction of liquid-liquid equilibria. In this work the binary interaction parameters have been estimated with and without closure equations for few ternary and quaternary systems. Parameters that satisfy the closure equations exhibit better root mean square deviation than those that do not satisfy the closure equations in most of the cases. A similar behavior is observed for NRTL model also.