Bridging fundamental approaches to model optimization for pharmacometricians, systems pharmacologists and statisticians is a critical issue. These fields rely primarily on Maximum Likelihood and Extended Least Squares metrics with iterative estimation of parameters. Our research combines adaptive chaos synchronization and grid search to estimate physiological and pharmacological systems with emergent properties by exploring deterministic methods that are more appropriate and have potentially superior performance than classical numerical approaches, which minimize the sum of squares or maximize the likelihood. We illustrate these issues with an established model of cortisol in human with nonlinear dynamics. The model describes cortisol kinetics over time, including its chaotic oscillations, by a delay differential equation. We demonstrate that chaos synchronization helps to avoid the tendency of the gradient-based optimization algorithms to end up in a local minimum. The subsequent analysis illustrates that the hybrid adaptive chaos synchronization for estimation of linear parameters with coarse-to-fine grid search for optimal values of non-linear parameters can be applied iteratively to accurately estimate parameters and effectively track trajectories for a wide class of noisy chaotic systems.