Let $K$ be a quadratic extension of $\mathbb{Q}$. Then we know that $K$ is always a subfield of the algebra $\mathrm{M}_2(\mathbb{Q})$.
A "clean" way to produce embeddings is to observe that since $\dim_\mathbb{Q}(K)=2$ multiplication by elements in $K$ defines an embedding $$ K\longrightarrow\mathrm{End}_\mathbb{Q}(K)\simeq\mathrm{M}_2(\mathbb{Q}), $$ i.e. the choice of a $\mathbb{Q}$-basis of $K$ defines a said embedding.
Now let $D$ be a quaternion algebra over $\mathbb{Q}$. We know that embeddings $K\hookrightarrow D$ exists if and only if $D\otimes_\mathbb{Q}K\simeq\mathrm{M}_2(K)$, and actually one can characterize the situation in terms of local conditions.
The actual construction of embeddings seems a bit messier: if one writes $K=\mathbb{Q}(\sqrt d)$, one can analyze the reduced norm form in $D$ to deduce the existence of an element $x\in D$ such that $x^2=d$ and then embed $K$ in $D$ via $\sqrt d\mapsto x$.
I wonder if there are "cleaner" ways to produce embeddings $K\to D$, more akin to the split case $D=\mathrm{M}_2(\mathbb{Q})$.
It should be clear that I stated the question for quadratic extensions of $\mathbb{Q}$ and quaternion algebras over $\mathbb{Q}$ for the sake of simplicity. The question remains exactly the same replacing $\mathbb{Q}$ with any number field $E$.