The enzyme carbon monoxide dehydrogenase is capable of efficiently converting CO 2 to CO and therefore can enable an affordable CO 2 recycling strategy. The reduction of CO 2 occurs at a peculiar nickel-ironsulfur cluster, following a mechanism that remains little understood. In this study, we have used ab initio molecular dynamics simulations to explore the free energy landscape of the reaction. We predict the existence of a COOH ligand that strongly interacts with the surrounding protein residues and favours a mechanism where a H 2 O molecule is eliminated before CO. We have taken advantages of the insights offered by our simulations to revisit the catalytic mechanism and the role of the residues surrounding the active centre in particular, thus assisting in the design of inorganic catalysts that mimic the enzyme.