I will use the notation of $\mathcal{D}(M)$ for global differential operators on $M$. This is a filtered $\mathbb{R}$-algebra, defined as follows:
$$\mathcal{D}(M)_0=C^\infty(M)$$
and for $n>0$, define $\mathcal{D}(M)_i$ inductively as:
$$\mathcal{D}(M)_i=\{T\in End_\mathbb{R}(C^\infty(M)):[T,f]\in\mathcal{D}(M)_{i-1}\text{ for all }f\in C^\infty(M)\}$$
where $[T,f]$ means commutator as linear operators on $C^\infty(M)$. Then set:
$$\mathcal{D}(M)=\bigcup\limits_i\mathcal{D}(M)_i$$
A good exercise is to check that $\mathcal{D}_1(M)=C^\infty(M)\oplus Vec(M)$, where $Vec(M)$ is the Lie algebra of smooth vector fields on $M$.
Another, equivalent way of defining it (although these aren't the only ones) is to define $\mathcal{D}(M)$ to be the set of linear operators $T\in End_\mathbb{R}(C^\infty(M))$ such that in local coordinates $(t_1,\dots,t_n)$ on $M$, $T$ is given by a finite sum $\sum\limits_{\alpha}f_\alpha\partial^\alpha$, where the sum is over multi-indices $\alpha$ (I hope the notation is clear).
Anyhow, now suppose that a Lie group $G$ acts $M$. Then $G$ acts on $C^\infty(M)$ via pullback. This is now a representation of $G$ on an algebra (it preserves the multiplicative structure of $C^\infty(M)$), so we can differentiate the action and get an action of the Lie algebra $\mathfrak{g}$ on $C^\infty(M)$ acting by derivations, i.e. a map of Lie algebras $\mathfrak{g}\to Vec(M)$.
This map $\mathfrak{g}\to Vec(M)$ of Lie algebras induces, by the universal property, a map of $\mathbb{R}$-algebras $\mathfrak{U}(\mathfrak{g})\to\mathcal{D}(M)$. This, I believe, is the desired map.