Intro:
Great answer given by Kevin. I would like to contribute an additional perspective. My experience with and my understanding of the Risk Neutral measure is entirely based on "no arbitrage" and "replication / hedging" arguments.
The way I would like to explain this view is via the following three-step construction:
(i) First, I want to build the intuition with a one-period discrete model: only a single stock and a risk-free account, no derivatives. The aim is to show that even without trying to price derivatives, one can create a mathematical object called a "risk-neutral probability measure", just by assuming no arbitrage in the model.
(ii) I then want to show that pricing a derivative by replication of its payoff with the underlying instrument and the risk-free rate instrument is equivalent to taking the expectation of the derivative pay-off under the risk-neutral measure & discounting it.
(iii) I then want to highlight that the discrete model converges to the well-known continuous Black-Scholes model.
Part 1: Discrete single-period model:
I assume that today's stock price is $S_0$, and one period from now, the stock price can be $S_0 * u=S_u$ or $S_0 * d=S_d$, with $u$ and $d$ being "up" and "down" multiplicative factors. I assume that the risk-free rate is $r$.
Without imposing some conditions on $u$, $d$ and $r$, there might be some arbitrage opportunities. If for example $e^r>u$, I could short the stock at time $t_0$ and invest the proceeds $S_0$ into the risk-free account: in both future states at time $t_1$, I could then buy the stock back for less than my proceeds $S_0e^r$ from the risk free (because $S_u$ and $S_d$ would both be less than $S_0e^r$).
If, on the other hand, $e^r<d$, I could borrow exactly $S_0$ money at $t_0$ to buy the stock, and in both future states of the world at time $t_1$, the stock value would be higher than what I would have to repay ($S_0e^r$): so there would be arbitrage again.
Imposing $d \leq e^r \leq u$, will ensure no arbitrage in the one-period model.
Now I am going to perform the following algebraic manipulation:
$$ S_0 = \frac{S_0(u-d)}{(u-d)}= \\= \frac{1}{e^r}\frac{S_0(u-d)e^r}{(u-d)}= \\ =\frac{1}{e^r}\frac{S_0(u-d)e^r+(S_0ud - S_0ud)}{(u-d)}=\\= \frac{1}{e^r}\left( \frac{S_0ue^r -(S_0ud)}{u-d} + \frac{-S_0de^r+(S_0ud)}{u-d} \right)=\\=\frac{1}{e^r}\left(S_0u \left( \frac{e^r -d}{u-d} \right) + S_0d \left(\frac{u-e^r}{u-d} \right) \right)$$
The no arbitrage condition $d \leq e^r \leq u$ will result in the following bounds:
$$0 \leq \frac{e^r -d}{u-d} \leq 1$$
$$0 \leq \frac{u-e^r}{u-d} \leq 1$$
Furthermore:
$$ \frac{e^r -d}{u-d} + \frac{u-e^r}{u-d} = 1 $$
Let's call $\frac{e^r -d}{u-d}:=p_u$ and $\frac{u-e^r}{u-d}:=p_d$. In the one period model, the stock going up and the stock going down are two different states of the world, i.e. there is no "intersection" between these states in the probabilistic sense. Therefore $p_u$ and $p_d$ are additive over disjoint sets and they are within the zero-one range, therefore mathematically, these parameters qualify as a probability measure.
Rewriting the algebraic manipulation above in terms of $p_u$ & $p_d$ yields the following:
$$ S_0 = \frac{S_u p_u + S_d p_d}{e^r} = \frac{1}{e^r}\mathbb{E} [S_1] $$
Also notice that in the entire construction above, we did not talk about the probabilities of the stock going up or down. Every market participant might have his or her Bayesian view of the world with different probabilities assigned to the stock going up or down. But the risk-neutral measure is agreed upon by the market as a whole as a consequence of no arbitrage.
This also brings up an interesting point: in my view, the risk neutral probabilities are probabilities only in the "mathematical object" sense. They do not actually represent "likelihoods", in the sense that we human beings like to interpret probabilistic events with.
Part 2: Pricing derivatives:
Let's assume we want to price a derivative on the stock with pay-off function $V(S_t)$ (could be a forward, option, whatever). The derivative pay-off in the two states will trivially be $V(S_u)$ and $V(S_d)$. We have two states, two underlying instruments: let's try to replicate the derivative pay-off in both states ($x$ is the number of stocks and $y$ is the amount invested in the risk-free account: I want to replicate the derivative pay-off in both states with $x$ stocks and $y$ risk-free investment):
$$ (i) x S_u + ye^r = V(S_u) $$
$$ (ii) x S_d + ye^r = V(S_d) $$
Solving gives:
$$ x = \frac{V(S_u)-V(S_d)}{S_0(u-d)} $$
$$ y = \frac{uV(S_d)-dV(S_u)}{(u-d)} \frac{1}{e^r} $$
Therefore the derivative price at time $t_0$ is the $x$ amount of the stock + $y$ amount invested in the risk-free account:
$$ V(S_0,t_0) = x*S_0 + y*1 = \\ = \frac{V(S_u)-V(S_d)}{S_0(u-d)}*S_0 + \frac{uV(S_d)-dV(S_u)}{(u-d)} \frac{1}{e^r}*1$$.
The above evaluates to:
$$\frac{1}{e^r}\left(V(S_u) \left( \frac{e^r -d}{u-d} \right) + V(S_d) \left(\frac{u-e^r}{u-d} \right) \right) $$
Notice that again we can write $\frac{e^r -d}{u-d}:=p_u$ and $\frac{u-e^r}{u-d}:=p_d$, where notably $p_u$ and $p_d$ are the same as in Part 1 above, Therefore, instead of having to compute the replication portfolio weights $x$ and $y$, the derivative can be priced as:
$$ V(S_0,t_0) = \frac{1}{e^r}\left(V(S_u) p_u + V(S_d) p_d \right) = \\ = \frac{1}{e^r} \mathbb{E}[V(S_1,t_1)]$$
Hopefully, by now you can see where I am going with this: the risk-neutral measure pricing technique has the following features:
(A) Is a consequence of no-arbitrage assumptions in the model
(B) Taking the expectation of a derivative pay-off and discounting it to today is the equivalent of: computing "replication portfolio" weights at each time-step, and pricing the derivative using these replicating weights at time $t_0$.
Part 3: Continous-time models:
Extending the one-period model leads to a multi-period "binomial tree" discrete model. Pricing a derivative on a multi-period tree would require working "backwards" from the terminal pay-off and computing the replicating portfolio pay-off at each node. Alternatively, the more convenient way is to use the risk-neutral expectation of the terminal pay-off and discounting it to "today": as that will produce the same result (as shown above) and will save us having to worry about the replicating portfolio weights.
There are multiple papers online showing how the binomial tree model converges to the Black-Scholes formula when the number of steps tends to infinity as $\delta t$ tends to zero (for example a great paper by John Hull here or alternatively this paper here).
As a sketch, we can show that the multi-period Binomial model for the stock converges to the well-known continuous Geometric Brownian Motion (GBM) model (which in turn can be used to derive the Black-Scholes formula directly by applying the risk-neutral expectation to the option pay-off at maturity where the stock process is simulated with GBM).
For the multi-period Binomial model, we have:
$$S_n=S_0u^kd^{n-k}$$
with $k\sim Bin(n,p)$. The "Cox-Ross-Rubinstein" parameters for "up" and "down" are then set as follows: $u:=e^{\eta\frac{T}{n}+\sqrt{\frac{T}{n}}\sigma}$, $d:=e^{\eta\frac{T}{n}-\sqrt{\frac{T}{n}}\sigma}$.
It is well known that the Binomial distribution converges to Normal, in fact, using CLT, we get:
$$\lim_{n\to\infty}k\xrightarrow{d}N(np,\sqrt{np(1-p)})$$
Going back to the equation for $S_n$, taking a log and substituting for $u$ and $d$, we get:
$$ln \left( \frac{S_n}{S_0} \right) = k \left( 2 \sqrt{\frac{T}{n}} \sigma \right) + n\left(\eta \frac{T}{n}-\sqrt{\frac{T}{n}}\sigma\right)$$
Which gives:
$$\mathbb{E}\left[ln \left( \frac{S_n}{S_0} \right)\right]=\sqrt{Tn}\sigma(2p-1)+\eta T$$
$$Var\left(ln \left( \frac{X_n}{X_0} \right) \right)=4T\sigma^2p(1-p))$$
So we see that:
$$\lim_{n\to\infty}\left( \frac{X_n}{X_0} \right)\xrightarrow{d}N(\sqrt{Tn}\sigma(2p-1)+\eta T,2\sqrt{T}\sigma\sqrt{p(1-p))}$$
Taking $p=0.5$, we get:
$$\lim_{n\to\infty}ln\left(\frac{X_n}{X_0}\right)\xrightarrow{d} N\left(\eta T,\sigma \sqrt{T}\right)$$
If we now set $\eta:=\mu -0.5\sigma^2$, we see that the discrete multiperiod Binomial model for the stock price converges exactly to the continuous Geometric-Brownian-Motion model in distribution.
Another interesting fact to note is that the replicating weight of the stock in part 1, i.e. $x$, converges to $N(d_1)$, i.e. the instantaneous option Delta.
I will conclude by producing the same summary as Kevin, but with the following additional comments:
Summary
Risk-neutral probability measures are artificial measures (agreed) made up of risk-aversion (SDF) and real-world probabilities (disagree here: don't think risk-aversion comes into it. I see it as an artificial measure entirely created by assuming the existence of no-arbitrage and completeness).
Derivatives can be priced relative to underlying assets. This hedging price can be computed as expectation with respect to the risk-neutral probability measure (agreed). Equivalent martingale measures are deeply related to the absence of arbitrage and completeness (agreed: I would say they are not just deeply related to these, they are the consequence of these).
The risk-neutral density can be estimated from observed market data (agreed: i.e. twice differentiating the Implied Vol surface with respect to strike). The risk-neutral framework connects many different approaches to derivatives pricing