The b1-adrenergic signaling system plays an important role in the functioning of cardiac cells. Experimental data shows that the activation of this system produces inotropy, lusitropy, and chronotropy in the heart, such as increased magnitude and relaxation rates of [Ca2+]i transients and contraction force, and increased heart rhythm. However, excessive stimulation of b1-adrenergic receptors leads to heart dysfunction and heart failure. In this paper, a comprehensive, experimentally based mathematical model of the b1-adrenergic signaling system for mouse ventricular myocytes is developed, which includes major subcellular functional compartments (caveolae, extracaveolae, and cytosol). The model describes biochemical reactions that occur during stimulation of b1-adrenoceptors, changes in ionic currents, and modifications of Ca2+ handling system. Simulations describe the dynamics of major signaling molecules, such as cyclic AMP and protein kinase A, in different subcellular compartments; the effects of inhibition of phosphodiesterases on cAMP production; kinetics and magnitudes of phosphorylation of ion channels, transporters, and Ca2+ handling proteins; modifications of action potential shape and duration; magnitudes and relaxation rates of [Ca2+]i transients; changes in intracellular and transmembrane Ca2+ fluxes; and [Na+]i fluxes and dynamics. The model elucidates complex interactions of ionic currents upon activation of b1- adrenoceptors at different stimulation frequencies, which ultimately lead to a relatively modest increase in action potential duration and significant increase in [Ca2+]i transients. In particular, the model includes two subpopulations of the L-type Ca2+ channels, in caveolae and extracaveolae compartments, and their effects on the action potential and [Ca2+]i transients are investigated. The presented model can be used by researchers for the interpretation of experimental data and for the developments of mathematical models for other species or for pathological conditions.