Comparative analysis of inner cavities and ligand migration in non-symbiotic AHb1 and AHb2

This study reports a comparative analysis of the topological properties of inner cavities and the intrinsic dynamics of non-symbiotic hemoglobins AHb1 and AHb2 from Arabidopsis thaliana. The two proteins belong to the 3/3 globin fold and have a sequence identity of about 60%. However, it is widely assumed that they have distinct physiological roles. In order to investigate the structure–function relationships in these proteins, we have examined the bis-histidyl and ligand-bound hexacoordinated states by atomistic simulations using in silico structural models. The results allow us to identify two main pathways to the distal cavity in the bis-histidyl hexacoordinated proteins. Nevertheless, a larger accessibility to small gaseous molecules is found in AHb2. This effect can be attributed to three factors: the mutation Leu35(AHb1) → Phe32(AHb2), the enhanced flexibility of helix B, and the more favorable energetic profile for ligand migration to the distal cavity. The net effect of these factors would be to facilitate the access of ligands, thus compensating the preference for the fully hexacoordination of AHb2, in contrast to the equilibrium between hexa- and pentacoordinated species in AHb1. On the other hand, binding of the exogenous ligand introduces distinct structural changes in the two proteins. A well-defined tunnel is formed in AHb1, which might be relevant to accomplish the proposed NO detoxification reaction. In contrast, no similar tunnel is found in AHb2, which can be ascribed to the reduced flexibility of helix E imposed by the larger number of salt bridges compared to AHb1. This feature would thus support the storage and transport functions proposed for AHb2. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.