Ischaemic Hepatitis (IH) or Hypoxic Hepatitis (HH) also known as centrilobular liver cell necrosis is an acute liver injury characterized by a rapid increase in serum aminotransferase. The liver injury typically results from different underlying medical conditions such as cardiac failure, respiratory failure and septic shock in which the liver becomes damaged due to deprivation of either blood or oxygen. IH is a potentially lethal condition that is often preventable if diagnosed timely. The role of mechanisms that cause IH is often not well understood, making it difficult to diagnose or accurately quantify the patterns of related biomarkers. In most patients, currently, the only way to determine a case of IH is to rule out all other possible conditions for liver injuries. A better understanding of the liver's response to IH is necessary to aid in its diagnosis, measurement, and improve outcomes. The goal of this study is to identify mechanisms that can alter associated biomarkers for reducing the density of damaged hepatocytes, and thus reduce the chances of IH. We develop a mathematical model capturing dynamics of hepatocytes in the liver through the rise and fall of associated liver enzymes aspartate transaminase (AST), alanine transaminase (ALT) and lactate dehydrogenase (LDH) related to the condition of IH. The model analysis provides a novel approach to predict the level of biomarkers given variations in the systemic oxygen in the body. Using IH patient data in the US, novel model parameters are described and then estimated for the first time to capture real-time dynamics of hepatocytes in the presence and absence of IH condition. The results may allow physicians to estimate the extent of liver damage in an IH patient based on their enzyme levels and receive faster treatment on a real-time basis.