Fourier Surrogate Models of Dilated Fitness Landscapes
A novel approach combining dilation functions and Fourier surrogate modeling to smooth noisy, multi-modal fitness landscapes in systems biology.
"Or how we learned to torture optimization problems until they confess."
Parameter Estimation (PE) is one of the most notoriously complex problems in Systems Biology. When inferring the kinetic parameters of biochemical systems—especially when species exist in very low amounts—the intrinsic noise of molecular collisions makes the system highly stochastic. This project introduces a novel computational approach to tame these wild, rugged fitness landscapes.
By combining dilation functions with Fourier surrogate modeling and filtering, we effectively "squash and smooth" the optimization landscape, allowing metaheuristics to find global optima that would otherwise remain hidden in the noise.
The Parameter Estimation Bottleneck
Finding an accurate parameterization to reproduce observed experimental behavior is critical for building predictive biochemical models. However, the Parameter Estimation (PE) problem is fraught with computational hurdles:
- Stochastic Noise: Random fluctuations in low-molecule-count systems mean that evaluating the exact same candidate parameterization can yield radically different fitness values.
- Rugged & Multi-modal: The fitness landscape is non-convex and non-separable, filled with local optima that cause standard optimization algorithms to converge prematurely.
- Log-uniform Distributions: Kinetic parameters naturally follow a log-uniform distribution, meaning the global optima tend to hide in the lowest orders of magnitude within the search space, making them incredibly difficult to locate.
Torturing the Landscape
To simultaneously tackle all these issues, we developed a methodology to manipulate the fitness landscape itself before attempting to optimize it.
- Dilation Functions: We apply a dilation function to “squash” the search space, stretching out the lower orders of magnitude where the global optima are most likely to reside.
- Fourier Surrogate Modeling & Filtering: To deal with the high-frequency noise introduced by stochastic simulations, we utilize a Fourier surrogate model. By applying a low-pass filter in the frequency domain, we smooth out the chaotic, high-frequency stochastic fluctuations, leaving behind a clean, navigable gradient.
This combination of spatial distortion and frequency filtering significantly simplifies low-dimensional PE instances, allowing optimization algorithms to efficiently lock onto the true parameters.
Project Highlights
Space Dilation
Strategically distorts the parameter search space to expose hidden global optima in the lowest orders of magnitude.
Fourier Filtering
Builds a surrogate model in the frequency domain to filter out the high-frequency noise inherent in stochastic biochemical simulations.
Robust Optimization
Transforms non-convex, multi-modal, and non-separable problems into smoothed landscapes that computational intelligence methods can easily solve.
Core Publication
2020
- Fourier Surrogate Models of Dilated Fitness Landscapes in Systems Biology: Or how we learned to torture optimization problems until they confessIn 2020 IEEE Conference on Computational Intelligence in Bioinformatics and Computational Biology (CIBCB), 2020