Adsorption isotherm-based statistical thermodynamic models can be used to determine solute concentration and solute and solvent activities in aqueous solutions. Recently, the number of adjustable parameters in the isotherm model of Dutcher et al. J. Phys. Chem. A/C 2011, 2012, 2013 were reduced for neutral solutes as well as symmetric 1:1 electrolytes by using a Coulombic model to describe the solute-solvent energy interactions (Ohm et al. J. Phys. Chem. A 2015, Nandy et al. J. Phys. Chem. A 2016). Here, the Coulombic treatment for symmetric electrolytes is extended to establish improved isotherm model equations for asymmetric 1-2 and 1-3 electrolyte systems. The Coulombic model developed here results in prediction of activities and other thermodynamic properties in multicomponent systems containing ions of arbitrary charge. The model is found to accurately calculate the osmotic coefficient over the entire solute concentration range with two model parameters, related to intermolecular solute-solute and solute-solvent spacing. The inorganic salts and acids treated here are generally considered to be fully dissociated. However, there are certain weak acids that do not dissociate completely, such as the bisulfate ion. In this work, partial dissociation of the bisulfate ion from sulfuric acid is treated as a mixture, with an additional model parameter that accounts for the dissociation ratio of the dissociated ions to nondissociated ions.