Open Access

Our model of the bovine estrous cycle is a set of ordinary differential equations which generates hormone profiles of successive estrous cycles with several follicular waves per cycle. It describes the growth and decay of the follicles and the corpus luteum, as well as the change of the key substances over time. In this work we describe recent improvements of this model, including the introduction of new components, and elimination of time delays. We validate our model by showing that the simulations agree with observations from synchronization studies and with measured progesterone data after a single dose administration of synthetic prostaglandin F2alpha.

High-performing dairy cows require a particular composition of nutritional ingredients, adapted to their individual requirements and depending on their production status. The optimal dimensioning of minerals in the diet, one being potassium, is indispensable for the prevention of imbalances. Potassium balance in cows is the result of potassium intake, distribution in the organism, and excretion, and it is closely related to glucose and electrolyte metabolism. In this paper, we present a dynamical model for potassium balance in lactating and nonlactating dairy cows based on ordinary differential equations. Parameter values were obtained from clinical trial data and from the literature. To verify the consistency of the model, we present simulation outcomes for 3 different scenarios: potassium balance in (1) nonlactating cows with varying feed intake, (2) nonlactating cows with varying potassium fraction in the diet, and (3) lactating cows with varying milk production levels. The results give insights into the short- and long-term potassium metabolism, providing an important step toward the understanding of the potassium network, the design of prophylactic feed additives, and possible treatment strategies.

The estrous cycle of mono-ovulatory species such as cows or humans, is known to show two or more waves of follicular growth and decline between two successive ovulations. Within each wave, there is one dominant follicle escorted by subordinate follicles of varying number. Under the surge of the luteinizing hormone a growing dominant follicle ovulates. Rarely the number of ovulating follicles exceeds one. In the biological literature, the change of hormonal concentrations and individually varying numbers of follicular receptors are made responsible for the selection of exactly one dominant follicle, yet a clear cause has not been identified. In this paper, we suggest a synergistic explanation based on competition, formulated by a parsimoniously defined system of ordinary differential equations (ODEs) that quantifies the time evolution of multiple follicles and their competitive interaction during one wave. Not discriminating between follicles, growth and decline are given by fixed rates. Competition is introduced via a growth-suppressing term, equally supported by all follicles. We prove that the number of dominant follicles is determined exclusively by the ratio of follicular growth and competition. This number turns out to be independent of the number of subordinate follicles. The asymptotic behavior of the corresponding dynamical system is investigated rigorously, where we demonstrate that the omega-limit set only contains fixed points. When also including follicular decline, our ODEs perfectly resemble ultrasound data of bovine follicles. Implications for the involved but not explicitly modeled hormones are discussed.