Discovery of diarylpyrimidine derivatives bearing piperazine sulfonyl as potent HIV-1 nonnucleoside reverse transcriptase inhibitors

HIV-1 reverse transcriptase is one of the most attractive targets for the treatment of AIDS. However, the rapid emergence of drug-resistant strains and unsatisfactory drug-like properties seriously limit the clinical application of HIV-1 non-nucleoside reverse transcriptase inhibitors (NNRTIs). Here we show that a series of piperazine sulfonyl-bearing diarylpyrimidine-based NNRTIs were designed to improve the potency against wild-type and NNRTI-resistant strains by enhancing backbone-binding interactions. Among them, compound 18b1 demonstrates single-digit nanomolar potency against the wild-type and five mutant HIV-1 strains, which is significantly better than the approved drug etravirine. The co-crystal structure analysis and molecular dynamics simulation studies were conducted to explain the broad-spectrum inhibitory activity of 18b1 against reverse transcriptase variants. Besides, compound 18b1 demonstrates improved water solubility, cytochrome P450 liability, and other pharmacokinetic properties compared to the currently approved diarylpyrimidine (DAPY) NNRTIs. Therefore, we consider compound 18b1 a potential lead compound worthy of further study.

Halo Library, a Tool for Rapid Identification of Ligand Binding Sites on Proteins Using Crystallographic Fragment Screening

X-ray crystallographic fragment screening (XCFS) uses fragment-sized molecules (∼60 to 300 Da) to access binding sites on proteins that may be inaccessible to larger drug-like molecules (>300 Da). Previous studies have shown that fragments containing halogen atoms bind more often to proteins than non-halogenated fragments. Here, we designed the Halo Library containing 46 halogenated fragments (including the "universal fragment" 4-bromopyrazole), a majority of which have been reported to bind to or inhibit one or more targets. The library was screened against the crystals of HIV-1 reverse transcriptase with the drug rilpivirine, yielding an overall hit rate of 26%. Two new binding sites were discovered, and several hot spots were identified. This small library may thus provide a convenient tool for rapidly assessing the feasibility of a target for XCFS, mapping hot spots and cryptic sites, as well as finding fragment binders that can be useful for developing drug leads.

Targeting HIV-1 Reverse Transcriptase Using a Fragment-Based Approach

Human immunodeficiency virus type I (HIV-1) is a retrovirus that infects cells of the host's immune system leading to acquired immunodeficiency syndrome and potentially death. Although treatments are available to prevent its progression, HIV-1 remains a major burden on health resources worldwide. Continued emergence of drug-resistance mutations drives the need for novel drugs that can inhibit HIV-1 replication through new pathways. The viral protein reverse transcriptase (RT) plays a fundamental role in the HIV-1 replication cycle, and multiple approved medications target this enzyme. In this study, fragment-based drug discovery was used to optimize a previously identified hit fragment (compound B-1), which bound RT at a novel site. Three series of compounds were synthesized and evaluated for their HIV-1 RT binding and inhibition. These series were designed to investigate different vectors around the initial hit in an attempt to improve inhibitory activity against RT. Our results show that the 4-position of the core scaffold is important for binding of the fragment to RT, and a lead compound with a cyclopropyl substitution was selected and further investigated. Requirements for binding to the NNRTI-binding pocket (NNIBP) and a novel adjacent site were investigated, with lead compound 27-a minimal but efficient NNRTI-offering a starting site for the development of novel dual NNIBP-Adjacent site inhibitors.

Biochemical and structural insights into SARS-CoV-2 polyprotein processing by Mpro

SARS-CoV-2, a human coronavirus, is the causative agent of the COVID-19 pandemic. Its genome is translated into two large polyproteins subsequently cleaved by viral papain-like protease and main protease (Mpro). Polyprotein processing is essential yet incompletely understood. We studied Mpro-mediated processing of the nsp7-11 polyprotein, whose mature products include cofactors of the viral replicase, and identified the order of cleavages. Integrative modeling based on mass spectrometry (including hydrogen-deuterium exchange and cross-linking) and x-ray scattering yielded a nsp7-11 structural ensemble, demonstrating shared secondary structural elements with individual nsps. The pattern of cross-links and HDX footprint of the C145A Mpro and nsp7-11 complex demonstrate preferential binding of the enzyme active site to the polyprotein junction sites and additional transient contacts to help orient the enzyme on its substrate for cleavage. Last, proteolysis assays were used to characterize the effect of inhibitors/binders on Mpro processing/inhibition using the nsp7-11 polyprotein as substrate.

Evolution of the SARS-CoV-2 proteome in three dimensions (3D) during the first 6 months of the COVID-19 pandemic

Understanding the molecular evolution of the SARS-CoV-2 virus as it continues to spread in communities around the globe is important for mitigation and future pandemic preparedness. Three-dimensional structures of SARS-CoV-2 proteins and those of other coronavirusess archived in the Protein Data Bank were used to analyze viral proteome evolution during the first 6 months of the COVID-19 pandemic. Analyses of spatial locations, chemical properties, and structural and energetic impacts of the observed amino acid changes in >48 000 viral isolates revealed how each one of 29 viral proteins have undergone amino acid changes. Catalytic residues in active sites and binding residues in protein-protein interfaces showed modest, but significant, numbers of substitutions, highlighting the mutational robustness of the viral proteome. Energetics calculations showed that the impact of substitutions on the thermodynamic stability of the proteome follows a universal bi-Gaussian distribution. Detailed results are presented for potential drug discovery targets and the four structural proteins that comprise the virion, highlighting substitutions with the potential to impact protein structure, enzyme activity, and protein-protein and protein-nucleic acid interfaces. Characterizing the evolution of the virus in three dimensions provides testable insights into viral protein function and should aid in structure-based drug discovery efforts as well as the prospective identification of amino acid substitutions with potential for drug resistance.

Structural basis of HIV inhibition by L-nucleosides: opportunities for drug development and repurposing

Infection with HIV can cripple the immune system and lead to AIDS. Hepatitis B virus (HBV) is a hepadnavirus that causes human liver diseases. Both pathogens are major public health problems affecting millions of people worldwide. The polymerases from both viruses are the most common drug target for viral inhibition, sharing common architecture at their active sites. The L-nucleoside drugs emtricitabine and lamivudine are widely used HIV reverse transcriptase (RT) and HBV polymerase (Pol) inhibitors. Nevertheless, structural details of their binding to RT(Pol)/nucleic acid remained unknown until recently. Here, we discuss the implications of these structures, alongside related complexes with L-dNTPs, for the development of novel L-nucleos(t)ide drugs, and prospects for repurposing them.

HIV-1 gp120 Antagonists Also Inhibit HIV-1 Reverse Transcriptase by Bridging the NNRTI and NRTI Sites

HIV-1 infection is typically treated using ≥2 drugs, including at least one HIV-1 reverse transcriptase (RT) inhibitor. Drugs targeting RT comprise nucleos(t)ide RT inhibitors (NRTIs) and non-nucleoside RT inhibitors (NNRTIs). NRTI-triphosphates bind at the polymerase active site and, following incorporation, inhibit DNA elongation. NNRTIs bind at an allosteric pocket ∼10 Å away from the polymerase active site. This study focuses on compounds ("NBD derivatives") originally developed to bind to HIV-1 gp120, some of which inhibit RT. We have determined crystal structures of three NBD compounds in complex with HIV-1 RT, correlating with RT enzyme inhibition and antiviral activity, to develop structure-activity relationships. Intriguingly, these compounds bridge the dNTP and NNRTI-binding sites and inhibit the polymerase activity of RT in the enzymatic assays (IC50 < 5 μM). Two of the lead compounds, NBD-14189 and NBD-14270, show potent antiviral activity (EC50 < 200 nM), and NBD-14270 shows low cytotoxicity (CC50 > 100 μM).

Crystal Structure of a Retroviral Polyprotein: Prototype Foamy Virus Protease-Reverse Transcriptase (PR-RT)

In most cases, proteolytic processing of the retroviral Pol portion of the Gag-Pol polyprotein precursor produces protease (PR), reverse transcriptase (RT), and integrase (IN). However, foamy viruses (FVs) express Pol separately from Gag and, when Pol is processed, only the IN domain is released. Here, we report a 2.9 Å resolution crystal structure of the mature PR-RT from prototype FV (PFV) that can carry out both proteolytic processing and reverse transcription but is in a configuration not competent for proteolytic or polymerase activity. PFV PR-RT is monomeric and the architecture of PFV PR is similar to one of the subunits of HIV-1 PR, which is a dimer. There is a C-terminal extension of PFV PR (101-145) that consists of two helices which are adjacent to the base of the RT palm subdomain, and anchors PR to RT. The polymerase domain of PFV RT consists of fingers, palm, thumb, and connection subdomains whose spatial arrangements are similar to the p51 subunit of HIV-1 RT. The RNase H and polymerase domains of PFV RT are connected by flexible linkers. Significant spatial and conformational (sub)domain rearrangements are therefore required for nucleic acid binding. The structure of PFV PR-RT provides insights into the conformational maturation of retroviral Pol polyproteins.

Revealing the Structural Plasticity of SARS-CoV-2 nsp7 and nsp8 Using Structural Proteomics

Coronavirus (CoV) nonstructural proteins (nsps) assemble to form the replication-transcription complex (RTC) responsible for viral RNA synthesis. nsp7 and nsp8 are important cofactors of the RTC, as they interact and regulate the activity of RNA-dependent RNA polymerase and other nsps. To date, no structure of the full-length SARS-CoV-2 nsp7:nsp8 complex has been published. The current understanding of this complex is based on structures from truncated constructs, with missing electron densities, or from related CoV species where SARS-CoV-2 nsp7 and nsp8 share upward of 90% sequence identity. Despite available structures solved using crystallography and cryo-EM representing detailed static snapshots of the nsp7:nsp8 complex, it is evident that the complex has a high degree of structural plasticity. However, relatively little is known about the conformational dynamics of the individual proteins and how they complex to interact with other nsps. Here, the solution-based structural proteomic techniques, hydrogen-deuterium exchange mass spectrometry (HDX-MS) and cross-linking mass spectrometry (XL-MS), illuminate the dynamics of SARS-CoV-2 full-length nsp7 and nsp8 proteins and the nsp7:nsp8 protein complex. Results presented from the two techniques are complementary and validate the interaction surfaces identified from the published three-dimensional heterotetrameric crystal structure of the SARS-CoV-2 truncated nsp7:nsp8 complex. Furthermore, mapping of XL-MS data onto higher-order complexes suggests that SARS-CoV-2 nsp7 and nsp8 do not assemble into a hexadecameric structure as implied by the SARS-CoV full-length nsp7:nsp8 crystal structure. Instead, our results suggest that the nsp7:nsp8 heterotetramer can dissociate into a stable dimeric unit that might bind to nsp12 in the RTC without significantly altering nsp7-nsp8 interactions.

2,4,5-Trisubstituted Pyrimidines as Potent HIV-1 NNRTIs: Rational Design, Synthesis, Activity Evaluation, and Crystallographic Studies

There is an urgent unmet medical need for novel human immunodeficiency virus type 1 (HIV-1) inhibitors that are effective against a variety of NNRTI-resistance mutations. We report our research efforts aimed at discovering a novel chemotype of anti-HIV-1 agents with improved potency against a variety of NNRTI-resistance mutations in this paper. Structural modifications of the lead K-5a2 led to the identification of a potent inhibitor 16c. 16c yielded highly potent anti-HIV-1 activities and improved resistance profiles compared with the approved drug etravirine. The co-crystal structure revealed the key role of the water networks surrounding the NNIBP for binding and for resilience against resistance mutations, while suggesting further extension of 16c toward the NNRTI-adjacent site as a lead development strategy. Furthermore, 16c demonstrated favorable pharmacokinetic and safety properties, suggesting the potential of 16c as a promising anti-HIV-1 drug candidate.

Resolving the Dynamic Motions of SARS-CoV-2 nsp7 and nsp8 Proteins Using Structural Proteomics

Coronavirus (CoV) non-structural proteins (nsps) assemble to form the replication-transcription complex (RTC) responsible for viral RNA synthesis. nsp7 and nsp8 are important cofactors of the RTC, as they interact and regulate the activity of RNA-dependent RNA polymerase (RdRp) and other nsps. To date, no structure of full-length SARS-CoV-2 nsp7:nsp8 complex has been published. Current understanding of this complex is based on structures from truncated constructs or with missing electron densities and complexes from related CoV species with which SARS-CoV-2 nsp7 and nsp8 share upwards of 90% sequence identity. Despite available structures being solved using crystallography and cryo-EM representing detailed snapshots of the nsp7:nsp8 complex, it is evident that the complex has a high degree of structural plasticity. However, relatively little is known about the conformational dynamics of the complex and how it assembles to interact with other nsps. Here, the solution-based structural proteomic techniques, hydrogen-deuterium exchange mass spectrometry (HDX-MS) and crosslinking mass spectrometry (XL-MS), illuminate the structural dynamics of the SARS-CoV-2 full-length nsp7:nsp8 complex. The results presented from the two techniques are complementary and validate the interaction surfaces identified from the published three-dimensional heterotetrameric crystal structure of SARS-CoV-2 truncated nsp7:nsp8 complex. Furthermore, mapping of XL-MS data onto higher order complexes suggests that SARS-CoV-2 nsp7 and nsp8 do not assemble into a hexadecameric structure as implied by the SARS-CoV full-length nsp7:nsp8 crystal structure. Instead our results suggest that the nsp7:nsp8 heterotetramer can dissociate into a stable dimeric unit that might bind to nsp12 in the RTC without altering nsp7-nsp8 interactions.

Evolution of the SARS-CoV-2 proteome in three dimensions (3D) during the first six months of the COVID-19 pandemic

Three-dimensional structures of SARS-CoV-2 and other coronaviral proteins archived in the Protein Data Bank were used to analyze viral proteome evolution during the first six months of the COVID-19 pandemic. Analyses of spatial locations, chemical properties, and structural and energetic impacts of the observed amino acid changes in >48,000 viral proteome sequences showed how each one of the 29 viral study proteins have undergone amino acid changes. Structural models computed for every unique sequence variant revealed that most substitutions map to protein surfaces and boundary layers with a minority affecting hydrophobic cores. Conservative changes were observed more frequently in cores versus boundary layers/surfaces. Active sites and protein-protein interfaces showed modest numbers of substitutions. Energetics calculations showed that the impact of substitutions on the thermodynamic stability of the proteome follows a universal bi-Gaussian distribution. Detailed results are presented for six drug discovery targets and four structural proteins comprising the virion, highlighting substitutions with the potential to impact protein structure, enzyme activity, and functional interfaces. Characterizing the evolution of the virus in three dimensions provides testable insights into viral protein function and should aid in structure-based drug discovery efforts as well as the prospective identification of amino acid substitutions with potential for drug resistance.

Integrative structural biology studies of HIV-1 reverse transcriptase binding to a high-affinity DNA aptamer

The high-resolution crystal structure of HIV-1 reverse transcriptase (RT) bound to a 38-mer DNA hairpin aptamer with low pM affinity was previously described. The high-affinity binding aptamer contained 2'-O-methyl modifications and a seven base-pair GC-rich tract and the structure of the RT-aptamer complex revealed specific contacts between RT and the template strand of the aptamer. Similar to all crystal structures of RT bound to nucleic acid template-primers, the aptamer bound RT with a bend in the duplex DNA. To understand the structural basis for the ultra-high-affinity aptamer binding, an integrative structural biology approach was used. Hydrogen-deuterium exchange coupled to liquid chromatography-mass spectrometry (HDX-MS) was used to examine the structural dynamics of RT alone and in the presence of the DNA aptamer. RT was selectively labeled with 15N to unambiguously identify peptides from each subunit. HDX of unliganded RT shows a mostly stable core. The p66 fingers and thumb subdomains, and the RNase H domain are relatively dynamic. HDX indicates that both the aptamer and a scrambled version significantly stabilize regions of RT that are dynamic in the absence of DNA. No substantial differences in RT dynamics are observed between aptamer and scrambled aptamer binding, despite a large difference in binding affinity. Small-angle X-ray scattering and circular dichroism spectroscopy were used to investigate the aptamer conformation in solution and revealed a pre-bent DNA that possesses both A- and B-form helical character. Both the 2'-O-methyl modifications and the GC tract appear to contribute to an energetically favorable conformation for binding to RT that contributes to the aptamer's ultra-high affinity for RT. The X-ray structure of RT with an RNA/DNA version of the aptamer at 2.8 Å resolution revealed a potential role of the hairpin positioning in affinity. Together, the data suggest that both the 2'-O-methyl modifications and the GC tract contribute to an energetically favorable conformation for high-affinity binding to RT.

IDD388 Polyhalogenated Derivatives as Probes for an Improved Structure-Based Selectivity of AKR1B10 Inhibitors

Human enzyme aldo-keto reductase family member 1B10 (AKR1B10) has evolved as a tumor marker and promising antineoplastic target. It shares high structural similarity with the diabetes target enzyme aldose reductase (AR). Starting from the potent AR inhibitor IDD388, we have synthesized a series of derivatives bearing the same halophenoxyacetic acid moiety with an increasing number of bromine (Br) atoms on its aryl moiety. Next, by means of IC₅₀ measurements, X-ray crystallography, WaterMap analysis, and advanced binding free energy calculations with a quantum-mechanical (QM) approach, we have studied their structure-activity relationship (SAR) against both enzymes. The introduction of Br substituents decreases AR inhibition potency but improves it in the case of AKR1B10. Indeed, the Br atoms in ortho position may impede these drugs to fit into the AR prototypical specificity pocket. For AKR1B10, the smaller aryl moieties of MK181 and IDD388 can bind into the external loop A subpocket. Instead, the bulkier MK184, MK319, and MK204 open an inner specificity pocket in AKR1B10 characterized by a π-π stacking interaction of their aryl moieties and Trp112 side chain in the native conformation (not possible in AR). Among the three compounds, only MK204 can make a strong halogen bond with the protein (-4.4 kcal/mol, using QM calculations), while presenting the lowest desolvation cost among all the series, translated into the most selective and inhibitory potency AKR1B10 (IC₅₀ = 80 nM). Overall, SAR of these IDD388 polyhalogenated derivatives have unveiled several distinctive AKR1B10 features (shape, flexibility, hydration) that can be exploited to design novel types of AKR1B10 selective drugs.

X-Ray Crystal Structure of the Full Length Human Chitotriosidase (CHIT1) Reveals Features of Its Chitin Binding Domain

Chitinases are enzymes that catalyze the hydrolysis of chitin. Human chitotriosidase (CHIT1) is one of the two active human chitinases, involved in the innate immune response and highly expressed in a variety of diseases. CHIT1 is composed of a catalytic domain linked by a hinge to its chitin binding domain (ChBD). This latter domain belongs to the carbohydrate-binding module family 14 (CBM14 family) and facilitates binding to chitin. So far, the available crystal structures of the human chitinase CHIT1 and the Acidic Mammalian Chitinase (AMCase) comprise only their catalytic domain. Here, we report a crystallization strategy combining cross-seeding and micro-seeding cycles which allowed us to obtain the first crystal structure of the full length CHIT1 (CHIT1-FL) at 1.95 Å resolution. The CHIT1 chitin binding domain (ChBD_CHIT1) structure shows a distorted β-sandwich 3D fold, typical of CBM14 family members. Accordingly, ChBD_CHIT1 presents six conserved cysteine residues forming three disulfide bridges and several exposed aromatic residues that probably are involved in chitin binding, including the highly conserved Trp465 in a surface- exposed conformation. Furthermore, ChBD_CHIT1 presents a positively charged surface which may be involved in electrostatic interactions. Our data highlight the strong structural conservation of CBM14 family members and uncover the structural similarity between the human ChBD_CHIT1, tachycitin and house mite dust allergens. Overall, our new CHIT1-FL structure, determined with an adapted crystallization approach, is one of the few complete bi-modular chitinase structures available and reveals the structural features of a human CBM14 domain.

Structural Determinants of the Selectivity of 3-Benzyluracil-1-acetic Acids toward Human Enzymes Aldose Reductase and AKR1B10

The human enzymes aldose reductase (AR) and AKR1B10 have been thoroughly explored in terms of their roles in diabetes, inflammatory disorders, and cancer. In this study we identified two new lead compounds, 2-(3-(4-chloro-3-nitrobenzyl)-2,4-dioxo-3,4-dihydropyrimidin-1(2H)-yl)acetic acid (JF0048, 3) and 2-(2,4-dioxo-3-(2,3,4,5-tetrabromo-6-methoxybenzyl)-3,4-dihydropyrimidin-1(2H)-yl)acetic acid (JF0049, 4), which selectively target these enzymes. Although 3 and 4 share the 3-benzyluracil-1-acetic acid scaffold, they have different substituents in their aryl moieties. Inhibition studies along with thermodynamic and structural characterizations of both enzymes revealed that the chloronitrobenzyl moiety of compound 3 can open the AR specificity pocket but not that of the AKR1B10 cognate. In contrast, the larger atoms at the ortho and/or meta positions of compound 4 prevent the AR specificity pocket from opening due to steric hindrance and provide a tighter fit to the AKR1B10 inhibitor binding pocket, probably enhanced by the displacement of a disordered water molecule trapped in a hydrophobic subpocket, creating an enthalpic signature. Furthermore, this selectivity also occurs in the cell, which enables the development of a more efficient drug design strategy: compound 3 prevents sorbitol accumulation in human retinal ARPE-19 cells, whereas 4 stops proliferation in human lung cancer NCI-H460 cells.

Substrate Specificity, Inhibitor Selectivity and Structure-Function Relationships of Aldo-Keto Reductase 1B15: A Novel Human Retinaldehyde Reductase

Human aldo-keto reductase 1B15 (AKR1B15) is a newly discovered enzyme which shares 92% amino acid sequence identity with AKR1B10. While AKR1B10 is a well characterized enzyme with high retinaldehyde reductase activity, involved in the development of several cancer types, the enzymatic activity and physiological role of AKR1B15 are still poorly known. Here, the purified recombinant enzyme has been subjected to substrate specificity characterization, kinetic analysis and inhibitor screening, combined with structural modeling. AKR1B15 is active towards a variety of carbonyl substrates, including retinoids, with lower kcat and Km values than AKR1B10. In contrast to AKR1B10, which strongly prefers all-trans-retinaldehyde, AKR1B15 exhibits superior catalytic efficiency with 9-cis-retinaldehyde, the best substrate found for this enzyme. With ketone and dicarbonyl substrates, AKR1B15 also shows higher catalytic activity than AKR1B10. Several typical AKR inhibitors do not significantly affect AKR1B15 activity. Amino acid substitutions clustered in loops A and C result in a smaller, more hydrophobic and more rigid active site in AKR1B15 compared with the AKR1B10 pocket, consistent with distinct substrate specificity and narrower inhibitor selectivity for AKR1B15.

The Effect of Halogen-to-Hydrogen Bond Substitution on Human Aldose Reductase Inhibition

The effect of halogen-to-hydrogen bond substitution on the binding energetics and biological activity of a human aldose reductase inhibitor has been studied using X-ray crystallography, IC50 measurements, advanced binding free energy calculations, and simulations. The replacement of Br or I atoms by an amine (NH2) group has not induced changes in the original geometry of the complex, which made it possible to study the isolated features of selected noncovalent interactions in a biomolecular complex.

New insights into the enzymatic mechanism of human chitotriosidase (CHIT1) catalytic domain by atomic resolution X-ray diffraction and hybrid QM/MM

Chitotriosidase (CHIT1) is a human chitinase belonging to the highly conserved glycosyl hydrolase family 18 (GH18). GH18 enzymes hydrolyze chitin, an N-acetylglucosamine polymer synthesized by lower organisms for structural purposes. Recently, CHIT1 has attracted attention owing to its upregulation in immune-system disorders and as a marker of Gaucher disease. The 39 kDa catalytic domain shows a conserved cluster of three acidic residues, Glu140, Asp138 and Asp136, involved in the hydrolysis reaction. Under an excess concentration of substrate, CHIT1 and other homologues perform an additional activity, transglycosylation. To understand the catalytic mechanism of GH18 chitinases and the dual enzymatic activity, the structure and mechanism of CHIT1 were analyzed in detail. The resolution of the crystals of the catalytic domain was improved from 1.65 Å (PDB entry 1waw) to 0.95-1.10 Å for the apo and pseudo-apo forms and the complex with chitobiose, allowing the determination of the protonation states within the active site. This information was extended by hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. The results suggest a new mechanism involving changes in the conformation and protonation state of the catalytic triad, as well as a new role for Tyr27, providing new insights into the hydrolysis and transglycosylation activities.

Structural analysis of sulindac as an inhibitor of aldose reductase and AKR1B10

Aldose reductase (AR, AKR1B1) and AKR1B10 are enzymes implicated in important pathologies (diabetes and cancer) and therefore they have been proposed as suitable targets for drug development. Sulindac is the metabolic precursor of the potent non-steroidal anti-inflammatory drug (NSAID) sulindac sulfide, which suppresses prostaglandin production by inhibition of cyclooxygenases (COX). In addition, sulindac has been found to be one of the NSAIDs with higher antitumoral activity, presumably through COX inhibition. However, sulindac anticancer activity could be partially mediated through COX-independent mechanisms, including the participation of AR and AKR1B10. Previously, it had been shown that sulindac and sulindac sulfone were good AR inhibitors and the structure of the ternary complex with NADP(+) and sulindac was described (PDB ID 3U2C). In this work, we determined the three-dimensional structure of AKR1B10 with sulindac and established structure-activity relationships (SAR) of sulindac and their derivatives with AR and AKR1B10. The difference in the IC50 values for sulindac between AR (0.36 μM) and AKR1B10 (2.7 μM) might be explained by the different positioning and stacking interaction given by Phe122/Phe123, and by the presence of two buried and ordered water molecules in AKR1B10 but not in AR. Moreover, SAR analysis shows that the substitution of the sulfinyl group is structurally allowed in sulindac derivatives. Hence, sulindac and its derivatives emerge as lead compounds for the design of more potent and selective AR and AKR1B10 inhibitors.

Identification of a novel polyfluorinated compound as a lead to inhibit the human enzymes aldose reductase and AKR1B10: structure determination of both ternary complexes and implications for drug design

Aldo-keto reductases (AKRs) are mostly monomeric enzymes which fold into a highly conserved (α/β)8 barrel, while their substrate specificity and inhibitor selectivity are determined by interaction with residues located in three highly variable external loops. The closely related human enzymes aldose reductase (AR or AKR1B1) and AKR1B10 are of biomedical interest because of their involvement in secondary diabetic complications (AR) and in cancer, e.g. hepatocellular carcinoma and smoking-related lung cancer (AKR1B10). After characterization of the IC50 values of both AKRs with a series of polyhalogenated compounds, 2,2',3,3',5,5',6,6'-octafluoro-4,4'-biphenyldiol (JF0064) was identified as a lead inhibitor of both enzymes with a new scaffold (a 1,1'-biphenyl-4,4'-diol). An ultrahigh-resolution X-ray structure of the AR-NADP(+)-JF0064 complex has been determined at 0.85 Å resolution, allowing it to be observed that JF0064 interacts with the catalytic residue Tyr48 through a negatively charged hydroxyl group (i.e. the acidic phenol). The non-competitive inhibition pattern observed for JF0064 with both enzymes suggests that this acidic hydroxyl group is also present in the case of AKR1B10. Moreover, the combination of surface lysine methylation and the introduction of K125R and V301L mutations enabled the determination of the X-ray crystallographic structure of the corresponding AKR1B10-NADP(+)-JF0064 complex. Comparison of the two structures has unveiled some important hints for subsequent structure-based drug-design efforts.

Modulation of aldose reductase inhibition by halogen bond tuning

In this paper, we studied a designed series of aldose reductase (AR) inhibitors. The series was derived from a known AR binder, which had previously been shown to form a halogen bond between its bromine atom and the oxygen atom of the Thr-113 side chain of AR. In the series, the strength of the halogen bond was modulated by two factors, namely bromine-iodine substitution and the fluorination of the aromatic ring in several positions. The role of the single halogen bond in AR-ligand binding was elucidated by advanced binding free energy calculations involving the semiempirical quantum chemical Hamiltonian. The results were complemented with ultrahigh-resolution X-ray crystallography and IC50 measurements. All of the AR inhibitors studied were shown by X-ray crystallography to bind in an identical manner. Further, it was demonstrated that it was possible to decrease the IC50 value by about 1 order of magnitude by tuning the strength of the halogen bond by a monoatomic substitution. The calculations revealed that the protein-ligand interaction energy increased upon the substitution of iodine for bromine or upon the addition of electron-withdrawing fluorine atoms to the ring. However, the effect on the binding affinity was found to be more complex due to the change of the solvation/desolvation properties within the ligand series. The study shows that it is possible to modulate the strength of a halogen bond in a protein-ligand complex as was designed based on the previous studies of low-molecular-weight complexes.

A missense mutation in ALDH1A3 causes isolated microphthalmia/anophthalmia in nine individuals from an inbred Muslim kindred

Nine affected individuals with isolated anophthalmia/microphthalmia from a large Muslim-inbred kindred were investigated. Assuming autosomal-recessive mode of inheritance, whole-genome linkage analysis, on DNA samples from four affected individuals, was undertaken. Homozygosity mapping techniques were employed and a 1.5-Mbp region, homozygous in all affected individuals, was delineated. The region contained nine genes, one of which, aldehyde dehydrogenase 1 (ALDH1A3), was a clear candidate. This gene seems to encode a key enzyme in the formation of a retinoic-acid gradient along the dorsoventral axis during an early eye development and the development of the olfactory system. Sanger sequence analysis revealed a missense mutation, causing a substitution of valine (Val) to methionine (Met) at position 71. Analyzing the p.Val71Met missense mutation using standard open access software (MutationTaster online, PolyPhen, SIFT/PROVEAN) predicts this variant to be damaging. Enzymatic activity, studied in vitro, showed no changes between the mutated and the wild-type ALDH1A3 protein.

Crystal packing modifies ligand binding affinity: the case of aldose reductase

The relationship between the structures of protein-ligand complexes existing in the crystal and in solution, essential in the case of fragment-based screening by X-ray crystallography (FBS-X), has been often an object of controversy. To address this question, simultaneous co-crystallization and soaking of two inhibitors with different ratios, Fidarestat (FID; K(d) = 6.5 nM) and IDD594 (594; K(d) = 61 nM), which bind to h-aldose reductase (AR), have been performed. The subatomic resolution of the crystal structures allows the differentiation of both inhibitors, even when the structures are almost superposed. We have determined the occupation ratio in solution by mass spectrometry (MS) Occ(FID)/Occ(594) = 2.7 and by X-ray crystallography Occ(FID)/Occ(594) = 0.6. The occupancies in the crystal and in solution differ 4.6 times, implying that ligand binding potency is influenced by crystal contacts. A structural analysis shows that the Loop A (residues 122-130), which is exposed to the solvent, is flexible in solution, and is involved in packing contacts within the crystal. Furthermore, inhibitor 594 contacts the base of Loop A, stabilizing it, while inhibitor FID does not. This is shown by the difference in B-factors of the Loop A between the AR-594 and AR-FID complexes. A stable loop diminishes the entropic energy barrier to binding, favoring 594 versus FID. Therefore, the effect of the crystal environment should be taken into consideration in the X-ray diffraction analysis of ligand binding to proteins. This conclusion highlights the need for additional methodologies in the case of FBS-X to validate this powerful screening technique, which is widely used.