Distribution and favorable binding sites of pyrroloquinoline and its analogues in a lipid bilayer studied by molecular dynamics simulations

The distribution of 1H-pyrrolo[3,2-h]quinoline (PQ), 11H-dipyrido[2,3-a]carbazole (PC) and 7-azaindole (7AI) at a water/membrane interface has been investigated by molecular dynamics (MD) simulations. The MD study focused on favorable binding sites of the azaaromatic probes across a dipalmitoylphosphatidylcholine (DPPC) bilayer. Our simulations show that PQ and PC are preferably accommodated at the hydrocarbon core of the bilayer below the glycerol moiety. In addition, it is found that the hydrophobic aromatic parts of the probes are located inside a more ordered region of DPPC, consisting of hydrophobic lipid chains. In contrast to PQ and PC, 7AI is characterized by a broad distribution between a DPPC interface and water, so that the three preferable binding sites are found across a water/membrane interface. It is found that in the sequence 7AI-PQ-PC, due to the increase of the number of aromatic rings and, hence, the hydrophobic character of the probes, the depth of the probe localization is gradually shifted deeper inside the hydrocarbon core of the bilayer. We found that the probe-lipid hydrogen-bonding contributes weakly to the favorable localizations of the azaaromatic probes inside the DPPC bilayer, so that the probe localization is mainly driven by electrostatic dipole-dipole and van der Waals interactions.