Over the last few years a significant part of the JET experimental programme has been dedicated to the development of the hybrid scenario, characterized, with respect to the standard H-mode scenario, by a flatter q profile in the plasma core, a value of q0 greater than or approximately equal to 1, higher £] and improved performances in terms of H98(y,2) factor. In support of the experiments an extensive modelling activity has been carried out to interpret the results. The main focus of this activity has been the investigation of the evolution of the current profile and the study of thermal transport. This was achieved by means of semipredictive transport simulations where the density behaviour was prescribed. However, since the shape of the density profile affects important plasma parameters, like the bootstrap current or the beam deposition, density should be modelled self-consistently to understand present day experiments and to extrapolate the results to future devices. In this paper we present the first results of fully predictive simulations of JET hybrid scenarios, performed in the framework of the ITM-TF ISM group. The behaviour of the density, temperature and current profiles is simulated self-consistently, although the rotation is not taken into account. The codes deployed to perform the simulations were JINTRAC, ASTRA and CRONOS, all equipped with the physics based GLF23 and the semi-empirical Bohm/gyro-Bohm transport models.