CryptoSite: Expanding the Druggable Proteome by Characterization and Prediction of Cryptic Binding Sites

Exploring the Application of SiteMap and Site Finder for Focused Cryptic Pocket Identification

The characterization of cryptic pockets has been elusive, despite substantial efforts. Computational modeling approaches, such as molecular dynamics (MD) simulations, can provide atomic-level details of binding site motions and binding pathways. However, the time scale that MD can achieve at a reasonable cost often limits its application for cryptic pocket identification. Enhanced sampling techniques can improve the efficiency of MD simulations by focused sampling of important regions of the protein, but prior knowledge of the simulated system is required to define the appropriate coordinates. In the case of a novel, unknown cryptic pocket, such information is not available, limiting the application of enhanced sampling techniques for cryptic pocket identification. In this work, we explore the ability of SiteMap and Site Finder, widely used commercial packages for pocket identification, to detect focus points on the protein and further apply other advanced computational methods. The information gained from this analysis enables the use of computational modeling, including enhanced MD sampling techniques, to explore potential cryptic binding pockets suggested by SiteMap and Site Finder. Here, we examined SiteMap and Site Finder results on 136 known cryptic pockets from a combination of the PocketMiner dataset (a recently curated set of cryptic pockets), the Cryptosite Set (a classic set of cryptic pockets), and Natural killer group 2D (NKG2D, a protein target where a cryptic pocket is confirmed). Our findings demonstrate the application of existing, well-studied tools in efficiently mapping potential regions harboring cryptic pockets.

The Malaria Box molecules: a source for targeting the RBD and NTD cryptic pocket of the spike glycoprotein in SARS-CoV-2

CONTEXT

SARS-CoV-2, responsible for COVID-19, has led to over 500 million infections and more than 6 million deaths globally. There have been limited effective treatments available. The study aims to find a drug that can prevent the virus from entering host cells by targeting specific sites on the virus's spike protein.

METHOD

We examined 13,397 compounds from the Malaria Box library against two specific sites on the spike protein: the receptor-binding domain (RBD) and a predicted cryptic pocket. Using virtual screening, molecular docking, molecular dynamics, and MMPBSA techniques, they evaluated the stability of two compounds. TCMDC-124223 showed high stability and binding energy in the RBD, while TCMDC-133766 had better binding energy in the cryptic pocket. The study also identified that the interacting residues are conserved, which is crucial for addressing various virus variants. The findings provide insights into the potential of small molecules as drugs against the spike protein.

Application of artificial intelligence and machine learning techniques to the analysis of dynamic protein sequences

We apply methods of Artificial Intelligence and Machine Learning to protein dynamic bioinformatics. We rewrite the sequences of a large protein data set, containing both folded and intrinsically disordered molecules, using a representation developed previously, which encodes the intrinsic dynamic properties of the naturally occurring amino acids. We Fourier analyze the resulting sequences. It is demonstrated that classification models built using several different supervised learning methods are able to successfully distinguish folded from intrinsically disordered proteins from sequence alone. It is further shown that the most important sequence property for this discrimination is the sequence mobility, which is the sequence averaged value of the residue-specific average alpha carbon B factor. This is in agreement with previous work, in which we have demonstrated the central role played by the sequence mobility in protein dynamic bioinformatics and biophysics. This finding opens a path to the application of dynamic bioinformatics, in combination with machine learning algorithms, to a range of significant biomedical problems.

Secondary Sites of the C-type Lectin-Like Fold

C-type lectins are a large superfamily of proteins involved in a multitude of biological processes. In particular, their involvement in immunity and homeostasis has rendered them attractive targets for diverse therapeutic interventions. They share a characteristic C-type lectin-like domain whose adaptability enables them to bind a broad spectrum of ligands beyond the originally defined canonical Ca2+-dependent carbohydrate binding. Together with variable domain architecture and high-level conformational plasticity, this enables C-type lectins to meet diverse functional demands. Secondary sites provide another layer of regulation and are often intricately linked to functional diversity. Located remote from the canonical primary binding site, secondary sites can accommodate ligands with other physicochemical properties and alter protein dynamics, thus enhancing selectivity and enabling fine-tuning of the biological response. In this review, we outline the structural determinants allowing C-type lectins to perform a large variety of tasks and to accommodate the ligands associated with it. Using the six well-characterized Ca2+-dependent and Ca2+-independent C-type lectin receptors DC-SIGN, langerin, MGL, dectin-1, CLEC-2 and NKG2D as examples, we focus on the characteristics of non-canonical interactions and secondary sites and their potential use in drug discovery endeavors.

SWISH-X, an Expanded Approach to Detect Cryptic Pockets in Proteins and at Protein-Protein Interfaces

Protein-protein interactions mediate most molecular processes in the cell, offering a significant opportunity to expand the set of known druggable targets. Unfortunately, targeting these interactions can be challenging due to their typically flat and featureless interaction surfaces, which often change as the complex forms. Such surface changes may reveal hidden (cryptic) druggable pockets. Here, we analyze a set of well-characterized protein-protein interactions harboring cryptic pockets and investigate the predictive power of current computational methods. Based on our observations, we developed a new computational strategy, SWISH-X (SWISH Expanded), which combines the established cryptic pocket identification capabilities of SWISH with the rapid temperature range exploration of OPES MultiThermal. SWISH-X is able to reliably identify cryptic pockets at protein-protein interfaces while retaining its predictive power for revealing cryptic pockets in isolated proteins, such as TEM-1 β-lactamase.

Protein Ensemble Generation Through Variational Autoencoder Latent Space Sampling

Mapping the ensemble of protein conformations that contribute to function and can be targeted by small molecule drugs remains an outstanding challenge. Here, we explore the use of variational autoencoders for reducing the challenge of dimensionality in the protein structure ensemble generation problem. We convert high-dimensional protein structural data into a continuous, low-dimensional representation, carry out a search in this space guided by a structure quality metric, and then use RoseTTAFold guided by the sampled structural information to generate 3D structures. We use this approach to generate ensembles for the cancer relevant protein K-Ras, train the VAE on a subset of the available K-Ras crystal structures and MD simulation snapshots, and assess the extent of sampling close to crystal structures withheld from training. We find that our latent space sampling procedure rapidly generates ensembles with high structural quality and is able to sample within 1 Å of held-out crystal structures, with a consistency higher than that of MD simulation or AlphaFold2 prediction. The sampled structures sufficiently recapitulate the cryptic pockets in the held-out K-Ras structures to allow for small molecule docking.

An allosteric mechanism for potent inhibition of SARS-CoV-2 main proteinase

The main proteinase (Mpro) of SARS-CoV-2 plays a critical role in cleaving viral polyproteins into functional proteins required for viral replication and assembly, making it a prime drug target for COVID-19. It is well known that noncompetitive inhibition offers potential therapeutic options for treating COVID-19, which can effectively reduce the likelihood of cross-reactivity with other proteins and increase the selectivity of the drug. Therefore, the discovery of allosteric sites of Mpro has both scientific and practical significance. In this study, we explored the binding characteristics and inhibiting process of Mpro activity by two recently reported allosteric inhibitors, pelitinib and AT7519 which were obtained by the X-ray screening experiments, to probe the allosteric mechanism via molecular dynamic (MD) simulations. We found that pelitinib and AT7519 can stably bind to Mpro far from the active site. The binding affinity is estimated to be -24.37 ± 4.14 and - 26.96 ± 4.05 kcal/mol for pelitinib and AT7519, respectively, which is considerably stable compared with orthosteric drugs. Furthermore, the strong binding caused clear changes in the catalytic site of Mpro, thus decreasing the substrate accessibility. The community network analysis also validated that pelitinib and AT7519 strengthened intra- and inter-domain communication of Mpro dimer, resulting in a rigid Mpro, which could negatively impact substrate binding. In summary, our findings provide the detailed working mechanism for the two experimentally observed allosteric sites of Mpro. These allosteric sites greatly enhance the 'druggability' of Mpro and represent attractive targets for the development of new Mpro inhibitors.

AHoJ-DB: A PDB-wide assignment of apo & holo relationships based on individual protein-ligand interactions

A single protein structure is rarely sufficient to capture the conformational variability of a protein. Both bound and unbound (holo and apo) forms of a protein are essential for understanding its geometry and making meaningful comparisons. Nevertheless, docking or drug design studies often still consider only single protein structures in their holo form, which are for the most part rigid. With the recent explosion in the field of structural biology, large, curated datasets are urgently needed. Here, we use a previously developed application (AHoJ) to perform a comprehensive search for apo-holo pairs for 468,293 biologically relevant protein-ligand interactions across 27,983 proteins. In each search, the binding pocket is captured and mapped across existing structures within the same UniProt, and the mapped pockets are annotated as apo or holo, based on the presence or absence of ligands. We assemble the results into a database, AHoJ-DB (www.apoholo.cz/db), that captures the variability of proteins with identical sequences, thereby exposing the agents responsible for the observed differences in geometry. We report several metrics for each annotated pocket, and we also include binding pockets that form at the interface of multiple chains. Analysis of the database shows that about 24% of the binding sites occur at the interface of two or more chains and that less than 50% of the total binding sites processed have an apo form in the PDB. These results can be used to train and evaluate predictors, discover potentially druggable proteins, and reveal protein- and ligand-specific relationships that were previously obscured by intermittent or partial data. Availability: http://apoholo.cz/db.

Toward physics-based precision medicine: Exploiting protein dynamics to design new therapeutics and interpret variants

The goal of precision medicine is to utilize our knowledge of the molecular causes of disease to better diagnose and treat patients. However, there is a substantial mismatch between the small number of food and drug administration (FDA)-approved drugs and annotated coding variants compared to the needs of precision medicine. This review introduces the concept of physics-based precision medicine, a scalable framework that promises to improve our understanding of sequence-function relationships and accelerate drug discovery. We show that accounting for the ensemble of structures a protein adopts in solution with computer simulations overcomes many of the limitations imposed by assuming a single protein structure. We highlight studies of protein dynamics and recent methods for the analysis of structural ensembles. These studies demonstrate that differences in conformational distributions predict functional differences within protein families and between variants. Thanks to new computational tools that are providing unprecedented access to protein structural ensembles, this insight may enable accurate predictions of variant pathogenicity for entire libraries of variants. We further show that explicitly accounting for protein ensembles, with methods like alchemical free energy calculations or docking to Markov state models, can uncover novel lead compounds. To conclude, we demonstrate that cryptic pockets, or cavities absent in experimental structures, provide an avenue to target proteins that are currently considered undruggable. Taken together, our review provides a roadmap for the field of protein science to accelerate precision medicine.

Synergy of Mutation-Induced Effects in Human Vitamin K Epoxide Reductase: Perspectives and Challenges for Allo-Network Modulator Design

The human Vitamin K Epoxide Reductase Complex (hVKORC1), a key enzyme transforming vitamin K into the form necessary for blood clotting, requires for its activation the reducing equivalents delivered by its redox partner through thiol-disulfide exchange reactions. The luminal loop (L-loop) is the principal mediator of hVKORC1 activation, and it is a region frequently harbouring numerous missense mutations. Four L-loop hVKORC1 mutants, suggested in vitro as either resistant (A41S, H68Y) or completely inactive (S52W, W59R), were studied in the oxidised state by numerical approaches (in silico). The DYNASOME and POCKETOME of each mutant were characterised and compared to the native protein, recently described as a modular protein composed of the structurally stable transmembrane domain (TMD) and the intrinsically disordered L-loop, exhibiting quasi-independent dynamics. The DYNASOME of mutants revealed that L-loop missense point mutations impact not only its folding and dynamics, but also those of the TMD, highlighting a strong mutation-specific interdependence between these domains. Another consequence of the mutation-induced effects manifests in the global changes (geometric, topological, and probabilistic) of the newly detected cryptic pockets and the alternation of the recognition properties of the L-loop with its redox protein. Based on our results, we postulate that (i) intra-protein allosteric regulation and (ii) the inherent allosteric regulation and cryptic pockets of each mutant depend on its DYNASOME; and (iii) the recognition of the redox protein by hVKORC1 (INTERACTOME) depend on their DYNASOME. This multifaceted description of proteins produces "omics" data sets, crucial for understanding the physiological processes of proteins and the pathologies caused by alteration of the protein properties at various "omics" levels. Additionally, such characterisation opens novel perspectives for the development of "allo-network drugs" essential for the treatment of blood disorders.

DynamicBind: predicting ligand-specific protein-ligand complex structure with a deep equivariant generative model

While significant advances have been made in predicting static protein structures, the inherent dynamics of proteins, modulated by ligands, are crucial for understanding protein function and facilitating drug discovery. Traditional docking methods, frequently used in studying protein-ligand interactions, typically treat proteins as rigid. While molecular dynamics simulations can propose appropriate protein conformations, they're computationally demanding due to rare transitions between biologically relevant equilibrium states. In this study, we present DynamicBind, a deep learning method that employs equivariant geometric diffusion networks to construct a smooth energy landscape, promoting efficient transitions between different equilibrium states. DynamicBind accurately recovers ligand-specific conformations from unbound protein structures without the need for holo-structures or extensive sampling. Remarkably, it demonstrates state-of-the-art performance in docking and virtual screening benchmarks. Our experiments reveal that DynamicBind can accommodate a wide range of large protein conformational changes and identify cryptic pockets in unseen protein targets. As a result, DynamicBind shows potential in accelerating the development of small molecules for previously undruggable targets and expanding the horizons of computational drug discovery.

Binding Mechanism of Riboswitch to Natural Ligand Elucidated by McMD-Based Dynamic Docking Simulations

Flavin mononucleotide riboswitches are common among many pathogenic bacteria and are therefore considered to be an attractive target for antibiotics development. The riboswitch binds riboflavin (RBF, also known as vitamin B2), and although an experimental structure of their complex has been solved with the ligand bound deep inside the RNA molecule in a seemingly unreachable state, the binding mechanism between these molecules is not yet known. We have therefore used our Multicanonical Molecular Dynamics (McMD)-based dynamic docking protocol to analyze their binding mechanism by simulating the binding process between the riboswitch aptamer domain and the RBF, starting from the apo state of the riboswitch. Here, the refinement stage was crucial to identify the native binding configuration, as several other binding configurations were also found by McMD-based docking simulations. RBF initially binds the interface between P4 and P6 including U61 and G62, which forms a gateway where the ligand lingers until this gateway opens sufficiently to allow the ligand to pass through and slip into the hidden binding site including A48, A49, and A85.

Targeting phosphatases: From molecule design to clinical trials

Phosphatase is a kind of enzyme that can dephosphorylate target proteins, which can be divided into serine/threonine phosphatase and tyrosine phosphatase according to its mode of action. Current evidence showed multiple phosphatases were highly correlated with diseases including various cancers, demonstrating them as potential targets. However, currently, targeting phosphatases with small molecules faces many challenges, resulting in no drug approved. In this case, phosphatases are even regarded as "undruggable" targets for a long time. Recently, a variety of strategies have been adopted in the design of small molecule inhibitors targeting phosphatases, leading many of them to enter into the clinical trials. In this review, we classified these inhibitors into 4 types, including (1) molecular glues, (2) small molecules targeting catalytic sites, (3) allosteric inhibition, and (4) bifunctional molecules (proteolysis targeting chimeras, PROTACs). These molecules with diverse strategies prove the feasibility of phosphatases as drug targets. In addition, the combination therapy of phosphatase inhibitors with other drugs has also entered clinical trials, which suggests a broad prospect. Thus, targeting phosphatases with small molecules by different strategies is emerging as a promising way in the modulation of pathogenetic phosphorylation.

A Structure-Based Allosteric Modulator Design Paradigm

Importance: Allosteric drugs bound to topologically distal allosteric sites hold a substantial promise in modulating therapeutic targets deemed undruggable at their orthosteric sites. Traditionally, allosteric modulator discovery has predominantly relied on serendipitous high-throughput screening. Nevertheless, the landscape has undergone a transformative shift due to recent advancements in our understanding of allosteric modulation mechanisms, coupled with a significant increase in the accessibility of allosteric structural data. These factors have extensively promoted the development of various computational methodologies, especially for machine-learning approaches, to guide the rational design of structure-based allosteric modulators. Highlights: We here presented a comprehensive structure-based allosteric modulator design paradigm encompassing 3 critical stages: drug target acquisition, allosteric binding site, and modulator discovery. The recent advances in computational methods in each stage are encapsulated. Furthermore, we delve into analyzing the successes and obstacles encountered in the rational design of allosteric modulators. Conclusion: The structure-based allosteric modulator design paradigm holds immense potential for the rational design of allosteric modulators. We hope that this review would heighten awareness of the use of structure-based computational methodologies in advancing the field of allosteric drug discovery.

Integrated Molecular Modeling and Machine Learning for Drug Design

Modern therapeutic development often involves several stages that are interconnected, and multiple iterations are usually required to bring a new drug to the market. Computational approaches have increasingly become an indispensable part of helping reduce the time and cost of the research and development of new drugs. In this Perspective, we summarize our recent efforts on integrating molecular modeling and machine learning to develop computational tools for modulator design, including a pocket-guided rational design approach based on AlphaSpace to target protein-protein interactions, delta machine learning scoring functions for protein-ligand docking as well as virtual screening, and state-of-the-art deep learning models to predict calculated and experimental molecular properties based on molecular mechanics optimized geometries. Meanwhile, we discuss remaining challenges and promising directions for further development and use a retrospective example of FDA approved kinase inhibitor Erlotinib to demonstrate the use of these newly developed computational tools.

In Silico Exploration of the Trypanothione Reductase (TryR) of L. mexicana

Human leishmaniasis is a neglected tropical disease which affects nearly 1.5 million people every year, with Mexico being an important endemic region. One of the major defense mechanisms of these parasites is based in the polyamine metabolic pathway, as it provides the necessary compounds for its survival. Among the enzymes in this route, trypanothione reductase (TryR), an oxidoreductase enzyme, is crucial for the Leishmania genus' survival against oxidative stress. Thus, it poses as an attractive drug target, yet due to the size and features of its catalytic pocket, modeling techniques such as molecular docking focusing on that region is not convenient. Herein, we present a computational study using several structure-based approaches to assess the druggability of TryR from L. mexicana, the predominant Leishmania species in Mexico, beyond its catalytic site. Using this consensus methodology, three relevant pockets were found, of which the one we call σ-site promises to be the most favorable one. These findings may help the design of new drugs of trypanothione-related diseases.

Mining gene expression data for rational identification of novel drug targets and vaccine candidates against the cattle tick, Rhipicephalus microplus

Control of complex parasites via vaccination remains challenging, with the current combination of vaccines and small drugs remaining the choice for an integrated control strategy. Studies conducted to date, are providing evidence that multicomponent vaccines will be needed for the development of protective vaccines against endo- and ectoparasites, though multicomponent vaccines require an in-depth understanding of parasite biology which remains insufficient for ticks. With the rapid development and spread of acaricide resistance in ticks, new targets for acaricide development also remains to be identified, along with novel targets that can be exploited for the design of lead compounds. In this study, we analysed the differential gene expression of Rhipicephalus microplus ticks that were fed on cattle vaccinated with a multi-component vaccine (Bm86 and 3 putative Bm86-binding proteins). The data was scrutinised for the identification of vaccine targets, small drug targets and novel pathways that can be evaluated in future studies. Limitations associated with targeting novel proteins for vaccine and/or drug design is also discussed and placed into the context of challenges arising when targeting large protein families and intracellular localised proteins. Lastly, this study provide insight into how Bm86-based vaccines may reduce successful uptake and digestion of the bloodmeal and overall tick fecundity.

Exploring Conformational Landscapes and Cryptic Binding Pockets in Distinct Functional States of the SARS-CoV-2 Omicron BA.1 and BA.2 Trimers: Mutation-Induced Modulation of Protein Dynamics and Network-Guided Prediction of Variant-Specific Allosteric Binding Sites

A significant body of experimental structures of SARS-CoV-2 spike trimers for the BA.1 and BA.2 variants revealed a considerable plasticity of the spike protein and the emergence of druggable binding pockets. Understanding the interplay of conformational dynamics changes induced by the Omicron variants and the identification of cryptic dynamic binding pockets in the S protein is of paramount importance as exploring broad-spectrum antiviral agents to combat the emerging variants is imperative. In the current study, we explore conformational landscapes and characterize the universe of binding pockets in multiple open and closed functional spike states of the BA.1 and BA.2 Omicron variants. By using a combination of atomistic simulations, a dynamics network analysis, and an allostery-guided network screening of binding pockets in the conformational ensembles of the BA.1 and BA.2 spike conformations, we identified all experimentally known allosteric sites and discovered significant variant-specific differences in the distribution of binding sites in the BA.1 and BA.2 trimers. This study provided a structural characterization of the predicted cryptic pockets and captured the experimentally known allosteric sites, revealing the critical role of conformational plasticity in modulating the distribution and cross-talk between functional binding sites. We found that mutational and dynamic changes in the BA.1 variant can induce the remodeling and stabilization of a known druggable pocket in the N-terminal domain, while this pocket is drastically altered and may no longer be available for ligand binding in the BA.2 variant. Our results predicted the experimentally known allosteric site in the receptor-binding domain that remains stable and ranks as the most favorable site in the conformational ensembles of the BA.2 variant but could become fragmented and less probable in BA.1 conformations. We also uncovered several cryptic pockets formed at the inter-domain and inter-protomer interface, including functional regions of the S2 subunit and stem helix region, which are consistent with the known role of pocket residues in modulating conformational transitions and antibody recognition. The results of this study are particularly significant for understanding the dynamic and network features of the universe of available binding pockets in spike proteins, as well as the effects of the Omicron-variant-specific modulation of preferential druggable pockets. The exploration of predicted druggable sites can present a new and previously underappreciated opportunity for therapeutic interventions for Omicron variants through the conformation-selective and variant-specific targeting of functional sites involved in allosteric changes.

Prediction of the Peptide-TIM3 Binding Site in Inhibiting TIM3-Galectin 9 Binding Pathways

T-cell immunoglobulin and mucin domain-containing protein-3 (TIM3) is an important receptor protein that modulates the immune system. The binding of TIM3 with Galectin 9 (GAL9) triggers immune system suppression, but the TIM3-GAL9 binding can be inhibited by binding of the peptide P26 to TIM3. A fast and accurate prediction of the P26-TIM3 binding site is crucial and a prerequisite for the investigation of P26-TIM3 interactions and TIM3-GAL9 binding pathways. Here, we present a machine learning approach, which considers protein conformational changes, to quickly identify the ligand-binding site on TIM3. Our results show that the P26 binding site is located near the C″-D loop of TIM3. Further simulations show that the binding pose is stabilized by a variety of electrostatic and hydrophobic interactions. Binding of P26 can alter the conformations of nearby glycan side chains on TIM3, providing possible mechanisms of how P26 inhibits TIM3-GAL9 binding pathways. The insights from this work will facilitate the identification of other peptides or antibodies that may also inhibit the TIM3-GAL9 pathways and eventually lead to improved attempts in the modulation of the TIM3-GAL9 immunosuppression pathways. The strategies and machine learning method can be generalized to study ligand-receptor binding when the conformational changes during the binding are important.

Worth the Weight: Sub-Pocket EXplorer (SubPEx), a Weighted Ensemble Method to Enhance Binding-Pocket Conformational Sampling

Structure-based virtual screening (VS) is an effective method for identifying potential small-molecule ligands, but traditional VS approaches consider only a single binding-pocket conformation. Consequently, they struggle to identify ligands that bind to alternate conformations. Ensemble docking helps address this issue by incorporating multiple conformations into the docking process, but it depends on methods that can thoroughly explore pocket flexibility. We here introduce Sub-Pocket EXplorer (SubPEx), an approach that uses weighted ensemble (WE) path sampling to accelerate binding-pocket sampling. As proof of principle, we apply SubPEx to three proteins relevant to drug discovery: heat shock protein 90, influenza neuraminidase, and yeast hexokinase 2. SubPEx is available free of charge without registration under the terms of the open-source MIT license: http://durrantlab.com/subpex/.

Discovery of a selective and biologically active low-molecular weight antagonist of human interleukin-1β

Human interleukin-1β (hIL-1β) is a pro-inflammatory cytokine involved in many diseases. While hIL-1β directed antibodies have shown clinical benefit, an orally available low-molecular weight antagonist is still elusive, limiting the applications of hIL-1β-directed therapies. Here we describe the discovery of a low-molecular weight hIL-1β antagonist that blocks the interaction with the IL-1R1 receptor. Starting from a low affinity fragment-based screening hit 1, structure-based optimization resulted in a compound (S)-2 that binds and antagonizes hIL-1β with single-digit micromolar activity in biophysical, biochemical, and cellular assays. X-ray analysis reveals an allosteric mode of action that involves a hitherto unknown binding site in hIL-1β encompassing two loops involved in hIL-1R1/hIL-1β interactions. We show that residues of this binding site are part of a conformationally excited state of the mature cytokine. The compound antagonizes hIL-1β function in cells, including primary human fibroblasts, demonstrating the relevance of this discovery for future development of hIL-1β directed therapeutics.

Grasping cryptic binding sites to neutralize drug resistance in the field of anticancer

Drug resistance is a significant obstacle to successful cancer treatment. The utilization and development of cryptic binding sites (CBSs) in proteins involved in cancer-related drug-resistance (CRDR) could help to overcome that drug resistance. However, there is no comprehensive review of the successful use of CBSs in addressing CRDR. Here, we have systematically summarized and analyzed the opportunities and challenges of using CBSs in addressing CRDR and revealed the key role that CBSs have in targeting CRDR. First, we have identified the CRDR targets and the corresponding CBSs. Second, we discuss the mechanisms by which CBSs can overcome CRDR. Finally, we have provided examples of successful CBS applications in addressing CRDR. We hope that this approach will provide guidance to biologists and chemists in effectively utilizing CBSs for the development of new drugs to alleviate CRDR.

Recent Developments in Ultralarge and Structure-Based Virtual Screening Approaches

Drug development is a wide scientific field that faces many challenges these days. Among them are extremely high development costs, long development times, and a small number of new drugs that are approved each year. New and innovative technologies are needed to solve these problems that make the drug discovery process of small molecules more time and cost efficient, and that allow previously undruggable receptor classes to be targeted, such as protein-protein interactions. Structure-based virtual screenings (SBVSs) have become a leading contender in this context. In this review, we give an introduction to the foundations of SBVSs and survey their progress in the past few years with a focus on ultralarge virtual screenings (ULVSs). We outline key principles of SBVSs, recent success stories, new screening techniques, available deep learning-based docking methods, and promising future research directions. ULVSs have an enormous potential for the development of new small-molecule drugs and are already starting to transform early-stage drug discovery.

Worth the weight: Sub-Pocket EXplorer (SubPEx), a weighted-ensemble method to enhance binding-pocket conformational sampling

Structure-based virtual screening (VS) is an effective method for identifying potential small-molecule ligands, but traditional VS approaches consider only a single binding-pocket conformation. Consequently, they struggle to identify ligands that bind to alternate conformations. Ensemble docking helps address this issue by incorporating multiple conformations into the docking process, but it depends on methods that can thoroughly explore pocket flexibility. We here introduce Sub-Pocket EXplorer (SubPEx), an approach that uses weighted ensemble (WE) path sampling to accelerate binding-pocket sampling. As proof of principle, we apply SubPEx to three proteins relevant to drug discovery: heat shock protein 90, influenza neuraminidase, and yeast hexokinase 2. SubPEx is available free of charge without registration under the terms of the open-source MIT license: http://durrantlab.com/subpex/.

Binding Site Detection Remastered: Enabling Fast, Robust, and Reliable Binding Site Detection and Descriptor Calculation with DoGSite3

Binding site prediction on protein structures is a crucial step in early phase drug discovery whenever experimental or predicted structure models are involved. DoGSite belongs to the widely used tools for this task. It is a grid-based method that uses a Difference-of-Gaussian filter to detect cavities on the protein surface. We recently reimplemented the first version of this method, released in 2010, focusing on improved binding site detection in the presence of ligands and optimized parameters for more robust, reliable, and fast predictions and binding site descriptor calculations. Here, we introduce the new version, DoGSite3, compare it to its predecessor, and re-evaluate DoGSite on published data sets for a large-scale comparative performance evaluation.

Bioorthogonal Tethering Enhances Drug Fragment Affinity for G Protein-Coupled Receptors in Live Cells

G protein-coupled receptors (GPCRs) modulate diverse cellular signaling pathways and are important drug targets. Despite the availability of high-resolution structures, the discovery of allosteric modulators remains challenging due to the dynamic nature of GPCRs in native membranes. We developed a strategy to covalently tether drug fragments adjacent to allosteric sites in GPCRs to enhance their potency and enable fragment-based drug screening in cell-based systems. We employed genetic code expansion to site-specifically introduce noncanonical amino acids with reactive groups in C-C chemokine receptor 5 (CCR5) near an allosteric binding site for the drug maraviroc. We then used molecular dynamics simulations to design heterobifunctional maraviroc analogues consisting of a drug fragment connected by a flexible linker to a reactive moiety capable of undergoing a bioorthogonal coupling reaction. We synthesized a library of these analogues and employed the bioorthogonal inverse electron demand Diels-Alder reaction to couple the analogues to the engineered CCR5 in live cells, which were then assayed using cell-based signaling assays. Tetherable low-affinity maraviroc fragments displayed an increase in potency for CCR5 engineered with reactive unnatural amino acids that were adjacent to the maraviroc binding site. The strategy we describe to tether novel drug fragments to GPCRs should prove useful to probe allosteric or cryptic binding site functionality in fragment-based GPCR-targeted drug discovery.

Mutual induced-fit mechanism drives binding between intrinsically disordered Bim and cryptic binding site of Bcl-xL

The intrinsically disordered region (IDR) of Bim binds to the flexible cryptic site of Bcl-xL, a pro-survival protein involved in cancer progression that plays an important role in initiating apoptosis. However, their binding mechanism has not yet been elucidated. We have applied our dynamic docking protocol, which correctly reproduced both the IDR properties of Bim and the native bound configuration, as well as suggesting other stable/meta-stable binding configurations and revealed the binding pathway. Although the cryptic site of Bcl-xL is predominantly in a closed conformation, initial binding of Bim in an encounter configuration leads to mutual induced-fit binding, where both molecules adapt to each other; Bcl-xL transitions to an open state as Bim folds from a disordered to an α-helical conformation while the two molecules bind each other. Finally, our data provides new avenues to develop novel drugs by targeting newly discovered stable conformations of Bcl-xL.

Predicting allosteric sites using fast conformational sampling as guided by coarse-grained normal modes

To computationally identify cryptic binding sites for allosteric modulators, we have developed a fast and simple conformational sampling scheme guided by coarse-grained normal modes solved from the elastic network models followed by atomistic backbone and sidechain reconstruction. Despite the complexity of conformational changes associated with ligand binding, we previously showed that simply sampling along each of the lowest 30 modes can adequately restructure cryptic sites so they are detectable by pocket finding programs like Concavity. Here, we applied this method to study four classical examples of allosteric regulation (GluR2 receptor, GroEL chaperonin, GPCR, and myosin). Our method along with alternative methods has been utilized to locate known allosteric sites and predict new promising allosteric sites. Compared with other sampling methods based on extensive molecular dynamics simulation, our method is both faster (1-2 h for an average-size protein of ∼400 residues) and more flexible (it can be easily integrated with any structure-based pocket finding methods), so it is suitable for high-throughput screening of large datasets of protein structures at the genome scale.

Identification of a new class of proteasome inhibitors based on a naphthyl-azotricyclic-urea-phenyl scaffold

Proteasomes play an important role in protein degradation and regulation of many cellular pathways by maintaining protein balance. Inhibitors of proteasomes disrupt this balance affecting proteins that are key in malignancies and as such have found applications in the treatment of multiple myeloma and mantle cell lymphoma. However, resistance mechanisms have been reported for these proteasome inhibitors including mutations at the β5 site which necessitates the constant development of new inhibitors. In this work, we report the identification of a new class of proteasome inhibitors, polycyclic molecules bearing a naphthyl-azotricyclic-urea-phenyl scaffold, from screening of the ZINC library of natural products. The most potent of these compounds showed evidence of dose dependency through proteasome assays with IC50 values in the low micromolar range, and kinetic analysis revealed competitive binding at the β5c site with an estimated inhibition constant, K i, of 1.15 μM. Inhibition was also shown for the β5i site of the immunoproteasome at levels similar to those of the constitutive proteasome. Structure-activity relationship studies identified the naphthyl substituent to be crucial for activity and this was attributed to enhanced hydrophobic interactions within β5c. Further to this, halogen substitution within the naphthyl ring enhanced the activity and allowed for π-π interactions with Y169 in β5c and Y130 and F124 in β5i. The combined data highlight the importance of hydrophobic and halogen interactions in β5 binding and assist in the design of next generation inhibitors of proteasomes.

Predicting locations of cryptic pockets from single protein structures using the PocketMiner graph neural network

Cryptic pockets expand the scope of drug discovery by enabling targeting of proteins currently considered undruggable because they lack pockets in their ground state structures. However, identifying cryptic pockets is labor-intensive and slow. The ability to accurately and rapidly predict if and where cryptic pockets are likely to form from a structure would greatly accelerate the search for druggable pockets. Here, we present PocketMiner, a graph neural network trained to predict where pockets are likely to open in molecular dynamics simulations. Applying PocketMiner to single structures from a newly curated dataset of 39 experimentally confirmed cryptic pockets demonstrates that it accurately identifies cryptic pockets (ROC-AUC: 0.87) >1,000-fold faster than existing methods. We apply PocketMiner across the human proteome and show that predicted pockets open in simulations, suggesting that over half of proteins thought to lack pockets based on available structures likely contain cryptic pockets, vastly expanding the potentially druggable proteome.

Development of Novel siRNA Therapeutics: A Review with a Focus on Inclisiran for the Treatment of Hypercholesterolemia

Over the past two decades, it was discovered that introducing synthetic small interfering RNAs (siRNAs) into the cytoplasm facilitates effective gene-targeted silencing. This compromises gene expression and regulation by repressing transcription or stimulating sequence-specific RNA degradation. Substantial investments in developing RNA therapeutics for disease prevention and treatment have been made. We discuss the application to proprotein convertase subtilisin/kexin type 9 (PCSK9), which binds to and degrades the low-density lipoprotein cholesterol (LDL-C) receptor, interrupting the process of LDL-C uptake into hepatocytes. PCSK9 loss-of-function modifications show significant clinical importance by causing dominant hypocholesterolemia and lessening the risk of cardiovascular disease (CVD). Monoclonal antibodies and small interfering RNA (siRNA) drugs targeting PCSK9 are a significant new option for managing lipid disorders and improving CVD outcomes. In general, monoclonal antibodies are restricted to binding with cell surface receptors or circulating proteins. Similarly, overcoming the intracellular and extracellular defenses that prevent exogenous RNA from entering cells must be achieved for the clinical application of siRNAs. N-acetylgalactosamine (GalNAc) conjugates are a simple solution to the siRNA delivery problem that is especially suitable for treating a broad spectrum of diseases involving liver-expressed genes. Inclisiran is a GalNAc-conjugated siRNA molecule that inhibits the translation of PCSK9. The administration is only required every 3 to 6 months, which is a significant improvement over monoclonal antibodies for PCSK9. This review provides an overview of siRNA therapeutics with a focus on detailed profiles of inclisiran, mainly its delivery strategies. We discuss the mechanisms of action, its status in clinical trials, and its prospects.

Drug specificity and affinity are encoded in the probability of cryptic pocket opening in myosin motor domains

The design of compounds that can discriminate between closely related target proteins remains a central challenge in drug discovery. Specific therapeutics targeting the highly conserved myosin motor family are urgently needed as mutations in at least six of its members cause numerous diseases. Allosteric modulators, like the myosin-II inhibitor blebbistatin, are a promising means to achieve specificity. However, it remains unclear why blebbistatin inhibits myosin-II motors with different potencies given that it binds at a highly conserved pocket that is always closed in blebbistatin-free experimental structures. We hypothesized that the probability of pocket opening is an important determinant of the potency of compounds like blebbistatin. To test this hypothesis, we used Markov state models (MSMs) built from over 2 ms of aggregate molecular dynamics simulations with explicit solvent. We find that blebbistatin's binding pocket readily opens in simulations of blebbistatin-sensitive myosin isoforms. Comparing these conformational ensembles reveals that the probability of pocket opening correctly identifies which isoforms are most sensitive to blebbistatin inhibition and that docking against MSMs quantitatively predicts blebbistatin binding affinities (R²=0.82). In a blind prediction for an isoform (Myh7b) whose blebbistatin sensitivity was unknown, we find good agreement between predicted and measured IC50s (0.67 μM vs. 0.36 μM). Therefore, we expect this framework to be useful for the development of novel specific drugs across numerous protein targets.

Room-temperature serial synchrotron crystallography of the human phosphatase PTP1B

Room-temperature X-ray crystallography provides unique insights into protein conformational heterogeneity, but obtaining sufficiently large protein crystals is a common hurdle. Serial synchrotron crystallography (SSX) helps to address this hurdle by allowing the use of many medium- to small-sized crystals. Here, a recently introduced serial sample-support chip system has been used to obtain the first SSX structure of a human phosphatase, specifically protein tyrosine phosphatase 1B (PTP1B) in the unliganded (apo) state. In previous apo room-temperature structures, the active site and allosteric sites adopted alternate conformations, including open and closed conformations of the active-site WPD loop and of a distal allosteric site. By contrast, in our SSX structure the active site is best fitted with a single conformation, but the distal allosteric site is best fitted with alternate conformations. This observation argues for additional nuance in interpreting the nature of allosteric coupling in this protein. Overall, our results illustrate the promise of serial methods for room-temperature crystallography, as well as future avant-garde crystallography experiments, for PTP1B and other proteins.

Comprehensive Survey of Consensus Docking for High-Throughput Virtual Screening

The rapid advances of 3D techniques for the structural determination of proteins and the development of numerous computational methods and strategies have led to identifying highly active compounds in computer drug design. Molecular docking is a method widely used in high-throughput virtual screening campaigns to filter potential ligands targeted to proteins. A great variety of docking programs are currently available, which differ in the algorithms and approaches used to predict the binding mode and the affinity of the ligand. All programs heavily rely on scoring functions to accurately predict ligand binding affinity, and despite differences in performance, none of these docking programs is preferable to the others. To overcome this problem, consensus scoring methods improve the outcome of virtual screening by averaging the rank or score of individual molecules obtained from different docking programs. The successful application of consensus docking in high-throughput virtual screening highlights the need to optimize the predictive power of molecular docking methods.

A Guide to In Silico Drug Design

The drug discovery process is a rocky path that is full of challenges, with the result that very few candidates progress from hit compound to a commercially available product, often due to factors, such as poor binding affinity, off-target effects, or physicochemical properties, such as solubility or stability. This process is further complicated by high research and development costs and time requirements. It is thus important to optimise every step of the process in order to maximise the chances of success. As a result of the recent advancements in computer power and technology, computer-aided drug design (CADD) has become an integral part of modern drug discovery to guide and accelerate the process. In this review, we present an overview of the important CADD methods and applications, such as in silico structure prediction, refinement, modelling and target validation, that are commonly used in this area.

Computer simulation of molecular recognition in biomolecular system: from in silico screening to generalized ensembles

Prediction of ligand-receptor complex structure is important in both the basic science and the industry such as drug discovery. We report various computation molecular docking methods: fundamental in silico (virtual) screening, ensemble docking, enhanced sampling (generalized ensemble) methods, and other methods to improve the accuracy of the complex structure. We explain not only the merits of these methods but also their limits of application and discuss some interaction terms which are not considered in the in silico methods. In silico screening and ensemble docking are useful when one focuses on obtaining the native complex structure (the most thermodynamically stable complex). Generalized ensemble method provides a free-energy landscape, which shows the distribution of the most stable complex structure and semi-stable ones in a conformational space. Also, barriers separating those stable structures are identified. A researcher should select one of the methods according to the research aim and depending on complexity of the molecular system to be studied.

Data-driven analysis and druggability assessment methods to accelerate the identification of novel cancer targets

Over the past few decades, drug discovery has greatly improved the outcomes for patients, but several challenges continue to hinder the rapid development of novel drugs. Addressing unmet clinical needs requires the pursuit of drug targets that have a higher likelihood to lead to the development of successful drugs. Here we describe a bioinformatic approach for identifying novel cancer drug targets by performing statistical analysis to ascertain quantitative changes in expression levels between protein-coding genes, as well as co-expression networks to classify these genes into groups. Subsequently, we provide an overview of druggability assessment methodologies to prioritize and select the best targets to pursue.

AHoJ: rapid, tailored search and retrieval of apo and holo protein structures for user-defined ligands

SUMMARY

Understanding the mechanism of action of a protein or designing better ligands for it, often requires access to a bound (holo) and an unbound (apo) state of the protein. Resources for the quick and easy retrieval of such conformations are severely limited. Apo-Holo Juxtaposition (AHoJ), is a web application for retrieving apo-holo structure pairs for user-defined ligands. Given a query structure and one or more user-specified ligands, it retrieves all other structures of the same protein that feature the same binding site(s), aligns them, and examines the superimposed binding sites to determine whether each structure is apo or holo, in reference to the query. The resulting superimposed datasets of apo-holo pairs can be visualized and downloaded for further analysis. AHoJ accepts multiple input queries, allowing the creation of customized apo-holo datasets.

AVAILABILITY AND IMPLEMENTATION

Freely available for non-commercial use at http://apoholo.cz. Source code available at https://github.com/cusbg/AHoJ-project.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Allosteric Hotspots in the Main Protease of SARS-CoV-2

Inhibiting the main protease of SARS-CoV-2 is of great interest in tackling the COVID-19 pandemic caused by the virus. Most efforts have been centred on inhibiting the binding site of the enzyme. However, considering allosteric sites, distant from the active or orthosteric site, broadens the search space for drug candidates and confers the advantages of allosteric drug targeting. Here, we report the allosteric communication pathways in the main protease dimer by using two novel fully atomistic graph-theoretical methods: Bond-to-bond propensity, which has been previously successful in identifying allosteric sites in extensive benchmark data sets without a priori knowledge, and Markov transient analysis, which has previously aided in finding novel drug targets in catalytic protein families. Using statistical bootstrapping, we score the highest ranking sites against random sites at similar distances, and we identify four statistically significant putative allosteric sites as good candidates for alternative drug targeting.

FPocketWeb: protein pocket hunting in a web browser

Detecting macromolecular (e.g., protein) cavities where small molecules bind is an early step in computer-aided drug discovery. Multiple pocket-detection algorithms have been developed over the past several decades. Among them, fpocket, created by Schmidtke and Le Guilloux, is particularly popular. Like many programs used in computational-biology research, fpocket requires users to download and install an executable file. That file must also be run via a command-line interface, further complicating use. An existing fpocket server application effectively addresses these challenges, but it requires users to upload their possibly proprietary structures to a third-party server. The FPocketWeb web app builds on this prior work. It runs the fpocket3 executable entirely in a web browser without requiring installation. The pocket-finding calculations occur on the user's computer rather than on a remote server. A working version of the open-source FPocketWeb app can be accessed free of charge from http://durrantlab.com/fpocketweb .

Mapping the binding sites of challenging drug targets

An increasing number of medically important proteins are challenging drug targets because their binding sites are too shallow or too polar, are cryptic and thus not detectable without a bound ligand or located in a protein-protein interface. While such proteins may not bind druglike small molecules with sufficiently high affinity, they are frequently druggable using novel therapeutic modalities. The need for such modalities can be determined by experimental or computational fragment based methods. Computational mapping by mixed solvent molecular dynamics simulations or the FTMap server can be used to determine binding hot spots. The strength and location of the hot spots provide very useful information for selecting potentially successful approaches to drug discovery.

FTMove: A Web Server for Detection and Analysis of Cryptic and Allosteric Binding Sites by Mapping Multiple Protein Structures

Protein mapping distributes many copies of different molecular probes on the surface of a target protein in order to determine binding hot spots, regions that are highly preferable for ligand binding. While mapping of X-ray structures by the FTMap server is inherently static, this limitation can be overcome by the simultaneous analysis of multiple structures of the protein. FTMove is an automated web server that implements this approach. From the input of a target protein, by PDB code, the server identifies all structures of the protein available in the PDB, runs mapping on them, and combines the results to form binding hot spots and binding sites. The user may also upload their own protein structures, bypassing the PDB search for similar structures. Output of the server consists of the consensus binding sites and the individual mapping results for each structure - including the number of probes located in each binding site, for each structure. This level of detail allows the users to investigate how the strength of a binding site relates to the protein conformation, other binding sites, and the presence of ligands or mutations. In addition, the structures are clustered on the basis of their binding properties. The use of FTMove is demonstrated by application to 22 proteins with known allosteric binding sites; the orthosteric and allosteric binding sites were identified in all but one case, and the sites were typically ranked among the top five. The FTMove server is publicly available at https://ftmove.bu.edu.

Coarse-graining protein structures into their dynamic communities with DCI, a dynamic community identifier

SUMMARY

A new dynamic community identifier (DCI) is presented that relies upon protein residue dynamic cross-correlations generated by Gaussian elastic network models to identify those residue clusters exhibiting motions within a protein. A number of examples of communities are shown for diverse proteins, including GPCRs. It is a tool that can immediately simplify and clarify the most essential functional moving parts of any given protein. Proteins usually can be subdivided into groups of residues that move as communities. These are usually densely packed local sub-structures, but in some cases can be physically distant residues identified to be within the same community. The set of these communities for each protein are the moving parts. The ways in which these are organized overall can aid in understanding many aspects of functional dynamics and allostery. DCI enables a more direct understanding of functions including enzyme activity, action across membranes and changes in the community structure from mutations or ligand binding. The DCI server is freely available on a web site (https://dci.bb.iastate.edu/).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A cryptic pocket in Ebola VP35 allosterically controls RNA binding

Protein-protein and protein-nucleic acid interactions are often considered difficult drug targets because the surfaces involved lack obvious druggable pockets. Cryptic pockets could present opportunities for targeting these interactions, but identifying and exploiting these pockets remains challenging. Here, we apply a general pipeline for identifying cryptic pockets to the interferon inhibitory domain (IID) of Ebola virus viral protein 35 (VP35). VP35 plays multiple essential roles in Ebola's replication cycle but lacks pockets that present obvious utility for drug design. Using adaptive sampling simulations and machine learning algorithms, we predict VP35 harbors a cryptic pocket that is allosterically coupled to a key dsRNA-binding interface. Thiol labeling experiments corroborate the predicted pocket and mutating the predicted allosteric network supports our model of allostery. Finally, covalent modifications that mimic drug binding allosterically disrupt dsRNA binding that is essential for immune evasion. Based on these results, we expect this pipeline will be applicable to other proteins.

Glycosylation is key for enhancing drug recognition into spike glycoprotein of SARS-CoV-2

The emergence of COVID-19 caused by SARS-CoV-2 and its spread since 2019 represents the major public health problem worldwide nowadays, which has generated a high number of infections and deaths. The spike protein (S protein) is the most studied protein of SARS-CoV-2, and key to host-cell entry through ACE2 receptor. This protein presents a large pattern of glycosylations with important roles in immunity and infection mechanisms. Therefore, understanding key aspects of the molecular mechanisms of these structures, during drug recognition in SARS-CoV-2, may contribute to therapeutic alternatives. In this work, we explored the impact of glycosylations on the drug recognition on two domains of the S protein, the receptor-binding domain (RBD) and the N-terminal domain (NTD) through molecular dynamics simulations and computational biophysics analysis. Our results show that glycosylations in the S protein induce structural stability and changes in rigidity/flexibility related to the number of glycosylations in the structure. These structural changes are important for its biological activity as well as the correct interaction of ligands in the RBD and NTD regions. Additionally, we evidenced a roto-translation phenomenon in the interaction of the ligand with RBD in the absence of glycosylation, which disappears due to the influence of glycosylation and the convergence of metastable states in RBM. Similarly, glycosylations in NTD promote an induced fit phenomenon, which is not observed in the absence of glycosylations; this process is decisive for the activity of the ligand at the cryptic site. Altogether, these results provide an explanation of glycosylation relevance in biophysical properties and drug recognition to S protein of SARS-CoV-2, which must be considered in the rational drug development and virtual screening targeting S protein.

Accelerated Ligand-Mapping Molecular Dynamics Simulations for the Detection of Recalcitrant Cryptic Pockets and Occluded Binding Sites

The identification and characterization of binding sites is a critical component of structure-based drug design (SBDD). Probe-based/cosolvent molecular dynamics (MD) methods that allow for protein flexibility have been developed to predict ligand binding sites. However, cryptic pockets that appear only upon ligand binding and occluded binding sites with no access to the solvent pose significant challenges to these methods. Here, we report the development of accelerated ligand-mapping MD (aLMMD), which combines accelerated MD with LMMD, for the detection of these challenging binding sites. The method was validated on five proteins with what we term "recalcitrant" cryptic pockets, which are deeply buried pockets that require extensive movement of the protein backbone to expose, and three proteins with occluded binding sites. In all the cases, aLMMD was able to detect and sample the binding sites. Our results suggest that aLMMD could be used as a general approach for the detection of such elusive binding sites in protein targets, thus providing valuable information for SBDD.

Small-molecule modulators of TRMT2A decrease PolyQ aggregation and PolyQ-induced cell death

Polyglutamine (polyQ) diseases are characterized by an expansion of cytosine-adenine-guanine (CAG) trinucleotide repeats encoding for an uninterrupted prolonged polyQ tract. We previously identified TRMT2A as a strong modifier of polyQ-induced toxicity in an unbiased large-scale screen in Drosophila melanogaster. This work aimed at identifying and validating pharmacological TRMT2A inhibitors as treatment opportunities for polyQ diseases in humans. Computer-aided drug discovery was implemented to identify human TRMT2A inhibitors. Additionally, the crystal structure of one protein domain, the RNA recognition motif (RRM), was determined, and Biacore experiments with the RRM were performed. The identified molecules were validated for their potency to reduce polyQ aggregation and polyQ-induced cell death in human HEK293T cells and patient derived fibroblasts. Our work provides a first step towards pharmacological inhibition of this enzyme and indicates TRMT2A as a viable drug target for polyQ diseases.

What Features of Ligands Are Relevant to the Opening of Cryptic Pockets in Drug Targets?

Small-molecule drug design aims to identify inhibitors that can specifically bind to a functionally important region on the target, i.e., an active site of an enzyme. Identification of potential binding pockets is typically based on static three-dimensional structures. However, small molecules may induce and select a dynamic binding pocket that is not visible in the apo protein, which presents a well-recognized challenge for structure-based drug discovery. Here, we assessed whether it is possible to identify features in molecules, which we refer to as inducers, that can induce the opening of cryptic pockets. The volume change between apo and bound protein conformations was used as a metric to differentiate chemical features in inducers vs. non-inducers. Based on the dataset of holo–apo pairs, classification models were built to determine an optimum threshold. The model analysis suggested that inducers preferred to be more hydrophobic and aromatic. The impact of sulfur was ambiguous, while phosphorus and halogen atoms were overrepresented in inducers. The fragment analysis showed that small changes in the structures of molecules can strongly affect the potential to induce a cryptic pocket. This analysis and developed model can be used to design inducers that can potentially open cryptic pockets for undruggable proteins.

Big data and artificial intelligence (AI) methodologies for computer-aided drug design (CADD)

There have been numerous advances in the development of computational and statistical methods and applications of big data and artificial intelligence (AI) techniques for computer-aided drug design (CADD). Drug design is a costly and laborious process considering the biological complexity of diseases. To effectively and efficiently design and develop a new drug, CADD can be used to apply cutting-edge techniques to various limitations in the drug design field. Data pre-processing approaches, which clean the raw data for consistent and reproducible applications of big data and AI methods are introduced. We include the current status of the applicability of big data and AI methods to drug design areas such as the identification of binding sites in target proteins, structure-based virtual screening (SBVS), and absorption, distribution, metabolism, excretion and toxicity (ADMET) property prediction. Data pre-processing and applications of big data and AI methods enable the accurate and comprehensive analysis of massive biomedical data and the development of predictive models in the field of drug design. Understanding and analyzing biological, chemical, or pharmaceutical architectures of biomedical entities related to drug design will provide beneficial information in the biomedical big data era.