REvoLd: Ultra-Large Library Screening with an Evolutionary Algorithm in Rosetta.

Ultra-large make-on-demand compound libraries now contain billions of readily available compounds. This represents a golden opportunity for in-silico drug discovery. One challenge, however, is the time and computational cost of an exhaustive screen of such large libraries when receptor flexibility is taken into account. We propose an evolutionary algorithm to search combinatorial make-on-demand chemical space efficiently without enumerating all molecules. We exploit the feature of make-on-demand compound libraries, namely that they are constructed from lists of substrates and chemical reactions. Our novel algorithm RosettaEvolutionaryLigand (REvoLd) explores the vast search space of combinatorial libraries for protein-ligand docking with full ligand and receptor flexibility through RosettaLigand. A benchmark of REvoLd on five drug targets showed improvements in hit rates by factors between 869 and 1,622 compared to random selections. REvoLd is available as an application within the Rosetta software suite.

Discovery of Protease-Activated Receptor 4 (PAR4)-Tethered Ligand Antagonists Using Ultralarge Virtual Screening.

Here, we demonstrate a structure-based small molecule virtual screening and lead optimization pipeline using a homology model of a difficult-to-drug G-protein-coupled receptor (GPCR) target. Protease-activated receptor 4 (PAR4) is activated by thrombin cleavage, revealing a tethered ligand that activates the receptor, making PAR4 a challenging target. A virtual screen of a make-on-demand chemical library yielded a one-hit compound. From the single-hit compound, we developed a novel series of PAR4 antagonists. Subsequent lead optimization via simultaneous virtual library searches and structure-based rational design efforts led to potent antagonists of thrombin-induced activation. Interestingly, this series of antagonists was active against PAR4 activation by the native protease thrombin cleavage but not the synthetic PAR4 agonist peptide AYPGKF.

Combining Rosetta Sequence Design with Protein Language Model Predictions Using Evolutionary Scale Modeling (ESM) as Restraint.

Computational protein sequence design has the ambitious goal of modifying existing or creating new proteins; however, designing stable and functional proteins is challenging without predictability of protein dynamics and allostery. Informing protein design methods with evolutionary information limits the mutational space to more native-like sequences and results in increased stability while maintaining functions. Recently, language models, trained on millions of protein sequences, have shown impressive performance in predicting the effects of mutations. Assessing Rosetta-designed sequences with a language model showed scores that were worse than those of their original sequence. To inform Rosetta design protocols with language model predictions, we added a new metric to restrain the energy function during design using the Evolutionary Scale Modeling (ESM) model. The resulting sequences have better language model scores and similar sequence recovery, with only a minor decrease in the fitness as assessed by Rosetta energy. In conclusion, our work combines the strength of recent machine learning approaches with the Rosetta protein design toolbox.

MutationExplorer: a webserver for mutation of proteins and 3D visualization of energetic impacts.

The possible effects of mutations on stability and function of a protein can only be understood in the context of protein 3D structure. The MutationExplorer webserver maps sequence changes onto protein structures and allows users to study variation by inputting sequence changes. As the user enters variants, the 3D model evolves, and estimated changes in energy are highlighted. In addition to a basic per-residue input format, MutationExplorer can also upload an entire replacement sequence. Previously the purview of desktop applications, such an upload can back-mutate PDB structures to wildtype sequence in a single step. Another supported variation source is human single nucelotide polymorphisms (SNPs), genomic coordinates input in VCF format. Structures are flexibly colorable, not only by energetic differences, but also by hydrophobicity, sequence conservation, or other biochemical profiling. Coloring by interface score reveals mutation impacts on binding surfaces. MutationExplorer strives for efficiency in user experience. For example, we have prepared 45 000 PDB depositions for instant retrieval and initial display. All modeling steps are performed by Rosetta. Visualizations leverage MDsrv/Mol*. MutationExplorer is available at: http://proteinformatics.org/mutation_explorer/.

Backbone chemical shift and secondary structure assignments for mouse siderocalin.

The lipocalin protein family is a structurally conserved group of proteins with a variety of biological functions defined by their ability to bind small molecule ligands and interact with partner proteins. One member of this family is siderocalin, a protein found in mammals. Its role is discussed in inflammatory processes, iron trafficking, protection against bacterial infections and oxidative stress, cell migration, induction of apoptosis, and cancer. Though it seems to be involved in numerous essential pathways, the exact mechanisms are often not fully understood. The NMR backbone assignments for the human siderocalin and its rat ortholog have been published before. In this work we describe the backbone NMR assignments of siderocalin for another important model organism, the mouse - data that might become important for structure-based drug discovery. Secondary structure elements were predicted based on the assigned backbone chemical shifts using TALOS-N and CSI 3.0, revealing a high content of beta strands and one prominent alpha helical region. Our findings correlate well with the known crystal structure and the overall conserved fold of the lipocalin family.

IHMCIF: An Extension of the PDBx/mmCIF Data Standard for Integrative Structure Determination Methods.

IHMCIF (github.com/ihmwg/IHMCIF) is a data information framework that supports archiving and disseminating macromolecular structures determined by integrative or hybrid modeling (IHM), and making them Findable, Accessible, Interoperable, and Reusable (FAIR). IHMCIF is an extension of the Protein Data Bank Exchange/macromolecular Crystallographic Information Framework (PDBx/mmCIF) that serves as the framework for the Protein Data Bank (PDB) to archive experimentally determined atomic structures of biological macromolecules and their complexes with one another and small molecule ligands (e.g., enzyme cofactors and drugs). IHMCIF serves as the foundational data standard for the PDB-Dev prototype system, developed for archiving and disseminating integrative structures. It utilizes a flexible data representation to describe integrative structures that span multiple spatiotemporal scales and structural states with definitions for restraints from a variety of experimental methods contributing to integrative structural biology. The IHMCIF extension was created with the benefit of considerable community input and recommendations gathered by the Worldwide Protein Data Bank (wwPDB) Task Force for Integrative or Hybrid Methods (wwpdb.org/task/hybrid). Herein, we describe the development of IHMCIF to support evolving methodologies and ongoing advancements in integrative structural biology. Ultimately, IHMCIF will facilitate the unification of PDB-Dev data and tools with the PDB archive so that integrative structures can be archived and disseminated through PDB.

Recent Advances in Automated Structure-Based De Novo Drug Design.

As the number of determined and predicted protein structures and the size of druglike 'make-on-demand' libraries soar, the time-consuming nature of structure-based computer-aided drug design calls for innovative computational algorithms. De novo drug design introduces in silico heuristics to accelerate searching in the vast chemical space. This review focuses on recent advances in structure-based de novo drug design, ranging from conventional fragment-based methods, evolutionary algorithms, and Metropolis Monte Carlo methods to deep generative models. Due to the historical limitation of de novo drug design generating readily available drug-like molecules, we highlight the synthetic accessibility efforts in each category and the benchmarking strategies taken to validate the proposed framework.

Combining machine learning with structure-based protein design to predict and engineer post-translational modifications of proteins.

Post-translational modifications (PTMs) of proteins play a vital role in their function and stability. These modifications influence protein folding, signaling, protein-protein interactions, enzyme activity, binding affinity, aggregation, degradation, and much more. To date, over 400 types of PTMs have been described, representing chemical diversity well beyond the genetically encoded amino acids. Such modifications pose a challenge to the successful design of proteins, but also represent a major opportunity to diversify the protein engineering toolbox. To this end, we first trained artificial neural networks (ANNs) to predict eighteen of the most abundant PTMs, including protein glycosylation, phosphorylation, methylation, and deamidation. In a second step, these models were implemented inside the computational protein modeling suite Rosetta, which allows flexible combination with existing protocols to model the modified sites and understand their impact on protein stability as well as function. Lastly, we developed a new design protocol that either maximizes or minimizes the predicted probability of a particular site being modified. We find that this combination of ANN prediction and structure-based design can enable the modification of existing, as well as the introduction of novel, PTMs. The potential applications of our work include, but are not limited to, glycan masking of epitopes, strengthening protein-protein interactions through phosphorylation, as well as protecting proteins from deamidation liabilities. These applications are especially important for the design of new protein therapeutics where PTMs can drastically change the therapeutic properties of a protein. Our work adds novel tools to Rosetta's protein engineering toolbox that allow for the rational design of PTMs.

Sequence and structural insights of monoleucine-based sorting motifs contained within the cytoplasmic domains of basolateral proteins.

Delivery to the correct membrane domain in polarized epithelial cells is a critical regulatory mechanism for transmembrane proteins. The trafficking of these proteins is directed by short amino acid sequences known as sorting motifs. In six basolaterally-localized proteins lacking the canonical tyrosine- and dileucine-based basolateral sorting motifs, a monoleucine-based sorting motif has been identified. This review will discuss these proteins with an identified monoleucine-based sorting motif, their conserved structural features, as well as the future directions of study for this non-canonical basolateral sorting motif.

MUTATIONEXPLORER- A WEBSERVER FOR MUTATION OF PROTEINS AND 3D VISUALIZATION OF ENERGETIC IMPACTS.

The possible effects of mutations on stability and function of a protein can only be understood in the context of protein 3D structure. The MutationExplorer webserver maps sequence changes onto protein structures and allows users to study variation by inputting sequence changes. As the user enters variants, the 3D model evolves, and estimated changes in energy are highlighted. In addition to a basic per-residue input format, MutationExplorer can also upload an entire replacement sequence. Previously the purview of desktop applications, such an upload can back-mutate PDB structures to wildtype sequence in a single step. Another supported variation source is human single nucelotide polymorphisms (SNPs), genomic coordinates input in VCF format.

Modeling membrane geometries implicitly in Rosetta.

Interactions between membrane proteins (MPs) and lipid bilayers are critical for many cellular functions. In the Rosetta molecular modeling suite, the implicit membrane energy function is based on a "slab" model, which represent the membrane as a flat bilayer. However, in nature membranes often have a curvature that is important for function and/or stability. Even more prevalent, in structural biology research MPs are reconstituted in model membrane systems such as micelles, bicelles, nanodiscs, or liposomes. Thus, we have modified the existing membrane energy potentials within the RosettaMP framework to allow users to model MPs in different membrane geometries. We show that these modifications can be utilized in core applications within Rosetta such as structure refinement, protein-protein docking, and protein design. For MP structures found in curved membranes, refining these structures in curved, implicit membranes produces higher quality models with structures closer to experimentally determined structures. For MP systems embedded in multiple membranes, representing both membranes results in more favorable scores compared to only representing one of the membranes. Modeling MPs in geometries mimicking the membrane model system used in structure determination can improve model quality and model discrimination.

Groebke Blackburn Bienaymé-mediated multi-component synthesis of selective HDAC6 inhibitors with anti-inflammatory properties.

Histone deacetylases (HDACs) are a class of enzymes that cleave acyl groups from lysine residues of histone and non-histone proteins. There are 18 human HDAC isoforms with different cellular targets and functions. Among them, HDAC6 was found to be overexpressed in different types of cancer. However, when used in monotherapy, HDAC6 inhibition by selective inhibitors fails to show pronounced anti-cancer effects. The HDAC6 enzyme also addresses non-histone proteins like α-tubulin and cortactin, making it important for cell migration and angiogenesis. Recently, the NLRP3 inflammasome was identified as an important regulator of inflammation and immune responses and, importantly, HDAC6 is critically involved the activation of the inflammasome. We herein report the design, synthesis and biological evaluation of a library of selective HDAC6 inhibitors. Starting from the previously published crystal structure of MAIP-032 in complex with CD2 of zHDAC6, we performed docking studies to evaluate additional possible interactions of the cap group with the L1-loop pocket. Based on the results we synthesized 13 novel HDAC6 inhibitors via the Groebke-Blackburn-Bienaymé three component reaction as the key step. Compounds 8k (HDAC1 IC50: 5.87 µM; HDAC6 IC50: 0.024 µM; selectivity factor (SF1/6): 245) and 8m (HDAC1 IC50: 3.07 µM; HDAC6 IC50: 0.026 µM; SF1/6: 118) emerged as the most potent and selective inhibitors of HDAC6 and outperformed the lead structure MAIP-032 (HDAC1 IC50: 2.20 µM; HDAC6 IC50: 0.058 µM; SF1/6: 38) both in terms of inhibitory potency and selectivity. Subsequent immunoblot analysis confirmed the high selectivity of 8k and 8m for HDAC6 in a cellular environment. While neither 8k and 8m nor the selectivity HDAC6 inhibitor tubastatin A showed antiproliferative effects in the U-87 MG glioblastoma cell line, compound 8m attenuated cell migration significantly in wound healing assays in U-87 MG cells. Moreover, in macrophages compounds 8k and 8m demonstrated significant inhibition of LPS-induced IL1B mRNA expression and TNF release. These findings suggest that our imidazo[1,2-a]pyridine-capped HDAC6 inhibitors may serve as promising candidates for the development of drugs to effectively treat NLRP3 inflammasome-driven inflammatory diseases.

Approximating Projections of Conformational Boltzmann Distributions with AlphaFold2 Predictions: Opportunities and Limitations.

Protein thermodynamics is intimately tied to biological function and can enable processes such as signal transduction, enzyme catalysis, and molecular recognition. The relative free energies of conformations that contribute to these functional equilibria evolved for the physiology of the organism. Despite the importance of these equilibria for understanding biological function and developing treatments for disease, computational and experimental methods capable of quantifying the energetic determinants of these equilibria are limited to systems of modest size. Recently, it has been demonstrated that the artificial intelligence system AlphaFold2 can be manipulated to produce structurally valid protein conformational ensembles. Here, we extend these studies and explore the extent to which AlphaFold2 contact distance distributions can approximate projections of the conformational Boltzmann distributions. For this purpose, we examine the joint probability distributions of inter-residue contact distances along functionally relevant collective variables of several protein systems. Our studies suggest that AlphaFold2 normalized contact distance distributions can correlate with conformation probabilities obtained with other methods but that they suffer from peak broadening. We also find that the AlphaFold2 contact distance distributions can be sensitive to point mutations. Overall, we anticipate that our findings will be valuable as the community seeks to model the thermodynamics of conformational changes in large biomolecular systems.

Benchmarking AlphaMissense pathogenicity predictions against cystic fibrosis variants.

Variants in the cystic fibrosis transmembrane conductance regulator gene (CFTR) result in cystic fibrosis-a lethal autosomal recessive disorder. Missense variants that alter a single amino acid in the CFTR protein are among the most common cystic fibrosis variants, yet tools for accurately predicting molecular consequences of missense variants have been limited to date. AlphaMissense (AM) is a new technology that predicts the pathogenicity of missense variants based on dual learned protein structure and evolutionary features. Here, we evaluated the ability of AM to predict the pathogenicity of CFTR missense variants. AM predicted a high pathogenicity for CFTR residues overall, resulting in a high false positive rate and fair classification performance on CF variants from the CFTR2.org database. AM pathogenicity score correlated modestly with pathogenicity metrics from persons with CF including sweat chloride level, pancreatic insufficiency rate, and Pseudomonas aeruginosa infection rate. Correlation was also modest with CFTR trafficking and folding competency in vitro. By contrast, the AM score correlated well with CFTR channel function in vitro-demonstrating the dual structure and evolutionary training approach learns important functional information despite lacking such data during training. Different performance across metrics indicated AM may determine if polymorphisms in CFTR are recessive CF variants yet cannot differentiate mechanistic effects or the nature of pathophysiology. Finally, AM predictions offered limited utility to inform on the pharmacological response of CF variants i.e., theratype. Development of new approaches to differentiate the biochemical and pharmacological properties of CFTR variants is therefore still needed to refine the targeting of emerging precision CF therapeutics.

The novel non-hallucinogenic compound DM506 (3-methyl-1,2,3,4,5,6-hexahydroazepino[4,5-b]indole) induces sedative- and anxiolytic-like activity in mice by a mechanism involving 5-HT2A receptor activation.

The anxiolytic and sedative-like effects of 3-methyl-1,2,3,4,5,6-hexahydroazepino[4,5-b]indole (DM506), a non-hallucinogenic compound derived from ibogamine, were studied in mice. The behavioral effects were examined using Elevated O-maze and novelty suppressed feeding (NSFT) tests, open field test, and loss of righting reflex (LORR) test. The results showed that 15 mg/kg DM506 induced acute and long-lasting anxiolytic-like activity in naive and stressed/anxious mice, respectively. Repeated administration of 5 mg/kg DM506 did not cause cumulative anxiolytic activity or any side effects. Higher doses of DM506 (40 mg/kg) induced sedative-like activity, which was inhibited by a selective 5-HT2A receptor antagonist, volinanserin. Electroencephalography results showed that 15 mg/kg DM506 fumarate increased the transition from a highly alert state (fast γ wavelength) to a more synchronized deep-sleeping activity (δ wavelength), which is reflected in the sedative/anxiolytic activity in mice but without the head-twitch response observed in hallucinogens. The functional, radioligand binding, and molecular docking results showed that DM506 binds to the agonist sites of human 5-HT2A (Ki = 24 nM) and 5-HT2B (Ki = 16 nM) receptors and activates them with a potency (EC50) of 9 nM and 3 nM, respectively. DM506 was relatively less potent and behaved as a partial agonist (efficacy <80%) for both receptor subtypes compared to the full agonist DOI (2,5-dimethoxy-4-iodoamphetamine). Our study showed for the first time that the non-hallucinogenic compound DM506 induces anxiolytic- and sedative-like activities in naïve and stressed/anxious mice in a dose-, time-, and volinanserin-sensitive manner, likely through mechanisms involving 5-HT2A receptor activation.

Divergent Folding-Mediated Epistasis Among Unstable Membrane Protein Variants.

Many membrane proteins are prone to misfolding, which compromises their functional expression at the plasma membrane. This is particularly true for the mammalian gonadotropin-releasing hormone receptor GPCRs (GnRHR). We recently demonstrated that evolutionary GnRHR modifications appear to have coincided with adaptive changes in cotranslational folding efficiency. Though protein stability is known to shape evolution, it is unclear how cotranslational folding constraints modulate the synergistic, epistatic interactions between mutations. We therefore compared the pairwise interactions formed by mutations that disrupt the membrane topology (V276T) or tertiary structure (W107A) of GnRHR. Using deep mutational scanning, we evaluated how the plasma membrane expression of these variants is modified by hundreds of secondary mutations. An analysis of 251 mutants in three genetic backgrounds reveals that V276T and W107A form distinct epistatic interactions that depend on both the severity and the mechanism of destabilization. V276T forms predominantly negative epistatic interactions with destabilizing mutations in soluble loops. In contrast, W107A forms positive interactions with mutations in both loops and transmembrane domains that reflect the diminishing impacts of the destabilizing mutations in variants that are already unstable. These findings reveal how epistasis is remodeled by conformational defects in membrane proteins and in unstable proteins more generally.

Dark nanodiscs for evaluating membrane protein thermostability by differential scanning fluorimetry.

Measuring protein thermostability provides valuable information on the biophysical rules that govern the structure-energy relationships of proteins. However, such measurements remain a challenge for membrane proteins. Here, we introduce a new experimental system to evaluate membrane protein thermostability. This system leverages a recently developed nonfluorescent membrane scaffold protein to reconstitute proteins into nanodiscs and is coupled with a nano-format of differential scanning fluorimetry (nanoDSF). This approach offers a label-free and direct measurement of the intrinsic tryptophan fluorescence of the membrane protein as it unfolds in solution without signal interference from the "dark" nanodisc. In this work, we demonstrate the application of this method using the disulfide bond formation protein B (DsbB) as a test membrane protein. NanoDSF measurements of DsbB reconstituted in dark nanodiscs loaded with 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 1,2-dimyristoyl-sn-glycero-3-phosphorylglycerol (DMPG) lipids show a complex biphasic thermal unfolding pattern with a minor unfolding transition followed by a major transition. The inflection points of the thermal denaturation curve reveal two distinct unfolding midpoint melting temperatures (Tm) of 70.5°C and 77.5°C, consistent with a three-state unfolding model. Further, we show that the catalytically conserved disulfide bond between residues C41 and C130 drives the intermediate state of the unfolding pathway for DsbB in a DMPC and DMPG nanodisc. To extend the utility of this method, we evaluate and compare the thermostability of DsbB in different lipid environments. We introduce this method as a new tool that can be used to understand how compositionally and biophysically complex lipid environments drive membrane protein stability.

Benchmarking AlphaMissense Pathogenicity Predictions Against Cystic Fibrosis Variants.

Variants in the cystic fibrosis transmembrane conductance regulator gene (CFTR) result in cystic fibrosis - a lethal autosomal recessive disorder. Missense variants that alter a single amino acid in the CFTR protein are among the most common cystic fibrosis variants, yet tools for accurately predicting molecular consequences of missense variants have been limited to date. AlphaMissense (AM) is a new technology that predicts the pathogenicity of missense variants based on dual learned protein structure and evolutionary features. Here, we evaluated the ability of AM to predict the pathogenicity of CFTR missense variants. AM predicted a high pathogenicity for CFTR residues overall, resulting in a high false positive rate and fair classification performance on CF variants from the CFTR2.org database. AM pathogenicity score correlated modestly with pathogenicity metrics from persons with CF including sweat chloride level, pancreatic insufficiency rate, and Pseudomonas aeruginosa infection rate. Correlation was also modest with CFTR trafficking and folding competency in vitro. By contrast, the AM score correlated well with CFTR channel function in vitro - demonstrating the dual structure and evolutionary training approach learns important functional information despite lacking such data during training. Different performance across metrics indicated AM may determine if polymorphisms in CFTR are recessive CF variants yet cannot differentiate mechanistic effects or the nature of pathophysiology. Finally, AM predictions offered limited utility to inform on the pharmacological response of CF variants i.e., theratype. Development of new approaches to differentiate the biochemical and pharmacological properties of CFTR variants is therefore still needed to refine the targeting of emerging precision CF therapeutics.

Tissue factor binds to and inhibits interferon-α receptor 1 signaling.

Tissue factor (TF), which is a member of the cytokine receptor family, promotes coagulation and coagulation-dependent inflammation. TF also exerts protective effects through unknown mechanisms. Here, we showed that TF bound to interferon-α receptor 1 (IFNAR1) and antagonized its signaling, preventing spontaneous sterile inflammation and maintaining immune homeostasis. Structural modeling and direct binding studies revealed binding of the TF C-terminal fibronectin III domain to IFNAR1, which restricted the expression of interferon-stimulated genes (ISGs). Podocyte-specific loss of TF in mice (PodΔF3) resulted in sterile renal inflammation, characterized by JAK/STAT signaling, proinflammatory cytokine expression, disrupted immune homeostasis, and glomerulopathy. Inhibiting IFNAR1 signaling or loss of Ifnar1 expression in podocytes attenuated these effects in PodΔF3 mice. As a heteromer, TF and IFNAR1 were both inactive, while dissociation of the TF-IFNAR1 heteromer promoted TF activity and IFNAR1 signaling. These data suggest that the TF-IFNAR1 heteromer is a molecular switch that controls thrombo-inflammation.

Time-resolved live-cell spectroscopy reveals EphA2 multimeric assembly.

Ephrin type-A receptor 2 (EphA2) is a receptor tyrosine kinase that initiates both ligand-dependent tumor-suppressive and ligand-independent oncogenic signaling. We used time-resolved, live-cell fluorescence spectroscopy to show that the ligand-free EphA2 assembles into multimers driven by two types of intermolecular interactions in the ectodomain. The first type entails extended symmetric interactions required for ligand-induced receptor clustering and tumor-suppressive signaling that inhibits activity of the oncogenic extracellular signal-regulated kinase (ERK) and protein kinase B (AKT) protein kinases and suppresses cell migration. The second type is an asymmetric interaction between the amino terminus and the membrane proximal domain of the neighboring receptors, which supports oncogenic signaling and promotes migration in vitro and tumor invasiveness in vivo. Our results identify the molecular interactions that drive the formation of the EphA2 multimeric signaling clusters and reveal the pivotal role of EphA2 assembly in dictating its opposing functions in oncogenesis.