Enhancement and wise archiving of astronomical images require an accurate estimate of the observational Point Spread Function (PSF). Although modelling of the telescope and its optics is a well-understood problem, PSF reconstruction becomes challenging when the observations include adaptive optics (AO) correction. The approach presented in this paper consists in feeding an end-to-end (E2E) simulation of the telescope, the instrument and its environment with the characterized disturbances from the telemetry and AO loop data, in order to produce the estimated PSFs. This method benefits from the developments made in the last years with respect to the estimation of external disturbances during AO correction, such as turbulence profile and its dynamics as well as sensor noise and vibrations characteristics. In particular, characterization of the turbulence profile in terms of strength, C 2 n (h), and outer scale, L 0 (h), is considered with an example on on-sky recorded AO telemetry from the GALACSI AO system. Once identified, the internal and external parameters of the observing conditions are used as inputs to carry out E2E simulations of the optical propagation and estimate the PSF. The method can be regarded as a "brute force" approach, as it is highly intensive in computer power; particularly for the ELTs. However, its ability to integrate complex combination of effects from all disturbances and not relying on analytical approximations for the aliasing or fitting errors makes the approach worth of a deeper study. E2E simulations have been used before in PSF reconstruction, but limited to a theoretical modelling of the system. Here, the development of the E2E simulation part is an ongoing work. A simplified AO system similar to the GALACSI WFM is currently simulated to obtain the PSF estimates and illustrate how such approach allows to account for the anisoplanatic effects and for the influence of the outer scale values.