Molchan-Golosov fractional Lévy processes (MG-FLPs) are introduced by way of a multivariate componentwise Molchan-Golosov transformation based on an n-dimensional driving Lévy process. Using results of fractional calculus and infinitely divisible distributions, we are able to calculate the conditional characteristic function of integrals driven by MG-FLPs. This leads to important predictions concerning multivariate fractional Brownian motion, fractional subordinators, and general fractional stochastic differential equations. Examples are the fractional Lévy Ornstein-Uhlenbeck and Cox-Ingersoll-Ross models. As an application we present a fractional credit model with a long range dependent hazard rate and calculate bond prices.