Assessing Structural Determinants of Zn²⁺ Binding to Human H_V1 via Multiple MD Simulations

Voltage-gated proton channels (H_V1) are essential for various physiological tasks but are strongly inhibited by Zn²⁺ cations. Some determinants of Zn²⁺ binding have been elucidated experimentally and in computational studies. However, the results have always been interpreted under the assumption that Zn²⁺ binds to monomeric H_V1 despite evidence that H_V1 expresses as a dimer and that the dimer has a higher affinity for zinc than the monomer and experimental data that suggest coordination in the dimer interface. The results of former studies are also controversial, e.g., supporting either one single or two binding sites. Some structural determinants of the binding are still elusive. We performed a series of molecular dynamics simulations to address different structures of the human proton channel, the monomer and two plausible dimer conformations, to compare their respective potential to interact with and bind Zn²⁺ via the essential histidines. The series consisted of several copies of the system to generate independent trajectories and increase the significance compared to a single simulation. The amount of time simulated totals 29.9 μs for 126 simulations of systems comprising ∼59,000 to ∼187,000 atoms. Our approach confirms the existence of two binding sites in monomeric and dimeric human H_V1. The dimer interface is more efficient for attracting and binding Zn²⁺ via the essential histidines than the monomer or a dimer with the histidines in the periphery. The higher affinity is due to the residues in the dimer interface that create an attractive electrostatic potential funneling the zinc cations toward the binding sites.