We can start with remarkably simple differential equations to define two functions: $\sin$ and $\exp$. Then we can use calculus to show how those functions are related in the form of the unit circle.
$\sin 0 = 0$; $\sin^\prime 0 = 1$; $\sin^{\prime\prime} = -sin$
$\exp 0 = 1$; $\exp^\prime = \exp$
These equations represent some very nice intuitions. $\sin$ represents a harmonic oscillator (the farther it is from 0, the faster it accelerates toward 0); and $\exp$ represents exponential growth (the bigger it is, the faster it grows).
The recursion in these definitions allows us to use the fact that $f=g \Rightarrow f^\prime = g^\prime$ to generate a sequence of iterated derivatives for each function:
$\sin, \sin^\prime, -\sin, -\sin^\prime, \sin, \sin^\prime, -\sin, -\sin^\prime, \dots$
$\exp, \exp, \dots$
Since the starting values are given for $\sin 0$, $\sin^\prime 0$, and $\exp 0$, these sequences give us the Maclaurin series for their respective functions, allowing us to calculate the functions to an arbitrary precision.
But let's go deeper. Let's define $E$ as $\theta \mapsto \exp (i\theta)$ and $S$ as $\theta \mapsto \sin^\prime \theta + i \sin \theta$. With a little algebra and some application of the chain rule, we find that the iterated derivatives for these functions are:
$E, iE, -E, -iE, E, iE, -E, -iE, \dots$
$S, iS, -S, -iS, S, iS, -S, -iS, \dots$
And $E(0) = S(0) = 1$.
Since the sequences and starting conditions are the same, that means the functions are the same. That is, $\exp (i\theta) = \sin^\prime \theta + i \sin \theta$, which is Euler's formula.
So for a complex number $x+yi$ with magnitude $1$, its angle $\theta$ is related to its coordinates $(x, y)$ by the equation $(x, y) = (\sin^\prime \theta, \sin \theta)$.
There's an even more intuitive reason why this makes sense.
A point orbiting the origin at a distance of $r$ is a point whose:
position vector $\mathbf{d}$ of magnitude $r$ points from the origin to the orbiting point
velocity vector $\mathbf{v}$ of magnitude $r$ points perpendicular to the position vector
acceleration $\mathbf{a}$ vector of magnitude $r$ points opposite to the position vector
And we can directly apply our definition of $\sin$ from above to show that a point at position $(x, y) = (\sin^\prime t, \sin t)$ at time $t$ indeed has those properties:
$\mathbf{d} = (x, y)$
$\mathbf{v} = \frac{d\mathbf{d}}{dt} = (\sin^{\prime\prime} t, \sin^\prime t) = (-\sin t, \sin^\prime t) = (-y, x)$
$\mathbf{a} = \frac{d\mathbf{v}}{dt} = (-\sin^\prime t, \sin^{\prime\prime} t) = (-\sin^\prime t, -\sin t) = (-x, -y) = -\mathbf{d}$
And according to the starting conditions, $\mathbf{d}_0 = (1, 0)$; $\mathbf{v}_0 = (0, 1)$; $\mathbf{a}_0 = (-1, 0)$