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

Molecular Modelling Approaches for the Analysis of Histamine Receptors and Their Interaction with Ligands

Several experimental techniques to analyse histamine receptors are available, e.g. pharmacological characterisation of known or new compounds by different types of assays or mutagenesis studies. To obtain insights into the histamine receptors on a molecular and structural level, crystal structures have to be determined and molecular modelling studies have to be performed. It is widely accepted to generate homology models of the receptor of interest based on an appropriate crystal structure as a template and to refine the resulting models by molecular dynamic simulations. A lot of modelling techniques, e.g. docking, QSAR or interaction fingerprint methods, are used to predict binding modes of ligands and pharmacological data, e.g. affinity or even efficacy. However, within the last years, molecular dynamic simulations got more and more important: First of all, molecular dynamic simulations are very helpful to refine the binding mode of a ligand to a histamine receptor, obtained by docking studies. Furthermore, with increasing computational performance it got possible to simulate complete binding pathways of ions or ligands from the aqueous extracellular phase into the allosteric or orthosteric binding pocket of histamine receptors.

Pathways and mechanisms for product release in the engineered haloalkane dehalogenases explored using classical and random acceleration molecular dynamics simulations

Eight mutants of the DhaA haloalkane dehalogenase carrying mutations at the residues lining two tunnels, previously observed by protein X-ray crystallography, were constructed and biochemically characterized. The mutants showed distinct catalytic efficiencies with the halogenated substrate 1,2,3-trichloropropane. Release pathways for the two dehalogenation products, 2,3-dichloropropane-1-ol and the chloride ion, and exchange pathways for water molecules, were studied using classical and random acceleration molecular dynamics simulations. Five different pathways, denoted p1, p2a, p2b, p2c, and p3, were identified. The individual pathways showed differing selectivity for the products: the chloride ion releases solely through p1, whereas the alcohol releases through all five pathways. Water molecules play a crucial role for release of both products by breakage of their hydrogen-bonding interactions with the active-site residues and shielding the charged chloride ion during its passage through a hydrophobic tunnel. Exchange of the chloride ions, the alcohol product, and the waters between the buried active site and the bulk solvent can be realized by three different mechanisms: (i) passage through a permanent tunnel, (ii) passage through a transient tunnel, and (iii) migration through a protein matrix. We demonstrate that the accessibility of the pathways and the mechanisms of ligand exchange were modified by mutations. Insertion of bulky aromatic residues in the tunnel corresponding to pathway p1 leads to reduced accessibility to the ligands and a change in mechanism of opening from permanent to transient. We propose that engineering the accessibility of tunnels and the mechanisms of ligand exchange is a powerful strategy for modification of the functional properties of enzymes with buried active sites.

Analysis of the activation mechanism of the guinea-pig Histamine H1-receptor

The Histamine H(1)-receptor (H1R), belonging to the amine receptor-class of family A of the G-protein coupled receptors (GPCRs) gets activated by agonists. The consequence is a conformational change of the receptor, which may involve the binding-pocket. So, for a good prediction of the binding-mode of an agonist, it is necessary to have knowledge about these conformational changes. Meanwhile some experimental data about the structural changes of GPCRs during activation exist. Based on homology modeling of the guinea-pig H1R (gpH1R), using the crystal structure of bovine rhodopsin as template, we performed several MD simulations with distance restraints in order to get an inactive and an active structure of the gpH1R. The calculations led to a Phe6.44/Trp6.48/Phe6.52-switch and linearization of the proline kinked transmembrane helix VI during receptor activation. Our calculations showed that the Trp6.48/Phe6.52-switch induces a conformational change in Phe6.44, which slides between transmembrane helices III and VI. Additionally we observed a hydrogen bond interaction of Ser3.39 with Asn7.45 in the inactive gpH1R, but because of a counterclockwise rotation of transmembrane helix III Ser3.39 establishes a water-mediated hydrogen bond to Asp2.50 in the active gpH1R. Additionally we simulated a possible mechanism for receptor activation with a modified LigPath-algorithm.