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.

Optimization Fourier Modeling Systems Biology

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.

  1. 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.
  2. 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

  1. Fourier Surrogate Models of Dilated Fitness Landscapes in Systems Biology: Or how we learned to torture optimization problems until they confess
    Marco S Nobile, Paolo Cazzaniga, Simone Spolaor, and 2 more authors
    In 2020 IEEE Conference on Computational Intelligence in Bioinformatics and Computational Biology (CIBCB), 2020