Application of Free Energy Perturbation (FEP+) to Understanding Ligand Selectivity: A Case Study to Assess Selectivity Between Pairs of Phosphodiesterases (PDE's).

Cyclic nucleotide phosphodiesterases (PDE's) are metalloenzymes that play a key role in regulating the levels of the ubiquitous second messengers, cyclic adenosine monophosphate (cAMP) and cyclic guanosine monophosphate (cGMP). In humans, 11 PDE protein families mediate numerous biochemical pathways throughout the body and are effective drug targets for the treatment of diseases ranging from central nervous system disorders to heart and pulmonary diseases. PDE's also share a highly conserved catalytic site (about 50%), thus making the design of selective drug candidates very challenging with classical structure-based design approaches given also the lack of publicly available co-crystal structures of pairs of PDE's with consistent biological data to be compared, as we show in our work. In this retrospective study, we apply free energy perturbation (FEP+) to predict the selectivity of inhibitors that bind two pairs of closely related PDE families: PDE9/1 and PDE5/6 where only 1 co-crystal structure per pair is publicly available. As another challenge, the p Ka of the PDE5/6 inhibitor is close to the experimental pH, making unclear the exact protonation state that should be used in the computational workflow. We demonstrate that running FEP+ on homology models constructed for these metalloenzymes accurately reproduces experimentally observed selectivity profiles also addressing the unclear protonation state to be used during computation with our recently developed p Ka-correction method. Based on these data, we conclude that FEP+ is a robust method for prediction of selectivity for this target class and may be helpful to address related lead optimization challenges in drug discovery.

αIIbß3 variants defined by next-generation sequencing: predicting variants likely to cause Glanzmann thrombasthenia.

Next-generation sequencing is transforming our understanding of human genetic variation but assessing the functional impact of novel variants presents challenges. We analyzed missense variants in the integrin αIIbß3 receptor subunit genes ITGA2B and ITGB3 identified by whole-exome or -genome sequencing in the ThromboGenomics project, comprising ∼32,000 alleles from 16,108 individuals. We analyzed the results in comparison with 111 missense variants in these genes previously reported as being associated with Glanzmann thrombasthenia (GT), 20 associated with alloimmune thrombocytopenia, and 5 associated with aniso/macrothrombocytopenia. We identified 114 novel missense variants in ITGA2B (affecting ∼11% of the amino acids) and 68 novel missense variants in ITGB3 (affecting ∼9% of the amino acids). Of the variants, 96% had minor allele frequencies (MAF) < 0.1%, indicating their rarity. Based on sequence conservation, MAF, and location on a complete model of αIIbß3, we selected three novel variants that affect amino acids previously associated with GT for expression in HEK293 cells. αIIb P176H and ß3 C547G severely reduced αIIbß3 expression, whereas αIIb P943A partially reduced αIIbß3 expression and had no effect on fibrinogen binding. We used receiver operating characteristic curves of combined annotation-dependent depletion, Polyphen 2-HDIV, and sorting intolerant from tolerant to estimate the percentage of novel variants likely to be deleterious. At optimal cut-off values, which had 69-98% sensitivity in detecting GT mutations, between 27% and 71% of the novel αIIb or ß3 missense variants were predicted to be deleterious. Our data have implications for understanding the evolutionary pressure on αIIbß3 and highlight the challenges in predicting the clinical significance of novel missense variants.

Identification of a µ-δ opioid receptor heteromer-biased agonist with antinociceptive activity.

G protein-coupled receptors play a pivotal role in many physiological signaling pathways. Mounting evidence suggests that G protein-coupled receptors, including opioid receptors, form dimers, and dimerization is necessary for receptor maturation, signaling, and trafficking. However, the physiological role of dimerization in vivo has not been well-explored because of the lack of tools to study these dimers in endogenous systems. To address this problem, we previously generated antibodies to µ-δ opioid receptor (µOR-δOR) dimers and used them to study the pharmacology and signaling by this heteromer. We also showed that the heteromer exhibits restricted distribution in the brain and that its abundance is increased in response to chronic morphine administration. Thus, the µOR-δOR heteromer represents a potentially unique target for the development of therapeutics to treat pain. Here, we report the identification of compounds targeting µOR-δOR heteromers through high-throughput screening of a small-molecule library. These compounds exhibit activity in µOR-δOR cells but not µOR or δOR cells alone. Among them, CYM51010 was found to be a µOR-δOR-biased ligand, because its activity is blocked by the µOR-δOR heteromer antibody. Notably, systemic administration of CYM51010 induced antinociceptive activity similar to morphine, and chronic administration of CYM51010 resulted in lesser antinociceptive tolerance compared with morphine. Taken together, these results suggest that CYM51010, a µOR-δOR-biased ligand, could serve as a scaffold for the development of a unique type (heteromer-biased) of drug that is more potent and without the severe side effects associated with conventional clinical opioids.

RUC-4: a novel αIIbß3 antagonist for prehospital therapy of myocardial infarction.

OBJECTIVE: Treatment of myocardial infarction within the first 1 to 2 hours with a thrombolytic agent, percutaneous coronary intervention, or an αIIbß3 antagonist decreases mortality and the later development of heart failure. We previously reported on a novel small molecule αIIbß3 antagonist, RUC-2, that has a unique mechanism of action. We have now developed a more potent and more soluble congener of RUC-2, RUC-4, designed to be easily administered intramuscularly by autoinjector to facilitate its use in the prehospital setting. Here, we report the properties of RUC-4 and the antiplatelet and antithrombotic effects of RUC-2 and RUC-4 in animal models. APPROACH AND RESULTS: RUC-4 was ≈ 20% more potent than RUC-2 in inhibiting human ADP-induced platelet aggregation and much more soluble in aqueous solutions (60-80 mg/mL). It shared RUC-2's specificity for αIIbß3 versus αVß3, did not prime the receptor to bind fibrinogen, or induce changes in ß3 identified by a conformation-specific monoclonal antibody. Both RUC-2 and RUC-4 prevented FeCl3-induced thrombotic occlusion of the carotid artery in mice and decreased microvascular thrombi in response to laser injury produced by human platelets infused into transgenic mice containing a mutated von Willebrand factor that reacts with human but not mouse platelets. Intramuscular injection of RUC-4 in nonhuman primates at 1.9 and 3.85 mg/kg led to complete inhibition of platelet aggregation within 15 minutes, with dose-dependent return of platelet aggregation after 4.5 to 24 hours. CONCLUSIONS: RUC-4 has favorable biochemical, pharmacokinetic, pharmacodynamic, antithrombotic, and solubility properties as a prehospital therapy of myocardial infarction, but the possibility of increased bleeding with therapeutic doses remains to be evaluated.

Talin-driven inside-out activation mechanism of platelet αIIbß3 integrin probed by multimicrosecond, all-atom molecular dynamics simulations.

Platelet aggregation is the consequence of the binding of extracellular bivalent ligands such as fibrinogen and von Willebrand factor to the high affinity, active state of integrin αIIbß3. This state is achieved through a so-called "inside-out" mechanism characterized by the membrane-assisted formation of a complex between the F2 and F3 subdomains of intracellular protein talin and the integrin ß3 tail. Here, we present the results of multi-microsecond, all-atom molecular dynamics simulations carried on the complete transmembrane (TM) and C-terminal (CT) domains of αIIbß3 integrin in an explicit lipid-water environment, and in the presence or absence of the talin-1 F2 and F3 subdomains. These large-scale simulations provide unprecedented molecular-level insights into the talin-driven inside-out activation of αIIbß3 integrin. Specifically, they suggest a preferred conformation of the complete αIIbß3 TM/CT domains in a lipid-water environment, and testable hypotheses of key intermolecular interactions between αIIbß3 integrin and the F2/F3 domains of talin-1. Notably, not only do these simulations give support to a stable left-handed reverse turn conformation of the αIIb juxtamembrane motif rather than a helical turn, but they raise the question as to whether TM helix separation is required for talin-driven integrin activation.

Bioactive conformations of two seminal delta opioid receptor penta-peptides inferred from free-energy profiles.

Delta-opioid (DOP) receptors are members of the G protein-coupled receptor (GPCR) sub-family of opioid receptors, and are evolutionarily related, with homology exceeding 70%, to cognate mu-opioid (MOP), kappa-opioid (KOP), and nociceptin opioid (NOP) receptors. DOP receptors are considered attractive drug targets for pain management because agonists at these receptors are reported to exhibit strong antinociceptive activity with relatively few side effects. Among the most potent analgesics targeting the DOP receptor are the linear and cyclic enkephalin analogs known as DADLE (Tyr-D-Ala-Gly-Phe-D-Leu) and DPDPE (Tyr-D-Pen-Gly-Phe-D-Pen), respectively. Several computational and experimental studies have been carried out over the years to characterize the conformational profile of these penta-peptides with the ultimate goal of designing potent peptidomimetic agonists for the DOP receptor. The computational studies published to date, however, have investigated only a limited range of timescales and used over-simplified representations of the solvent environment. We provide here a thorough exploration of the conformational space of DADLE and DPDPE in an explicit solvent, using microsecond-scale molecular dynamics and bias-exchange metadynamics simulations. Free-energy profiles derived from these simulations point to a small number of DADLE and DPDPE conformational minima in solution, which are separated by relatively small energy barriers. Candidate bioactive forms of these peptides are selected from identified common spatial arrangements of key pharmacophoric points within all sampled conformations.

A novel class of ion displacement ligands as antagonists of the αIIbß3 receptor that limit conformational reorganization of the receptor.

A collection of αIIbß3 integrin receptor antagonists possessing a unique MIDAS metal ion displacement mechanism of action is presented. Insight into these agents' structure-activity relationships, binding modality, and pharmacokinetic and pharmacodynamic profiles highlight the potential of these small molecule ion displacement ligands as attractive candidates for clinical development.

Probing the dimerization interface of Leishmania infantum trypanothione reductase with site-directed mutagenesis and short peptides.

Binding at the interface: We tested the inhibitory activity of a set of peptide sequences derived from an α-helix of the dimeric trypanothione reductase from Leishmania infantum. Replacement of a glutamic acid residue with a lysine promoted monomer dissociation and enzyme inhibition.

Discovery of a novel selective kappa-opioid receptor agonist using crystal structure-based virtual screening.

Kappa-opioid (KOP) receptor agonists exhibit analgesic effects without activating reward pathways. In the search for nonaddictive opioid therapeutics and novel chemical tools to study physiological functions regulated by the KOP receptor, we screened in silico its recently released inactive crystal structure. A selective novel KOP receptor agonist emerged as a notable result and is proposed as a new chemotype for the study of the KOP receptor in the etiology of drug addiction, depression, and/or pain.

Closed headpiece of integrin αIIbß3 and its complex with an αIIbß3-specific antagonist that does not induce opening.

The platelet integrin α(IIb)ß(3) is essential for hemostasis and thrombosis through its binding of adhesive plasma proteins. We have determined crystal structures of the α(IIb)ß(3) headpiece in the absence of ligand and after soaking in RUC-1, a novel small molecule antagonist. In the absence of ligand, the α(IIb)ß(3) headpiece is in a closed conformation, distinct from the open conformation visualized in presence of Arg-Gly-Asp (RGD) antagonists. In contrast to RGD antagonists, RUC-1 binds only to the α(IIb) subunit. Molecular dynamics revealed nearly identical binding. Two species-specific residues, α(IIb) Y190 and α(IIb) D232, in the RUC-1 binding site were confirmed as important by mutagenesis. In sharp contrast to RGD-based antagonists, RUC-1 did not induce α(IIb)ß(3) to adopt an open conformation, as determined by gel filtration and dynamic light scattering. These studies provide insights into the factors that regulate integrin headpiece opening, and demonstrate the molecular basis for a novel mechanism of integrin antagonism.

Ligand-induced modulation of the free-energy landscape of G protein-coupled receptors explored by adaptive biasing techniques.

Extensive experimental information supports the formation of ligand-specific conformations of G protein-coupled receptors (GPCRs) as a possible molecular basis for their functional selectivity for signaling pathways. Taking advantage of the recently published inactive and active crystal structures of GPCRs, we have implemented an all-atom computational strategy that combines different adaptive biasing techniques to identify ligand-specific conformations along pre-determined activation pathways. Using the prototypic GPCR ß2-adrenergic receptor as a suitable test case for validation, we show that ligands with different efficacies (either inverse agonists, neutral antagonists, or agonists) modulate the free-energy landscape of the receptor by shifting the conformational equilibrium towards active or inactive conformations depending on their elicited physiological response. Notably, we provide for the first time a quantitative description of the thermodynamics of the receptor in an explicit atomistic environment, which accounts for the receptor basal activity and the stabilization of different active-like states by differently potent agonists. Structural inspection of these metastable states reveals unique conformations of the receptor that may have been difficult to retrieve experimentally.

Rationale for the opposite stereochemistry of the major monoadducts and interstrand crosslinks formed by mitomycin C and its decarbamoylated analogue at CpG steps in DNA and the effect of cytosine modification on reactivity.

Mitomycin C (MMC) is a potent antitumour agent that forms a covalent bond with the 2-amino group of selected guanines in the minor groove of double-stranded DNA following intracellular reduction of its quinone ring and opening of its aziridine moiety. At some 5'-CG-3' (CpG) steps the resulting monofunctional adduct can evolve towards a more deleterious bifunctional lesion, which is known as an interstrand crosslink (ICL). MMC reactivity is enhanced when the cytosine bases are methylated (5 MC) and decreased when they are replaced with 5-F-cytosine (5FC) whereas the stereochemical preference of alkylation changes upon decarbamoylation. We have studied three duplex oligonucleotides of general formula d(CGATAAXGCTAACG) in which X stands for C, 5MC or 5FC. Using a combination of molecular dynamics simulations in aqueous solution, quantum mechanics and continuum electrostatics, we have been able to (i) obtain a large series of snapshots that facilitate an understanding in atomic detail of the distinct stereochemistry of monoadduct and ICL formation by MMC and its decarbamoylated analogue, (ii) provide an explanation for the altered reactivity of MMC towards DNA molecules containing 5MC or 5FC, and (iii) show the distinct accommodation in the DNA minor groove of the different covalent modifications, particularly the most cytotoxic C1α and C1ß ICLs.

Structure-based virtual screening of small-molecule antagonists of platelet integrin αIIbß3 that do not prime the receptor to bind ligand.

Integrin αIIbß3 has emerged as an important therapeutic target for thrombotic vascular diseases owing to its pivotal role in mediating platelet aggregation through interaction with adhesive ligands. In the search for effective anti-thrombotic agents that can be administered orally without inducing the high-affinity ligand binding state, we recently discovered via high-throughput screening of 33,264 compounds a novel, αIIbß3-selective inhibitor (RUC-1) of adenosine-5'-diphosphate (ADP) -induced platelet aggregation that exhibits a different chemical scaffold and mode of binding with respect to classical Arg-Gly-Asp (RGD)-mimicking αIIbß3 antagonists. Most importantly, RUC-1 and its higher-affinity derivative, RUC-2, do not induce major conformational changes in the protein ß3 subunit or prime the receptor to bind ligand. To identify additional αIIbß3-selective chemotypes that inhibit platelet aggregation through similar mechanisms, we screened in silico over 2.5 million commercially available, 'lead-like' small molecules based on complementarity to the predicted binding mode of RUC-2 into the RUC-1-αIIbß3 crystal structure. This first reported structure-based virtual screening application to the αIIbß3 integrin led to the identification of 2 αIIbß3-selective antagonists out of 4 tested, which compares favorably with the 0.003 % "hit rate" of our previous high-throughput chemical screening study. The newly identified compounds, like RUC-1 and RUC-2, showed specificity for αIIbß3 compared to αVß3 and did not prime the receptor to bind ligand. They thus may hold promise as αIIbß3 antagonist therapeutic scaffolds.

Synthesis, modeling and evaluation of 3'-(1-aryl-1H-tetrazol-5-ylamino)-substituted 3'-deoxythymidine derivatives as potent and selective human mitochondrial thymidine kinase inhibitors.

Based on the presumed binding mode of an earlier identified inhibitor, we herein report new 3'-modified nucleosides as potent and selective inhibitors of mitochondrial thymidine kinase (TK2). A series of thirteen 3'-amino-, 3'-guanidino- and 3'-tetrazole-containing nucleosides were synthesized and evaluated for their TK2 inhibitory activity. Within the tetrazole series, compounds with nanomolar inhibitory activity were identified. A homology model of TK2 allowed to elucidate the observed activities. Introduction of a 2-bromovinyl group on C-5 of the pyrimidine base of the most promising 3'-derivative further improved the inhibitory activity, and caused a significant increase in the selectivity for TK2 versus TK1. Interestingly, for the current series of analogues, a strong correlation was observed between TK2 and Drosophila melanogaster dNK inhibition, further substantiating the phylogenetic relationship between these two nucleoside kinases.

Protein-protein interactions at an enzyme-substrate interface: characterization of transient reaction intermediates throughout a full catalytic cycle of Escherichia coli thioredoxin reductase.

A large collection of structural snapshots along a full catalytic cycle of Escherichia coli thioredoxin reductase (TrxR) has been generated and characterized using a combination of theoretical methods. Molecular models were built starting from the available X-ray crystallographic structures of dimeric wild-type TrxR in the flavin-oxidizing conformation and a C135S TrxR mutant enzyme in a flavin-reducing conformation "trapped" by a cross-link between Cys138 of TrxR and Cys32 of C35S mutant thioredoxin (Trx). The transition between these two extreme states, which is shown to be reproduced in a normal mode analysis, as well as natural cofactor binding and dissociation, were simulated for the wild-type species using unrestrained and targeted molecular dynamics following docking of oxidized Trx to reduced TrxR. The whole set of simulations provides a comprehensive structural framework for understanding the mechanism of disulfide reduction in atomic detail and identifying the most likely intermediates that facilitate entry of NADPH and exit of NADP(+). The crucial role assigned to Arg73 and Lys36 of Trx in substrate binding and complex stabilization was ascertained when R73G, R73D, and K36A site-directed mutants of Trx were shown to be impaired to different extents in their ability to be reduced by TrxR. On the basis of previous findings and the results reported herein, E. coli TrxR appears as a beautifully engineered molecular machine that is capable of synchronizing cofactor capture and ejection with substrate binding and redox activity through an interdomain twisting motion.

Intramolecular cation-pi interactions as the driving force to restrict the conformation of certain nucleosides.

Despite the well-established importance of intermolecular cation-pi interactions in molecular recognition, intramolecular cation-pi interactions have been less studied. Here we describe how the simultaneous presence of an aromatic ring at the 5'-position of an inosine derivative and a positively charged imidazolium ring in the purine base drive the conformation of the nucleoside toward a very major conformer in solution that is stabilized by an intramolecular cation-pi interaction. Therefore, the cation-pi interaction between imidazolium ions and aromatic rings can also be proposed in the design of small molecules where this type of interaction is desirable. The imidazolium ion can be obtained by a simple acidification of the pH of the media. So a simple change in pH can shift the conformational equilibrium from a random to a restricted conformation stabilized by an intramolecular cation-pi interaction. Thus the here described nucleosides can be considered as a new class of pH-dependent conformationally switchable molecules.

Human mitochondrial thymidine kinase is selectively inhibited by 3'-thiourea derivatives of beta-thymidine: identification of residues crucial for both inhibition and catalytic activity.

Substituted 3'-thiourea derivatives of beta-thymidine (dThd) and 5'-thiourea derivatives of alpha-dThd have been evaluated for their inhibitory activity against recombinant human cytosolic dThd kinase-1 (TK-1), human mitochondrial TK-2, herpes simplex virus type 1 (HSV-1) TK, and varicella-zoster virus TK. Several substituted 3'-thiourea derivatives of beta-dThd proved highly inhibitory to and selective for TK-2 (IC(50) value, 0.15-3.1 microM). The 3'-C-branched p-methylphenyl (compound 1) and 3-CF(3)-4-Cl-phenyl (compound 7) thiourea derivatives of beta-dThd showed competitive inhibition of TK-2 when dThd was used as the variable substrate (K(i) values, 0.40 and 0.05 microM, respectively), but uncompetitive inhibition in the presence of variable concentrations of ATP (K(i) values, 15 and 2.0 microM, respectively). These kinetic properties of compounds 1 and 7 against TK-2 could be accounted for by molecular modeling showing that two hydrogen bonds can be formed between the thiourea nitrogens of compound 7 and the oxygens of the gamma-phosphate of ATP. The importance of several active-site residues was assessed by site-directed mutagenesis experiments on TK-2 and the related HSV-1 TK. The low K(i)/K(m) ratios for compounds 1 and 7 (0.38 and 0.039 against dThd, and 0.75 and 0.12 against ATP, respectively) indicate that these compounds are among the most potent inhibitors of TK-2 described so far. In addition, a striking close correlation was found between the inhibitory activities of the test compounds against TK-2 and Mycobacterium tuberculosis thymidylate kinase that is strongly indicative of close structural and/or functional similarities between both enzymes in relation to their mode of interaction with these nucleoside analog inhibitors.

Structure, physiological role, and specific inhibitors of human thymidine kinase 2 (TK2): present and future.

Human mitochondrial thymidine kinase (TK2) is a pyrimidine deoxynucleoside kinase (dNK) that catalyzes the phosphorylation of pyrimidine deoxynucleosides to their corresponding deoxynucleoside 5'-monophosphates by gamma-phosphoryl transfer from ATP. In resting cells, TK2 is suggested to play a key role in the mitochondrial salvage pathway to provide pyrimidine nucleotides for mitochondrial DNA (mtDNA) synthesis and maintenance. However, recently the physiological role of TK2turned out to have direct clinical relevance as well. Point mutations in the gene encoding TK2 have been correlated to mtDNA disorders in a heterogeneous group of patients suffering from the so-called mtDNA depletion syndrome (MDS). TK2 activity could also be involved in mitochondrial toxicity associated to prolonged treatment with antiviral nucleoside analogues like AZT and FIAU. Therefore, TK2 inhibitors can be considered as valuable tools to unravel the role of TK2 in the maintenance and homeostasis of mitochondrial nucleotide pools and mtDNA, and to clarify the contribution of TK2 activity to mitochondrial toxicity of certain antivirals. Highly selective TK-2 inhibitors having an acyclic nucleoside structure and efficiently discriminating between TK-2 and the closely related TK-1 have already been reported. It is actually unclear whether these agents efficiently reach the inner mitochondrial compartment. In the present review article,structural features of TK2, MDS-related mutations observed in TK2 and their role in MDS will be discussed. Also, an update on novel and selective TK2 inhibitors will be provided.

Identification of a Mycothiol-Dependent Nitroreductase from Mycobacterium tuberculosis.

The success of Mycobacterium tuberculosis (Mtb) as a pathogen depends on the redundant and complex mechanisms it has evolved for resisting nitrosative and oxidative stresses inflicted by host immunity. Improving our understanding of these defense pathways can reveal vulnerable points in Mtb pathogenesis. In this study, we combined genetic, structural, computational, biochemical, and biophysical approaches to identify a novel enzyme class represented by Rv2466c. We show that Rv2466c is a mycothiol-dependent nitroreductase of Mtb and can reduce the nitro group of a novel mycobactericidal compound using mycothiol as a cofactor. In addition to its function as a nitroreductase, Rv2466c confers partial protection to menadione stress.

Antitumor activity, X-ray crystal structure, and DNA binding properties of thiocoraline A, a natural bisintercalating thiodepsipeptide.

The marine natural product thiocoraline A displayed approximately equal cytotoxic activity at nanomolar concentrations in a panel of 12 human cancer cell lines. X-ray diffraction analyses of orthorhombic crystals of this DNA-binding drug revealed arrays of docked pairs of staple-shaped molecules in which one pendent hydroxyquinoline chromophore from each cysteine-rich molecule appears intercalated between the two chromophores of a facing molecule. This arrangement is in contrast to the proposed mode of binding to DNA that shows the two drug chromophores clamping two stacked base pairs, in agreement with the nearest-neighbor exclusion principle. Proof of DNA sequence recognition was obtained from both classical DNase I footprinting experiments and determination of the melting temperatures of several custom-designed fluorescently labeled oligonucleotides. A rationale for the DNA-binding behavior was gained when models of thiocoraline clamping a central step embedded in several octanucleotides were built and studied by means of unrestrained molecular dynamics simulations in aqueous solution.