Predicting residue ionization of OmpF channel using Constant pH Molecular Dynamics as benchmarking

Electrostatic interactions play a key role in protein structure function. There is a large family of mesoscopic protein channels whose selectivity is mainly controlled by the protein electrostatic properties and ion specific channel interactions play a minor role. The knowledge of the charge state of the ionizable residues over a wide pH range, often summarized in their pKa, stands as the most valuable information for structure-function studies of many protein channels. However, experimental pKa determination is a difficult task, typically accomplished using Nuclear Magnetic Resonance only in a limited number of membrane proteins. Thus, the pKa calculation is the most frequently used alternative. Constant pH Molecular Dynamics (CpHMD) simulation provides arguably the most accurate pKa prediction method in proteins containing many charged residues since it captures the coupling between conformational dynamics and residue protonation. Here we study the charge state of a general diffusion porin, OmpF, in which protons exert a crucial regulation of the channel discrimination of small inorganic ions as well as antibiotic translocation. We examine the pKa prediction using different methods, with the CpHMD simulations as benchmarking, and discuss the somewhat unusual titration of several acidic residues. The most widely used pKa prediction methods, though useful for globular proteins, fail to capture the specificities of channel proteins embedded in biological membranes. This is the first attempt we know to use CpHMD to study the pH- dependent charge of a large multiionic channel (with over three hundred ionizable residues) embedded in a lipid membrane.