Three-dimensional lipophilicity characterization of molecular pores and channel-like cavities

Molecular lipophilicity is a useful property for assessing molecular similarity or complementarity within the context of computer-aided drug design. As well, local contributions to solvent affinity help us to understand both dynamics and conformational stability in biomolecules. In this work, we discuss an approach to characterize the local contributions to hydrophobicity by using one- and two-dimensional representations of molecular channel-like cavities. The method monitors how a phenomenological lipophilicity potential (based on fragmental atom contributions) changes over a continuum of "molecular tubes" used for modeling channels and pores. Our results convey a relatively detailed picture of the spatial distribution of water affinity. The procedure can then be used as a complement to the hydrophobicity scales based on averaging contributions from single amino acids. In addition, we can study how the water affinity changes for inner and outer regions of the pores. As an application, we compute the 3D distribution of lipophilicity in the "pore conformation" of gramicidin A. The qualitative trends indicated by our results are broadly consistent with computer simulations of the gramicidin channel in the presence of hydrated ions. The behavior revealed by the simulations can then be incorporated to produce an improved, simple 2D model for water-channel interactions.