Hostname: page-component-78c5997874-4rdpn Total loading time: 0 Render date: 2024-11-16T19:26:05.879Z Has data issue: false hasContentIssue false

Further solutions to an isotope dilution model for partitioning phenylalanine and tyrosine between milk protein synthesis and other metabolic fates by the mammary gland of lactating dairy cows

Published online by Cambridge University Press:  11 May 2022

L. A. Crompton
Affiliation:
Department of Animal Sciences, School of Agriculture, Policy and Development, University of Reading, Whiteknights, Reading RG6 6EU, UK
L. A. McKnight
Affiliation:
Centre for Nutrition Modelling, Department of Animal Biosciences, University of Guelph, Guelph, Ontario N1G 2W1, Canada Trouw Nutrition AgResearch Canada, 150 Research Lane, Guelph, ON N1G 4T2, Canada
C. K. Reynolds
Affiliation:
Department of Animal Sciences, School of Agriculture, Policy and Development, University of Reading, Whiteknights, Reading RG6 6EU, UK
J. L. Ellis
Affiliation:
Centre for Nutrition Modelling, Department of Animal Biosciences, University of Guelph, Guelph, Ontario N1G 2W1, Canada
J. Dijkstra
Affiliation:
Animal Nutrition Group, Wageningen University & Research, 6700 AH Wageningen, The Netherlands
J. France*
Affiliation:
Centre for Nutrition Modelling, Department of Animal Biosciences, University of Guelph, Guelph, Ontario N1G 2W1, Canada
*
Author for correspondence: J. France, E-mail: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Phenylalanine (PHE) and to a lesser extent TYR are two commonly used amino acid tracers for measuring protein metabolism in a variety of species and tissues. The model examined in this paper was developed to resolve trans-organ and stable isotope dilution data collected from experiments with lactating dairy cows using these tracers. Two methods of solving the model, i.e. as two four-pool submodels, one representing PHE and the other TYR, or as an integrated eight-pool model, are investigated and the alternative solutions are contrasted. Solving the model as the two four-pool submodels rather than the integrated 8-pool model is preferred as the equations are slightly simpler and their application less susceptible to any compounding of measurement errors. The data used to illustrate the model were taken from experiments conducted to investigate the effects of high and low protein diets on the partitioning of PHE and TYR between milk protein synthesis and other metabolic fates by the mammary gland.

Type
Animal Research Paper
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
Copyright © The Author(s), 2022. Published by Cambridge University Press

Introduction

The aim of dairy researchers and producers alike is to increase the conversion efficiency of dietary nutrients into milk by dairy animals. The efficiency of converting dietary nitrogen into milk protein output is poor at 25–30% (Lobley, Reference Lobley2003) and in recent years has become a focus of attention, due to the increasing problem of environmental pollution related to emissions of nitrogen (N) to the environment, principally ammonia and nitrous oxide to air, nitrate to groundwater and particulate N to surface waters from dairy production systems (Dijkstra et al., Reference Dijkstra, Oenema, van Groenigen, Spek, van Vuuren and Bannink2013). Milk protein content and yield can be increased by dietary protein supplementation or by gastrointestinal infusion of protein or amino acids. However, the responses attained cannot be predicted accurately using current requirement-based feeding schemes for dairy cows (e.g. Thomas, Reference Thomas2004), or by supply-driven response models (e.g. NASEM, 2021) in view of the large metabolic flexibility of lactating cattle to handle variation in supply of nutrients, in particular supply of amino acids (e.g. Nichols et al. Reference Nichols, Bannink and Dijkstra2019b). The synthesis of milk protein in the mammary gland requires high amounts of amino acids, either extracted from the circulation or synthesized de novo in the gland. To address such issues, research programmes have focused on identifying and understanding the factors and mechanisms regulating the partitioning of amino acids between milk and constitutive proteins in the mammary gland.

Much of the earlier knowledge on amino acid and protein metabolism in the lactating mammary gland was derived from the perfused mammary gland (e.g. Roets et al., Reference Roets, Massart-Leen, Verbeke and Peeters1983) and from in vivo measurements of tissue protein synthesis (e.g. Baracos et al., Reference Baracos, Brun-Bellut and Marie1991) using dairy goats. Such studies provided information on the metabolic pathways of milk synthesis but were generally limited by single measurements. Studies of tissue protein metabolism require estimates of the rate of both protein synthesis and degradation, which are routinely estimated in different animals over dissimilar periods. To overcome these difficulties, several laboratories have developed an alternative indirect approach to repeatedly measure the partition and metabolism of amino acids across tissue beds in vivo and estimate the rates of constitutive protein turnover. The technique involves sampling the blood supplying and draining a tissue bed and measuring the rate of blood flow across the tissue in combination with the unidirectional uptake and release of a tracer and tracee amino acid of choice. The procedures have been applied across the mammary gland of lactating goats (Oddy et al., Reference Oddy, Lindsay and Fleet1988; Bequette et al., Reference Bequette, Kyle, Crompton, Anderson and Hanigan2002) and dairy cows (France et al., Reference France, Bequette, Lobley, Metcalf, Wray-Cahen, Dhanoa, Backwell, Hanigan, MacRae and Beever1995; Raggio et al., Reference Raggio, Lemosquet, Lobley, Rulquin and Lapierre2006; Huang et al., Reference Huang, Yoder, Teixeira and Hanigan2021), generally using leucine as a single tracer amino acid.

Phenylalanine (PHE) and to a lesser extent tyrosine (TYR) have been widely used as tracer amino acids, for the study of muscle protein metabolism in particular, in several species, because they have minimal enrichment gradients between extracellular and intracellular compartments, low endogenous turnovers and neither amino acid is synthesized or oxidized in muscle. The major site of catabolism for PHE and TYR is the liver, where PHE can be converted to TYR by hydroxylation and both PHE and TYR are catabolized. For the mammary gland, early in vitro studies demonstrated significant conversion of PHE to TYR in bovine secretory tissue (Jorgensen and Larson, Reference Jorgensen and Larson1968). Later in vivo evidence suggests hydroxylation of PHE across the mammary gland in ruminants is negligible, even under conditions of increased PHE supply (e.g. Bequette et al., Reference Bequette, Backwell, Kyle, Calder, Buchan, Crompton, France and MacRae1999); other recent in vivo data with dairy cattle indicate the mammary gland uptake to output ratio of PHE and TYR does not differ from unity at low levels of metabolizable protein supply, but is significantly higher and lower than unity, respectively, at high levels of metabolizable protein supply, indicating PHE hydroxylation may depend on protein supply (e.g. Nichols et al., Reference Nichols, Bannink, Doelman and Dijkstra2019a). Therefore, in the mammary gland during established lactation, where the amount of secretory tissue is in equilibrium, the uptake of the sum of PHE and TYR reflects their output in milk protein. That PHE and TYR appear to be taken up by the mammary gland in sufficient quantities to match their requirement for milk protein synthesis has led to the two amino acids being extensively used to estimate blood flow rate across the mammary gland by application of the Fick principle (e.g. Cant et al., Reference Cant, DePeters and Baldwin1993).

The model examined herein was initially developed by Crompton et al. (Reference Crompton, France, Reynolds, Mills, Hanigan, Ellis and Dijkstra2014) to resolve trans-organ and isotope dilution data collected from experiments with lactating dairy cows. The experiments were undertaken to investigate the effects of high and low protein diets on the partitioning of PHE and TYR between milk protein synthesis and constitutive protein synthesis. Stable isotopes were used as the tracer in these studies. The present paper describes an evolution of the work reported by Crompton et al. (Reference Crompton, France, Reynolds, Mills, Hanigan, Ellis and Dijkstra2014) in which a second method of solving the model is proposed and the alternative solutions contrasted, and a larger data set is investigated. The model provides an effective means of generating information about the fates of phenylalanine and tyrosine in the mammary gland and could be used as part of a more complex system describing amino acid metabolism in the whole ruminant or other species.

The model

Overall scheme

The overall scheme is depicted in Fig. 1. It contains four intracellular and four extracellular pools. The intracellular pools are free PHE (pool 4), PHE in milk protein (pool 3), free TYR (pool 5) and TYR in milk protein (pool 6), while the extracellular ones represent arterial PHE and TYR (pools 1 and 2) and venous PHE and TYR (pools 7 and 8). The flows of PHE and TYR between pools and into and out of the system are shown as arrowed lines. The extracellular arterial PHE pool 1 has a single inflow: entry into the pool, F 10, and two outflows: uptake by the mammary gland, F 41, and release into the extracellular venous PHE pool, F 71. The extracellular arterial TYR pool 2 also has a single inflow: entry into the pool, F 20, and two outflows: uptake by the mammary gland, F 52, and release into the extracellular venous TYR pool, F 82. The milk protein-bound PHE pool 3 has a single inflow: from free PHE, F 34, and two outflows: secretion of protein in milk, F 03, and degradation, F 43. The intracellular free PHE pool 4 has three inflows: from the degradation of constitutive mammary gland protein, F 40, from the extracellular arterial pool, F 41, and from degradation of milk protein, F 43. The pool has five outflows: secretion in milk,$F_{04}^{( m ) }$, synthesis of constitutive mammary gland protein, $F_{04}^{( s ) }$, incorporation into milk protein, F 34, hydroxylation to the intracellular free TYR pool, F 54, and outflow to the extracellular venous PHE pool, F 74. The intracellular free TYR pool 5 has four inflows: from the degradation of constitutive mammary gland protein, F 50, from the extracellular arterial TYR pool, F 52, from the intracellular PHE pool, F 54, and from the degradation of milk protein, F 56. The pool has five outflows: secretion in milk, $F_{05}^{( m ) }$, oxidation and TYR degradation products, $F_{05}^{( o ) }$, synthesis of constitutive mammary gland protein, $F_{05}^{( s ) }$, incorporation into milk protein, F 65, and outflow to the extracellular venous TYR pool, F 85. The milk protein-bound TYR pool 6 has one inflow: from the intracellular free TYR pool, F 65, and two outflows: secretion of protein in milk, F 06, and degradation, F 56. The extracellular venous PHE pool 7 has two inflows: bypass from the arterial PHE pool, F 71, and release from the intracellular PHE pool, F 74, and one outflow out of the system, F 07. The extracellular venous TYR pool 8 also has two inflows: bypass from the arterial TYR pool, F 82, and release from the intracellular TYR pool, F 85, and one outflow from the system, F 08.

Fig. 1. Scheme for the uptake and utilization of PHE and TYR by the mammary gland of lactating dairy cows as described by Crompton et al. (Reference Crompton, France, Reynolds, Mills, Hanigan, Ellis and Dijkstra2014). The small circles indicate flows out of the system that need to be measured experimentally.

This scheme can be solved as an eight-pool model (Crompton et al., Reference Crompton, France, Reynolds, Mills, Hanigan, Ellis and Dijkstra2014). Alternatively, it can also be solved by decomposing it into two four-pool schemes (i.e. a PHE sub-model and a TYR sub-model), then linking the two schemes. The PHE and TYR sub-models are both similar structurally to the model of LEU kinetics presented by France et al. (Reference France, Bequette, Lobley, Metcalf, Wray-Cahen, Dhanoa, Backwell, Hanigan, MacRae and Beever1995).

PHE sub-model

The schemes adopted for the movement of total and labelled PHE in the PHE sub-model are shown in Figs 2(a) and (b) respectively. The fundamental equations are (mathematical notation is defined in Table 1):

(1)$$\displaystyle{{{\rm d}Q_1} \over {{\rm d}t}} = F_{10}-F_{41}-F_{71}$$
(2)$$\displaystyle{{{\rm d}Q_3} \over {{\rm d}t}} = F_{34}-F_{03}-F_{43}$$
(3)$$\displaystyle{{{\rm d}Q_4} \over {{\rm d}t}} = F_{40} + F_{41} + F_{43}-F_{04}^{( m ) } -F_{04}^{( s ) } -F_{34}-F_{54}-F_{74}$$
(4)$$\displaystyle{{{\rm d}Q_7} \over {{\rm d}t}} = F_{71} + F_{74}-F_{07}$$

and for [13C] labelled PHE:

(5)$$\displaystyle{{{\rm d}q_1} \over {{\rm d}t}} = I_1-e_1( {F_{41} + F_{71}} ) $$
(6)$$\displaystyle{{{\rm d}q_3} \over {{\rm d}t}} = e_4F_{34}-e_3( {F_{03} + F_{43}} ) $$
(7)$$\displaystyle{{{\rm d}q_4} \over {{\rm d}t}} = e_1F_{41} + e_3F_{43}-e_4( {F_{04}^{( m ) } + F_{04}^{( s ) } + F_{34} + F_{54} + F_{74}} ) $$
(8)$$\displaystyle{{{\rm d}q_7} \over {{\rm d}t}} = e_1F_{71} + e_4F_{74}-e_7F_{07}$$

Fig. 2. Scheme for the uptake and utilization of PHE by the mammary gland of lactating dairy cows: (a) total PHE and (b) [13C] labelled PHE. The small circles in Fig. 2(a) indicate flows out of the system which need to be measured experimentally.

Table 1. Principle symbols used for the kinetic model

a Total material (i.e. tracee + tracer).

When the system is in steady state with respect to both total and labelled PHE, the derivative terms in these 8 differential equations are zero. For the scheme assumed, the enrichment of the intracellular milk protein-bound pool equalizes with that of the free pool as steady state is approached (i.e. e 3 = e 4). After equating intracellular enrichments and eliminating redundant equations, the four differential equations for labelled PHE, Eqns (5) to (8), yield the following three identities:

(9)$$I_1-e_1( {F_{41} + F_{71}} ) = 0$$
(10)$$e_1F_{41}-e_3( {F_{04}^{( m ) } + F_{04}^{( s ) } + F_{34} + F_{54} + F_{74}-F_{43}} ) = 0$$
(11)$$e_1F_{71} + e_3F_{74}-e_7F_{07} = 0$$

To obtain steady state solutions to this sub-model, it is assumed that free PHE in milk, PHE secreted in milk protein and PHE removal from the venous pool (i.e. $F_{04}^{( m ) }$, F 03 and F 07, respectively) can be measured experimentally. Algebraic manipulation of Eqns (1) to (4) with the derivatives set to zero, together with Eqns (9) to (11) gives:

(12)$$F_{10} = I_1/e_1$$
(13)$$\overline {F_{34}-F_{43}} = F_{03}$$
(14)$$F_{71} = \left({\displaystyle{{e_7-e_3} \over {e_1-e_3}}} \right)F_{07}{\rm ; \;\ }e_1\ne e_3$$
(15)$$F_{41} = F_{10}-F_{71}$$
(16)$$F_{74} = F_{07}-F_{71}$$
(17)$$F_{40} = \left({\displaystyle{{e_1-e_3} \over {e_3}}} \right)F_{41}$$
(18)$$\overline {F_{04}^{( s ) } + F_{54}} = F_{40} + F_{41}-F_{04}^{( m ) } -\overline {F_{34}-F_{43}} -F_{74}$$

where for these equations the italics denote steady state values of flows and enrichments, and the over-lining indicates coupled flows (those which cannot be determined separately by the sub-model). The net flow $\overline {F_{34}-F_{43}}$ may be uncoupled by assuming that a fixed proportion (~0.1) of the nascent milk protein is cleaved during the docking and secretory processes (Razooki Hasan et al., Reference Razooki Hasan, White and Mayer1982).

TYR sub-model

The schemes adopted for the movement of total and labelled TYR in the TYR sub-model are shown in Figs 3(a) and (b) respectively. The fundamental equations are:

(19)$$\displaystyle{{{\rm d}Q_2} \over {{\rm d}t}} = F_{20}-F_{52}-F_{82}$$
(20)$$\displaystyle{{{\rm d}Q_5} \over {{\rm d}t}} = F_{50} + F_{52} + F_{54} + F_{56}-F_{05}^{( m ) } -F_{05}^{( o ) } -F_{05}^{( s ) } -F_{65}-F_{85}$$
(21)$$\displaystyle{{{\rm d}Q_6} \over {{\rm d}t}} = F_{65}-F_{06}-F_{56}$$
(22)$$\displaystyle{{{\rm d}Q_8} \over {{\rm d}t}} = F_{82} + F_{85}-F_{08}$$

and for [2H] labelled TYR:

(23)$$\displaystyle{{{\rm d}\phi _2} \over {{\rm d}t}} = {\rm \Phi }_2-\varepsilon _2( {F_{52} + F_{82}} ) $$
(24)$$\displaystyle{{{\rm d}\phi _5} \over {{\rm d}t}} = \varepsilon _2F_{52} + \varepsilon _6F_{56}-\varepsilon _5( {F_{05}^{( m ) } + F_{05}^{( o ) } + F_{05}^{( s ) } + F_{65} + F_{85}} ) $$
(25)$$\displaystyle{{{\rm d}\phi _6} \over {{\rm d}t}} = \varepsilon _5F_{65}-\varepsilon _6( {F_{06} + F_{56}} ) $$
(26)$$\displaystyle{{{\rm d}\phi _8} \over {{\rm d}t}} = \varepsilon _ 2F_{82} + \varepsilon _5F_{85}-\varepsilon _8F_{08}$$

Fig. 3. Scheme for the uptake and utilization of TYR by the mammary gland of lactating dairy cows: (a) total TYR and (b) [2H] labelled TYR. The small circles in Fig. 3(a) indicate flows out of the system which need to be measured experimentally.

When the system is in steady state with respect to both total and labelled TYR, the derivative terms in these 8 differential equations are zero. For the scheme assumed, the enrichment of the intracellular milk protein-bound pool equalizes with that of the free pool as steady state is approached (i.e. ɛ 6 = ɛ 5). After equating intracellular enrichments and eliminating redundant equations, the four differential equations for labelled TYR, Eqns (23) to (26), yield the following three identities:

(27)$${\rm \Phi }_2-\varepsilon _2( {F_{52} + F_{82}} ) = 0$$
(28)$$\varepsilon _2F_{52}-\varepsilon _6( {F_{05}^{( m ) } + F_{05}^{( o ) } + F_{05}^{( s ) } + F_{65} + F_{85}-F_{56}} ) = 0$$
(29)$$\varepsilon _ 2F_{82} + \varepsilon _6F_{85}-\varepsilon _8F_{08} = 0$$

To obtain steady state solutions to the sub-model, it is assumed that free TYR in milk, CO2 production, TYR secreted in milk protein and TYR removal from the venous pool (i.e. $F_{05}^{( m ) }$, $F_{05}^{( o ) }$, F 06 and F 08, respectively) can be measured experimentally. Algebraic manipulation of Eqns (19) to (22) with the derivatives set to zero, together with Eqns (27) to (29), gives:

(30)$$F_{20} = {\rm \Phi }_2/\varepsilon _2$$
(31)$$\overline {F_{65}-F_{56}} = F_{06}$$
(32)$$F_{82} = \left({\displaystyle{{\varepsilon_8-\varepsilon_6} \over {\varepsilon_2-\varepsilon_6}}} \right)F_{08}{\rm ; \;\ }\varepsilon _2\ne \varepsilon _6$$
(33)$$F_{52} = F_{20}-F_{82}$$
(34)$$F_{85} = F_{08}-F_{82}$$
(35)$$\overline {F_{50} + F_{54}} = \left({\displaystyle{{\varepsilon_2-\varepsilon_6} \over {\varepsilon_6}}} \right)F_{52}$$
(36)$$F_{05}^{( s ) } = \overline {F_{50} + F_{54}} + F_{52}-F_{05}^{( m ) } -F_{05}^{( o ) } -\overline {F_{65}-F_{56}} -F_{85}$$

where the italics denote steady state values of flows and enrichments for these equations, and over-lining indicates coupled flows (which cannot be separately estimated by the sub-model). The net flow $\overline {F_{65}-F_{56}}$ may be uncoupled as described for the uncoupling of $\overline {F_{34}-F_{43}}$.

Linking the PHE and TYR sub-models

The two sub-models can be linked by considering constitutive mammary protein. Assuming a fixed protein composition for constitutive mammary tissue, then the ratios of TYR to PHE in protein synthesized and protein degraded (μmol TYR/μmol PHE) are equal:

(37)$$\displaystyle{{F_{50}} \over {F_{40}}} = \displaystyle{{F_{05}^{( s) } } \over {F_{04}^{( s) } }}$$

This assumption allows Eqns (18) and (35) to be uncoupled. Differencing these coupled flows:

(38)$$\overline {F_{04}^{( s) } + F_{54}} -\overline {F_{50} + F_{54}} = F_{04}^{( s) } -F_{50} = b$$

Using Eqn (37) to substitute for F 50 in the above equation:

$$F_{04}^{( s) } -\displaystyle{{F_{05}^{( s) } F_{40}} \over {F_{04}^{( s) } }} = b$$
$$[ {F_{04}^{( s) } } ] ^2-bF_{04}^{( s) } -F_{05}^{( s) } F_{40} = 0$$

Solving this quadratic:

(39)$$F_{04}^{( s) } = \displaystyle{{b + \sqrt {b^2 + 4F_{05}^{( s) } F_{40}} } \over 2}$$

Note that only positive roots of this quadratic are permissible, so any negative roots must be discarded. Therefore:

(40)$$F_{50} = F_{04}^{( s) } -b$$
(41)$$F_{54} = \overline {F_{04}^{( s) } + F_{54}} -F_{04}^{( s) } $$

The overall scheme can now be solved by computing Eqns (12) to (18), (30) to (36) and (38) to (41) sequentially.

Application

To illustrate application of the model, 4 datasets obtained with 3 cows were taken from an experiment conducted at our laboratories in the United Kingdom with multi-catheterized and mid-lactation cows. The data were taken from a trial using multiparous Holstein-Friesian dairy cows (average body weight 622 kg) 21 weeks into lactation with an average milk yield of 23.2 kg/day. The cows were fed hourly and ad libitum by auto-feeders at two levels of dietary crude protein (CP), based on a diet consisting of chopped Lucerne hay and grass silage [50% of diet dry matter (DM)] with the remaining 50% of diet DM provided as either a low (L, 108 g/kg) or high (H, 206 g/kg) protein concentrate. Dietary CP levels were 117 and 168 g/kg DM for L and H respectively, and so provided different levels of PHE and TYR supply to the small intestine for absorption. The average daily intakes were 20.1 kg DM. Diets were fed for 6 weeks before the cows were given constant jugular vein infusions of sterile saline for 3 day, followed by a buffered mixture of essential amino acids for a further 3 day. The essential amino acids were administered at a daily rate equivalent to the essential amino acids in 600 g milk protein (316 g essential amino acids/day). On the final day of each 3 day infusion, the animals received a primed, constant jugular vein infusion of [1-13C]PHE (350 mg/h) and [2,3,5,6-2H]TYR (100 mg/h) in sterile saline for 6 h, and 6 hourly blood sample sets were taken simultaneously from catheters in the carotid artery and subcutaneous abdominal vein for the measurement of blood flow rate (by PAH dilution) and nutrient metabolism by the mammary gland. Blood samples were centrifuged at 4°C for 20 min at 2000 g and the plasma stored at −20°C for subsequent analysis.

The relevant experimental measurements are given in Table 2. They are reported for two animals during the saline infusion (1 low protein diet; 1 high protein) and two animals during the amino acid infusion (both low protein diet). Values are based on plasma rather than whole blood. Phenylalanine and TYR measurements are based on free rather than total (i.e. free plus bound) plasma PHE and TYR. The output of PHE and TYR in milk protein were calculated using the protein content of milk and the amino acid composition of milk proteins (Maas et al., Reference Maas, France and McBride1997). The effective isotope infusion rates to the mammary gland, I 1 and Φ2 were calculated from the arterial concentrations and enrichments of PHE and TYR and plasma flow rate across the gland. Flows F 07 and F 08 were determined from venous PHE and TYR concentration and plasma flow rate. Flows F 04(m) and F 05(m) were assigned a value of zero (for further justification see Mehaia and Al-Kanhal, Reference Mehaia and Al-Kanhal1992). As the enrichment of the intracellular free pool was not measured directly, e4 was assumed to equal e3 and ɛ 6 was assumed to equal ɛ 5 as steady state is approached, in line with other reports (e.g. Huang et al., Reference Huang, Yoder, Teixeira and Hanigan2021). There was no detectable appearance of labelled 13CO2 across the mammary gland, indicating zero oxidation of TYR (F 05(o) = 0); an observation supported by the study of Lemosquet et al. (Reference Lemosquet, Lobley, Koopman, Van Loon, Kies and Lapierre2010). The solutions to the split model described herein are shown in Table 3. Combining input data reported herein and those reported by Crompton et al. (Reference Crompton, France, Reynolds, Mills, Hanigan, Ellis and Dijkstra2014) enables the solutions from the two 4-pool schemes to be contrasted with corresponding solutions obtained using the integrated 8-pool model (Crompton et al., Reference Crompton, France, Reynolds, Mills, Hanigan, Ellis and Dijkstra2014). The averaged flows from both models are shown in Table 4 and highlight the unanimity of calculated flows between the two schemes. Linking the PHE and TYR sub-models affects flows F 04(s), F 50 and F 54 representing constitutive protein synthesis and degradation and PHE hydroxylation, compared to the 8-pool model.

Table 2. Experimental and other inputs (symbols are defined in the text and Table 1)

a Essential amino acid infusion.

b Average over the 3-day infusion period.

Table 3. Phenylalanine and tyrosine uptake and partition by the mammary gland of four lactating dairy cows obtained using the two four-pool models (symbols are defined in the text and Table 1)

a Essential amino acid infusion.

Table 4. Phenylalanine and tyrosine uptake and partition by the mammary gland of lactating dairy cows, obtained using the two four-pool models (symbols are defined in the text and Table 1) and the corresponding solutions obtained using the eight-pool model of Crompton et al. (Reference Crompton, France, Reynolds, Mills, Hanigan, Ellis and Dijkstra2014)

Values are means across datasets (both those reported here and those reported by Crompton et al. (Reference Crompton, France, Reynolds, Mills, Hanigan, Ellis and Dijkstra2014)). Figures in parentheses are standard errors of the means.

An analysis of measurement errors in experimental enrichments and infusion rates on model solutions was undertaken. Input datasets from the integrated model of Crompton et al. (Reference Crompton, France, Reynolds, Mills, Hanigan, Ellis and Dijkstra2014) and the split model described herein were averaged to provide the initial unperturbed values for e 1, e 3, e 7, ɛ2, ɛ6, ɛ8, I 1, Φ2, F 03, F 06, F 07, F 08. Inputs to the split model were then perturbed sequentially by 0, ±10% and ±20%. Each calculated flow (y, μmol/min) was plotted against the perturbation (x, %), and a five-point linear regression of y on x was performed to determine the slope of the line produced. The average slope was subsequently scaled by its corresponding unperturbed average flow value, giving the dimensions of the scaled slope of % change in y per % change in x. Results of the error assessment are presented in Table 5. In general, errors in infusion rates and prescribed flows had less impact on the sensitivity of model solutions than errors in the measurement of isotopic enrichment. Perturbing all inputs caused marked changes in the flow representing PHE hydroxylation (F 54).

Table 5. Average slope (%) for each of the flows calculated using the model described herein, obtained by perturbing each input in turn (symbols are defined in the text and Table 1)a

a The slope for each flow is expressed relative to the value of the flow obtained when no perturbation is made. Only slopes which differ from zero are shown.

b Model solved by perturbing each input in turn by 0%, ±10% and ±20%.

c Values calculated from the mean of inputs reported in Table 2 and inputs reported by Crompton et al. (Reference Crompton, France, Reynolds, Mills, Hanigan, Ellis and Dijkstra2014).

Discussion

Increasing the efficiency of conversion of feed N into milk and meat N in ruminant production is an integral part of the effort to increase global food production while decreasing agriculture's environmental impact. Improving N utilization in the ruminant is dependent on a clear understanding of post-absorptive amino acid metabolism. The present model describes the partitioning of the indispensable amino acid PHE (and TYR) in the bovine mammary gland. It was constructed to interpret isotope dilution data from in vivo trans-organ studies with dairy cows undertaken at our laboratories. The model gave estimates of PHE flow across the mammary gland, rates of PHE and TYR incorporation into constitutive and export protein synthesis, and the rate of hydroxylation of PHE to TYR. The four datasets gave biologically feasible solutions, i.e. the computed flows were all non-negative solutions.

Hepatic removal of group 1 amino acids (amino acids transferred from the mammary arterial blood supply into milk in a 1:1 ratio, including PHE + TYR) affects the efficiency of conversion of absorbed amino acids into milk amino acids and its associated changes in milk N secretion and urine N excretion (Nichols et al., Reference Nichols, Bannink and Dijkstra2019b). Phenylalanine and TYR net uptake across the mammary gland were highest for the cows receiving the essential amino acid infusion and lowest for the cow receiving saline (range 66 to 96 μmol PHE/min; 62 to 95 μmol TYR/min). The ratio of mammary net uptake to milk output varied from 0.85 to 1.10 for PHE and 0.82 to 0.87 for TYR. The ratio of TYR to PHE in synthesized constitutive protein and degraded constitutive protein was the same in each animal (range 0.70 to 1.72) and averaged 1.17. The ratio of constitutive protein synthesis to degradation was 0.84 for both PHE and TYR. However, in mid to late lactation dairy cows, the ratio of PHE and TYR synthesis to degradation should be equal. Model estimates of intracellular PHE and TYR partitioning must be interpreted with caution due to methodological limitations and imposed assumptions. The rate of hydroxylation of PHE to TYR (F 54) was small, on average representing 4.2% of PHE inflow. When the integrated model (Crompton et al., Reference Crompton, France, Reynolds, Mills, Hanigan, Ellis and Dijkstra2014) was used to solve the combined data, three out of the 8 datasets gave small negative values (i.e. infeasible solutions) for the rate of PHE hydroxylation. This indicates that the compounding of measurement errors is likely to be of greater concern when using the integrated scheme as opposed to the split model.

The model is considered simple because it requires relatively few measurements. Due to this simplicity, measurements need to be accurate to avoid errors in the flow calculations. In our model, plasma is the only entity exchanging labelled and unlabelled amino acids with the mammary glands. It might be argued that whole blood would give more reliable measurements regarding amino acid uptake due to the presence of packed cells. However, studies have demonstrated that erythrocytes make only a minor contribution to amino acid uptake by the mammary gland (Mackle et al., Reference Mackle, Dwyer, Ingvartsen, Chouinard, Ross and Bauman2000), and therefore plasma isotopic enrichment should be a suitable indicator of the labelling of the mammary intracellular pool. There is evidence for ruminants that a proportion of total amino acids in circulating blood and plasma is bound in the form of peptides, and that there appears to be a removal of these by the mammary gland (Hanigan et al., Reference Hanigan, Calvert, De Peters, Reis and Baldwin1991). In the present model, any hydrolysis of peptides at the intracellular level would be masked in the estimate of the constitutive protein degradation (flows F 40 and F 50). The model relies on the assumption that enrichments have reached their steady state, and 6 h constant infusions of isotopically labelled amino acids have been used by others as well (e.g. Huang et al., Reference Huang, Yoder, Teixeira and Hanigan2021). However, the 6 h of infusion used in trials reported here might not have been sufficient to enable PHE and TYR present in milk protein to reach their true plateaux. This last hypothesis is supported by Bequette et al. (Reference Bequette, Backwell, Kyle, Calder, Buchan, Crompton, France and MacRae1999) who reported for goats that ~20 h of infusion were required for [1-13C] PHE in casein to effectively reach a true plateau.

Secreted milk contains a heterogeneous mixture of proteins. In the dairy cow, casein proteins comprise approximately 80% of total milk protein whilst the remainder is made up of various whey proteins (Miller et al., Reference Miller, Witherly and Clark1990). All casein proteins and some 70% of whey proteins are synthesized in the mammary gland. The remaining whey proteins are synthesized in the liver, transported to the mammary gland and then secreted in milk (Larson, Reference Larson1979). The present model omits this influx of preformed proteins to the mammary gland and assumes that all protein secreted in milk is synthesized in the mammary gland. Lapierre et al. (Reference Lapierre, Lobley, Doepel, Raggio, Rulquin and Lemosquet2012) estimated that some 3.4% of amino acids in milk protein may be contributed by blood-derived proteins. Following the synthesis of most milk proteins, a signal sequence on the newly synthesized protein is recognized by a specific recognition protein. The signal sequence and the specific recognition protein are cleaved during this process and presumably degraded intracellularly since they do not appear in secreted milk. As stated earlier, a value of 0.1 of total protein synthesis is ascribed in the present model calculations to allow for this retention and re-entry process. The process of milk protein synthesis from intracellular amino acids, and the reverse process of milk protein degradation, differ from transport of intracellular amino acids into milk to result in free amino acids in milk. The associated flows were therefore considered separately in constructing the model. However, given that amount of free amino acid in milk is generally low, for model application it is not strictly necessary to determine free PHE and TYR output in milk experimentally.

Despite its limitations, the model described provides a useful vehicle for obtaining information on the uptake and partitioning of PHE and TYR by the bovine mammary gland, indicating aspects of regulation that could be manipulated to direct more of the amino acid towards milk protein synthesis. Solving the model as the two four-pool scheme rather than in integrated 8-pool form is preferred as the equations are simpler and their application less susceptible to any compounding of measurement errors.

Concluding remarks

Solving the mammary model as two four-pool submodels rather than the integrated 8-pool scheme is preferred as the equations are slightly simpler and their application less susceptible to the compounding of measurement errors. The mammary model per se can be applied to other amino acids with similar metabolic fates within the tissue of study, thus used in research employing stable isotope techniques. In terms of practical usage, model solutions permit calculation of PHE and TYR flows for milk protein synthesis and constitutive protein turnover. The efficient use of amino acid N by dairy cows is of major practical importance given that amino acid surpluses are voided in urinary N compounds, in particular urea, which is highly susceptible to volatilization and associated environmental issues such as ammonia deposition and nitrous oxide (a greenhouse gas) formation. Furthermore, linking models describing the dynamics of several amino acids allows quantitative description of inter-organ amino acid metabolism, suggesting aspects of regulation that could be manipulated to direct more amino acids towards milk protein and therefore reduce urinary N excretion.

Financial support

This work was funded, in part, through DEFRA project LS3656 and the Canada Research Chairs Program. The experimental work was funded by a consortium of DEFRA, BBSRC, the Milk Development Council, Purina Mills LLC and NUTRECO Inc.

Authors' contributions

Conceptualization LAC, CKR and JF; methodology LAC, CKR, JD and JF; formal analysis LAC, LAM, JLE and JF; writing – original draft preparation LAC and JF; writing – review and editing LAC, LAM, CKR, JLE, JD and JF. All authors have read and agreed to the published version of the manuscript.

Conflict of interest

The authors declare there are no conflicts of interest.

Ethical standards

All experimental procedures used were licensed, regulated, and inspected by the UK Home Office under the Animals (Scientific Procedures) Act, 1986.

References

Baracos, V, Brun-Bellut, J and Marie, M (1991) Tissue protein synthesis in lactating and dry goats. British Journal of Nutrition 66, 451465.CrossRefGoogle ScholarPubMed
Bequette, BJ, Backwell, FRC, Kyle, CE, Calder, AG, Buchan, V, Crompton, LA, France, J and MacRae, JC (1999) Vascular sources of phenylalanine, tyrosine, lysine and methionine for casein synthesis in lactating goats. Journal of Dairy Science 82, 362377.CrossRefGoogle ScholarPubMed
Bequette, BJ, Kyle, CE, Crompton, LA, Anderson, SE and Hanigan, MD (2002) Protein metabolism in lactating goats subjected to the insulin clamp. Journal of Dairy Science 85, 15461555.CrossRefGoogle Scholar
Cant, JP, DePeters, EJ and Baldwin, RL (1993) Mammary amino acid utilization in dairy cows fed fat and its relationship to milk protein depression. Journal of Dairy Science 76, 762774.CrossRefGoogle ScholarPubMed
Crompton, LA, France, J, Reynolds, CK, Mills, JAN, Hanigan, MD, Ellis, JL and Dijkstra, J (2014) An isotope dilution model for partitioning phenylalanine and tyrosine uptake by the mammary gland of lactating dairy cows. Journal of Theoretical Biology 359, 5460.CrossRefGoogle ScholarPubMed
Dijkstra, J, Oenema, O, van Groenigen, JW, Spek, JW, van Vuuren, AM and Bannink, A (2013) Diet effects on urine composition of cattle and N2O emissions. Animal 7(Suppl. 2), 292302.CrossRefGoogle ScholarPubMed
France, J, Bequette, BJ, Lobley, GE, Metcalf, JA, Wray-Cahen, D, Dhanoa, MS, Backwell, FRC, Hanigan, MD, MacRae, JC and Beever, DE (1995) An isotope dilution for partitioning leucine uptake by the bovine mammary gland. Journal of Theoretical Biology 172, 369377.CrossRefGoogle Scholar
Hanigan, MD, Calvert, CC, De Peters, EJ, Reis, BL and Baldwin, RL (1991) Whole blood and plasma amino acid uptakes by lactating bovine mammary glands. Journal of Dairy Science 74, 24842490.CrossRefGoogle ScholarPubMed
Huang, X, Yoder, PS, Teixeira, IAMA and Hanigan, MD (2021) Assessing amino acid uptake and metabolism in mammary glands of lactating dairy cows intravenously infused with methionine, lysine, and histidine or with leucine and isoleucine. Journal of Dairy Science 104, 30323051.CrossRefGoogle ScholarPubMed
Jorgensen, GN and Larson, BL (1968) Conversion of phenylalanine to tyrosine in the bovine mammary secretory cell. Biochimica et Biophysica Acta 165, 121126.CrossRefGoogle ScholarPubMed
Lapierre, H, Lobley, GE, Doepel, L, Raggio, G, Rulquin, H and Lemosquet, S (2012) Mammary metabolism of amino acids in dairy cows. Journal of Animal Science 90, 17081721.CrossRefGoogle ScholarPubMed
Larson, BL (1979) Biosynthesis and secretion of milk proteins: a review. Journal of Dairy Research 46, 161174.CrossRefGoogle ScholarPubMed
Lemosquet, S, Lobley, GE, Koopman, R, Van Loon, LJC, Kies, AK and Lapierre, H (2010) A large supply of phenylalanine is not oxidized by the mammary gland of dairy cows. Energy and Protein Metabolism and Nutrition in Sustainable Animal Production (Ed GM Crovetto), EAAP Publication No. 134. Wageningen, the Netherlands: Wageningen Academic Publishers, pp. 137138.Google Scholar
Lobley, GE (2003) Protein turnover - what does it mean for animal production? Canadian Journal of Animal Science 83, 327340.CrossRefGoogle Scholar
Maas, JA, France, J and McBride, BW (1997) A mechanistic model of milk protein synthesis in the lactating bovine mammary gland. Journal of Theoretical Biology 187, 363378.CrossRefGoogle ScholarPubMed
Mackle, TR, Dwyer, DA, Ingvartsen, KL, Chouinard, PY, Ross, DA and Bauman, DE (2000) Effects of insulin and postruminal supply of protein on use of amino acids by the mammary gland for milk protein synthesis. Journal of Dairy Science 83, 93105.CrossRefGoogle ScholarPubMed
Mehaia, MA and Al-Kanhal, MA (1992) Taurine and other free amino acids in milk of camel, goat, cow and man. Milchwissenschaft 47, 351353.Google Scholar
Miller, MJS, Witherly, SA and Clark, DA (1990) Casein: a milk protein with diverse biologic consequences. Proceedings of the Society for Experimental Biology 195, 143159.CrossRefGoogle ScholarPubMed
National Academies of Science, Engineering and Medicine (NASEM) (2021) Nutrient Requirements of Dairy Cattle, 8th ed. Washington, DC, USA: National Academy Press, 502 pp.Google Scholar
Nichols, K, Bannink, A, Doelman, J and Dijkstra, J (2019a) Mammary gland metabolite utilization in response to exogenous glucose or long-chain fatty acids at low and high metabolizable protein levels. Journal of Dairy Science 102, 71507167.CrossRefGoogle Scholar
Nichols, K, Bannink, A and Dijkstra, J (2019b) Energy and nitrogen balance of dairy cattle as affected by provision of different essential amino acid profiles at the same metabolizable protein supply. Journal of Dairy Science 102, 89638976.CrossRefGoogle Scholar
Oddy, VH, Lindsay, DB and Fleet, LR (1988) Protein synthesis and degradation in the mammary gland of goats. Journal of Dairy Research 55, 143154.CrossRefGoogle ScholarPubMed
Raggio, G, Lemosquet, S, Lobley, GE, Rulquin, H and Lapierre, H (2006) Effect of casein and propionate supply on mammary protein metabolism in lactating dairy cows. Journal of Dairy Science 89, 43404351.CrossRefGoogle ScholarPubMed
Razooki Hasan, H, White, DA and Mayer, RJ (1982) Extensive destruction of newly synthesized casein in mammary explants in organ culture. Biochemical Journal 202, 133138.CrossRefGoogle ScholarPubMed
Roets, E, Massart-Leen, A-M, Verbeke, R and Peeters, G (1983) Metabolism of leucine by the isolated perfused goat udder. Journal of Dairy Research 50, 413424.CrossRefGoogle ScholarPubMed
Thomas, C (2004) Feed into Milk - A New Applied Feeding System for Dairy Cows. Nottingham, UK: Nottingham University Press, 68 pp.Google Scholar
Figure 0

Fig. 1. Scheme for the uptake and utilization of PHE and TYR by the mammary gland of lactating dairy cows as described by Crompton et al. (2014). The small circles indicate flows out of the system that need to be measured experimentally.

Figure 1

Fig. 2. Scheme for the uptake and utilization of PHE by the mammary gland of lactating dairy cows: (a) total PHE and (b) [13C] labelled PHE. The small circles in Fig. 2(a) indicate flows out of the system which need to be measured experimentally.

Figure 2

Table 1. Principle symbols used for the kinetic model

Figure 3

Fig. 3. Scheme for the uptake and utilization of TYR by the mammary gland of lactating dairy cows: (a) total TYR and (b) [2H] labelled TYR. The small circles in Fig. 3(a) indicate flows out of the system which need to be measured experimentally.

Figure 4

Table 2. Experimental and other inputs (symbols are defined in the text and Table 1)

Figure 5

Table 3. Phenylalanine and tyrosine uptake and partition by the mammary gland of four lactating dairy cows obtained using the two four-pool models (symbols are defined in the text and Table 1)

Figure 6

Table 4. Phenylalanine and tyrosine uptake and partition by the mammary gland of lactating dairy cows, obtained using the two four-pool models (symbols are defined in the text and Table 1) and the corresponding solutions obtained using the eight-pool model of Crompton et al. (2014)

Figure 7

Table 5. Average slope (%) for each of the flows calculated using the model described herein, obtained by perturbing each input in turn (symbols are defined in the text and Table 1)a