A new dynamic model for predicting road cycling individual time trials with optimal control was created. The model included both lateral and longitudinal bicycle dynamics, 3D road geometry, and anaerobic source depletion. The prediction of the individual time trial performance was formulated as an optimal control problem and solved with an indirect approach to find the pacing and cornering strategies in the respect of the physical/physiological limits of the system. The model was tested against the velocity and power output data collected by professional cyclists in two individual time trial Giro d'Italia data sets: the first data set (Rovereto, n = 15) was used to adjust the parameters of the model and the second data set (Verona, n = 13) was used to test the predictive ability of the model. The simulated velocity fell in the CI 95% of the experimental data for 32 and 18% of the duration of the course for Rovereto and Verona stages, respectively. The simulated power output fell in the CI 95% of the experimental data for 50 and 25% of the duration of the course for Rovereto and Verona stages respectively. This framework can be used to input rider's physical/physiological characteristics, 3D road geometry, and conditions to generate realistic velocity and power output predictions in individual time trials. It, therefore, constitutes a tool that could be used by coaches and athletes to plan the pacing and cornering strategies before the race.