KineticLaw.py 文件源码

def setHill(self, parameters):

        t_reactants = self.getReactantsFormula()
        t_modifiers = self.getModifiersFormula()

        t_kcat = parameters[0].symbol.getInternalMathFormula()
        t_kd = parameters[1].symbol.getInternalMathFormula()
        t_n = parameters[2].symbol.getInternalMathFormula()

        t_r_pow = t_reactants**t_n
        t_formula = t_kcat*t_modifiers*t_r_pow/(t_r_pow + t_kd)

        if len(self.reaction.listOfReactants) > 0:
            first_reactant = self.reaction.listOfReactants[0].getSpecies()
            if not first_reactant.hasOnlySubstanceUnits:
                t_formula *= first_reactant.getCompartment().symbol.getInternalMathFormula()

        elif len(self.reaction.listOfModifiers) > 0:
            first_modifier = self.reaction.listOfModifiers[0].getSpecies()
            if not first_modifier.hasOnlySubstanceUnits:
                t_formula *= first_modifier.getCompartment().symbol.getInternalMathFormula()

        self.__definition.setInternalMathFormula(simplify(t_formula))