Discovery of a novel androgen receptor antagonist, MEL-6, with stereoselective activity and optimization of its metabolic stability

A new chemical scaffold with antagonistic activity towards the androgen receptor (AR) was identified. The parent compound, (3-Methoxy-N-[1-methyl-2-(4-phenyl-1-piperazinyl)-2-(2-thienyl)ethyl]benzamide) referred to as MEL-6, binds in the ligand binding pocket of AR and induces an antagonistic conformation of the ligand binding domain, even in presence of the antagonist-to-agonist switch mutations W741C, T877A and F876L-T877A. MEL-6 has antiproliferative effects on several AR positive prostate cancer cell lines. We further identified AR as the specific target of MEL-6 since it demonstrates little effect on other steroid receptors. In LNCaP cells it also inhibits the androgen-regulated transcriptome. These findings identify MEL-6 as a promising candidate for treatment of patients with prostate tumors that have become resistant to current clinically used AR antagonists. Analytical studies on the chemical composition of MEL-6 identified the presence of four isomers (two enantiomeric pairs), among which one isomer is responsible for the antiandrogenic activity. We therefore developed a synthetic route towards the selective preparation of the active enantiomeric pair. Various MEL-6-like analogues had improved metabolic stability while maintaining antiandrogenic activity. Metabolite identification of MEL-6 derivatives pinpointed N-dealkylation of the piperazine as the main mode for inactivation by liver enzymes. For further structural optimization, MEL-6 derivatives were purchased or synthesized having alterations on the N-phenyl group of the piperazine, the benzoyl group and additionally substituting the thiophen-2-yl ring of MEL-6 to a phenyl ring. This optimization process resulted in compound 12b with sustained AR inhibition and a 4-fold increased half-life due to the 1-(5-chloro-2-methylphenyl)-piperazine substitution, thienyl-to-phenyl substitution and chloro in para-position of the benzoyl group.

Structural mechanism underlying variations in DNA binding by the androgen receptor

The androgen receptor (AR) is a steroid activated transcription factor which recognizes DNA motifs resembling inverted repeats of a conserved 5'-AGAACA-3'-like hexanucleotides separated by a three-nucleotide spacer from a similar, but less conserved hexanucleotide. Here, we report the structures of the human AR DNA binding domain (DBD) bound to two natural AREs (C3 and MTV) in head-to-head dimer conformations, diffracting at 2.05 Å and 2.25 Å, respectively. These structures help to explain the impact of androgen insensitivity mutations on the structure integrity, DNA binding and DBD dimerization. The binding affinity of the AR DBD to different DNA motifs were measured by the BioLayer Interferometry (BLI) and further validated by Molecular Dynamics (MD) simulations. This shows that the high binding affinity of the first DBD to the upstream 5'-AGAACA-3' motif induces the cooperative binding of the second DBD to the second hexanucleotide. Our data indicate identical interaction of the DBDs to the upstream hexanucleotides, while forming an induced closer contact of the second DBD on the non-canonical hexanucleotides. The variation in binding between the DBD monomers are the result of differences in DNA occupancy, protein-protein interactions, DNA binding affinity, and DNA binding energy profiles. We propose this has functional consequences.

Rational evolution for altering the ligand preference of estrogen receptor alpha

Estrogen receptor α is commonly used in synthetic biology to control the activity of genome editing tools. The activating ligands, estrogens, however, interfere with various cellular processes, thereby limiting the applicability of this receptor. Altering its ligand preference to chemicals of choice solves this hurdle but requires adaptation of unspecified ligand-interacting residues. Here, we provide a solution by combining rational protein design with multi-site-directed mutagenesis and directed evolution of stably integrated variants in Saccharomyces cerevisiae. This method yielded an estrogen receptor variant, named TERRA, that lost its estrogen responsiveness and became activated by tamoxifen, an anti-estrogenic drug used for breast cancer treatment. This tamoxifen preference of TERRA was maintained in mammalian cells and mice, even when fused to Cre recombinase, expanding the mammalian synthetic biology toolbox. Not only is our platform transferable to engineer ligand preference of any steroid receptor, it can also profile drug-resistance landscapes for steroid receptor-targeted therapies.

Insights in inulin binding and inulin oligosaccharide formation by novel multi domain endo-inulinases from Botrytis cinerea

World-wide, pathogenic fungi such as Botrytis cinerea cause tremendous yield losses in terms of food production and post-harvest food decay. Many fungi produce inulin-type oligosaccharides (IOSs) from inulin through endo-inulinases which typically show a two domain structure. B.cinerea lacks a two domain endo-inulinase but contains a three domain structure instead. Genome mining revealed three and four domain (d4) enzymes in the fungal kingdom. Here, three and two domain enzymes were compared in their capacity to produce IOSs from inulin. Hill kinetics were observed in three domain enzymes as compared to Michaelis-Menten kinetics in two domain enzymes, suggesting that the N-terminal extension functions as a carbohydrate binding module. Analysis of the IOS product profiles generated from purified GF6, GF12, GF16 and GF18 inulins and extensive sugar docking approaches led to enhanced insights in the active site functioning, revealing subtle differences between the endo-inulinases from Aspergillus niger and B. cinerea. Improved insights in structure-function relationships in fungal endo-inulinases offer opportunities to develop superior enzymes for the production of specific IOS formulations to improve plant and animal health (priming agents, prebiotics).

Gollop-Wolfgang Complex Is Associated with a Monoallelic Variation in WNT11

Gollop-Wolfgang complex (GWC) is a rare congenital limb anomaly characterized by tibial aplasia with femur bifurcation, ipsilateral bifurcation of the thigh bone, and split hand and monodactyly of the feet, resulting in severe and complex limb deformities. The genetic basis of GWC, however, has remained elusive. We studied a three-generation family with four GWC-affected family members. An analysis of whole-genome sequencing results using a custom pipeline identified the WNT11 c.1015G>A missense variant associated with the phenotype. In silico modelling and an in vitro reporter assay further supported the link between the variant and GWC. This finding further contributes to mapping the genetic heterogeneity underlying split hand/foot malformations in general and in GWC specifically.

Heterozygous mutations in the C-terminal domain of COPA underlie a complex autoinflammatory syndrome

Mutations in the N-terminal WD40 domain of coatomer protein complex subunit α (COPA) cause a type I interferonopathy, typically characterized by alveolar hemorrhage, arthritis, and nephritis. We described 3 heterozygous mutations in the C-terminal domain (CTD) of COPA (p.C1013S, p.R1058C, and p.R1142X) in 6 children from 3 unrelated families with a similar syndrome of autoinflammation and autoimmunity. We showed that these CTD COPA mutations disrupt the integrity and the function of coat protein complex I (COPI). In COPAR1142X and COPAR1058C fibroblasts, we demonstrated that COPI dysfunction causes both an anterograde ER-to-Golgi and a retrograde Golgi-to-ER trafficking defect. The disturbed intracellular trafficking resulted in a cGAS/STING-dependent upregulation of the type I IFN signaling in patients and patient-derived cell lines, albeit through a distinct molecular mechanism in comparison with mutations in the WD40 domain of COPA. We showed that CTD COPA mutations induce an activation of ER stress and NF-κB signaling in patient-derived primary cell lines. These results demonstrate the importance of the integrity of the CTD of COPA for COPI function and homeostatic intracellular trafficking, essential to ER homeostasis. CTD COPA mutations result in disease by increased ER stress, disturbed intracellular transport, and increased proinflammatory signaling.

Exploration of 1,2,3-triazolo fused triterpenoids as inhibitors of human coronavirus 229E targeting the viral nsp15 protein

The coronavirus disease-19 (COVID-19) pandemic has raised major interest in innovative drug concepts to suppress human coronavirus (HCoV) infections. We previously reported on a class of 1,2,3-triazolo fused betulonic acid derivatives causing strong inhibition of HCoV-229E replication via the viral nsp15 protein, which is proposedly related to compound binding at an intermonomer interface in hexameric nsp15. In the present study, we further explored the structure-activity relationship (SAR), by varying the substituent at the 1,2,3-triazolo ring as well as the triterpenoid skeleton. The 1,2,3-triazolo fused triterpenoids were synthesized by a multicomponent triazolization reaction, which has been developed in-house. Several analogs possessing a betulin, oleanolic acid, or ursolic acid core displayed favorable activity and selectivity (EC50 values for HCoV-229E: 1.6-3.5 μM), but neither of them proved as effective as the lead compound containing betulonic acid. The 18β-glycyrrhetinic acid-containing analogs had low selectivity. The antiviral findings were rationalized by in silico docking in the available structure of the HCoV-229E nsp15 protein. The new SAR insights will aid the further development of these 1,2,3-triazolo fused triterpenoid compounds as a unique type of coronavirus inhibitors.

Human Autosomal Recessive DNA Polymerase Delta 3 Deficiency Presenting as Omenn Syndrome

The DNA polymerase δ complex (PolD), comprising catalytic subunit POLD1 and accessory subunits POLD2, POLD3, and POLD4, is essential for DNA synthesis and is central to genome integrity. We identified, by whole exome sequencing, a homozygous missense mutation (c.1118A > C; p.K373T) in POLD3 in a patient with Omenn syndrome. The patient exhibited severely decreased numbers of naïve T cells associated with a restricted T-cell receptor repertoire and a defect in the early stages of TCR recombination. The patient received hematopoietic stem cell transplantation at age 6 months. He manifested progressive neurological regression and ultimately died at age 4 years. We performed molecular and functional analysis of the mutant POLD3 and assessed cell cycle progression as well as replication-associated DNA damage. Patient fibroblasts showed a marked defect in S-phase entry and an enhanced number of double-stranded DNA break-associated foci despite normal expression levels of PolD components. The cell cycle defect was rescued by transduction with WT POLD3. This study validates autosomal recessive POLD3 deficiency as a novel cause of profound T-cell deficiency and Omenn syndrome.

Intracellular BAPTA directly inhibits PFKFB3, thereby impeding mTORC1-driven Mcl-1 translation and killing MCL-1-addicted cancer cells

Intracellular Ca2+ signals control several physiological and pathophysiological processes. The main tool to chelate intracellular Ca2+ is intracellular BAPTA (BAPTAi), usually introduced into cells as a membrane-permeant acetoxymethyl ester (BAPTA-AM). Previously, we demonstrated that BAPTAi enhanced apoptosis induced by venetoclax, a BCL-2 antagonist, in diffuse large B-cell lymphoma (DLBCL). This finding implied a novel interplay between intracellular Ca2+ signaling and anti-apoptotic BCL-2 function. Hence, we set out to identify the underlying mechanisms by which BAPTAi enhances cell death in B-cell cancers. In this study, we discovered that BAPTAi alone induced apoptosis in hematological cancer cell lines that were highly sensitive to S63845, an MCL-1 antagonist. BAPTAi provoked a rapid decline in MCL-1-protein levels by inhibiting mTORC1-driven Mcl-1 translation. These events were not a consequence of cell death, as BAX/BAK-deficient cancer cells exhibited similar downregulation of mTORC1 activity and MCL-1-protein levels. Next, we investigated how BAPTAi diminished mTORC1 activity and identified its ability to impair glycolysis by directly inhibiting 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase 3 (PFKFB3) activity, a previously unknown effect of BAPTAi. Notably, these effects were also induced by a BAPTAi analog with low affinity for Ca2+. Consequently, our findings uncover PFKFB3 inhibition as an Ca2+-independent mechanism through which BAPTAi impairs cellular metabolism and ultimately compromises the survival of MCL-1-dependent cancer cells. These findings hold two important implications. Firstly, the direct inhibition of PFKFB3 emerges as a key regulator of mTORC1 activity and a promising target in MCL-1-dependent cancers. Secondly, cellular effects caused by BAPTAi are not necessarily related to Ca2+ signaling. Our data support the need for a reassessment of the role of Ca2+ in cellular processes when findings were based on the use of BAPTAi.

Polyketide Synthase-Mediated O-Methyloxime Formation in the Biosynthesis of the Oximidine Anticancer Agents

Bacterial trans-acyltransferase polyketide synthases (trans-AT PKSs) are modular megaenzymes that employ unusual catalytic domains to assemble diverse bioactive natural products. One such PKS is responsible for the biosynthesis of the oximidine anticancer agents, oxime-substituted benzolactone enamides that inhibit vacuolar H+ -ATPases. Here, we describe the identification of the oximidine gene cluster in Pseudomonas baetica and the characterization of four novel oximidine variants, including a structurally simpler intermediate that retains potent anticancer activity. Using a combination of in vivo, in vitro and computational approaches, we experimentally elucidate the oximidine biosynthetic pathway and reveal an unprecedented mechanism for O-methyloxime formation. We show that this process involves a specialized monooxygenase and methyltransferase domain and provide insight into their activity, mechanism and specificity. Our findings expand the catalytic capabilities of trans-AT PKSs and identify potential strategies for the production of novel oximidine analogues.

A widely applicable and cost-effective method for specific RNA-protein complex isolation

Although methodological advances have been made over the past years, a widely applicable, easily scalable and cost-effective procedure that can be routinely used to isolate specific ribonucleoprotein complexes (RNPs) remains elusive. We describe the "Silica-based Acidic Phase Separation (SAPS)-capture" workflow. This versatile method combines previously described techniques in a cost-effective, optimal and widely applicable protocol. The specific RNP isolation procedure is performed on a pre-purified RNP sample instead of cell lysate. This combination of protocols results in an increased RNP/bead ratio and by consequence a reduced experimental cost. To validate the method, the 18S rRNP of S. cerevisiae was captured and to illustrate its applicability we isolated the complete repertoire of RNPs in A. thaliana. The procedure we describe can provide the community with a powerful tool to advance the study of the ribonome of a specific RNA molecule in any organism or tissue type.

Study of novel androgen receptor V770 variant in androgen insensitivity syndrome patients reveals the transitional state of the androgen receptor ligand binding domain homodimer

We report the discovery of the androgen receptor missense mutation V770D, that was found in two sisters suffering from complete androgen insensitivity. Experimental validation of AR V770 variants demonstrated that AR V770D was transcriptionally inactive due to the inability to dimerize and a reduced ligand binding affinity. The more conservative AR V770A variant showed a dimerization defect at low levels of DHT with a partial recovery of the transcriptional activity and of the receptor's ability to dimerize when increasing the DHT levels. With V770 located outside of the proposed LBD dimerization interface of the AR LBD homodimer crystal structure, the effects of the V770A mutation on AR dimerization were unexpected. We therefore explored whether the AR LBD dimerization interface would be better described by an alternative dimerization mode based on available human homodimeric LBD crystal structures of other nuclear receptors. Superposition of the monomeric AR LBD in the homodimeric crystal structures of GR, PR, ER, CAR, TRβ, and HNF-4α showed that the GR-like LBD dimer model was energetically the most stable. Moreover, V770 was a key energy residue in the GR-like LBD dimer while it was not involved in the stabilization of the AR LBD homodimer according to the crystal structure. Additionally, the observation that 4 AIS mutations impacted the stability of the AR LBD dimer while 16 mutations affected the GR-like LBD dimer, suggested that the AR LBD dimer crystal is a snapshot of a breathing AR LBD homodimer that can transition into the GR-like LBD dimer model.

In vitro characterization of a novel Arg102 mutation in the ADAMTS13 metalloprotease domain

BACKGROUND

Congenital thrombotic thrombocytopenic purpura is caused by defects in the ADAMTS13 gene. ADAMTS13 is normally preactivated by conformational changes of the Metalloprotease (M) domain. Studying a novel congenital thrombotic thrombocytopenic purpura p.R102S mutation in the M domain, which results in undetectable ADAMTS13 activity in the patient, could help to explain the patients' phenotype and to elucidate the currently unclear mechanism of allosteric preactivation.

OBJECTIVES

To investigate the in vitro effect of p.R102S mutation on ADAMTS13 secretion, activity, and allosteric preactivation.

METHODS

Molecular modeling was used to study the effect of the mutation on the stability of ADAMTS13. Recombinant mutant ADAMTS13 was generated by transient and stable transfection of, respectively, CHO K1 and HEK293-T cells. ADAMTS13 antigen was measured in enzyme-linked immunosorbent assay. ADAMTS13 activity was measured in a FRETS-VWF73 assay. Allosteric preactivation was assessed in FRETS-VWF73 assay, using monoclonal antibody (mAb) 17G2 that normally induces a ∼2-fold increase in activity, and in enzyme-linked immunosorbent assay using mAb 6A6 recognizing a cryptic epitope in the M domain that becomes exposed after binding of 17G2.

RESULTS

p.R102S mutation destabilizes the interactions between the M and Disintegrin-like (D) domain. p.R102S mutant secretion was impaired (35% of wild type) and activity was severely reduced (12% of wild type). p.R102S mutant could still be activated and the cryptic epitope of 6A6 was still fully exposed by 17G2 addition.

CONCLUSION

p.R102S mutation destabilizes the M-D domain interactions, causing impaired ADAMTS13 secretion and activity, which explains the patients' phenotype. Allosteric preactivation of ADAMTS13 remains conserved in the presence of the p.R102S mutation.

Exploiting Ligand-binding Domain Dimerization for Development of Novel Androgen Receptor Inhibitors

Currently, all clinically used androgen receptor (AR) antagonists target the AR ligand-binding pocket and inhibit T and dihydrotestosterone (DHT) binding. Resistance to these inhibitors in prostate cancer frequently involves AR-dependent mechanisms resulting in a retained AR dependence of the tumor. More effective or alternative AR inhibitors are therefore required to limit progression in these resistant stages. Here, we applied the structural information of the ligand-binding domain (LBD) dimerization interface to screen in silico for inhibitors. A completely new binding site, the Dimerisation Inhibiting Molecules (DIM) pocket, was identified at the LBD dimerization interface. Selection of compounds that fit the DIM pocket via virtual screening identified the DIM20 family of compounds which inhibit AR transactivation and dimerization of the full-length AR as well as the isolated LBDs. Via biolayer interferometry, reversible dose-dependent binding to the LBD was confirmed. While DIM20 does not compete with 3H-DHT for binding in the LBP, it limits the maximal activity of the AR indicative of a noncompetitive binding to the LBD. DIM20 and DIM20.39 specifically inhibit proliferation of AR-positive prostate cancer cell lines, with only marginal effects on AR-negative cell lines such as HEK 293 and PC3. Moreover, combination treatment of DIM compounds with enzalutamide results in synergistic antiproliferative effects which underline the specific mechanism of action of the DIM compounds.

RF03 | PMON286 Engineering Steroid Receptors to Respond to Chemicals of Choice

In synthetic biology, small molecules are commonly used for altering gene expression through the (in)activation of gene switches. Current switches, however, have a narrow species range, are non-titratable, or use non-inert inducers which evoke important side effects. Universal gene switches (XRs) that have dynamic ranges without interference with the host cell physiology will overcome these limitations. Ideal starting points for developing such XRs are the steroid receptors (SRs). The SRs are frequently applied as regulator domain or gene switch in organisms in which they are not naturally expressed, like microbes and plants. Several SRs, such as the estrogen receptor alpha (ERα), have been redesigned to respond to synthetic effectors but these variants require inducer amounts that are detrimental to the host or affect endogenous SRs.

By combining in silico modelling with directed evolution in the yeast Saccharomyces cerevisiae, we developed a platform and strategy for evolving SRs towards a chemical of choice. Validation of this method was achieved by first evolving the ERα to employ tamoxifen as an agonist, yielding the preTERRA variant. Subsequently, activation by estradiol was removed, rendering a tamoxifen receptor (TERRA). We show that TERRA gets activated by subtoxic tamoxifen amounts while remaining insensitive for physiological estradiol levels. Further evaluation demonstrated TERRA to be transferable to mammalian cells. Ongoing research aims at converting the TERRA receptor into a regulatory domain for controlling nuclear translocation of Cre recombinases and Cas9 proteins in both cell lines and mouse models.

The same platform was employed for adapting ERα towards the plasticizer Bisphenol A. This yielded an ERα variant that is no longer inducible by estradiol but is responsive to lower concentrations of Bisphenol in yeast. The transferal to cell lines will provide us with an additional regulatory domain. The ability to engineer any SR towards a synthetic ligand of choice not only expands the synthetic biology toolbox but can find applications in landscaping SR-related drug resistance and ligand-inducibility.

Presentation: Saturday, June 11, 2022 1:48 p.m. - 1:53 p.m., Monday, June 13, 2022 12:30 p.m. - 2:30 p.m.

RF15 | PMON291 Crystal Structures of Androgen DNA Binding Domain Interacting with its Response Elements (C3(1) ARE and MMTV-GRE)

Prostate cancer (PCa) is the second leading cancer in the world in men. The major therapeutic target in advanced PCa is the ligand-binding pocket (LBP) of the androgen receptor (AR). The anti-androgens used to treat cancer are binding to this domain and compete with the natural agonist. Over time the cancer grows resistant towards this therapy. Nevertheless, the cancer still typically depends on the AR, indicating other sites of the AR may still serve as drug targets. The DNA binding domain (DBD) is a likely good therapeutic target as it is involved in essential processes such AR dimerization, translocation to the nucleus and binding to the androgen response elements (AREs). However the structure and the DNA binding interactions with AREs are not fully explored yet. Only one crystal structure (pdb: 1R4I) is available. In this structure, the DBD is bound to an artificial direct repeat DNA sequences (ADR3), resulting in an unexpected dimeric head-to-head arrangement. A better understanding on how the AR recognizes and binds to different native DNA sequences which are mainly organized as semi inverted repeats is crucial for rational drug design. Here, we performed protein crystallography of the DBD in complex with various natural ARE sequences, including the mouse mammary tumour virus (MMTV-GRE) motif. Surprisingly, this latter element is still able to be activated by a dimerization deficient AR.

Our first crystal structures were hallmarked by a high degree of disorder and appeared to bind both in a head-to-head as a head-to-tail fashion. Rational engineering of the DNA duplexes used for crystallography yielded well diffracting crystals for two different complexes. These two structures were diffracted to a range of 2 to 2.7 Å. Both structures with natural sequences appeared to bind in the head-to-head arrangement including the MMTV-GRE which consist out of semi inverted repeat DNA, whereas the C3(1)ARE contains the canonical fully inverted repeats. The structure reveals a slightly bend arrangement from the previously reported structure. The displacement is most visible at the DNA recognition helix bound to the half site and the so-called lever arm, which positions this helix in this region with a RMSD difference of 1.3 A. This indicates that the structural conformation changes of the protein is driven towards improved DNA-base pair recognition. Indeed the protein is able to make an additional set of non-covalent interactions with the base pairs of MMTV-GRE, leading to higher affinity. Both these extra interactions and conformational change lead to a higher affinity binding of protein-DNA, which seems to overcome the barrier introduced by the dimerization breaking mutations. This interaction surface could be a druggable site for structure-based drug design in the future, targeting AR-DBD dimerization.

Presentation: Sunday, June 12, 2022 12:30 p.m. - 12:35 p.m., Monday, June 13, 2022 12:30 p.m. - 2:30 p.m.

A 1-aminocyclopropane-1-carboxylic-acid (ACC) dipeptide elicits ethylene responses through ACC-oxidase mediated substrate promiscuity

Plants produce the volatile hormone ethylene to regulate many developmental processes and to deal with (a)biotic stressors. In seed plants, ethylene is synthesized from 1-aminocyclopropane-1-carboxylic acid (ACC) by the dedicated enzyme ACC oxidase (ACO). Ethylene biosynthesis is tightly regulated at the level of ACC through ACC synthesis, conjugation and transport. ACC is a non-proteinogenic amino acid, which also has signaling roles independent from ethylene. In this work, we investigated the biological function of an uncharacterized ACC dipeptide. The custom-synthesized di-ACC molecule can be taken up by Arabidopsis in a similar way as ACC, in part via Lysine Histidine Transporters (e.g., LHT1). Using Nano-Particle Assisted Laser Desoprtion/Ionization (Nano-PALDI) mass-spectrometry imaging, we revealed that externally fed di-ACC predominantly localizes to the vasculature tissue, despite it not being detectable in control hypocotyl segments. Once taken up, the ACC dimer can evoke a triple response phenotype in dark-grown seedlings, reminiscent of ethylene responses induced by ACC itself, albeit less efficiently compared to ACC. Di-ACC does not act via ACC-signaling, but operates via the known ethylene signaling pathway. In vitro ACO activity and molecular docking showed that di-ACC can be used as an alternative substrate by ACO to form ethylene. The promiscuous nature of ACO for the ACC dimer also explains the higher ethylene production rates observed in planta, although this reaction occurred less efficiently compared to ACC. Overall, the ACC dipeptide seems to be transported and converted into ethylene in a similar way as ACC, and is able to augment ethylene production levels and induce subsequent ethylene responses in Arabidopsis.

A disease-associated missense mutation in CYP4F3 affects the metabolism of leukotriene B4 via disruption of electron transfer

BACKGROUND

Cytochrome P450 4F3 (CYP4F3) is an ω-hydroxylase that oxidizes leukotriene B4 (LTB4), prostaglandins, and fatty acid epoxides. LTB4 is synthesized by leukocytes and acts as a chemoattractant for neutrophils, making it an essential component of the innate immune system. Recently, involvement of the LTB4 pathway was reported in various immunological disorders such as asthma, arthritis, and inflammatory bowel disease. We report a 26-year-old female with a complex immune phenotype, mainly marked by exhaustion, muscle weakness, and inflammation-related conditions. The molecular cause is unknown, and symptoms have been aggravating over the years.

METHODS

Whole exome sequencing was performed and validated; flow cytometry and enzyme-linked immunosorbent assay were used to describe patient's phenotype. Function and impact of the mutation were investigated using molecular analysis: co-immunoprecipitation, western blot, and enzyme-linked immunosorbent assay. Capillary electrophoresis with ultraviolet detection was used to detect LTB4 and its metabolite and in silico modelling provided structural information.

RESULTS

We present the first report of a patient with a heterozygous de novo missense mutation c.C1123 > G;p.L375V in CYP4F3 that severely impairs its activity by 50% (P < 0.0001), leading to reduced metabolization of the pro-inflammatory LTB4. Systemic LTB4 levels (1034.0 ± 75.9 pg/mL) are significantly increased compared with healthy subjects (305.6 ± 57.0 pg/mL, P < 0.001), and immune phenotyping shows increased total CD19+ CD27- naive B cells (25%) and decreased total CD19+ CD27+ IgD- switched memory B cells (19%). The mutant CYP4F3 protein is stable and binding with its electron donors POR and Cytb5 is unaffected (P > 0.9 for both co-immunoprecipitation with POR and Cytb5). In silico modelling of CYP4F3 in complex with POR and Cytb5 suggests that the loss of catalytic activity of the mutant CYP4F3 is explained by a disruption of an α-helix that is crucial for the electron shuffling between the electron carriers and CYP4F3. Interestingly, zileuton still inhibits ex vivo LTB4 production in patient's whole blood to 2% of control (P < 0.0001), while montelukast and fluticasone do not (99% and 114% of control, respectively).

CONCLUSIONS

A point mutation in the catalytic domain of CYP4F3 is associated with high leukotriene B4 plasma levels and features of a more naive adaptive immune response. Our data provide evidence for the pathogenicity of the CYP4F3 variant as a cause for the observed clinical features in the patient. Inhibitors of the LTB4 pathway such as zileuton show promising effects in blocking LTB4 production and might be used as a future treatment strategy.

Small-molecule profiling for steroid receptor activity using a universal steroid receptor reporter assay

A critical step in the development of novel drug candidates for the treatment of steroid related diseases is ensuring the absence of crosstalk with steroid receptors (SRs). Establishing this SR cross-reactivity profile requires multiple reporter assays as each SR associates with its unique enhancer region, a labor intensive and time-consuming approach. To overcome this need for multi-reporter assays, we established a steroid receptor inducible luciferase reporter assay (SRi-Luc) that allows side-by-side examination of agonistic and antagonistic properties of small-molecules on all steroid receptors. This state-of-the-art SRi-Luc consists of a unique alteration of four distinct keto-steroid- and estrogen response elements. As proof of principle, the SRi-Luc assay was used to profile a set of novel designed steroidal 1,2,3-triazoles. These triazolized steroidal compounds were developed via our in-house triazolization methodology, in which an enolizable ketone is converted into a triazolo-fused or -linked analog by treatment with a primary amine or ammonium salt in the presence of 4-nitrophenyl azide. From these designed steroidal 1,2,3-triazoles, six successfully reduced androgen receptor activity by 40 %. Although opted as antiandrogens, their cross-reactivity with other SRs was apparent in our SRi-Luc assay and rendered them unsuited for further antagonist development and clinical use. Overall, the SRi-Luc overcomes the need of multi-reporter assays for the profiling of small-molecules on all SRs. This not only reduces the risk of introducing biases, it as well accelerates early-stage drug discovery when designing particular SR selective (ant)agonists or characterizing off-target effects of lead molecules acting on any drug target.

The T850D Phosphomimetic Mutation in the Androgen Receptor Ligand Binding Domain Enhances Recruitment at Activation Function 2

Several key functions of the androgen receptor (AR) such as hormone recognition and co-regulator recruitment converge in the ligand binding domain (LBD). Loss- or gain-of-function of the AR contributes to pathologies such as the androgen insensitivity syndrome and prostate cancer. Here, we describe a gain-of-function mutation of the surface-exposed threonine at position 850, located at the amino-terminus of Helix 10 (H10) in the AR LBD. Since T850 phosphorylation was reported to affect AR function, we created the phosphomimetic mutation T850D. The AR T850D variant has a 1.5- to 2-fold increased transcriptional activity with no effect on ligand affinity. In the androgen responsive LNCaP cell line grown in medium with low androgen levels, we observed a growth advantage for cells in which the endogenous AR was replaced by AR T850D. Despite the distance to the AF2 site, the AR T850D LBD displayed an increased affinity for coactivator peptides as well as the ²³FQNLF²⁷ motif of AR itself. Molecular Dynamics simulations confirm allosteric transmission of the T850D mutation towards the AF2 site via extended hydrogen bond formation between coactivator peptide and AF2 site. This mechanistic study thus confirms the gain-of-function character of T850D and T850 phosphorylation for AR activity and reveals details of the allosteric communications within the LBD.

The androgen receptor depends on ligand-binding domain dimerization for transcriptional activation

Whereas dimerization of the DNA-binding domain of the androgen receptor (AR) plays an evident role in recognizing bipartite response elements, the contribution of the dimerization of the ligand-binding domain (LBD) to the correct functioning of the AR remains unclear. Here, we describe a mouse model with disrupted dimerization of the AR LBD (AR^Lmon/Y ). The disruptive effect of the mutation is demonstrated by the feminized phenotype, absence of male accessory sex glands, and strongly affected spermatogenesis, despite high circulating levels of testosterone. Testosterone replacement studies in orchidectomized mice demonstrate that androgen-regulated transcriptomes in AR^Lmon/Y mice are completely lost. The mutated AR still translocates to the nucleus and binds chromatin, but does not bind to specific AR binding sites. In vitro studies reveal that the mutation in the LBD dimer interface also affects other AR functions such as DNA binding, ligand binding, and co-regulator binding. In conclusion, LBD dimerization is crucial for the development of AR-dependent tissues through its role in transcriptional regulation in vivo. Our findings identify AR LBD dimerization as a possible target for AR inhibition.

Betulonic Acid Derivatives Interfering with Human Coronavirus 229E Replication via the nsp15 Endoribonuclease

To develop antiviral therapeutics against human coronavirus (HCoV) infections, suitable coronavirus drug targets and corresponding lead molecules must be urgently identified. Here, we describe the discovery of a class of HCoV inhibitors acting on nsp15, a hexameric protein component of the viral replication-transcription complexes, endowed with immune evasion-associated endoribonuclease activity. Structure-activity relationship exploration of these 1,2,3-triazolo-fused betulonic acid derivatives yielded lead molecule 5h as a strong inhibitor (antiviral EC50: 0.6 μM) of HCoV-229E replication. An nsp15 endoribonuclease active site mutant virus was markedly less sensitive to 5h, and selected resistance to the compound mapped to mutations in the N-terminal part of HCoV-229E nsp15, at an interface between two nsp15 monomers. The biological findings were substantiated by the nsp15 binding mode for 5h, predicted by docking. Hence, besides delivering a distinct class of inhibitors, our study revealed a druggable pocket in the nsp15 hexamer with relevance for anti-coronavirus drug development.

Repurposing the Antidepressant Sertraline as SHMT Inhibitor to Suppress Serine/Glycine Synthesis-Addicted Breast Tumor Growth

Metabolic rewiring is a hallmark of cancer that supports tumor growth, survival, and chemotherapy resistance. Although normal cells often rely on extracellular serine and glycine supply, a significant subset of cancers becomes addicted to intracellular serine/glycine synthesis, offering an attractive drug target. Previously developed inhibitors of serine/glycine synthesis enzymes did not reach clinical trials due to unfavorable pharmacokinetic profiles, implying that further efforts to identify clinically applicable drugs targeting this pathway are required. In this study, we aimed to develop therapies that can rapidly enter the clinical practice by focusing on drug repurposing, as their safety and cost-effectiveness have been optimized before. Using a yeast model system, we repurposed two compounds, sertraline and thimerosal, for their selective toxicity against serine/glycine synthesis-addicted breast cancer and T-cell acute lymphoblastic leukemia cell lines. Isotope tracer metabolomics, computational docking, enzymatic assays, and drug-target interaction studies revealed that sertraline and thimerosal inhibit serine/glycine synthesis enzymes serine hydroxymethyltransferase and phosphoglycerate dehydrogenase, respectively. In addition, we demonstrated that sertraline's antiproliferative activity was further aggravated by mitochondrial inhibitors, such as the antimalarial artemether, by causing G1-S cell-cycle arrest. Most notably, this combination also resulted in serine-selective antitumor activity in breast cancer mouse xenografts. Collectively, this study provides molecular insights into the repurposed mode-of-action of the antidepressant sertraline and allows to delineate a hitherto unidentified group of cancers being particularly sensitive to treatment with sertraline. Furthermore, we highlight the simultaneous inhibition of serine/glycine synthesis and mitochondrial metabolism as a novel treatment strategy for serine/glycine synthesis-addicted cancers.

FraC nanopores with adjustable diameter identify the mass of opposite-charge peptides with 44 dalton resolution

A high throughput single-molecule method for identifying peptides and sequencing proteins based on nanopores could reduce costs and increase speeds of sequencing, allow the fabrication of portable home-diagnostic devices, and permit the characterization of low abundance proteins and heterogeneity in post-translational modifications. Here we engineer the size of Fragaceatoxin C (FraC) biological nanopore to allow the analysis of a wide range of peptide lengths. Ionic blockades through engineered nanopores distinguish a variety of peptides, including two peptides differing only by the substitution of alanine with glutamate. We also find that at pH 3.8 the depth of the peptide current blockades scales with the mass of the peptides irrespectively of the chemical composition of the analyte. Hence, this work shows that FraC nanopores allow direct readout of the mass of single peptide in solution, which is a crucial step towards the developing of a real-time and single-molecule protein sequencing device.

Design of tryptophan-containing mutants of the symmetrical Pizza protein for biophysical studies

β-propeller proteins are highly symmetrical, being composed of a repeated motif with four anti-parallel β-sheets arranged around a central axis. Recently we designed the first completely symmetrical β-propeller protein, Pizza6, consisting of six identical tandem repeats. Pizza6 is expected to prove a useful building block for bionanotechnology, and also a tool to investigate the folding and evolution of β-propeller proteins. Folding studies are made difficult by the high stability and the lack of buried Trp residues to act as monitor fluorophores, so we have designed and characterized several Trp-containing Pizza6 derivatives. In total four proteins were designed, of which three could be purified and characterized. Crystal structures confirm these mutant proteins maintain the expected structure, and a clear redshift of Trp fluorescence emission could be observed upon denaturation. Among the derivative proteins, Pizza6-AYW appears to be the most suitable model protein for future folding/unfolding kinetics studies as it has a comparable stability as natural β-propeller proteins.

Fragger: a protein fragment picker for structural queries

Protein modeling and design activities often require querying the Protein Data Bank (PDB) with a structural fragment, possibly containing gaps. For some applications, it is preferable to work on a specific subset of the PDB or with unpublished structures. These requirements, along with specific user needs, motivated the creation of a new software to manage and

query 3D protein fragments. Fragger is a protein fragment picker that allows protein fragment databases to be created and queried. All fragment lengths are supported and any set of PDB files can be used to create a database. Fragger can efficiently search a fragment database with a query fragment and a distance threshold. Matching fragments are ranked by distance to the query. The query fragment can have structural gaps and the allowed amino acid sequences matching a query can be constrained via a regular expression of one-letter amino acid codes. Fragger also incorporates a tool to compute the backbone RMSD of one versus many fragments in high throughput. Fragger should be useful for protein design, loop grafting and related structural

bioinformatics tasks.

Computational Studies of the Active and Inactive Regulatory Domains of Response Regulator PhoP Using Molecular Dynamics Simulations

The response regulator PhoP is part of the PhoP/PhoQ two-component system, which is responsible for regulating the expression of multiple genes involved in controlling virulence, biofilm formation, and resistance to antimicrobial peptides. Therefore, modulating the transcriptional function of the PhoP protein is a promising strategy for developing new antimicrobial agents. There is evidence suggesting that phosphorylation-mediated dimerization in the regulatory domain of PhoP is essential for its transcriptional function. Disruption or stabilization of protein-protein interactions at the dimerization interface may inhibit or enhance the expression of PhoP-dependent genes. In this study, we performed molecular dynamics simulations on the active and inactive dimers and monomers of the PhoP regulatory domains, followed by pocket-detecting screenings and a quantitative hot-spot analysis in order to assess the druggability of the protein. Consistent with prior hypothesis, the calculation of the binding free energy shows that phosphorylation enhances dimerization of PhoP. Furthermore, we have identified two different putative binding sites at the dimerization active site (the α4-β5-α5 face) with energetic "hot-spot" areas, which could be used to search for modulators of protein-protein interactions. This study delivers insight into the dynamics and druggability of the dimerization interface of the PhoP regulatory domain, and may serve as a basis for the rational identification of new antimicrobial drugs.

An integrated map of HIV genome-wide variation from a population perspective

Background

The HIV pandemic is characterized by extensive genetic variability, which has challenged the development of HIV drugs and vaccines. Although HIV genomes have been classified into different types, groups, subtypes and recombinants, a comprehensive study that maps HIV genome-wide diversity at the population level is still lacking to date. This study aims to characterize HIV genomic diversity in large-scale sequence populations, and to identify driving factors that shape HIV genome diversity.

Results

A total of 2996 full-length genomic sequences from 1705 patients infected with 16 major HIV groups, subtypes and circulating recombinant forms (CRFs) were analyzed along with structural, immunological and peptide inhibitor information. Average nucleotide diversity of HIV genomes was almost 50% between HIV-1 and HIV-2 types, 37.5% between HIV-1 groups, 14.7% between HIV-1 subtypes, 8.2% within individual HIV-1 subtypes and less than 1% within single patients. Along the HIV genome, diversity patterns and compositions of nucleotides and amino acids were highly similar across different groups, subtypes and CRFs. Current HIV-derived peptide inhibitors were predominantly derived from conserved, solvent accessible and intrinsically ordered structures in the HIV-1 subtype B genome. We identified these conserved regions in Capsid, Nucleocapsid, Protease, Integrase, Reverse transcriptase, Vpr and the GP41 N terminus as potential drug targets. In the analysis of factors that impact HIV-1 genomic diversity, we focused on protein multimerization, immunological constraints and HIV-human protein interactions. We found that amino acid diversity in monomeric proteins was higher than in multimeric proteins, and diversified positions were preferably located within human CD4 T cell and antibody epitopes. Moreover, intrinsic disorder regions in HIV-1 proteins coincided with high levels of amino acid diversity, facilitating a large number of interactions between HIV-1 and human proteins.

Conclusions

This first large-scale analysis provided a detailed mapping of HIV genomic diversity and highlighted drug-target regions conserved across different groups, subtypes and CRFs. Our findings suggest that, in addition to the impact of protein multimerization and immune selective pressure on HIV-1 diversity, HIV-human protein interactions are facilitated by high variability within intrinsically disordered structures.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-015-0148-6) contains supplementary material, which is available to authorized users.

Duplication of a promiscuous transcription factor drives the emergence of a new regulatory network

The emergence of new genes throughout evolution requires rewiring and extension of regulatory networks. However, the molecular details of how the transcriptional regulation of new gene copies evolves remain largely unexplored. Here we show how duplication of a transcription factor gene allowed the emergence of two independent regulatory circuits. Interestingly, the ancestral transcription factor was promiscuous and could bind different motifs in its target promoters. After duplication, one paralogue evolved increased binding specificity so that it only binds one type of motif, whereas the other copy evolved a decreased activity so that it only activates promoters that contain multiple binding sites. Interestingly, only a few mutations in both the DNA-binding domains and in the promoter binding sites were required to gradually disentangle the two networks. These results reveal how duplication of a promiscuous transcription factor followed by concerted cis and trans mutations allows expansion of a regulatory network.

Androgen receptor antagonists for prostate cancer therapy

Androgen deprivation is the mainstay therapy for metastatic prostate cancer (PCa). Another way of suppressing androgen receptor (AR) signaling is via AR antagonists or antiandrogens. Despite being frequently prescribed in clinical practice, there is conflicting evidence concerning the role of AR antagonists in the management of PCa. In the castration-resistant settings of PCa, docetaxel has been the only treatment option for decades. With recent evidence that castration-resistant PCa is far from AR-independent, there has been an increasing interest in developing new AR antagonists. This review gives a concise overview of the clinically available antiandrogens and the experimental AR antagonists that tackle androgen action with a different approach.

In silico, in vitro and cellular analysis with a kinome-wide inhibitor panel correlates cellular LRRK2 dephosphorylation to inhibitor activity on LRRK2

Leucine-rich repeat kinase 2 (LRRK2) is a complex, multidomain protein which is considered a valuable target for potential disease-modifying therapeutic strategies for Parkinson's disease (PD). In mammalian cells and brain, LRRK2 is phosphorylated and treatment of cells with inhibitors of LRRK2 kinase activity can induce LRRK2 dephosphorylation at a cluster of serines including Ser910/935/955/973. It has been suggested that phosphorylation levels at these sites reflect LRRK2 kinase activity, however kinase-dead variants of LRRK2 or kinase activating variants do not display altered Ser935 phosphorylation levels compared to wild type. Furthermore, Ser910/935/955/973 are not autophosphorylation sites, therefore, it is unclear if inhibitor induced dephosphorylation depends on the activity of compounds on LRRK2 or on yet to be identified upstream kinases. Here we used a panel of 160 ATP competitive and cell permeable kinase inhibitors directed against all branches of the kinome and tested their activity on LRRK2 in vitro using a peptide-substrate-based kinase assay. In neuronal SH-SY5Y cells overexpressing LRRK2 we used compound-induced dephosphorylation of Ser935 as readout. In silico docking of selected compounds was performed using a modeled LRRK2 kinase structure. Receiver operating characteristic plots demonstrated that the obtained docking scores to the LRRK2 ATP binding site correlated with in vitro and cellular compound activity. We also found that in vitro potency showed a high degree of correlation to cellular compound induced LRRK2 dephosphorylation activity across multiple compound classes. Therefore, acute LRRK2 dephosphorylation at Ser935 in inhibitor treated cells involves a strong component of inhibitor activity on LRRK2 itself, without excluding a role for upstream kinases. Understanding the regulation of LRRK2 phosphorylation by kinase inhibitors aids our understanding of LRRK2 signaling and may lead to development of new classes of LRRK2 kinase inhibitors.

A rotation-translation invariant molecular descriptor of partial charges and its use in ligand-based virtual screening

Background

Measures of similarity for chemical molecules have been developed since the dawn of chemoinformatics. Molecular similarity has been measured by a variety of methods including molecular descriptor based similarity, common molecular fragments, graph matching and 3D methods such as shape matching. Similarity measures are widespread in practice and have proven to be useful in drug discovery. Because of our interest in electrostatics and high throughput ligand-based virtual screening, we sought to exploit the information contained in atomic coordinates and partial charges of a molecule.

Results

A new molecular descriptor based on partial charges is proposed. It uses the autocorrelation function and linear binning to encode all atoms of a molecule into two rotation-translation invariant vectors. Combined with a scoring function, the descriptor allows to rank-order a database of compounds versus a query molecule. The proposed implementation is called ACPC (AutoCorrelation of Partial Charges) and released in open source. Extensive retrospective ligand-based virtual screening experiments were performed and other methods were compared with in order to validate the method and associated protocol.

Conclusions

While it is a simple method, it performed remarkably well in experiments. At an average speed of 1649 molecules per second, it reached an average median area under the curve of 0.81 on 40 different targets; hence validating the proposed protocol and implementation.

Functional conservation of HIV-1 Gag: implications for rational drug design

Background

HIV-1 replication can be successfully blocked by targeting gag gene products, offering a promising strategy for new drug classes that complement current HIV-1 treatment options. However, naturally occurring polymorphisms at drug binding sites can severely compromise HIV-1 susceptibility to gag inhibitors in clinical and experimental studies. Therefore, a comprehensive understanding of gag natural diversity is needed.

Findings

We analyzed the degree of functional conservation in 10862 full-length gag sequences across 8 major HIV-1 subtypes and identified the impact of natural variation on known drug binding positions targeted by more than 20 gag inhibitors published to date. Complete conservation across all subtypes was detected in 147 (29%) out of 500 gag positions, with the highest level of conservation observed in capsid protein. Almost half (41%) of the 136 known drug binding positions were completely conserved, but all inhibitors were confronted with naturally occurring polymorphisms in their binding sites, some of which correlated with HIV-1 subtype. Integration of sequence and structural information revealed one drug binding pocket with minimal genetic variability, which is situated at the N-terminal domain of the capsid protein.

Conclusions

This first large-scale analysis of full-length HIV-1 gag provided a detailed mapping of natural diversity across major subtypes and highlighted the considerable variation in current drug binding sites. Our results contribute to the optimization of gag inhibitors in rational drug design, given that drug binding sites should ideally be conserved across all HIV-1 subtypes.

Fragment based drug design: from experimental to computational approaches

Fragment based drug design has emerged as an effective alternative to high throughput screening for the identification of lead compounds in drug discovery in the past fifteen years. Fragment based screening and optimization methods have achieved credible success in many drug discovery projects with one approved drug and many more compounds in clinical trials. The fragment based drug design starts with the identification of fragments or low molecular weight compounds that generally bind with weak affinity to the target of interest. The fragments that form high quality interactions are then optimized to lead compounds with high affinity and selectivity. The weak affinity of fragments for their target requires the use of biophysical techniques such as nuclear magnetic resonance, X-ray crystallography or surface plasmon resonance to identify hits. These techniques are very sensitive and some of them provide detailed protein fragment interaction information that is important for fragment to lead optimization. Despite the huge advances in technology in the past years, experimental methods of fragment screening suffer several challenges such as low throughput, high cost of instruments and experiments, high protein and fragment concentration requirements. To address challenges posed by experimental screening approaches, computational methods were developed that play an important role in fragment library design, fragment screening and optimization of initial fragment hits. The computational approaches of fragment screening and optimization are most useful when they are used in combination with experimental approaches. The use of virtual fragment based screening in combination with experimental methods has fostered the application of fragment based drug design to important biological targets including protein-protein interactions and membrane proteins such as GPCRs. This review provides an overview of experimental and computational screening approaches used in fragment based drug discovery with an emphasis on recent successes achieved in discovering potent lead molecules using these approaches.

Pharmacophore modelling as a virtual screening tool for the discovery of small molecule protein-protein interaction inhibitors

Pharmacophore searches have become a popular tool for virtual screening of libraries to identify novel active substances that can be potentially developed into drugs. While they have been applied for years on common drug targets, their application in the discovery of protein-protein interaction inhibitors remains limited. This review describes current pharmacophore modelling methods applied in the discovery of novel inhibitors targeting protein-protein interactions. We first address the mimicry of protein-protein interactions with their respective inhibitors as observed in crystal structure complexes. This mimicry can be exploited to derive a pharmacophore query from protein-protein complex structures. We then discuss several cases where pharmacophore queries were utilized for the discovery of first-in-class inhibitors of their respective protein-protein interaction targets. These examples have demonstrated the usefulness of pharmacophore modelling in the quest for protein-protein interaction inhibitors.

Mapping the binding interface between an HIV-1 inhibiting intrabody and the viral protein Rev

HIV-1 Rev is the key protein in the nucleocytoplasmic export and expression of the late viral mRNAs. An important aspect for its function is its ability to multimerize on these mRNAs. We have recently identified a llama single-domain antibody (Nb190) as the first inhibitor targeting the Rev multimerization function in cells. This nanobody is a potent intracellular antibody that efficiently inhibits HIV-1 viral production. In order to gain insight into the Nb190-Rev interaction interface, we performed mutational and docking studies to map the interface between the nanobody paratope and the Rev epitope. Alanine mutants of the hyper-variable domains of Nb190 and the Rev multimerization domains were evaluated in different assays measuring Nb190-Rev interaction or viral production. Seven residues within Nb190 and five Rev residues are demonstrated to be crucial for epitope recognition. These experimental data were used to perform docking experiments and map the Nb190-Rev structural interface. This Nb190-Rev interaction model can guide further studies of the Nb190 effect on HIV-1 Rev function and could serve as starting point for the rational development of smaller entities binding to the Nb190 epitope, aimed at interfering with protein-protein interactions of the Rev N-terminal domain.

An integrated fragment based screening approach for the discovery of small molecule modulators of the VWF-GPIbα interaction

An integrated approach comprising STD NMR screening, pharmacophore based analogue selection and a bioassay is presented for the discovery of a stabilizer of the clinically relevant VWF-GPIbα protein-protein interaction.

Reconstruction of ancestral metabolic enzymes reveals molecular mechanisms underlying evolutionary innovation through gene duplication

Resurrection of ancient fungal maltase enzymes uncovers the molecular details of how repeated gene duplications allow the evolution of protein variants with different functions.

Gene duplications are believed to facilitate evolutionary innovation. However, the mechanisms shaping the fate of duplicated genes remain heavily debated because the molecular processes and evolutionary forces involved are difficult to reconstruct. Here, we study a large family of fungal glucosidase genes that underwent several duplication events. We reconstruct all key ancestral enzymes and show that the very first preduplication enzyme was primarily active on maltose-like substrates, with trace activity for isomaltose-like sugars. Structural analysis and activity measurements on resurrected and present-day enzymes suggest that both activities cannot be fully optimized in a single enzyme. However, gene duplications repeatedly spawned daughter genes in which mutations optimized either isomaltase or maltase activity. Interestingly, similar shifts in enzyme activity were reached multiple times via different evolutionary routes. Together, our results provide a detailed picture of the molecular mechanisms that drove divergence of these duplicated enzymes and show that whereas the classic models of dosage, sub-, and neofunctionalization are helpful to conceptualize the implications of gene duplication, the three mechanisms co-occur and intertwine.

Darwin's theory of evolution is one of gradual change, yet evolution sometimes takes remarkable leaps. Such evolutionary innovations are often linked to gene duplication through one of three basic scenarios: an extra copy can increase protein levels, different ancestral subfunctions can be split over the copies and evolve distinct regulation, or one of the duplicates can develop a novel function. Although there are numerous examples for all these trajectories, the underlying molecular mechanisms remain obscure, mostly because the preduplication genes and proteins no longer exist. Here, we study a family of fungal metabolic enzymes that hydrolyze disaccharides, and that all originated from the same ancestral gene through repeated duplications. By resurrecting the ancient genes and proteins using high-confidence predictions from many fungal genome sequences available, we show that the very first preduplication enzyme was promiscuous, preferring maltose-like substrates but also showing trace activity towards isomaltose-like sugars. After duplication, specific mutations near the active site of one copy optimized the minor activity at the expense of the major ancestral activity, while the other copy further specialized in maltose and lost the minor activity. Together, our results reveal how the three basic trajectories for gene duplicates cannot be separated easily, but instead intertwine into a complex evolutionary path that leads to innovation.

Identification and characterization of MEL-3, a novel AR antagonist that suppresses prostate cancer cell growth

Antiandrogens are an important component of prostate cancer therapy as the androgen receptor (AR) is the key regulator of prostate cancer growth and survival. Current AR antagonists, such as bicalutamide and hydroxyflutamide, have a low affinity for the AR and as a result block AR signaling insufficiently. Moreover, many patients develop a resistance for bicalutamide or hydroxyflutamide during therapy or show a clinical improvement after withdrawal of the antiandrogen. New and more effective AR antagonists are needed to ensure follow-up of these patients. We therefore developed a screening system to identify novel AR antagonists from a collection of compounds. MEL-3 [8-(propan-2-yl)-5,6-dihydro-4H-pyrazino[3,2,1-jk]carbazole] was selected as potent inhibitor of the AR and was further characterized in vitro. On different prostate cancer cell lines MEL-3 displayed an improved therapeutic profile compared with bicalutamide. Not only cell growth was inhibited but also the expression of androgen-regulated genes: PSA and FKBP5. Prostate cancer is often associated with mutated ARs that respond to a broadened spectrum of ligands including the current antiandrogens used in the clinic, hydroxyflutamide and bicalutamide. The activity of two mutant receptors (AR T877A and AR W741C) was shown to be reduced in presence of MEL-3, providing evidence that MEL-3 can potentially be a follow-up treatment for bicalutamide- and hydroxyflutamide-resistant patients. The mechanism of action of MEL-3 on the molecular level was further explored by comparing the structure-activity relationship of different chemical derivatives of MEL-3 with the in silico docking of MEL-3 derivatives in the binding pocket of the AR.

Computational investigation of the HIV-1 Rev multimerization using molecular dynamics simulations and binding free energy calculations

The HIV Rev protein mediates the nuclear export of viral mRNA, and is thereby essential for the production of late viral proteins in the replication cycle. Rev forms a large organized multimeric protein-protein complex for proper functioning. Recently, the three-dimensional structures of a Rev dimer and tetramer have been resolved and provide the basis for a thorough structural analysis of the binding interaction. Here, molecular dynamics (MD) and binding free energy calculations were performed to elucidate the forces thriving dimerization and higher order multimerization of the Rev protein. It is found that despite the structural differences between each crystal structure, both display a similar behavior according to our calculations. Our analysis based on a molecular mechanics-generalized Born surface area (MM/GBSA) and a configurational entropy approach demonstrates that the higher order multimerization site is much weaker than the dimerization site. In addition, a quantitative hot spot analysis combined with a mutational analysis reveals the most contributing amino acid residues for protein interactions in agreement with experimental results. Additional residues were found in each interface, which are important for the protein interaction. The investigation of the thermodynamics of the Rev multimerization interactions performed here could be a further step in the development of novel antiretrovirals using structure based drug design. Moreover, the variability of the angle between each Rev monomer as measured during the MD simulations suggests a role of the Rev protein in allowing flexibility of the arginine rich domain (ARM) to accommodate RNA binding.