Scheff JD, Mavroudis PD, Calvano SE, Lowry SF, Androulakis IP. Modeling autonomic regulation of cardiac function and heart rate variability in human endotoxemia. Physiol Genomics 43: 951-964, 2011. First published June 14, 2011 doi:10.1152/physiolgenomics.00040.2011.-Heart rate variability (HRV), the quantification of beat-to-beat variability, has been studied as a potential prognostic marker in inflammatory diseases such as sepsis. HRV normally reflects significant levels of variability in homeostasis, which can be lost under stress. Much effort has been placed in interpreting HRV from the perspective of quantitatively understanding how stressors alter HRV dynamics, but the molecular and cellular mechanisms that give rise to both homeostatic HRV and changes in HRV have received less focus. Here, we develop a mathematical model of human endotoxemia that incorporates the oscillatory signals giving rise to HRV and their signal transduction to the heart. Connections between processes at the cellular, molecular, and neural levels are quantitatively linked to HRV. Rhythmic signals representing autonomic oscillations and circadian rhythms converge to modulate the pattern of heartbeats, and the effects of these oscillators are diminished in the acute endotoxemia response. Based on the semimechanistic model developed herein, homeostatic and acute stress responses of HRV are studied in terms of these oscillatory signals. Understanding the loss of HRV in endotoxemia serves as a step toward understanding changes in HRV observed clinically through translational applications of systems biology based on the relationship between biological processes and clinical outcomes. systems biology; inflammation; circadian HEART RATE VARIABILITY (HRV) is generally defined as the quantification of the distribution of time intervals between successive heartbeats. Reduction in HRV, a manifestation of altered autonomic function under stress, is potentially a useful predictor of outcome in myocardial infarction (42), congestive heart failure (63), diabetic neuropathy (58), and neonatal sepsis (46). Diminished HRV has also been observed in critically ill patients in intensive care units (ICU) (52), which motivates interest in HRV as a critical variable in the recovery from critical illness (50). Due to this clinical relevance, dynamic characteristics of HRV have been assessed by time domain, frequency domain (76), and nonlinear metrics (46, 60). The majority of HRV research has thus far focused on the interpretation of the patterns of HRV (45) rather than linking cellular-level mechanisms to patterns (11). The realization that health may be characterized by a certain degree of variability of human heart signals motivates the hypothesis that appropriate physiological variability is the manifestation of robust dynamics of control signals whose fluctuations equip the host with the ability to anticipate external and internal disturbances. We hypothesize that these variable dynamics are driven by the convergence of rhythmic physiological signals on...