Development of a dynamic mathematical model for investigating mammary gland metabolism in lactating cows

SUMMARY A dynamic mathematical model of a closed mammary system in lactating cows was developed to incorporate the setpoint concept of tissue activity, using equations where nutrient supply and absorption are locally regulated so as to maintain a given rate of milk protein yield. The model consists of 12 differential equations, 11 of which are concerned with intracellular biochemical compartments and one describes the volume of tissue actively perfused by blood (AP). The intracellular compartments are: amino acids (AAs), acetate, fatty acids (FAs), β-hydroxybutyrate, glucose-6-phosphate, fructose-6-phosphate, phospho-glyceraldehyde, pyruvate, mitochondrial acetyl-CoA, adenosine triphosphate (ATP) and adenosine diphosphate (ADP). The model simulates mechanisms which are aimed at reproducing and, thereby, explain variations in mammary plasma flow (MPF) observed experimentally. The AP changes according to variations in the metabolic status or in the metabolic requirements of the gland. Should the tissue energy charge (i.e. ATP/ADP ratio) exceed a baseline ratio, then AP decreases and consequently MPF declines. Conversely, when milk protein yield increases, AP increases and MPF rises. In the present model, AA uptake by the mammary gland is inhibited by intracellular AAs. It is also assumed that, when milk protein yield diminishes, the respiratory chain and ATP synthesis become uncoupled and consequently ATP yield is reduced. Model evaluation included behavioural analysis and sensitivity analysis. Behaviour analysis was conducted to test whether the model mechanisms reproduced the scenarios from which the model hypotheses were developed, and took into consideration: an increase in arterial glucose concentration (HIGLC), increases in arterial concentrations of non-esterified FAs, triacylglycerol and β-hydroxybutyrate (HIFAT), a 50% reduction of arterial histidine concentration (LOHIS), and a hyperinsulinaemic euglycaemic clamp (HIINS). Both HIGLC and HIFAT resulted in a decrease in MPF and in milk protein yield; moreover, the scenario HIGLC also produced a notable decrease in the extraction of glucose. The scenario LOHIS resulted in increased MPF and extraction of His from plasma. However these responses were not sufficiently large to prevent a severe reduction of milk protein yield which was accompanied by a reduction in the extraction of other essential AAs. The scenario HIINS resulted in an increase of MPF and of milk protein yield, in the extraction of His and of other essential AAs. Model sensitivity analysis focused on variation of both affinity and inhibition constants of some of the Michaelis–Menten equations. Improvements in model structure and directions for future research suggested by the modelling analysis are discussed.