LigPath: a module for predictive calculation of a ligand's pathway into a receptor-application to the gpH1-receptor.

Until now, the access of ligands into the binding pocket of a G-protein coupled receptor has scarcely been studied using molecular-modeling techniques because of the lack of sufficient algorithms. Neither with Monte-Carlo- nor with Molecular Dynamics Simulations can the penetration of a ligand into the binding pocket of a receptor be calculated because of the excessive amount of computing time needed. Therefore, a new algorithm LigPath for approximate calculation of a ligand's pathway into the binding pocket has been developed. This new algorithm is based on a linkage of directional guiding of the ligand, Monte-Carlo-Search and minimization. In order to evaluate the performance of the algorithm, the guinea-pig histamine H(1) receptor was investigated in combination with one of its potent agonists, histaprodifen, which is proposed to bind in a pocket deep between the transmembrane helices of the receptor. Our calculations show that the amino acids Tyr194, Phe193, Phe436 and Phe433 guide the positively charged histaprodifen from the extracellular part of the receptor into the binding pocket.