“…Initial CO positions in the active site of the wild type protein and the protein mutations were created manually in accordance with previous experimental and theoretical studies of NgbCO (see Figure 1). 30,34,109,114 As suggested in the literature, the dissociated CO preferentially adapts two different orientations, where either CO is located perpendicular to the plane defined by the distal histidine (form A) or parallel to that plane (form B). Furthermore, in each of the two docking forms A and B CO can be placed initially in two opposite orientations before the geometry optimization, one where the C atom is closer to the heme center (denoted in our study as 1) and one where the O atom is closer to the heme center (denoted as 2) which led to 36 protein systems denoted in this work as: F106A δ 1, F106A δ 2, F106A ϵ 1, F106A ϵ 2, F28L δ 1, F28L δ 2, F28L ϵ 1, F28L ϵ 2, F28W δ 1, F28W δ 2, F28W ϵ 1, F28W ϵ 2, H64A1, H64A2, H64Q1, H64Q2, H64V1, H64V2, H64 δ 1, H64 δ 2, H64 ϵ 1, H64 ϵ 2, H64 δϵ 1, H64 δϵ 2, K67A δ 1, K67A δ 2, K67A ϵ 1, K67A ϵ 2, K67R δ 1, K67R δ 2, K67R ϵ 1, K67R ϵ 2, K67T δ 1, K67T δ 2, K67T ϵ 1, and K67T ϵ2 , where 1 and 2 denote complexes obtained from different starting orientations.…”