A mathematical model is proposed for predicting the copolymer composition of the microbially synthesized polyhydroxyalkanoate (PHA) copolymers. Based on the biochemical reactions involved in the precursor formation and polymerization pathways, the model correlates the copolymer composition with the cultivation conditions, the enzyme levels and selectivity, and the metabolic pathways. It suggests the following points: (1) in the case of a sole carbon source, the copolymer composition depends mainly on the topology of the metabolic pathways and the selectivity of both the enzymes involved in the precursor formation and the polymerization route; (2) the copolymer composition can be varied in a wide range via alteration of the flux ratio of different types of monomers channeled from two or more independent and simultaneous pathways; (3) the enzymes which should be over-expressed or inhibited to obtain the desired copolymer composition can be predicted. For example, inhibition of the b-oxidation pathway will increase the content of themonomerunits with longer chain length. To test the model, various experiments were envisaged by varying cultivation time, concentration and chain length of the sole carbon source, and molar ratio of the cosubstrates. The predictions from the model agree well with the experimental results. Therefore, the proposed model will be useful in predicting the PHA copolymer composition under different biochemical reaction conditions. In other words, it can provide a guide for the synthesis of desired PHA copolymers.