*These authors contributed equally to this workTruncated hemoglobins build one of the three branches of the globin protein superfamily. They display a characteristic two-on-two a-helical sandwich fold and are clustered into three groups (I, II and III) based on distinct structural features. Truncated hemoglobins are present in eubacteria, cyanobacteria, protozoa and plants. Here we present a structural, spectroscopic and molecular dynamics characterization of a group-II truncated hemoglobin, encoded by the PSHAa0030 gene from Pseudoalteromonas haloplanktis TAC125 (Ph-2/2HbO), a cold-adapted Antarctic marine bacterium hosting one flavohemoglobin and three distinct truncated hemoglobins. The Ph-2/2HbO aquo-met crystal structure (at 2.21A resolution) shows typical features of group-II truncated hemoglobins, namely the twoon-two a-helical sandwich fold, a helix Φ preceding the proximal helix F, and a heme distal-site hydrogen-bonded network that includes water molecules and several distal-site residues, including His(58)CD1. Analysis of Ph-2/2HbO by electron paramagnetic resonance, resonance Raman and electronic absorption spectra, under varied solution conditions, shows that Ph-2/2HbO can access diverse heme ligation states. Among these, detection of a low-spin heme hexa-coordinated species suggests that residue Tyr(42)B10 can undergo large conformational changes in order to act as the sixth heme-Fe ligand. Altogether, the results show that Ph-2/2HbO maintains the general structural features of group-II truncated hemoglobins but displays enhanced conformational flexibility in the proximity of the heme cavity, a property probably related to the functional challenges, such as low Abbreviations 5c, penta-coordinated heme state; 6c, hexa-coordinated heme state; At-2/2HbO, TrHbII from Agrobacterium tumefaciens; Ath-2/2HbO,

Predicting function from sequence is an important goal in current biological research, and although, broad functional assignment is possible when a protein is assigned to a family, predicting functional specificity with accuracy is not straightforward. If function is provided by key structural properties and the relevant properties can be computed using the sequence as the starting point, it should in principle be possible to predict function in detail. The truncated hemoglobin family presents an interesting benchmark study due to their ubiquity, sequence diversity in the context of a conserved fold and the number of characterized members. Their functions are tightly related to O2 affinity and reactivity, as determined by the association and dissociation rate constants, both of which can be predicted and analyzed using in-silico based tools. In the present work we have applied a strategy, which combines homology modeling with molecular based energy calculations, to predict and analyze function of all known truncated hemoglobins in an evolutionary context. Our results show that truncated hemoglobins present conserved family features, but that its structure is flexible enough to allow the switch from high to low affinity in a few evolutionary steps. Most proteins display moderate to high oxygen affinities and multiple ligand migration paths, which, besides some minor trends, show heterogeneous distributions throughout the phylogenetic tree, again suggesting fast functional adaptation. Our data not only deepens our comprehension of the structural basis governing ligand affinity, but they also highlight some interesting functional evolutionary trends.