Reshaping Echinocandin Antifungal Drugs To Circumvent Glucan Synthase Point-Mutation-Mediated Resistance

Echinocandins are a class of antifungal drugs that inhibit the activity of the β-(1,3)-glucan synthase complex, which synthesizes fungal cell wall β-(1,3)-glucan. Echinocandin resistance is linked to mutations in the FKS gene, which encodes the catalytic subunit of the glucan synthase complex. We present a molecular-docking-based model that provides insight into how echinocandins interact with the target Fks protein: echinocandins form a ternary complex with both Fks and membrane lipids. We used reductive dehydration of alcohols to generate dehydroxylated echinocandin derivatives and evaluated their potency against a panel of Candida pathogens constructed by introducing resistance-conferring mutations in the FKS gene. We found that removing the hemiaminal alcohol, which drives significant conformational alterations in the modified echinocandins, reduced their efficacy. Conversely, eliminating the benzylic alcohol of echinocandins enhanced potency by up to two orders of magnitude, in a manner dependent upon the resistance-conferring mutation. Strains that have developed resistance to either rezafungin, the most recently clinically approved echinocandin, or its dehydroxylated derivative RZF-1, exhibit high resistance to rezafungin while demonstrating moderate resistance to RZF-1. These findings provide valuable insight for combating echinocandin resistance through chemical modifications.

Identification and characterization of the key enzyme in the biosynthesis of the neurotoxin β-ODAP in grass pea

Grass pea (Lathyrus sativus L.) is a grain legume commonly grown in Asia and Africa for food and forage. It is a highly nutritious and robust crop, capable of surviving both droughts and floods. However, it produces a neurotoxic compound, β-N-oxalyl-L-α,β-diaminopropionic acid (β-ODAP), which can cause a severe neurological disorder when consumed as a primary diet component. While the catalytic activity associated with β-ODAP formation was demonstrated more than 50 years ago, the enzyme responsible for this activity has not been identified. Here, we report on the identity, activity, 3D structure, and phylogenesis of this enzyme—β-ODAP synthase (BOS). We show that BOS belongs to the benzylalcohol O-acetyltransferase, anthocyanin O-hydroxycinnamoyltransferase, anthranilate N-hydroxycinnamoyl/benzoyltransferase, deacetylvindoline 4-O-acetyltransferase superfamily of acyltransferases and is structurally similar to hydroxycinnamoyl transferase. Using molecular docking, we propose a mechanism for its catalytic activity, and using heterologous expression in tobacco leaves (Nicotiana benthamiana), we demonstrate that expression of BOS in the presence of its substrates is sufficient for β-ODAP production in vivo. The identification of BOS may pave the way toward engineering β-ODAP–free grass pea cultivars, which are safe for human and animal consumption.

Phylogeny and abiotic conditions shape the diel floral emission patterns of desert Brassicaceae species

A key facet of floral scent is diel fluctuations in emission, often studied in the context of plant-pollinator interactions, while contributions of environment and phylogeny remain overlooked. Here, we ask if these factors are involved in shaping temporal variations in scent emission. To that end, we coupled light/dark floral emission measurements of 17 desert Brassicaceae species with environmental and phylogenetic data to explore the individual/combined impacts of these predictors on diel emission patterns. We further investigated these patterns by conducting high-resolution emission measurements in a subset of genetically distant species with contrasting temporal dynamics. While diel shifts in magnitude and richness of emission were strongly affected by genetic relatedness, they also reflect the environmental conditions under which the species grow. Specifically, light/dark emission ratios were negatively affected by an increase in winter temperatures, known to impact both plant physiology and insect locomotion, and sandy soil fractions, previously shown to exert stress that tempers with diel metabolic rhythms. Additionally, the biosynthetic origins of the compounds were associated with their corresponding production patterns, possibly to maximize emission efficacy. Using a multidisciplinary chemical/ecological approach, we uncover and differentiate the main factors shaping floral scent diel fluctuations, highlighting their consequences under changing global climate.

Conserved interactions required for inhibition of the main protease of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)

The COVID-19 pandemic caused by the SARS-CoV-2 requires a fast development of antiviral drugs. SARS-CoV-2 viral main protease (Mpro, also called 3C-like protease, 3CLpro) is a potential target for drug design. Crystal and co-crystal structures of the SARS-CoV-2 Mpro have been solved, enabling the rational design of inhibitory compounds. In this study we analyzed the available SARS-CoV-2 and the highly similar SARS-CoV-1 crystal structures. We identified within the active site of the Mpro, in addition to the inhibitory ligands' interaction with the catalytic C145, two key H-bond interactions with the conserved H163 and E166 residues. Both H-bond interactions are present in almost all co-crystals and are likely to occur also during the viral polypeptide cleavage process as suggested from docking of the Mpro cleavage recognition sequence. We screened in silico a library of 6900 FDA-approved drugs (ChEMBL) and filtered using these key interactions and selected 29 non-covalent compounds predicted to bind to the protease. Additional screen, using DOCKovalent was carried out on DrugBank library (11,414 experimental and approved drugs) and resulted in 6 covalent compounds. The selected compounds from both screens were tested in vitro by a protease activity inhibition assay. Two compounds showed activity at the 50 µM concentration range. Our analysis and findings can facilitate and focus the development of highly potent inhibitors against SARS-CoV-2 infection.

The mitochondrial carrier Citrin plays a role in regulating cellular energy during carcinogenesis

Citrin, encoded by SLC25A13 gene, is an inner mitochondrial transporter that is part of the malate-aspartate shuttle, which regulates the NAD+/NADH ratio between the cytosol and mitochondria. Citrullinemia type II (CTLN-II) is an inherited disorder caused by germline mutations in SLC25A13, manifesting clinically in growth failure that can be alleviated by dietary restriction of carbohydrates. The association of citrin with glycolysis and NAD+/NADH ratio led us to hypothesize that it may play a role in carcinogenesis. Indeed, we find that citrin is upregulated in multiple cancer types and is essential for supplementing NAD+ for glycolysis and NADH for oxidative phosphorylation. Consequently, citrin deficiency associates with autophagy, whereas its overexpression in cancer cells increases energy production and cancer invasion. Furthermore, based on the human deleterious mutations in citrin, we found a potential inhibitor of citrin that restricts cancerous phenotypes in cells. Collectively, our findings suggest that targeting citrin may be of benefit for cancer therapy.

A computational combinatorial approach identifies a protein inhibitor of superoxide dismutase 1 misfolding, aggregation, and cytotoxicity

Molecular agents that specifically bind and neutralize misfolded and toxic superoxide dismutase 1 (SOD1) mutant proteins may find application in attenuating the disease progression of familial amyotrophic lateral sclerosis. However, high structural similarities between the wild-type and mutant SOD1 proteins limit the utility of this approach. Here we addressed this challenge by converting a promiscuous natural human IgG-binding domain, the hyperthermophilic variant of protein G (HTB1), into a highly specific aggregation inhibitor (designated HTB1_M) of two familial amyotrophic lateral sclerosis-linked SOD1 mutants, SOD1^G93A and SOD1^G85R We utilized a computational algorithm for mapping protein surfaces predisposed to HTB1 intermolecular interactions to construct a focused HTB1 library, complemented with an experimental platform based on yeast surface display for affinity and specificity screening. HTB1_M displayed high binding specificity toward SOD1 mutants, inhibited their amyloid aggregation in vitro, prevented the accumulation of misfolded proteins in living cells, and reduced the cytotoxicity of SOD1^G93A expressed in motor neuron-like cells. Competition assays and molecular docking simulations suggested that HTB1_M binds to SOD1 via both its α-helical and β-sheet domains at the native dimer interface that becomes exposed upon mutated SOD1 misfolding and monomerization. Our results demonstrate the utility of computational mapping of the protein-protein interaction potential for designing focused protein libraries to be used in directed evolution. They also provide new insight into the mechanism of conversion of broad-spectrum immunoglobulin-binding proteins, such as HTB1, into target-specific proteins, thereby paving the way for the development of new selective drugs targeting the amyloidogenic proteins implicated in a variety of human diseases.

Prediction of homoprotein and heteroprotein complexes by protein docking and template-based modeling: A CASP-CAPRI experiment

We present the results for CAPRI Round 30, the first joint CASP-CAPRI experiment, which brought together experts from the protein structure prediction and protein-protein docking communities. The Round comprised 25 targets from amongst those submitted for the CASP11 prediction experiment of 2014. The targets included mostly homodimers, a few homotetramers, and two heterodimers, and comprised protein chains that could readily be modeled using templates from the Protein Data Bank. On average 24 CAPRI groups and 7 CASP groups submitted docking predictions for each target, and 12 CAPRI groups per target participated in the CAPRI scoring experiment. In total more than 9500 models were assessed against the 3D structures of the corresponding target complexes. Results show that the prediction of homodimer assemblies by homology modeling techniques and docking calculations is quite successful for targets featuring large enough subunit interfaces to represent stable associations. Targets with ambiguous or inaccurate oligomeric state assignments, often featuring crystal contact-sized interfaces, represented a confounding factor. For those, a much poorer prediction performance was achieved, while nonetheless often providing helpful clues on the correct oligomeric state of the protein. The prediction performance was very poor for genuine tetrameric targets, where the inaccuracy of the homology-built subunit models and the smaller pair-wise interfaces severely limited the ability to derive the correct assembly mode. Our analysis also shows that docking procedures tend to perform better than standard homology modeling techniques and that highly accurate models of the protein components are not always required to identify their association modes with acceptable accuracy. Proteins 2016; 84(Suppl 1):323-348. © 2016 Wiley Periodicals, Inc.

1-(sulfonyl)-5-(arylsulfonyl)indoline as activators of the tumor cell specific M2 isoform of pyruvate kinase

Cancer cells preferentially use glycolysis rather than oxidative phosphorylation for their rapid growth. They consume large amount of glucose to produce lactate even when oxygen is abundant, a phenomenon known as the Warburg effect. This metabolic change originates from a shift in the expression of alternative spliced isoforms of the glycolytic enzyme pyruvate kinase (PK), from PKM1 to PKM2. While PKM1 is constitutively active, PKM2 is switched from an inactive dimer form to an active tetramer form by small molecule activators. The prevalence of PKM2 in cancer cells relative to the prevalence of PKM1 in many normal cells, suggests a therapeutic strategy whereby activation of PKM2 may counter the abnormal cellular metabolism in cancer cells, and consequently decreased cellular proliferation. Herein we describe the discovery and optimization of a series of PKM2 activators derived from the 2-((2,3-dihydrobenzo[b][1,4] dioxin-6-yl)thio)-1-(2-methyl-1-(methylsulfonyl)indolin-5-yl) ethanone scaffold. The synthesis, SAR analysis, enzyme active site docking, enzymatic reaction kinetics, selectivity and pharmaceutical properties are discussed.

Abstract 3219: Changing the metabolism of cancer cells with PKM2 activators - a path to a cancer metabolism drug

Pyruvate kinase (PK) M2 which catalyzes the last step in glycolysis, is the alternative spliced isoform expressed in cancer cells, and a key player in exerting the Warburg effect. One of the mechanisms by which PKM2 modulate cancer cell metabolism is by switching between the low activity monomer and the high activity tetramer forms. This process is controlled by the varying concentration of an upstream glycolytic intermediate, FBP. These changes enable the cancer cell to manage its usage of glucose carbon backbones, whether for ATP production or for biomass generation, according to its changing demands. Further, it has been recently shown that the elevated levels of ROS in cancer cells contribute to the decreased activity of PKM2 to support NADPH production to increase cellular anti-oxidation capacity to sustain proliferation. Our goal is to disrupt the metabolic adaptation of cancer cells with small molecule PKM2 modulators. We hypothesized that an activator will redirect the consumption of nutrients, especially of glucose, away from biomass production and ultimately send the cells to die. Using our proprietary structure-based technology, we identified several series of novel allosteric PKM2 activators. Chemical optimization resulted in potent compounds with AC50 as low as 10nM, which were selective against the other PK isoforms. These compounds were proven to stabilize the active tetramer form of PKM2 in cancer cells. Bioenergetic experiments in several cell lines demonstrated that not only do these agents reduce lactate production; they also reduce the oxygen consumption rate. Analysis of cell cycle showed that treatment with PKM2 activators causes the cells to arrest at the G1 phase. In outcome-based assays, these compounds significantly reduced the proliferation rate of various cancer cell lines, and this effect was sensitive to media conditions, such as glucose levels. Taken together, our data supports the hypothesis that activation of PKM2 effectively deprives the cancer cell of building blocks and reduces the detoxification capacity that are required to support growth and proliferation. An in vivo colorectal cancer HT29 cell line mouse xenograft model with a modestly active compound (IC50=0.9uM) demonstrated tumor growth inhibition greater than 50% (100 mg/kg Q2D and 200 and 400 mg/kg IP QD). The compound was very safe in mouse, even at the highest exposure levels (200 and 400 mg/kg IP QD), indicating that these efficacious doses are significantly lower than the MTD. Additional xenograft models are ongoing. Taken together, there is strong support for the effect of potent PKM2 small molecule activators on the cellular metabolism of cancer cells, demonstrating statistically significant anti-cancer effect in an animal model of colorectal cancer. The favorable DMPK profile of these compounds further supports their development as anti-proliferative agents, both as a single agent and in combination therapy.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 3219. doi:1538-7445.AM2012-3219

Small molecule correctors of F508del-CFTR discovered by structure-based virtual screening

Folding correctors of F508del-CFTR were discovered by in silico structure-based screening utilizing homology models of CFTR. The intracellular segment of CFTR was modeled and three cavities were identified at inter-domain interfaces: (1) Interface between the two Nucleotide Binding Domains (NBDs); (2) Interface between NBD1 and Intracellular Loop (ICL) 4, in the region of the F508 deletion; (3) multi-domain interface between NBD1:2:ICL1:2:4. We hypothesized that compounds binding at these interfaces may improve the stability of the protein, potentially affecting the folding yield or surface stability. In silico structure-based screening was performed at the putative binding-sites and a total of 496 candidate compounds from all three sites were tested in functional assays. A total of 15 compounds, representing diverse chemotypes, were identified as F508del folding correctors. This corresponds to a 3% hit rate, ~tenfold higher than hit rates obtained in corresponding high-throughput screening campaigns. The same binding sites also yielded potentiators and, most notably, compounds with a dual corrector-potentiator activity (dual-acting). Compounds harboring both activity types may prove to be better leads for the development of CF therapeutics than either pure correctors or pure potentiators. To the best of our knowledge this is the first report of structure-based discovery of CFTR modulators.

Docking to single-domain and multiple-domain proteins: Old and new challenges

The diverse selection of targets in the CAPRI experiments provides grounds for determining the limits of our rigid-body docking program MolFit, and for extending it. We find that the sensitivity of MolFit is high, enabling it to produce reasonably accurate docking solutions when the structures undergo moderate local conformation changes upon complex formation or when the docked molecules are modeled. Yet the ranks of these solutions are sometimes too low to meet the requirements of CAPRI assessment. This indicates that the selectivity of MolFit, which was optimized for docking of unbound X-ray structures, and which relies on the availability of external data from biochemical and bioinformatic sources, needs readjustment in order to meet the challenges presented by NMR or modeled structures. A different challenge is presented by large global conformation changes such as movements of domains. We show that such changes can be accommodated within the rigid-body approximation by employing rigid multibody multistage docking procedures. We also address the difficulty of ranking results from 2-body and multibody docking scans in cases in which there are no external data favoring one option over the other.

Experimental and Computational Characterization of the Dimerization of the PTS-regulation Domains of BglG from Escherichia coli

BglG and LicT are transcriptional antiterminators from Escherichia coli and Bacillus subtilis, respectively, that control the expression of genes and operons involved in transport and catabolism of carbohydrates. Both proteins contain a duplicate conserved domain, the PTS-regulation domain (PRD), and they are regulated by phosphorylation on specific, highly conserved histidine residues located in the PRDs. However, despite their similar function and the high sequence identity, experimental evidence implies different modes of regulation. Thus, BglG must be de-phosphorylated on PRD2 in order to form active dimers, whereas activation of LicT requires de-phosphorylation on PRD1 and phosphorylation on PRD2. Here we address two goals. First, we test in vivo and in silico the effect of point mutations in the PRDs of BglG on the PRD-PRD dimerization. Second, we explore computationally the effect of histidine phosphorylation on PRD dimerization in BglG and LicT. We find excellent correspondence between the experimental and computational measures of the influence of mutations on PRD dimerization in BglG. This establishes that the geometric-electrostatic complementarity scores computed with the program MolFit provide a good measure of the effects of mutations in this system. In addition, it indicates that the dimerization mode of the separately expressed PRDs of BglG is similar to the dimers formed by activated LicT. The computations also show that phosphorylation of the histidine residues in PRD1 of either BglG or LicT leads to a strong electrostatic repulsion. Conversely, the phosphorylation of one histidine residue in PRD2 of LicT leads to improved electrostatic complementarity at the PRD2-PRD2 interface, whereas the corresponding phosphorylation in BglG has negligible contribution. This different conduct may be attributed to a single replacement in the sequence of PRD2 in BglG compared to LicT, Ala262 versus Asp261, respectively.

Weighted geometric docking: incorporating external information in the rotation-translation scan

Weighted geometric docking is a prediction algorithm that matches weighted molecular surfaces. Each molecule is represented by a grid of complex numbers, storing information about the shape of the molecule in the real part and weight information in the imaginary part. The weights are based on experimental biochemical and biophysical data or on theoretical analyses of amino acid conservation or correlation patterns in multiple-sequence alignments of homologous proteins. Only a few surface residues on either one or both molecules are weighted. In contrast to methods that use postscan filtering based on biochemical information, our method incorporates the external data in the rotation-translation search, producing a different set of docking solutions biased toward solutions in which the up-weighted residues are at the interface. Similarly, interactions involving specified residues can be impeded. The weighted geometric algorithm was applied to five systems for which regular geometric docking of the unbound molecules gave poor results. We obtained much better ranking of the nearly correct prediction and higher statistical significance when weighted geometric docking was used. The method was successful even when the weighted portion of the surface corresponded only partially and approximately to the binding site.

Prediction of the unknown: inspiring experience with the CAPRI experiment

We submitted predictions for all seven targets in the CAPRI experiment. For four targets, our submitted models included acceptable, medium accuracy predictions of the structures of the complexes, and for a fifth target we identified the location of the binding site of one of the molecules. We used a weighted-geometric docking algorithm in which contacts involving specified parts of the surfaces of either one or both molecules were up-weighted or down-weighted. The weights were based on available structural and biochemical data or on sequence analyses. The weighted-geometric docking proved very useful for five targets, improving the complementarity scores and the ranks of the nearly correct solutions, as well as their statistical significance. In addition, the weighted-geometric docking promoted formation of clusters of similar solutions, which include more accurate predictions.

Prediction of the structure of the complex between the 30S ribosomal subunit and colicin E3 via weighted-geometric docking

Colicin E3 kills Escherichia coli cells by ribonucleolytic cleavage in the 16S rRNA. The cleavage occurs at the ribosomal decoding A-site between nucleotides A1493 and G1494. The breaking of this single phosphodiester bond results in a complete termination of protein biosynthesis leading to cell death. A model structure of the complex of the ribosomal subunit 30S and colicin E3 was constructed by means of a new weighted-geometric docking algorithm, in which interactions involving specified parts of the molecular surface can be up-weighted, allowing incorporation of experimental data in the docking search. Our model, together with available experimental data, predicts the role of the catalytic residues of colicin E3. In addition, it suggests that bound acidic immunity protein inhibits the enzymatic activity of colicin E3 by electrostatic repulsion of the negatively charged substrate.

The RNA helicase DbpA exhibits a markedly different conformation in the ADP-bound state when compared with the ATP- or RNA-bound states

The motor enzymes that belong to the family of RNA helicases catalyze the strand separation of duplex RNA via ATP hydrolysis. Among these enzymes, Escherichia coli DbpA is a unique RNA helicase because it possesses ATPase-specific activity toward the peptidyl transferase center in 23 S ribosomal RNA. For this reason, it has been the subject of numerous biochemical and structure-function studies. The ATP-stimulated unwinding activity of DbpA toward specific and nonspecific RNA duplexes has been demonstrated. However, the underlying molecular and structural basis, which facilitates its helicase activities, is presently not known. We combined time-dependent limited proteolysis digestion, fluorescence spectroscopy, and three-dimensional structural homology modeling techniques to study the structural conformations of DbpA with respect to its binding to stoichiometric ratios of RNA and cofactors. We show that the conformational state of DbpA is markedly different in the ADP-bound state than in any other state (ATP- or RNA-bound). These results, together with structural homology studies, suggest that a hinge region located in the core domain of DbpA mediates such conformational changes.

On the interaction of colicin E3 with the ribosome

Colicin E3 is a protein that kills Escherichia coli cells by a process that involves binding to a surface receptor, entering the cell and inactivating its protein biosynthetic machinery. Colicin E3 kills cells by a catalytic mechanism of a specific ribonucleolytic cleavage in 16S rRNA at the ribosomal decoding A-site between A1493 and G1494 (E. coli numbering system). The breaking of this single phosphodiester bond results in a complete cessation of protein biosynthesis and cell death. The inactive E517Q mutant of the catalytic domain of colicin E3 binds to 30S ribosomal subunits of Thermus thermophilus, as demonstrated by an immunoblotting assay. A model structure of the complex of the ribosomal subunit 30S and colicin E3, obtained via docking, explains the role of the catalytic residues, suggests a catalytic mechanism and provides insight into the specificity of the reaction. Furthermore, the model structure suggests that the inhibitory action of bound immunity is due to charge repulsion of this acidic protein by the negatively charged rRNA backbone

Phage display of a cellulose binding domain from Clostridium thermocellum and its application as a tool for antibody engineering

Phage display of antibody fragments has proved to be a powerful tool for the isolation and in vitro evolution of these biologically important molecules. However, the general usefulness of this technology is still limited by some technical difficulties. One of the most debilitating obstacles to the widespread application of the technology is the accumulation of "insert loss" clones in the libraries; phagemid clones from which the DNA encoding part or all of the cloned antibody fragment had been deleted. Another difficulty arises when phage technology is applied for cloning hybridoma-derived antibody genes, where myeloma derived light chains, irrelevant to the hybridoma's antibody specificity may be fortuitously cloned. Here, we report the construction of a novel phage-display system designed to address these problems. In our system a single-chain Fv (scFv) is expressed as an in-frame fusion protein with a cellulose-binding domain (CBD) derived from the Clostridium thermocellum cellulosome. The CBD domain serves as an affinity tag allowing rapid phage capture and concentration from crude culture supernatants, and immunological detection of both displaying phage and soluble scFv produced thereof. We demonstrate the utility of our system in solving the technical difficulties described above, and in speeding up the process of scFv isolation from combinatorial antibody repertoires.