Reverse engineering genetic networks using nonlinear saturation kinetics.

Area of Science:

• Systems Biology
• Computational Biology
• Genomics

Background:

• Gene regulatory networks (GRNs) govern cellular functions through complex genetic interactions.
• Nonlinear dynamics in biological systems are often modeled using Michaelis-Menten (MM) and Hill equations.
• MM kinetics are widely used in biochemistry but less explored for GRN reverse engineering.

Purpose of the Study:

• To develop a novel framework for GRN inference utilizing Michaelis-Menten kinetics.
• To create a biologically relevant and computationally efficient model for GRN analysis.
• To improve parameter estimation accuracy and speed in GRN modeling.

Main Methods:

• A coupled system of equations based on MM kinetics was proposed for GRN modeling.
• Parameter estimation was performed using trigonometric differential evolution (TDE) and a hill-climbing local search.
• A novel mutation operation was introduced to prevent population stagnation and premature convergence.
• The model was validated using both in vivo benchmark datasets and in silico data from GeneNetWeaver.

Main Results:

• The initial model assumed invariant Michaelis-Menten constants, which was later refined to allow distinct constants for each gene regulation.
• This refinement resulted in a decoupled GRN model with significantly faster parameter estimation.
• The proposed model demonstrated competitive performance compared to existing GRN inference methods.
• The inclusion of a novel mutation operation enhanced search exploitation and avoided convergence issues.

Conclusions:

• The developed MM kinetics-based framework offers a promising approach for GRN inference.
• The improved model enhances biological relevance and computational efficiency in GRN analysis.
• The method shows potential for accurate reconstruction of gene regulatory interactions from biological data.