Suppose a stock $S$ follows $$dS(t) = \alpha(t)S(t)dt + \sigma(t)S(t)dW(t),$$ where $W(t)$ is a Brownian motion under $P$. Also suppose there is a short rate process $r(t)$. My question would be is it possible to price a stock using the risk-neutral framework, i.e. can I say $$S(t) = E^{Q}[e^{-\int_t^Tr(s)ds}S(T) \mid \mathcal{F}_t]$$ for some $T$? More specifically, say I am currently at time $t=0$ and I simulate $S$ under $Q$ (basically change the drift from $\alpha(t)$ to $r(t)$) $N$ times up to time $T$ and I want to compute what would be a price $S(t)$ for some integer $t > 0$. Can I just average $S(T)$ over $N$ and discount up to time $t$?
If this is a valid approach furthermore assume that I computed $S(t+1)$ using the same method and I am about to decide in which security to invest for a $t+1$ horizon. Then, the rate of return, $S(t+1)/S(t)-1,$ is $r(t)$. For any other stock, say $\tilde S$ with different drift but the same diffusion, the rate of return under risk neutral measure would be again $r(t)$ and obviously this simulation would not give me useful information for my investment decision. I could however model both of them under $P$ and then pick the one that has higher expected return. Could you elaborate why risk neutral modelling does not work for portfolio choice problem?