Well, in the symplectic case, the situation is somehow much simpler as in the general Poisson case where you only can speak about coisotropic (there is no good meaning of minimal coisotropic as the rank may vary). In the symplectic case you have a theorem of Weinstein which states that a there is a tubular neighbourhood of $L$ which is symplectomorphic to a neighbourhood of the zero section of $T^*L$. Thus the question of a module structure is reduced to the case of a cotangent bundle since star products are local. For cotangent bundles there is a good understanding whether you can have a module structure on the functions on the configuration space $L$: the characteristic class of $\star$ has to be trivial. In fact, together with Martin Bordemann and Nikolai Neumaier we constructed such module structures in a series of papers in the end of the nineties. Also Markus Pflaum has some papers on this. Thus the global statement is that on $L$ you have a module structure for $\star$ iff the characteristic class of $\star$ is trivial in an tubular neighbourhood of $L$.


The module structures have (for particular star products) a very nice interpretation as global symbol calculus for differential operators on $L$. Moreover, if the char. class is not trivial but at least integral (up to some $2\pi$'s) then there is a module structure on the sections of some line bundle over $L$, coming quite close to the functions on $L$. Physically, this is important for the quantization of Dirac's magnetic monopole.


As DamienC already said, the situation in the general Poisson case or even in the general coisotropic case on symplectic manifolds is much more involved. Here my answer to 
http://mathoverflow.net/questions/64452/in-the-dictionary-between-poisson-and-quantum-what-corresponds-to-coisotropic/64478#64478
might also be of interest for you.


Oh, I forgot: the first order term can be obtained as in the flat case, at least morally. On the configuration space $L$ you chose your favorite connection then the first order term of the module structure is something like "half the Poisson bracket" which means that \begin{equation}
f \bullet \psi =  
\iota^* f \psi
+  i \hbar \iota^* 
\frac{\partial f}{\partial p_i} \frac{\partial \psi}{\partial q^i}
 + \cdots
\end{equation}
(modulo some constants I forgot) where $q^i$ are coordinates on $L$ and $p_i$ are the corresponding fiber coordinates on the cotangent bundle. This has an intrinsic meaning. Here $\iota: L \longrightarrow T^*L$ is the zero section...