We present an investigation of the molecular basis of the modulation of oxygen affinity in heme proteins using computer simulation. QM-MM calculations are applied to explore distal and proximal effects on O(2) binding to the heme, while classical molecular dynamics simulations are employed to investigate ligand migration across the polypeptide to the active site. Trends in binding energies and in the kinetic constants are illustrated through a number of selected examples highlighting the virtues and the limitations of the applied methodologies. These examples cover a wide range of O(2)-affinities, and include: the truncated-N and truncated-O hemoglobins from Mycobacterium tuberculosis, the mammalian muscular O(2) storage protein: myoglobin, the hemoglobin from the parasitic nematode Ascaris lumbricoides, the oxygen transporter in the root of leguminous plants: leghemoglobin, the Cerebratulus lacteus nerve tissue hemoglobin, and the Alcaligenes xyloxidans cytochrome c'.