High-Intensity Focused Ultrasound (HIFU) is a popular non-invasive therapeutic tool and widely used in many clinical settings. The simulation models used for HIFU responses are computationally expensive and time-consuming. Among many numerical HIFU simulation methods, the Rayleigh-Sommerfeld approach is considered to be a reliable one. However, Rayleigh-Sommerfeld is suitable for homogeneous medium, and for a heterogeneous media, many approximations should be made in order to reduce the calculation time. In this study, we propose a fast methodology for estimating focused ultrasound pressure-temperature field responses through layered tissue media. A computationally efficient nonlinear angular spectrum-based method that can address the effects of varying attenuations, reflections and refractions from tissue layers is implemented to calculate reference datasets. From the simulation datasets, a profile function coupled with a GUI code is constructed for estimating the pressure-temperature response by using a Gaussian function and a Genetic Algorithm. The HIFU response model illustrated in this study can be advantageous and time-effective when multiple simulations are required on a similar complex patient model, and it can be used to guide and reduce the sets of simulations required for planning HIFU treatment.