As you probably know, one of the observations about the universe is that there exist frames of references where particles move in straight lines when they seem to be far away in isolation from all the other particles.
But in general, the goal of physics is to predict the trajectories of interacting particles that are not in isolation (non-straight lines in general).
Then it follows that the equations of physics should be able to tell us the acceleration of a particle, if we provide the current state of the system as input (positions and velocities of all the particles in the system).
Acceleration is the quantity we care about in the end, as the goal of physics is to predict trajectories.
Now, how do we actually make inferences about how acceleration behaves? Obviously: by making observations about acceleration.
Observation 1: In an isolated system of two particles (say, billiard balls), their accelerations at any point of time during the collision is observed to have the property:
$$-\frac{a_1(t)}{a_2(t)}=c$$
The minus sign indicates the opposite directions of $a_1$ and $a_2$. $c$ is always the same positive constant. This observation allows you to say "Particle 2 is $c$ times as massive as Particle 1". We're defining relative mass right now. In any interaction involving the two particles, particle 1 "suffers" $c$ times as much change in velocity as particle 2, in the same amount of time.
We've defined mass as the constant measuring measuring resistance to change in the state of motion. We can further choose a reference mass, call it 1 unit, and relative to it, ascribe mass values ($m_i$) to all the particles of the universe. Then the equation would become:
$$-\frac{a_2}{a_1}=\frac{m_1}{m_2}$$
or $$m_1a_1+m_2a_2=0$$
or $$m_1v_1+m_2v_2=constant$$
Moving on, so far we've only observed a relationship between two accelerations $a_1$ and $a_2$ and that happen in pair. There's still no formula to measure either $a_1$ or $a_2$, given the initial state of the particles (positions, velocities).
Enter gravitation. The acceleration formula that has been observed to explain the motion of planets is $$\frac{Gm}{r^2}$$, $m$ being the mass of the sun and $r$ being the distance between Sun and Earth.
For the first time, we actually are able to calculate the acceleration in a system of two interacting bodies.
But, hold on, this acceleration formula coupled with the previous observation is enough to calculate the trajectory of the system. Then, who cares about Force?? Why introduce a middle man called Force??
I mean..you could introduce a quantity with the formula $F_1=m_1a_1=\frac{Gm_1m_2}{r^2}$. Then the calculation of acceleration ($a_1=\frac{F_1}{m_1})$ would go one step down the road for no reason. Why introduce Force when acceleration itself has a much simpler formula?
For one, you can't deny Force has a convenient property. From the first observation, $m_1a_1=-m_2a_2$. Unlike acceleration, this guy occurs in equal an opposite pairs.
For another, gravity isnt't the only interaction that we've observed:
Enter electromagnetism: Now we're in a realm where the observed laws are best described in terms of the quantity $ma$. The observed laws are $F=k\frac{q_1q_2}{r^2}$ and $F=qvB$. $m$ no longer simply cancels out in the division. In fact, it sticks around and acceleration ends up having an uglier formula.
This, combined with the "equal and opposite property", now makes Force the natural quantity to work with. Acceleration becomes the derived quantity by $a=\frac{F}{m}$
Another reason Force is a convenient quantity: Work energy theorem
With just gravity, you could simply write down the work energy theorem as:
$$\frac{1}{2}(v_2^2-v_1^2)=\int \frac{GM}{x^2}dx$$
No force involved in the above formula. You could forcefully involve it by multiplying both sides by $m$.
This again changes in electromagnetism:
$$\frac{1}{2}m(v_2^2-v_1^2)=\int k \frac{q_1q_2}{x^2}dx$$
This is a fundamental equation in which the RHS is naturally an integral of Force. You can no longer cancel out the mass