Control and inhibition analysis of complex formation processes
1 Division of Mathematical Science, Department of Systems Innovation, Graduate School of Engineering Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka, 560-8581, Japan
2 Division of Cancer Cell Research, Institute of Medical Science, University of Tokyo, 4-6-1 Shirokanedai, Minato-ku, Tokyo, 108-8639, Japan
3 Mathematical Oncology, Institute of Medical Science, University of Tokyo, 4-6-1 Shirokanedai, Minato-ku, Tokyo, 108-8639, Japan
Theoretical Biology and Medical Modelling 2012, 9:33 doi:10.1186/1742-4682-9-33Published: 3 August 2012
Proteolytic degradation of the extracellular matrix (ECM) is a key event in tumour metastasis and invasion. Matrix metalloproteinases (MMPs) are a family of endopeptidases that degrade most of the components of the ECM. Several broad-spectrum MMP inhibitors (MMPIs) have been developed, but have had little success due to side effects. Thus, it is important to develop mathematical methods to provide new drug treatment strategies. Matrix metalloproteinase 2 (MMP2) activation occurs via a mechanism involving complex formation that consists of membrane type 1 MMP (MT1-MMP), tissue inhibitor of matrix metalloproteinase 2 (TIMP2) and MMP2. Here, we focus on developing a method for analysing the complex formation process.
We used control analysis to investigate inhibitor responses in complex formation processes. The essence of the analysis is to define the response coefficient which measures the inhibitory efficiency, a small fractional change of concentration of a targeting molecule in response to a small fractional change of concentration of an inhibitor. First, by using the response coefficient, we investigated models for general classes of complex formation processes: chain reaction systems composed of ordered steps, and chain reaction systems and site-binding reaction systems composed of unordered multi-branched steps. By analysing the ordered step models, we showed that parameter-independent inequalities between the response coefficients held. For the unordered multi-branched step models, we showed that independence of the response coefficients with respect to equilibrium constants held. As an application of our analysis, we discuss a mathematical model for the MMP2 activation process. By putting the experimentally derived parameter values into the model, we were able to conclude that the TIMP2 and MMP2 interaction is the most efficient interaction to consider in selecting inhibitors.
Our result identifies a new drug target in the process of the MMP2 activation. Thus, our analysis will provide new insight into the design of more efficient drug strategies for cancer treatment.