CryoEM structure and Alphafold molecular modelling of a novel molluscan hemocyanin.

Hemocyanins are multimeric oxygen transport proteins present in the blood of arthropods and molluscs, containing up to 8 oxygen-binding functional units per monomer. In molluscs, hemocyanins are assembled in decamer 'building blocks' formed of 5 dimer 'plates', routinely forming didecamer or higher-order assemblies with d5 or c5 symmetry. Here we describe the cryoEM structures of the didecamer (20-mer) and tridecamer (30-mer) forms of a novel hemocyanin from the slipper limpet Crepidula fornicata (SLH) at 7.0 and 4.7 Å resolution respectively. We show that two decamers assemble in a 'tail-tail' configuration, forming a partially capped cylinder, with an additional decamer adding on in 'head-tail' configuration to make the tridecamer. Analysis of SLH samples shows substantial heterogeneity, suggesting the presence of many higher-order multimers including tetra- and pentadecamers, formed by successive addition of decamers in head-tail configuration. Retrieval of sequence data for a full-length isoform of SLH enabled the use of Alphafold to produce a molecular model of SLH, which indicated the formation of dimer slabs with high similarity to those found in keyhole limpet hemocyanin. The fit of the molecular model to the cryoEM density was excellent, showing an overall structure where the final two functional units of the subunit (FU-g and FU-h) form the partial cap at one end of the decamer, and permitting analysis of the subunit interfaces governing the assembly of tail-tail and head-tail decamer interactions as well as potential sites for N-glycosylation. Our work contributes to the understanding of higher-order oligomer formation in molluscan hemocyanins and demonstrates the utility of Alphafold for building accurate structural models of large oligomeric proteins.

Identification of the molecular determinants of antagonist potency in the allosteric binding pocket of human P2X4.

P2X receptors are a family of ATP-gated cation channels comprising seven subtypes in mammals, which play key roles in nerve transmission, pain sensation and inflammation. The P2X4 receptor in particular has attracted significant interest from pharmaceutical companies due to its physiological roles in neuropathic pain and modulation of vascular tone. A number of potent small-molecule P2X4 receptor antagonists have been developed, including the allosteric P2X4 receptor antagonist BX430, which is approximately 30-fold more potent at human P2X4 compared with the rat isoform. A single amino-acid difference between human and rat P2X4 (I312T), located in an allosteric pocket, has previously been identified as critical for BX430 sensitivity, implying that BX430 binds in this pocket. Using a combination of mutagenesis, functional assay in mammalian cells and in silico docking we confirmed these findings. Induced-fit docking, permitting the sidechains of the amino-acids of P2X4 to move, showed that BX430 could access a deeper portion of the allosteric pocket, and that the sidechain of Lys-298 was important for shaping the cavity. We then performed blind docking of 12 additional P2X4 antagonists into the receptor extracellular domain, finding that many of these compounds favored the same pocket as BX430 from their calculated binding energies. Induced-fit docking of these compounds in the allosteric pocket enabled us to show that antagonists with high potency (IC50 ≤ 100 nM) bind deep in the allosteric pocket, disrupting a network of interacting amino acids including Asp-85, Ala-87, Asp-88, and Ala-297, which are vital for transmitting the conformational change following ATP binding to channel gating. Our work confirms the importance of Ile-312 for BX430 sensitivity, demonstrates that the allosteric pocket where BX430 binds is a plausible binding pocket for a series of P2X4 antagonists, and suggests a mode of action for these allosteric antagonists involving disruption of a key structural motif required for the conformational change induced in P2X4 when ATP binds.

Identification of Small Molecular Chaperones Binding P23H Mutant Opsin through an In Silico Structure-Based Approach.

N-terminal P23H opsin mutation accounts for most of retinitis pigmentosa (RP) cases. P23H functions and folding can be rescued by small chaperone ligands, which contributes to validate mutant opsin as a suitable target for pharmacological treatment of RP. However, the lack of structural details on P23H mutant opsin strongly impairs drug design, and new chemotypes of effective chaperones of P23H opsin are in high demand. Here, a computational-boosted workflow combining homology modeling with molecular dynamics (MD) simulations and virtual screening was used to select putative P23H opsin chaperones among different libraries through a structure-based approach. In vitro studies corroborated the reliability of the structural model generated in this work and identified a number of novel chemotypes of safe and effective chaperones able to promote P23H opsin trafficking to the outer cell membrane.

Identification of a novel P2X7 antagonist using structure-based virtual screening.

P2X4 and P2X7 receptors are ATP-gated ion channels, which play important roles in neuropathic and inflammatory pain, and as such they are important drug targets in diseases of inflammatory origin. While several compounds targeting P2X4 and P2X7 receptors have been developed using traditional high-throughput screening approaches, relatively few compounds have been developed using structure-based design. We initially set out to develop compounds targeting human P2X4, by performing virtual screening on the orthosteric (ATP-binding) pocket of a molecular model of human P2X4 based on the crystal structure of the Danio rerio receptor. The screening of a library of approximately 300,000 commercially available drug-like compounds led to the initial selection of 17 compounds; however, none of these compounds displayed a significant antagonist effect at P2X4 in a Fluo-4 ATP-induced calcium influx assay. When the same set of compounds was tested against human P2X7 in an ATP-stimulated Yo-Pro1 dye uptake assay, one compound (an indeno(1,2-b)pyridine derivative; GP-25) reduced the response by greater than 50% when applied at a concentration of 30 µM. GP-25 displayed an IC50 value of 8.7 µM at human P2X7 and 24.4 µM at rat P2X7, and was confirmed to be active using whole-cell patch clamp electrophysiology and not cytotoxic. Schild analysis suggested that mode of action of GP-25 was orthosteric. Screening of a further 16 commercially available analogues of GP-25 led to the discovery of five additional compounds with antagonist activity at human P2X7, enabling us to investigate the structure-activity relationship. Finally, docking of the R- and S-enantiomers of GP-25 into the orthosteric pocket of molecular models of human P2X4 and human P2X7 revealed that, while both enantiomers were able to make multiple interactions between their carboxyl moieties and conserved positively charged amino-acids in human P2X7, only the S-enantiomer of GP-25 was able to do this in human P2X4, potentially explaining the lack of activity of GP-25 at this receptor.

Ligand-based rational design, synthesis and evaluation of novel potential chemical chaperones for opsin.

Inherited blinding diseases retinitis pigmentosa (RP) and a subset of Leber's congenital amaurosis (LCA) are caused by the misfolding and mistrafficking of rhodopsin molecules, which aggregate and accumulate in the endoplasmic reticulum (ER), leading to photoreceptor cell death. One potential therapeutic strategy to prevent the loss of photoreceptors in these conditions is to identify opsin-binding compounds that act as chemical chaperones for opsin, aiding its proper folding and trafficking to the outer cell membrane. Aiming to identify novel compounds with such effect, a rational ligand-based approach was applied to the structure of the visual pigment chromophore, 11-cis-retinal, and its locked analogue 11-cis-6mr-retinal. Following molecular docking studies on the main chromophore binding site of rhodopsin, 49 novel compounds were synthesized according to optimized one-to seven-step synthetic routes. These agents were evaluated for their ability to compete for the chromophore binding site of opsin, and their capacity to increase the trafficking of the P23H opsin mutant from the ER to the cell membrane. Different new molecules displayed an effect in at least one assay, acting either as chemical chaperones or as stabilizers of the 9-cis-retinal-rhodopsin complex. These compounds could provide the basis to develop novel therapeutics for RP and LCA.

Computational Studies towards the Identification of Novel Rhodopsin-Binding Compounds as Chemical Chaperones for Misfolded Opsins.

Accumulation of misfolded and mistrafficked rhodopsin on the endoplasmic reticulum of photoreceptor cells has a pivotal role in the pathogenesis of retinitis pigmentosa and a subset of Leber's congenital amaurosis. One potential strategy to reduce rhodopsin misfolding and aggregation in these conditions is to use opsin-binding compounds as chemical chaperones for opsin. Such molecules have previously shown the ability to aid rhodopsin folding and proper trafficking to the outer cell membranes of photoreceptors. As means to identify novel chemical chaperones for rhodopsin, a structure-based virtual screening of commercially available drug-like compounds (300,000) was performed on the main binding site of the visual pigment chromophore, the 11-cis-retinal. The best 24 virtual hits were examined for their ability to compete for the chromophore-binding site of opsin. Among these, four small molecules demonstrated the ability to reduce the rate constant for the formation of the 9-cis-retinal-rhodopsin complex, while five molecules surprisingly enhanced the formation of this complex. Compound 7, 13, 20 and 23 showed a weak but detectable increase in the trafficking of the P23H mutant, widely used as a model for both retinitis pigmentosa and Leber's congenital amaurosis, from the ER to the cell membrane. The compounds did not show any relevant cytotoxicity in two different human cell lines, with the only exception of 13. Based on the structures of these active compounds, a series of in silico studies gave important insights on the potential structural features required for a molecule to act either as chemical chaperone or as stabiliser of the 11-cis-retinal-rhodopsin complex. Thus, this study revealed a series of small molecules that represent a solid foundation for the future development of novel therapeutics against these severe inherited blinding diseases.

Drosophila taste neurons as an agonist-screening platform for P2X receptors.

The P2X receptor family of ATP-gated cation channels are attractive drug targets for pain and inflammatory disease, but no subtype-selective agonists, and few partially selective agonists have been described to date. As proof-of-concept for the discovery of novel P2X receptor agonists, here we demonstrate the use of Drosophila taste neurons heterologously expressing rat P2X2 receptors as a screening platform. We demonstrate that wild-type rat P2X2 expressed in Drosophila is fully functional (ATP EC50 8.7 µM), and that screening of small (2 µl) volumes of a library of 80 adenosine nucleotide analogues is rapid and straightforward. We have determined agonist potency and specificity profiles for rat P2X2 receptors; triphosphate-bearing analogues display broad activity, tolerating a number of substitutions, and diphosphate and monophosphate analogues display very little activity. While several ATP analogues gave responses of similar magnitude to ATP, including the previously identified agonists ATPγS and ATPαS, we were also able to identify a novel agonist, the synthetic analogue 2-fluoro-ATP, and to confirm its agonist activity on rat P2X2 receptors expressed in human cells. These data validate our Drosophila platform as a useful tool for the analysis of agonist structure-activity relationships, and for the screening and discovery of novel P2X receptor agonists.

A new antiviral scaffold for human norovirus identified with computer-aided approaches on the viral polymerase.

Human norovirus is the leading cause of acute gastroenteritis worldwide, affecting every year 685 million people. In about one third of cases, this virus affects children under five years of age, causing each year up to 200,000 child deaths, mainly in the developing countries. Norovirus outbreaks are associated with very significant economic losses, with an estimated societal cost of 60 billion dollars per year. Despite the marked socio-economic consequences associated, no therapeutic options or vaccines are currently available to treat or prevent this infection. One promising target to identify new antiviral agents for norovirus is the viral polymerase, which has a pivotal role for the viral replication and lacks closely homologous structures in the host. Starting from the scaffold of a novel class of norovirus polymerase inhibitors recently discovered in our research group with a computer-aided method, different new chemical modifications were designed and carried out, with the aim to identify improved agents effective against norovirus replication in cell-based assays. While different new inhibitors of the viral polymerase were found, a further computer-aided ligand optimisation approach led to the identification of a new antiviral scaffold for norovirus, which inhibits human norovirus replication at low-micromolar concentrations.

The Molecular Determinants of Small-Molecule Ligand Binding at P2X Receptors.

P2X receptors are trimeric eukaryotic ATP-gated cation channels. Extracellular ATP-their physiological ligand-is released as a neurotransmitter and in conditions of cell damage such as inflammation, and substantial evidence implicates P2X receptors in diseases including neuropathic pain, cancer, and arthritis. In 2009, the first P2X crystal structure, Danio rerio P2X4 in the apo- state, was published, and this was followed in 2012 by the ATP-bound structure. These structures transformed our understanding of the conformational changes induced by ATP binding and the mechanism of ligand specificity, and enabled homology modeling of mammalian P2X receptors for ligand docking and rational design of receptor modulators. P2X receptors are attractive drug targets, and a wide array of potent, subtype-selective modulators (mostly antagonists) have been developed. In 2016, crystal structures of human P2X3 in complex with the competitive antagonists TNP-ATP and A-317491, and Ailuropoda melanoleuca P2X7 in complex with a series of allosteric antagonists were published, giving fascinating insights into the mechanism of channel antagonism. In this article we not only summarize current understanding of small-molecule modulator binding at P2X receptors, but also use this information in combination with previously published structure-function data and molecular docking experiments, to hypothesize a role for the dorsal fin loop region in differential ATP potency, and describe novel, testable binding conformations for both the semi-selective synthetic P2X7 agonist 2'-(3')-O-(4-benzoyl)benzoyl ATP (BzATP), and the P2X4-selective positive allosteric modulator ivermectin. We find that the distal benzoyl group of BzATP lies in close proximity to Lys-127, a residue previously implicated in BzATP binding to P2X7, potentially explaining the increased potency of BzATP at rat P2X7 receptors. We also present molecular docking of ivermectin to rat P2X4 receptors, illustrating a plausible binding conformation between the first and second transmembrane domains which not only tallies with previous mutagenesis studies, but would also likely have the effect of stabilizing the open channel structure, consistent with the mode of action of this positive allosteric modulator. From our docking simulations and analysis of sequence homology we propose a series of mutations likely to confer ivermectin sensitivity to human P2X1.

Haptic-driven, interactive drug design: implementing a GPU-based approach to evaluate the induced fit effect.

In this paper, we describe a hybrid meta-heuristics method of energy minimization and conformational sampling and its application into our haptic-driven molecular modelling simulator. The proposed method has been designed to suit real-time molecular docking simulations, where the time-lapse between two successive ligand poses is relatively small. In these situations, the energy minimization problem becomes increasingly complex and chaotic. The algorithm is tuned to take advantage of recent advances in GPU computing with asynchronous kernel execution, which has allowed us to include full protein flexibility in the real-time interactive, haptic-driven simulations. Finally, in this paper, we will also discuss the implementation of such high-performance computing approaches in our software, discussing the results of our initial validation studies, highlighting the advantages and limitations of such interactive methodology.