Universal stabilization of the influenza hemagglutinin by structure-based redesign of the pH switch regions.

For an efficacious vaccine immunogen, influenza hemagglutinin (HA) needs to maintain a stable quaternary structure, which is contrary to the inherently dynamic and metastable nature of class I fusion proteins. In this study, we stabilized HA with three substitutions within its pH-sensitive regions where the refolding starts. An X-ray structure reveals how these substitutions stabilize the intersubunit ß-sheet in the base and form an interprotomeric aliphatic layer across the stem while the native prefusion HA fold is retained. The identification of the stabilizing substitutions increases our understanding of how the pH sensitivity is structurally accomplished in HA and possibly other pH-sensitive class I fusion proteins. Our stabilization approach in combination with the occasional back mutation of rare amino acids to consensus results in well-expressing stable trimeric HAs. This repair and stabilization approach, which proves broadly applicable to all tested influenza A HAs of group 1 and 2, will improve the developability of influenza vaccines based on different types of platforms and formats and can potentially improve efficacy.

Structural basis for recognition of the central conserved region of RSV G by neutralizing human antibodies.

Respiratory syncytial virus (RSV) is a major cause of severe lower respiratory tract infections in infants and the elderly, and yet there remains no effective treatment or vaccine. The surface of the virion is decorated with the fusion glycoprotein (RSV F) and the attachment glycoprotein (RSV G), which binds to CX3CR1 on human airway epithelial cells to mediate viral attachment and subsequent infection. RSV G is a major target of the humoral immune response, and antibodies that target the central conserved region of G have been shown to neutralize both subtypes of RSV and to protect against severe RSV disease in animal models. However, the molecular underpinnings for antibody recognition of this region have remained unknown. Therefore, we isolated two human antibodies directed against the central conserved region of RSV G and demonstrated that they neutralize RSV infection of human bronchial epithelial cell cultures in the absence of complement. Moreover, the antibodies protected cotton rats from severe RSV disease. Both antibodies bound with high affinity to a secreted form of RSV G as well as to a peptide corresponding to the unglycosylated central conserved region. High-resolution crystal structures of each antibody in complex with the G peptide revealed two distinct conformational epitopes that require proper folding of the cystine noose located in the C-terminal part of the central conserved region. Comparison of these structures with the structure of fractalkine (CX3CL1) alone or in complex with a viral homolog of CX3CR1 (US28) suggests that RSV G would bind to CX3CR1 in a mode that is distinct from that of fractalkine. Collectively, these results build on recent studies demonstrating the importance of RSV G in antibody-mediated protection from severe RSV disease, and the structural information presented here should guide the development of new vaccines and antibody-based therapies for RSV.

3D-e-Chem-VM: Structural Cheminformatics Research Infrastructure in a Freely Available Virtual Machine.

3D-e-Chem-VM is an open source, freely available Virtual Machine ( http://3d-e-chem.github.io/3D-e-Chem-VM/ ) that integrates cheminformatics and bioinformatics tools for the analysis of protein-ligand interaction data. 3D-e-Chem-VM consists of software libraries, and database and workflow tools that can analyze and combine small molecule and protein structural information in a graphical programming environment. New chemical and biological data analytics tools and workflows have been developed for the efficient exploitation of structural and pharmacological protein-ligand interaction data from proteomewide databases (e.g., ChEMBLdb and PDB), as well as customized information systems focused on, e.g., G protein-coupled receptors (GPCRdb) and protein kinases (KLIFS). The integrated structural cheminformatics research infrastructure compiled in the 3D-e-Chem-VM enables the design of new approaches in virtual ligand screening (Chemdb4VS), ligand-based metabolism prediction (SyGMa), and structure-based protein binding site comparison and bioisosteric replacement for ligand design (KRIPOdb).

Steering Siglec-Sialic Acid Interactions on Living Cells using Bioorthogonal Chemistry.

Sialic acid sugars that terminate cell-surface glycans form the ligands for the sialic acid binding immunoglobulin-like lectin (Siglec) family, which are immunomodulatory receptors expressed by immune cells. Interactions between sialic acid and Siglecs regulate the immune system, and aberrations contribute to pathologies like autoimmunity and cancer. Sialic acid/Siglec interactions between living cells are difficult to study owing to a lack of specific tools. Here, we report a glycoengineering approach to remodel the sialic acids of living cells and their binding to Siglecs. Using bioorthogonal chemistry, a library of cells with more than sixty different sialic acid modifications was generated that showed dramatically increased binding toward the different Siglec family members. Rational design reduced cross-reactivity and led to the discovery of three selective Siglec-5/14 ligands. Furthermore, glycoengineered cells carrying sialic acid ligands for Siglec-3 dampened the activation of Siglec-3+ monocytic cells through the NF-κB and IRF pathways.

Involvement of a Hydrophobic Pocket and Helix†11 in Determining the Modes of Action of Prenylated Flavonoids and Isoflavonoids in the Human Estrogen Receptor.

Six prenylated (iso)flavonoids were purified from a licorice root extract and subjected to competition experiments with six commercially available (iso)flavonoids. The agonistic and antagonistic activities of these compounds towards both hERα (human estrogen receptor alpha) and hERß were determined. Differences in the modes of action (agonist or antagonist) were observed for the various compounds tested. In general, each compound had the same mode of action towards both ERs. In silico modeling was performed in order to study the differences in estrogenicity observed between the compounds. It is suggested that prenyl chains fit into a hydrophobic pocket present in the hER, resulting in an increased agonistic activity. In addition, it was shown that an increase in length (≈1.7 Å) of pyran prenylated isoflavonoids resulted in an antagonistic mode of action. This might be caused by collision of the pyran ring with helix 11 in the ligand binding cavity of the hER.

Statin Lactonization by Uridine 5'-Diphospho-glucuronosyltransferases (UGTs).

Statins are cholesterol-lowering drugs that have proven to be effective in lowering the risk of major cardiovascular events. Although well tolerated, statin-induced myopathies are the most common side effects. Compared to their pharmacologically active acid form, statin lactones are more potent inducers of toxicity. They can be formed by glucuronidation mediated by uridine 5'-diphospho-glucuronosyltransferases (UGTs), but a systematic characterization of subtype specificity and kinetics of lactonization is lacking. Here, we demonstrate for six clinically relevant statins that only UGT1A1, 1A3, and 2B7 contribute significantly to their lactonization. UGT1A3 appeared to have the highest lactonization capacity with marked differences in statin conversion rates: pitavastatin ≫ atorvastatin > cerivastatin > lovastatin > rosuvastatin (simvastatin not converted). Using in silico modeling we could identify a probable statin interaction region in the UGT binding pocket. Polymorphisms in these regions of UGT1A1, 1A3, and 2B7 may be a contributing factor in statin-induced myopathies, which could be used in personalization of statin therapy with improved safety.

Discovery of small molecule CD40-TRAF6 inhibitors.

The CD154-CD40 receptor complex plays a pivotal role in several inflammatory pathways. Attempts to inhibit the formation of this complex have resulted in systemic side effects. Downstream inhibition of the CD40 signaling pathway therefore seems a better way to ameliorate inflammatory disease. To relay a signal, the CD40 receptor recruits adapter proteins called tumor necrosis factor receptor-associated factors (TRAFs). CD40-TRAF6 interactions are known to play an essential role in several inflammatory diseases. We used in silico, in vitro, and in vivo experiments to identify and characterize compounds that block CD40-TRAF6 interactions. We present in detail our drug docking and optimization pipeline and show how we used it to find lead compounds that reduce inflammation in models of peritonitis and sepsis. These compounds appear to be good leads for drug development, given the observed absence of side effects and their demonstrated efficacy for peritonitis and sepsis in mouse models.

In silico identification and in vitro validation of potential cholestatic compounds through 3D ligand-based pharmacophore modeling of BSEP inhibitors.

Drug-induced cholestasis is a frequently observed side effect of drugs and is often caused by an unexpected interaction with the bile salt export pump (BSEP/ABCB11). BSEP is the key membrane transporter responsible for the transport of bile acids from hepatocytes into bile. Here, we developed a pharmacophore model that describes the molecular features of compounds associated with BSEP inhibitory activity. To generate input and validation data sets, in vitro experiments with membrane vesicles overexpressing human BSEP were used to assess the effect of compounds (50 µM) on BSEP-mediated (3)H-taurocholic acid transport. The model contains two hydrogen bond acceptor/anionic features, two hydrogen bond acceptor vector features, four hydrophobic/aromatic features, and exclusion volumes. The pharmacophore was validated against a set of 59 compounds, including registered drugs. The model recognized 9 out of 12 inhibitors (75%), which could not be identified based on general parameters, such as molecular weight or SlogP, alone. Finally, the model was used to screen a virtual compound database. A number of compounds found via virtual screening were tested and displayed statistically significant BSEP inhibition, ranging from 13 ± 1% to 67 ± 7% of control (P < 0.05). In conclusion, we developed and validated a pharmacophore model that describes molecular features found in BSEP inhibitors. The model may be used as an in silico screening tool to identify potentially harmful drug candidates at an early stage in drug development.

Influence of azide incorporation on binding affinity by small papain inhibitors.

In order to develop affinity-based biosensor platforms, appropriate ligands with a functional handle for immobilization onto a biosensor surface are required. To this end, a library of papain inhibitors was designed and synthesized, containing different azide linkers for subsequent immobilization by 'click' chemistry, in this particular case by copper-free, strain-promoted azide-alkyne cycloaddition (SPAAC). Furthermore, a molecular docking study was performed to obtain a better insight as to at which position such azide handles could be tolerated without affecting binding affinity. Although the azide moiety is small, in some cases its introduction strongly influenced the binding affinity. For one class of inhibitors a swapped binding mode was proposed to explain the results. In addition, a specific site for linker introduction was identified, which did not significantly affect the binding affinity.

Efficacious human metapneumovirus vaccine based on AI-guided engineering of a closed prefusion trimer.

The prefusion conformation of human metapneumovirus fusion protein (hMPV Pre-F) is critical for eliciting the most potent neutralizing antibodies and is the preferred immunogen for an efficacious vaccine against hMPV respiratory infections. Here we show that an additional cleavage event in the F protein allows closure and correct folding of the trimer. We therefore engineered the F protein to undergo double cleavage, which enabled screening for Pre-F stabilizing substitutions at the natively folded protomer interfaces. To identify these substitutions, we developed an AI convolutional classifier that successfully predicts complex polar interactions often overlooked by physics-based methods and visual inspection. The combination of additional processing, stabilization of interface regions and stabilization of the membrane-proximal stem, resulted in a Pre-F protein vaccine candidate without the need for a heterologous trimerization domain that exhibited high expression yields and thermostability. Cryo-EM analysis shows the complete ectodomain structure, including the stem, and a specific interaction of the newly identified cleaved C-terminus with the adjacent protomer. Importantly, the protein induces high and cross-neutralizing antibody responses resulting in near complete protection against hMPV challenge in cotton rats, making the highly stable, double-cleaved hMPV Pre-F trimer an attractive vaccine candidate.

High-affinity inhibitors of Zymomonas mobilis tRNA-guanine transglycosylase through convergent optimization.

The tRNA-modifying enzyme tRNA-guanine transglycosylase (TGT) has been recognized as a drug target for the treatment of the foodborne illness shigellosis. The active site of TGT consists of three pockets: the central guanine/preQ1 recognition site and the ribose-33 and ribose-34 pockets. In previous work, lin-benzoguanines and lin-benzohypoxanthines, which differ by the presence of an exocyclic NH2 group in the former and its absence in the latter, were used as central scaffolds that bind to the guanine/preQ1 recognition site and allow suitable functionalization along exit vectors targeting the two ribose pockets. The substituents for both of these two pockets have been optimized individually. Here, a series of bifunctionalized inhibitors that occupy both ribose pockets are reported for the first time. Dissociation constants Kd down to the picomolar range were measured for the bifunctionalized lin-benzoguanine-based ligands and Kd values in the nanomolar range were measured for the corresponding lin-benzohypoxanthine-based ligands. The binding mode of all inhibitors was elucidated by X-ray crystal structure analysis. A remarkable influence of the crystallization protocol on the solvation pattern in the solid state and the residual mobility of the bound ligands was observed.

RSV A2-Based Prefusion F Vaccine Candidates Induce RSV A and RSV B Cross Binding and Neutralizing Antibodies and Provide Protection against RSV A and RSV B Challenge in Preclinical Models.

RSV is divided into two antigenic subtypes, RSV A and RSV B, which is largely based on the variation in the G protein, while the fusion protein F is more conserved and a target for antibody-mediated neutralization. Here we evaluate the breadth of the protective immune responses across RSV A and RSV B subtypes, induced by vaccines based on the RSV A-based fusion protein, stabilized in the prefusion conformation (preF) in preclinical models. Immunization of naïve cotton rats with preF subunit or preF encoded by a replication incompetent Adenoviral 26, induced antibodies capable of neutralizing recent RSV A and RSV B clinical isolates, as well as protective efficacy against a challenge with RSV A and RSV B strains. Similarly, induction of cross-neutralizing antibodies was observed after immunization with Ad26-encoded preF, preF protein or a mix of both (Ad26/preF protein) in RSV pre-exposed mice and African Green Monkeys. Transfer of serum of human subjects immunized with Ad26/preF protein into cotton rats provide protection against challenges with both RSV A and RSV B, with complete protection against both strains observed in the lower respiratory tract. In contrast, almost no protection against RSV A and B infection was observed after the transfer of a human serum pool isolated pre-vaccination. These results collectively show that the RSV A-based monovalent Ad26/preF protein vaccine induced neutralizing antibodies, as well as protection against both RSV A and RSV B subtypes in animals, including by passive transfer of human antibodies alone, suggesting that clinical efficacy against both subtypes can be achieved.

Enhancing breadth and durability of humoral immune responses in non-human primates with an adjuvanted group 1 influenza hemagglutinin stem antigen.

Seasonal influenza vaccines must be updated annually and suboptimally protect against strains mismatched to the selected vaccine strains. We previously developed a subunit vaccine antigen consisting of a stabilized trimeric influenza A group 1 hemagglutinin (H1) stem protein that elicits broadly neutralizing antibodies. Here, we further optimized the stability and manufacturability of the H1 stem antigen (H1 stem v2, also known as INFLUENZA G1 mHA) and characterized its formulation and potency with different adjuvants in vitro and in animal models. The recombinant H1 stem antigen (50 µg) was administered to influenza-naïve non-human primates either with aluminum hydroxide [Al(OH)3] + NaCl, AS01B, or SLA-LSQ formulations at week 0, 8 and 34. These SLA-LSQ formulations comprised of varying ratios of the synthetic TLR4 agonist 'second generation synthetic lipid adjuvant' (SLA) with liposomal QS-21 (LSQ). A vaccine formulation with aluminum hydroxide or SLA-LSQ (starting at a 10:25 µg ratio) induced HA-specific antibodies and breadth of neutralization against a panel of influenza A group 1 pseudoviruses, comparable with vaccine formulated with AS01B, four weeks after the second immunization. A formulation with SLA-LSQ in a 5:2 µg ratio contained larger fused or aggregated liposomes and induced significantly lower humoral responses. Broadly HA stem-binding antibodies were detectable for the entire period after the second vaccine dose up to week 34, after which they were boosted by a third vaccine dose. These findings inform about potential adjuvant formulations in clinical trials with an H1 stem-based vaccine candidate.

From lin-benzoguanines to lin-benzohypoxanthines as ligands for Zymomonas mobilis tRNA-guanine transglycosylase: replacement of protein-ligand hydrogen bonding by importing water clusters.

The foodborne illness shigellosis is caused by Shigella bacteria that secrete the highly cytotoxic Shiga toxin, which is also formed by the closely related enterohemorrhagic Escherichia coli (EHEC). It has been shown that tRNA-guanine transglycosylase (TGT) is essential for the pathogenicity of Shigella flexneri. Herein, the molecular recognition properties of a guanine binding pocket in Zymomonas mobilis TGT are investigated with a series of lin-benzohypoxanthine- and lin-benzoguanine-based inhibitors that bear substituents to occupy either the ribose-33 or the ribose-34 pocket. The three inhibitor scaffolds differ by the substituent at C(6) being H, NH(2), or NH-alkyl. These differences lead to major changes in the inhibition constants, pK(a) values, and binding modes. Compared to the lin-benzoguanines, with an exocyclic NH(2) at C(6), the lin-benzohypoxanthines without an exocyclic NH(2) group have a weaker affinity as several ionic protein-ligand hydrogen bonds are lost. X-ray cocrystal structure analysis reveals that a new water cluster is imported into the space vacated by the lacking NH(2) group and by a conformational shift of the side chain of catalytic Asp102. In the presence of an N-alkyl group at C(6) in lin-benzoguanine ligands, this water cluster is largely maintained but replacement of one of the water molecules in the cluster leads to a substantial loss in binding affinity. This study provides new insight into the role of water clusters at enzyme active sites and their challenging substitution by ligand parts, a topic of general interest in contemporary structure-based drug design.

Pharmacophore fingerprint-based approach to binding site subpocket similarity and its application to bioisostere replacement.

Bioisosteres have been defined as structurally different molecules or substructures that can form comparable intermolecular interactions, and therefore, fragments that bind to similar protein structures exhibit a degree of bioisosterism. We present KRIPO (Key Representation of Interaction in POckets): a new method for quantifying the similarities of binding site subpockets based on pharmacophore fingerprints. The binding site fingerprints have been optimized to improve their performance for both intra- and interprotein family comparisons. A range of attributes of the fingerprints was considered in the optimization, including the placement of pharmacophore features, whether or not the fingerprints are fuzzified, and the resolution and complexity of the pharmacophore fingerprints (2-, 3-, and 4-point fingerprints). Fuzzy 3-point pharmacophore fingerprints were found to represent the optimal balance between computational resource requirements and the identification of potential replacements. The complete PDB was converted into a database comprising almost 300,000 optimized fingerprints of local binding sites together with their associated ligand fragments. The value of the approach is demonstrated by application to two crystal structures from the Protein Data Bank: (1) a MAP kinase P38 structure in complex with a pyridinylimidazole inhibitor (1A9U) and (2) a complex of thrombin with melagatran (1K22). Potentially valuable bioisosteric replacements for all subpockets of the two studied protein are identified.

Stabilizing the closed SARS-CoV-2 spike trimer.

The trimeric spike (S) protein of SARS-CoV-2 is the primary focus of most vaccine design and development efforts. Due to intrinsic instability typical of class I fusion proteins, S tends to prematurely refold to the post-fusion conformation, compromising immunogenic properties and prefusion trimer yields. To support ongoing vaccine development efforts, we report the structure-based design of soluble S trimers with increased yields and stabilities, based on introduction of single point mutations and disulfide-bridges. We identify regions critical for stability: the heptad repeat region 1, the SD1 domain and position 614 in SD2. We combine a minimal selection of mostly interprotomeric mutations to create a stable S-closed variant with a 6.4-fold higher expression than the parental construct while no longer containing a heterologous trimerization domain. The cryo-EM structure reveals a correctly folded, predominantly closed pre-fusion conformation. Highly stable and well producing S protein and the increased understanding of S protein structure will support vaccine development and serological diagnostics.

Crystal structure analysis and in silico pKa calculations suggest strong pKa shifts of ligands as driving force for high-affinity binding to TGT.

A novel ligand series is presented to inhibit tRNA-guanine transglycosylase (TGT), a protein with a significant role in the pathogenicity mechanism of Shigella flexneri, the causative agent of Shigellosis. The enzyme exchanges guanine in the wobble position of tRNA(Asn,Asp,His,Tyr) against a modified base. To prevent the base-exchange reaction, several series of inhibitors have already been designed, synthesized, and tested. One aim of previous studies was to address a hydrophobic pocket with different side chains attached to the parent skeletons. Disappointingly, no significant increase in binding affinity could be observed that could be explained by the disruption of a conserved water cluster. The ligand series examined in this study are based on the known scaffold lin-benzoguanine. Different side chains were introduced leading to 2-amino-lin-benzoguanines, which address a different pocket of the protein and avoid disruption of the water cluster. With the introduction of an amino group in the 2-position, a dramatic increase in binding affinity can be experienced. To explain this significant gain in binding affinity, Poisson-Boltzmann calculations were performed to explore pK(a) changes of ligand functional groups upon protein binding, they can differ significantly on going from aqueous solution to protein environment. For all complexes, a permanent protonation of the newly designed ligands is suggested, leading to a charge-assisted hydrogen bond in the protein-ligand complex. This increased strength in hydrogen bonding takes beneficial effect on binding affinity of the ligands, resulting in low-nanomolar binders. Crystal structures and docking emphasize the importance of the newly created charge-assisted hydrogen bond. A detailed analysis of the crystal structures in complex with substituted 2-amino-lin-benzoguanines indicate pronounced disorder of the attached side chains addressing the ribose 33 binding pocket. Docking suggests multiple orientations of these side chains. Obviously, an entropic advantage of the residual mobility experienced by these ligands in the bound state is beneficial and reveals an overall improved protein binding.

High-affinity inhibitors of tRNA-guanine transglycosylase replacing the function of a structural water cluster.

The tRNA-modifying enzyme tRNA-guanine transglycosylase (TGT) is essential for the pathogenic mechanism of Shigella flexneri, the causing agent of the bacterial diarrheal disease shigellosis. Herein, the synthesis of a new class of rationally designed 6-amino-imidazo[4,5-g]quinazolin-8(7H)-one- (lin-benzoguanine) based inhibitors of TGT are reported. In order to accommodate a small hydrophobic crevice opening near the binding site of ribose-34, 2-aminoethyl substituents were introduced in position 4 of the heterocyclic scaffold. For this purpose, a synthetic sequence consisting of iodination, Suzuki cross-coupling, hydroboration, Mitsunobu reaction, and Gabriel synthesis was employed, furnishing a primary amine that served as a common intermediate for the preparation of a series of derivatives. The resulting ligands displayed very low inhibition constants, down to K(i)=2 nM. Substantial additional inhibitory potency is gained by interaction of terminal lipophilic groups attached to the substituent at position 4 with the hydrophobic crevice shaped by Val45 and Leu68. At the same time, the secondary ammonium center in the substituent displaces a cluster of water molecules, solvating the catalytic residues Asp102 and Asp280, without loss in binding affinity. In addition, a synthetic intermediate with an unusual 3,6,7,8,9,10-hexahydroimidazo[4,5-g][1,3]benzodiazepine core, as confirmed by X-ray analysis, is reported.

Drug Discovery and Repurposing Inhibits a Major Gut Pathogen-Derived Oncogenic Toxin.

Objective: The human intestinal microbiome plays an important role in inflammatory bowel disease (IBD) and colorectal cancer (CRC) development. One of the first discovered bacterial mediators involves Bacteroides fragilis toxin (BFT, also named as fragilysin), a metalloprotease encoded by enterotoxigenic Bacteroides fragilis (ETBF) that causes barrier disruption and inflammation of the colon, leads to tumorigenesis in susceptible mice, and is enriched in the mucosa of IBD and CRC patients. Thus, targeted inhibition of BFT may benefit ETBF carrying patients. Design: By applying two complementary in silico drug design techniques, drug repositioning and molecular docking, we predicted potential BFT inhibitory compounds. Top candidates were tested in vitro on the CRC epithelial cell line HT29/c1 for their potential to inhibit key aspects of BFT activity, being epithelial morphology changes, E-cadherin cleavage (a marker for barrier function) and increased IL-8 secretion. Results: The primary bile acid and existing drug chenodeoxycholic acid (CDCA), currently used for treating gallstones, cerebrotendinous xanthomatosis, and constipation, was found to significantly inhibit all evaluated cell responses to BFT exposure. The inhibition of BFT resulted from a direct interaction between CDCA and BFT, as confirmed by an increase in the melting temperature of the BFT protein in the presence of CDCA. Conclusion: Together, our results show the potential of in silico drug discovery to combat harmful human and microbiome-derived proteins and more specifically suggests a potential for retargeting CDCA to inhibit the pro-oncogenic toxin BFT.

Molecular interaction fingerprint approaches for GPCR drug discovery.

Protein-ligand interaction fingerprints (IFPs) are binary 1D representations of the 3D structure of protein-ligand complexes encoding the presence or absence of specific interactions between the binding pocket amino acids and the ligand. Various implementations of IFPs have been developed and successfully applied for post-processing molecular docking results for G Protein-Coupled Receptor (GPCR) ligand binding mode prediction and virtual ligand screening. Novel interaction fingerprint methods enable structural chemogenomics and polypharmacology predictions by complementing the increasing amount of GPCR structural data. Machine learning methods are increasingly used to derive relationships between bioactivity data and fingerprint descriptors of chemical and structural information of binding sites, ligands, and protein-ligand interactions. Factors that influence the application of IFPs include structure preparation, binding site definition, fingerprint similarity assessment, and data processing and these factors pose challenges as well possibilities to optimize interaction fingerprint methods for GPCR drug discovery.