Analytical expressions are derived for the evaluation of energy gradients in the zeroth order regular approximation (ZORA) to the Dirac equation. The electrostatic shift approximation is used to avoid gauge dependence problems. Comparison is made to the quasirelativistic Pauli method, the limitations of which are highlighted. The structures and first metal-carbonyl bond dissociation energies for the transition metal complexes W(CO)(6), Os(CO)(5), and Pt(CO)(4) are calculated, and basis set effects are investigated.