Context. Mass-loss by radiatively line-driven winds is central to our understanding of massive star evolution either single or in multiple systems. It for instance plays a key role in making massive star evolution different at different metallicities, specially in the case of very massive stars (M * ≥ 25M ).Aims. Here we present evolutionary models for a set of massive stars, introducing a new prescription for the mass-loss rate obtained from hydrodynamical calculations in which the wind velocity profile, (r), and the line-acceleration, g line , are obtained in a self consistently way. These new prescriptions cover the most part of the Main-Sequence phase of O-type stars. Methods. We perform a grid of self-consistent mass-loss rates Ṁsc , for a set of standard evolutionary tracks (i.e., using the old prescription for mass-loss rate) under different values for initial mass and metallicity. Based in this grid, we elaborate a statistical analysis to create a new simple formula predicting the values of Ṁsc just from the stellar parameters, without assuming any extra condition for the wind description. Therefore, replacing mass-loss rates at the Main Sequence stage from the standard Vink's formula by our new recipe, we generate a new set of evolutionary tracks for M ZAMS = 25, 40, 70 and 120 M and metallicities Z = 0.014 (Galactic), Z = 0.006 (LMC), and Z = 0.002 (SMC). Results. Our new derived formula for mass-loss rate predicts a dependence Ṁ ∝ Z a , where a is not longer constant but dependent on the stellar mass: ranging from a ∼ 0.53 when M * ∼ 120 M , to a ∼ 1.02 when M * ∼ 25 M . We found important differences between standard and our new self-consistent tracks. Models adopting the new recipe for Ṁ (which starts being around ∼ 3 times weaker than mass-loss rate from the old formulation) retain more mass during their evolution, which is expressed in larger radii and consequently more luminous tracks over the Hertzsprung-Russell diagram. These differences are more prominent for the cases of M ZAMS = 70 and 120 M at solar metallicity, where we found self-consistent tracks are ∼ 0.1 dex brighter and keep extra mass up to 20 M , compared with the classical models using the previous formulation for mass-loss rate. Later increments in the mass-loss rate for tracks when self-consistency is not longer used, attributed to the LBV stage, produce different final stellar radii and masses before the end of H-burning stage, which are analysed case to case. Moreover, we observed remarkable differences for the evolution of the radionuclide isotope 26 Al in the core and the surface of the star. Since Ṁsc are weaker than the commonly adopted values for evolutionary tracks, self-consistent tracks predict a later modification in the abundance number of 26 Al in the stellar winds. This new behaviour could provide useful information about the real contribution of this isotope from massive stars to the Galactic interstellar medium.