Sequential Application of Ligand and Structure Based Modeling Approaches to Index Chemicals for Their hH4R Antagonism

The human histamine H₄ receptor (hH₄R), a member of the G-protein coupled receptors (GPCR) family, is an increasingly attractive drug target. It plays a key role in many cell pathways and many hH₄R ligands are studied for the treatment of several inflammatory, allergic and autoimmune disorders, as well as for analgesic activity. Due to the challenging difficulties in the experimental elucidation of hH₄R structure, virtual screening campaigns are normally run on homology based models. However, a wealth of information about the chemical properties of GPCR ligands has also accumulated over the last few years and an appropriate combination of these ligand-based knowledge with structure-based molecular modeling studies emerges as a promising strategy for computer-assisted drug design. Here, two chemoinformatics techniques, the Intelligent Learning Engine (ILE) and Iterative Stochastic Elimination (ISE) approach, were used to index chemicals for their hH₄R bioactivity. An application of the prediction model on external test set composed of more than 160 hH₄R antagonists picked from the chEMBL database gave enrichment factor of 16.4. A virtual high throughput screening on ZINC database was carried out, picking ∼4000 chemicals highly indexed as H₄R antagonists'' candidates. Next, a series of 3D models of hH₄R were generated by molecular modeling and molecular dynamics simulations performed in fully atomistic lipid membranes. The efficacy of the hH₄R 3D models in discrimination between actives and non-actives were checked and the 3D model with the best performance was chosen for further docking studies performed on the focused library. The output of these docking studies was a consensus library of 11 highly active scored drug candidates. Our findings suggest that a sequential combination of ligand-based chemoinformatics approaches with structure-based ones has the potential to improve the success rate in discovering new biologically active GPCR drugs and increase the enrichment factors in a synergistic manner.     

Trio under threat: can we secure the future of rhinos,elephants and tigers in Malaysia?

Three of Malaysia’s endangered large mammal species are experiencing contrasting futures. Populations of the Sumatran rhino (Dicerorhinus sumatrensis) have dwindled to critically low numbers in Peninsular Malaysia (current estimates need to be revised) and the state of Sabah (less than 40 individuals estimated). In the latter region, a bold intervention involving the translocation of isolated rhinos is being developed to concentrate them into a protected area to improve reproduction success rates. For the Asian elephant (Elephas maximus), recently established baselines for Peninsular Malaysia (0.09 elephants/km² estimated from one site) and Sabah (between 0.56 and 2.15 elephants/km² estimated from four sites) seem to indicate globally significant populations based on dung count surveys. Similar surveys are required to monitor elephant population trends at these sites and to determine baselines elsewhere. The population status of the Malayan tiger (Panthera tigris jacksoni) in Peninsular Malaysia, however, remains uncertain as only a couple of scientifically defensible camera-trapping surveys (1.66 and 2.59 tigers/100 km² estimated from two sites) have been conducted to date. As conservation resources are limited, it may be prudent to focus tiger monitoring and protection efforts in priority areas identified by the National Tiger Action Plan for Malaysia. Apart from reviewing the conservation status of rhinos, elephants and tigers and threats facing them, we highlight existing and novel conservation initiatives, policies and frameworks that can help secure the long-term future of these iconic species in Malaysia.     

New vistas in GPCR 3D structure prediction

Human G-protein coupled receptors (hGPCRs) comprise the most prominent family of validated drug targets. More than 50% of approved drugs reveal their therapeutic effects by targeting this family. Accurate models would greatly facilitate the process of drug discovery and development. However, 3-D structure prediction of GPCRs remains a challenge due to limited availability of resolved structure. The X-ray structures have been solved for only four such proteins. The identity between hGPCRs and the potential templates is mostly less than 30%, well below the level at which sequence alignment can be done regularly. In this study, we analyze a large database of human G-protein coupled receptors that are members of family A in order to optimize usage of the available crystal structures for molecular modeling of hGPCRs. On the basis of our findings in this study, we propose to regard specific parts from the trans-membrane domains of the reference receptor helices as appropriate template for constructing models of other GPCRs, while other residues require other techniques for their remodeling and refinement. The proposed hypothesis in the current study has been tested by modeling human β2-adrenergic receptor based on crystal structures of bovine rhodopsin (1F88) and human A2A adenosine receptor (3EML). The results have shown some improvement in the quality of the predicted models compared to Modeller software.     

Reference-free detection of isolated SNPs

Detecting single nucleotide polymorphisms (SNPs) between genomes is becoming a routine task with next-generation sequencing. Generally, SNP detection methods use a reference genome. As non-model organisms are increasingly investigated, the need for reference-free methods has been amplified. Most of the existing reference-free methods have fundamental limitations: they can only call SNPs between exactly two datasets, and/or they require a prohibitive amount of computational resources. The method we propose, discoSnp, detects both heterozygous and homozygous isolated SNPs from any number of read datasets, without a reference genome, and with very low memory and time footprints (billions of reads can be analyzed with a standard desktop computer). To facilitate downstream genotyping analyses, discoSnp ranks predictions and outputs quality and coverage per allele. Compared to finding isolated SNPs using a state-of-the-art assembly and mapping approach, discoSnp requires significantly less computational resources, shows similar precision/recall values, and highly ranked predictions are less likely to be false positives. An experimental validation was conducted on an arthropod species (the tick Ixodes ricinus) on which de novo sequencing was performed. Among the predicted SNPs that were tested, 96% were successfully genotyped and truly exhibited polymorphism.     

Structure-toxicity relationships for different types of dinuclear platinum complexes

Nine structurally distinct dinuclear platinum complexes have been evaluated in a novel model system for the investigation of renal epithelial toxicity of platinum drugs. The results showed that these compounds are toxic when applied at the basolateral side of renal epithelia, whereas their toxic effects on the apical side are negligible. Such a difference in toxicity of the complexes has been found to result from their poor uptake through the apical membrane, as compared to the basolateral membrane. Toxicity of the compounds on the basolateral side varies depending on their structure. Structure-toxicity relationships for the group of complexes with rigid ligands and for the group of complexes with flexible ligands are discussed. Among the dinuclear complexes with rigid ligands, sterically hindered complexes are less toxic, due to their poor uptake and low reactivity towards glutathione. Within the group of complexes with flexible ligands, cis-configured isomers are more toxic than their trans-counterparts.     

A two-state model for Ca2+/CaM-dependent protein kinase II (alphaCaMKII) in response to persistent Ca2+ stimulation in hippocampal neurons

Persistent elevation of the intracellular free Ca(2+) concentration [Ca(2+)](i) is neurotoxic and therefore it is important to understand how it affects downstream components of the Ca(2+) signaling pathway. The response of calmodulin (CaM) and alphaCa(2+)/CaM-dependent protein kinase II (alphaCaMKII), to intracellular Ca(2+) overload in hippocampal neurons is studied by confocal imaging of fluorescently tagged proteins. Transient and persistent redistribution of CaM and alphaCaMKII together is seen from the cytosol to dendritic and somatic punctae. Typical persistent redistribution occurs following a lag of 138+/-(S.E.M.) 12 s and is complete at 460+/-(S.E.M.) 34 s (n=18), lack of Thr(286)-autophosphorylation of alphaCaMKII however promotes the formation of early transient punctae (peak at 40 s). In contrast, the T286D-mimick of phospho-Thr(286)-alphaCaMKII forms punctae with a delay >10 min, indicating that Thr(286)-autophosphorylation is antagonistic to CaMKII clustering. A two-state model is proposed in which phospho-Thr(286)-alphaCaMKII, formed immediately upon Ca(2+) stimulation, is primarily responsible for target interactions and memory functions of alphaCaMKII. However, a distinct clustering form denoted alphaCaMKII(c), generated upon persistent intracellular free Ca(2+) elevation, is deposited in the punctae which are made of self-interacting CaM/CaMKII complexes. Punctate deposition disables both the interactions and the activity of CaMKII.     

Structural basis for feedforward control in the PINK1/Parkin pathway

PINK1 and parkin constitute a mitochondrial quality control system mutated in Parkinson’s disease. PINK1, a kinase, phosphorylates ubiquitin to recruit parkin, an E3 ubiquitin ligase, to mitochondria. PINK1 controls both parkin localization and activity through phosphorylation of both ubiquitin and the ubiquitin‐like (Ubl) domain of parkin. Here, we observed that phospho‐ubiquitin can bind to two distinct sites on parkin, a high‐affinity site on RING1 that controls parkin localization and a low‐affinity site on RING0 that releases parkin autoinhibition. Surprisingly, ubiquitin vinyl sulfone assays, ITC, and NMR titrations showed that the RING0 site has higher affinity for phospho‐ubiquitin than phosphorylated Ubl in trans. We observed parkin activation by micromolar concentrations of tetra‐phospho‐ubiquitin chains that mimic mitochondria bearing multiple phosphorylated ubiquitins. A chimeric form of parkin with the Ubl domain replaced by ubiquitin was readily activated by PINK1 phosphorylation. In all cases, mutation of the binding site on RING0 abolished parkin activation. The feedforward mechanism of parkin activation confers robustness and rapidity to the PINK1‐parkin pathway and likely represents an intermediate step in its evolutionary development.