Metabolic control analysis (MCA) of pyruvate dehydrogenase multienzyme (PDH) complex of eucaryotic cells has been carried out using both in vitro and in vivo mechanistic models. Flux control coefficients (FCC) for the sensitivity of pyruvate decarboxylation rate to activities of various PDH complex reactions are determined. FCCs are shown to be strong functions of both pyruvate levels and various components of PDH complex. With the in vitro model, FCCs are shown to be sensitive to only the E1 component of the PDH complex at low pyruvate concentrations. At high pyruvate concentrations, the control is shared by all of the components, with E1 having a negative influence while the other three components, E2, X, and K, exert a positive control over the pyruvate decarboxylation rate. An unusual behavior of deactivation of the E1 component leading to higher net PDH activity is shown to be linked to the combined effect of protein X acylation and E1 deactivation. The steady-state analysis of the in vivo model reveals multiple steady state behavior of pyruvate metabolism with two stable and one unstable steady-states branches. FCCs also display multiplicity, showing completely different control distribution exerted by pyruvate and PDH components on three branches. At low pyruvate concentrations, pyruvate supply dominates the decarboxylation rate and PDH components do not exert any significant control. Reverse control distribution is observed at high pyruvate concentration. The effect of dilution due to cell growth on pyruvate metabolism is investigated in detail. While pyruvate dilution effects are shown to be negligible under all conditions, significant PDH complex dilution effects are observed under certain conditions. Comparison of in vitro and in vivo models shows that PDH components exert different degrees of control outside and inside the cells. At high pyruvate levels, PDH components are shown to exert a higher degree of control when reactions are taking place inside the cells as compared to the in vitro situation.