If you think instead of the prequantum line bundle (i.e. the complex line bundle associated to your prequantum circle bundle using the standard representation of the circle on $\mathbb{C}$) then the sections of this prequantum line bundle are the wavefunctions in quantum mechanics (so the circle bundle is capturing something about the phase). Of course, most of the time in quantum mechanics, your symplectic form is exact (e.g. a cotangent bundle) and you don't need to worry about these being sections of a bundle (they're just complex-valued functions).
However, when your symplectic manifold is not exact then you need to figure out whether you want wavefunctions to be functions or sections of some bundle. I believe the following reasoning explains why you should pick the prequantum line bundle.
If we take the wavefunctions to be sections of some complex line bundle (and pick a unitary connection on the bundle) then we can try to quantise the observable $F$ by associating the operator $\nabla_{V_F}+2\pi iF$ on the space of wavefunctions (where $V_F$ is the Hamiltonian vector field associated to $F$). Now the commutator of two such operators involves a curvature term (from the commutator of the covariant derivatives) and because Dirac tells you that commutators should agree with Poisson brackets, this tells us that the curvature of the bundle should be the symplectic form. This tells you which bundle to pick.
In the end, you want your space of wavefunctions to be something more like the space of functions of position (not of both position and momentum); you therefore pick a "polarisation" of your symplectic manifold (something like a Lagrangian foliation) and restrict attention to sections of the prequantum bundle which are covariantly constant along the polarisation. For example, in a cotangent bundle, you could polarise using the Lagrangian foliation by cotangent fibres, and your wavefunctions are precisely the functions on the base manifold; or if you're on Euclidean space, you could also polarise using a horizontal Lagrangian foliation and you'd get wavefunctions of momentum (related to the wavefunctions of position by Fourier transform).
I learned this stuff (and much more) from Tyurin's beautiful book "Quantization, classical and quantum field theory and theta functions"; it's on pages 1-2.