Algorithm for perturbing thermodynamically infeasible metabolic networks

Networks of reactions inside the cell are constrained by the laws of mass and energy balance. The second law of thermodynamics further constrains the topology of the metabolic network, by disallowing interactions that can sustain themselves when they are isolated from the rest of the network. We present an algorithm that identifies such subnetworks and makes them feasible by deleting reactions. This perturbs the network topology, but preserves the optimality of the mass balance solutions. This method has been applied to a full network of Escherichia coli.