The β- and β-adrenergic signaling systems play different roles in the functioning of cardiac cells. Experimental data show that the activation of the β-adrenergic signaling system produces significant inotropic, lusitropic, and chronotropic effects in the heart, whereas the effects of the β-adrenergic signaling system is less apparent. In this paper, a comprehensive compartmentalized experimentally based mathematical model of the combined β- and β-adrenergic signaling systems in mouse ventricular myocytes is developed to simulate the experimental findings and make testable predictions of the behavior of the cardiac cells under different physiological conditions. Simulations describe the dynamics of major signaling molecules in different subcellular compartments; kinetics and magnitudes of phosphorylation of ion channels, transporters, and Ca handling proteins; modifications of action potential shape and duration; and [Ca] and [Na] dynamics upon stimulation of β- and β-adrenergic receptors (β- and β-ARs). The model reveals physiological conditions when β-ARs do not produce significant physiological effects and when their effects can be measured experimentally. Simulations demonstrated that stimulation of β-ARs with isoproterenol caused a marked increase in the magnitude of the L-type Ca current, [Ca] transient, and phosphorylation of phospholamban only upon additional application of pertussis toxin or inhibition of phosphodiesterases of type 3 and 4. The model also made testable predictions of the changes in magnitudes of [Ca] and [Na] fluxes, the rate of decay of [Na] concentration upon both combined and separate stimulation of β- and β-ARs, and the contribution of phosphorylation of PKA targets to the changes in the action potential and [Ca] transient.