Investigating Cryptic Binding Sites by Molecular Dynamics Simulations

Dynamics-based drug discovery by time-resolved cryo-EM

Rational structure-based drug design (SBDD) depends on high-resolution structural models of target macromolecules or their complexes. However, the lack of atomic-level functional molecular dynamics hinders the applications of SBDD and limits their effective translation into clinically successful therapeutics. Time-resolved cryo-electron microscopy (cryo-EM) has emerged as a powerful tool in structural biology, capable of capturing high-resolution snapshots of biomolecular machines in action. Unlike molecular dynamics (MD) simulations, time-resolved cryo-EM can visualize rare intermediate states across a broader range of timescales, providing invaluable insights into drug-binding kinetics, dynamic protein-ligand interactions, and allosteric regulation. Integration of time-resolved cryo-EM with machine learning (ML) and artificial intelligence (AI) expands SBDD into a dynamics-based approach, allowing for more accurate pharmacological modeling of challenging drug targets that are beyond the reach of MD simulations. Time-resolved cryo-EM can help researchers to identify novel druggable conformations, overcome drug resistance, and reduce the time and cost of clinical translations. Despite current challenges, the future development of time-resolved cryo-EM with AI and in situ imaging strategy, such as cryo-electron tomography, holds the potential to revolutionize drug discovery by revealing in vivo molecular dynamics of drug actions at an unprecedented spatiotemporal scale.

Enhanced identification of small molecules binding to hnRNPA1 via cryptic pockets mapping coupled with X-Ray fragment screening

The human heterogeneous nuclear ribonucleoprotein (hnRNP) A1 is a prototypical RNA-binding protein essential for regulating a wide range of post-transcriptional events in cells. As a multifunctional protein with a key role in RNA metabolism, deregulation of its functions has been linked to neurodegenerative diseases, tumour aggressiveness and chemoresistance, which has fuelled efforts to develop novel therapeutics that modulate its RNA binding activities. Here, using a combination of Molecular Dynamics (MD) simulations and graph neural network pocket predictions, we showed that hnRNPA1 N-terminal RNA binding domain (UP1) contains several cryptic pockets capable of binding small molecules. To identify chemical entities for the development of potent drug candidates and experimentally validate identified druggable hotspots, we carried out a large fragment screening on UP1 protein crystals. Our screen identified 36 hits that extensively sample UP1 functional regions involved in RNA recognition and binding, as well as map hotspots onto novel protein interaction surfaces. We observed a wide range of ligand-induced conformational variation, by stabilisation of dynamic protein regions. Our high-resolution structures, the first of an hnRNP in complex with a fragment or small molecule, provide rapid routes for the rational development of a range of different inhibitors and chemical tools for studying molecular mechanisms of hnRNPA1-mediated splicing regulation.

New strategies to enhance the efficiency and precision of drug discovery

Drug discovery plays a crucial role in medicinal chemistry, serving as the cornerstone for developing new treatments to address a wide range of diseases. This review emphasizes the significance of advanced strategies, such as Click Chemistry, Targeted Protein Degradation (TPD), DNA-Encoded Libraries (DELs), and Computer-Aided Drug Design (CADD), in boosting the drug discovery process. Click Chemistry streamlines the synthesis of diverse compound libraries, facilitating efficient hit discovery and lead optimization. TPD harnesses natural degradation pathways to target previously undruggable proteins, while DELs enable high-throughput screening of millions of compounds. CADD employs computational methods to refine candidate selection and reduce resource expenditure. To demonstrate the utility of these methodologies, we highlight exemplary small molecules discovered in the past decade, along with a summary of marketed drugs and investigational new drugs that exemplify their clinical impact. These examples illustrate how these techniques directly contribute to advancing medicinal chemistry from the bench to bedside. Looking ahead, Artificial Intelligence (AI) technologies and interdisciplinary collaboration are poised to address the growing complexity of drug discovery. By fostering a deeper understanding of these transformative strategies, this review aims to inspire innovative research directions and further advance the field of medicinal chemistry.

Computational advances in discovering cryptic pockets for drug discovery

A number of promising therapeutic target proteins have been considered "undruggable" due to the lack of well-defined ligandable pockets. Substantial research in protein dynamics has elucidated the existence of "cryptic" pockets that only exist transiently and become favorable for binding in the presence of a ligand. These pockets provide an avenue to target challenging proteins, inspiring the development of multiple computational methods. This review highlights established cryptic pocket modeling approaches like mixed solvent molecular dynamics and recent applications of enhanced sampling and AI-based methods in therapeutically relevant proteins.

Promising tools for future drug discovery and development in antiarrhythmic therapy

Arrhythmia refers to irregularities in the rate and rhythm of the heart, with symptoms spanning from mild palpitations to life-threatening arrhythmias and sudden cardiac death. The complex molecular nature of arrhythmias complicates the selection of appropriate treatment. Current therapies involve the use of antiarrhythmic drugs (class I-IV) with limited efficacy and dangerous side effects and implantable pacemakers and cardioverter-defibrillators with hardware-related complications and inappropriate shocks. The number of novel antiarrhythmic drugs in the development pipeline has decreased substantially during the last decade and underscores uncertainties regarding future developments in this field. Consequently, arrhythmia treatment poses significant challenges, prompting the need for alternative approaches. Remarkably, innovative drug discovery and development technologies show promise in helping advance antiarrhythmic therapies. In this article, we review unique characteristics and the transformative potential of emerging technologies that offer unprecedented opportunities for transitioning from traditional antiarrhythmics to next-generation therapies. We assess stem cell technology, emphasizing the utility of innovative cell profiling using multiomics, high-throughput screening, and advanced computational modeling in developing treatments tailored precisely to individual genetic and physiological profiles. We offer insights into gene therapy, peptide, and peptibody approaches for drug delivery. We finally discuss potential strengths and weaknesses of such techniques in reducing adverse effects and enhancing overall treatment outcomes, leading to more effective, specific, and safer therapies. Altogether, this comprehensive overview introduces innovative avenues for personalized rhythm therapy, with particular emphasis on drug discovery, aiming to advance the arrhythmia treatment landscape and the prevention of sudden cardiac death. SIGNIFICANCE STATEMENT: Arrhythmias and sudden cardiac death account for 15%-20% of deaths worldwide. However, current antiarrhythmic therapies are ineffective and have dangerous side effects. Here, we review the field of arrhythmia treatment underscoring the slow progress in advancing the cardiac rhythm therapy pipeline and the uncertainties regarding evolution of this field. We provide information on how emerging technological and experimental tools can help accelerate progress and address the limitations of antiarrhythmic drug discovery.

CryptoBench: cryptic protein-ligand binding sites dataset and benchmark

MOTIVATION

Structure-based methods for detecting protein-ligand binding sites play a crucial role in various domains, from fundamental research to biomedical applications. However, current prediction methodologies often rely on holo (ligand-bound) protein conformations for training and evaluation, overlooking the significance of the apo (ligand-free) states. This oversight is particularly problematic in the case of cryptic binding sites (CBSs) where holo-based assessment yields unrealistic performance expectations.

RESULTS

To advance the development in this domain, we introduce CryptoBench, a benchmark dataset tailored for training and evaluating novel CBS prediction methodologies. CryptoBench is constructed upon a large collection of apo-holo protein pairs, grouped by UniProtID, clustered by sequence identity, and filtered to contain only structures with substantial structural change in the binding site. CryptoBench comprises 1107 structures with predefined cross-validation splits, making it the most extensive CBS dataset to date. To establish a performance baseline, we measured the predictive power of sequence- and structure-based CBS residue prediction methods using the benchmark. We selected PocketMiner as the state-of-the-art representative of the structure-based methods for CBS detection, and P2Rank, a widely-used structure-based method for general binding site prediction that is not specifically tailored for cryptic sites. For sequence-based approaches, we trained a neural network to classify binding residues using protein language model embeddings. Our sequence-based approach outperformed PocketMiner and P2Rank across key metrics, including area under the curve, area under the precision-recall curve, Matthew's correlation coefficient, and F1 scores. These results provide baseline benchmark results for future CBS and potentially also non-CBS prediction endeavors, leveraging CryptoBench as the foundational platform for further advancements in the field.

AVAILABILITY AND IMPLEMENTATION

The CryptoBench dataset, including the benchmark model, is available on Open Science Framework-https://osf.io/pz4a9/. The code and tutorial are available at the GitHub repository-https://github.com/skrhakv/CryptoBench/.

Enhanced identification of small molecules binding to hnRNPA1 via cryptic pockets mapping coupled with X-Ray fragment screening

The human heterogeneous nuclear ribonucleoprotein (hnRNP) A1 is a prototypical RNA-binding protein essential in regulating a wide range of post-transcriptional events in cells. As a multifunctional protein with a key role in RNA metabolism, deregulation of its functions has been linked to neurodegenerative diseases, tumour aggressiveness and chemoresistance, which has fuelled efforts to develop novel therapeutics that modulates its RNA binding activities. Here, using a combination of Molecular Dynamics (MD) simulations and graph neural network pockets predictions, we showed that hnRNPA1 N-terminal RNA binding domain (UP1) contains several cryptic pockets capable of binding small molecules. To identify chemical entities for development of potent drug candidates and experimentally validate identified druggable hotspots, we carried out a large fragment screening on UP1 protein crystals. Our screen identified 36 hits which extensively samples UP1 functional regions involved in RNA recognition and binding, as well as mapping hotspots onto novel protein interaction surfaces. We observed a wide range of ligand-induced conformational variation, by stabilisation of dynamic protein regions. Our high-resolution structures, the first of an hnRNP in complex with a fragment or small molecule, provides rapid routes for the rational development of a range of different inhibitors and chemical tools for studying molecular mechanisms of hnRNPA1 mediated splicing regulation.

Wetting of a Dynamically Patterned Surface Is a Time-Dependent Matter

In nature and many technological applications, aqueous solutions are in contact with patterned surfaces, which are dynamic over time scales spanning from ps to μs. For instance, in biology, exposed polar and apolar residues of biomolecules form a pattern, which fluctuates in time due to side chain and conformational motions. At metal/and oxide/water interfaces, the pattern is formed by surface topmost atoms, and fluctuations are due to, e.g., local surface polarization and rearrangements in the adsorbed water layer. All these dynamics have the potential to influence key processes such as wetting, energy relaxation, and biological function. Yet, their impact on the water H-bond network remains often elusive. Here, we leverage molecular dynamics to address this fundamental question at a self-assembled monolayer (SAM)/water interface, where ns dynamics is induced by frustrating SAM-water interactions via methylation of the terminal -OH groups of poly(ethylene glycol) (PEG) chains. We find that surface dynamics couples to the water H-bond network, inducing a response on the same ns time scale. This leads to time fluctuations of local wetting, oscillating from hydrophobic to hydrophilic environments. Our results suggest that rather than average properties, it is the local─ both in time and space─ solvation that determines the chemical-physical properties of dynamically patterned surfaces in water.

Revisiting the Plasmodium falciparum druggable genome using predicted structures and data mining

The identification of novel drug targets for the purpose of designing small molecule inhibitors is key component to modern drug discovery. In malaria parasites, discoveries of antimalarial targets have primarily occurred retroactively by investigating the mode of action of compounds found through phenotypic screens. Although this method has yielded many promising candidates, it is time- and resource-consuming and misses targets not captured by existing antimalarial compound libraries and phenotypic assay conditions. Leveraging recent advances in protein structure prediction and data mining, we systematically assessed the Plasmodium falciparum genome for proteins amenable to target-based drug discovery, identifying 867 candidate targets with evidence of small molecule binding and blood stage essentiality. Of these, 540 proteins showed strong essentiality evidence and lack inhibitors that have progressed to clinical trials. Expert review and rubric-based scoring of this subset based on additional criteria such as selectivity, structural information, and assay developability yielded 67 high priority candidates. This study also provides a genome-wide data resource and implements a generalizable framework for systematically evaluating and prioritizing novel pathogenic disease targets.

Artificial Intelligence-Driven Computational Approaches in the Development of Anticancer Drugs

The integration of AI has revolutionized cancer drug development, transforming the landscape of drug discovery through sophisticated computational techniques. AI-powered models and algorithms have enhanced computer-aided drug design (CADD), offering unprecedented precision in identifying potential anticancer compounds. Traditionally, cancer drug design has been a complex, resource-intensive process, but AI introduces new opportunities to accelerate discovery, reduce costs, and optimize efficiency. This manuscript delves into the transformative applications of AI-driven methodologies in predicting and developing anticancer drugs, critically evaluating their potential to reshape the future of cancer therapeutics while addressing their challenges and limitations.

Exploration of Cryptic Pockets Using Enhanced Sampling Along Normal Modes: A Case Study of KRAS G12D

Identification of cryptic pockets has the potential to open new therapeutic opportunities by discovering ligand binding sites that remain hidden in static apo structures of a target protein. Moreover, allosteric cryptic pockets can become valuable for designing target-selective ligands when the natural ligand binding sites are conserved in variants of a protein. For example, before an allosteric cryptic pocket was discovered, KRAS was considered undruggable due to its smooth surface and conservation of the GDP/GTP binding pocket across the wild type and oncogenic isoforms. Recent identification of the Switch-II cryptic pocket in the KRASG12C mutant and FDA approval of anticancer drugs targeting this site underscores the importance of cryptic pockets in solving pharmaceutical challenges. Here, we present a newly developed approach for the exploration of cryptic pockets using weighted ensemble molecular dynamics simulations with inherent normal modes as progress coordinates applied to the wild type KRAS and the G12D mutant. We performed extensive all-atomic simulations (>400 μs) with and without several cosolvents (xenon, ethanol, benzene), and analyzed trajectories using three distinct methods to search for potential binding pockets. These methods have been applied as a proof-of-concept to KRAS and have shown they can predict known cryptic binding sites. Furthermore, we performed ligand-binding simulations of a known inhibitor (MRTX1133) to shed light on the nature of cryptic pockets in KRASG12D and the role of conformational selection vs induced-fit mechanism in the formation of these cryptic pockets.

CosolvKit: a Versatile Tool for Cosolvent MD Preparation and Analysis

Cosolvent molecular dynamics (MDs) are an increasingly popular form of simulations where small molecule cosolvents are added to water-solvated protein systems. These simulations can perform diverse target characterization tasks, including cryptic and allosteric pocket identification and pharmacophore profiling and supplement suites of enhanced sampling methods to explore protein conformational landscapes. The behavior of these systems is tied to the cosolvents used, so the ability to define diverse and complex mixtures is critical in dictating the outcome of the simulations. However, existing methods for preparing cosolvent simulations only support a limited number of predefined cosolvents and concentrations. Here, we present CosolvKit, a tool for the preparation and analysis of systems composed of user-defined cosolvents and concentrations. This tool is modular, supporting the creation of files for multiple MD engines, as well as direct access to OpenMM simulations, and offering access to a variety of generalizable small-molecule force fields. To the best of our knowledge, CosolvKit represents the first generalized approach for the construction of these simulations.

Unbiased MD simulations identify lipid binding sites in lipid transfer proteins

The characterization of lipid binding to lipid transfer proteins (LTPs) is fundamental to understand their molecular mechanism. However, several structures of LTPs, and notably those proposed to act as bridges between membranes, do not provide the precise location of their endogenous lipid ligands. To address this limitation, computational approaches are a powerful alternative methodology, but they are often limited by the high flexibility of lipid substrates. Here, we develop a protocol based on unbiased coarse-grain molecular dynamics simulations in which lipids placed away from the protein can spontaneously bind to LTPs. This approach accurately determines binding pockets in LTPs and provides a working hypothesis for the lipid entry pathway. We apply this approach to characterize lipid binding to bridge LTPs of the Vps13-Atg2 family, for which the lipid localization inside the protein is currently unknown. Overall, our work paves the way to determine binding pockets and entry pathways for several LTPs in an inexpensive, fast, and accurate manner.

Integrative residue-intuitive machine learning and MD Approach to Unveil Allosteric Site and Mechanism for β2AR

Allosteric drugs offer a new avenue for modern drug design. However, the identification of cryptic allosteric sites presents a formidable challenge. Following the allostery nature of residue-driven conformation transition, we propose a state-of-the-art computational pipeline by developing a residue-intuitive hybrid machine learning (RHML) model coupled with molecular dynamics (MD) simulation, through which we can efficiently identify the allosteric site and allosteric modulator as well as reveal their regulation mechanism. For the clinical target β2-adrenoceptor (β2AR), we discover an additional allosteric site located around residues D792.50, F2826.44, N3187.45 and S3197.46 and one putative allosteric modulator ZINC5042. Using Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) and protein structure network (PSN), the allosteric potency and regulation mechanism are probed to further improve identification accuracy. Benefiting from sufficient computational evidence, the experimental assays then validate our predicted allosteric site, negative allosteric potency and regulation pathway, showcasing the effectiveness of the identification pipeline in practice. We expect that it will be applicable to other target proteins.

Insights into bacterial interactions: Comparing fluorine-containing 1,2,4-triazoles to antibiotics using molecular docking and molecular dynamics approaches

Understanding the interactions between drugs and enzymes is crucial for designing effective therapeutics. This study employed a combination of molecular docking and molecular dynamics (MD) simulations to evaluate the binding affinity, stability, and dynamic behavior of two new compounds (compound 1 and compound 2) compared to vancomycin and meropenem against Staphylococcus aureus and Serratia marcescens bacterial enzymes. Molecular docking studies provided insights into the binding interactions and affinities of these compounds, revealing that both compound 1 and compound 2 exhibit promising binding profiles. In particular, compound 1 demonstrated lower binding energies with key enzymes from Staphylococcus aureus compared to vancomycin, suggesting enhanced potential. MD simulations further elucidated the dynamic stability of these complexes. Results indicated that compound 1 maintains consistent binding modes with low RMSD and RMSF values, implying stable interactions. In contrast, vancomycin exhibited high RMSD and RMSF values in some enzyme complexes, reflecting potential instability. Compound 2 showed competitive stability and binding behavior compared to meropenem, with comparable RMSD and RMSF values across various enzyme targets. These findings highlight the potential of compound 1 and compound 2 as viable candidates for further development, offering insights into their stability and efficacy as new therapeutic agents.

AlphaFold-SFA: Accelerated sampling of cryptic pocket opening, protein-ligand binding and allostery by AlphaFold, slow feature analysis and metadynamics

Sampling rare events in proteins is crucial for comprehending complex phenomena like cryptic pocket opening, where transient structural changes expose new binding sites. Understanding these rare events also sheds light on protein-ligand binding and allosteric communications, where distant site interactions influence protein function. Traditional unbiased molecular dynamics simulations often fail to sample such rare events, as the free energy barrier between metastable states is large relative to the thermal energy. This renders these events inaccessible on the timescales typically simulated by unbiased molecular dynamics, limiting our understanding of these critical processes. In this paper, we proposed a novel unsupervised learning approach termed as slow feature analysis (SFA) which aims to extract slowly varying features from high-dimensional temporal data. SFA trained on small unbiased molecular dynamics simulations launched from AlphaFold generated conformational ensembles manages to capture rare events governing cryptic pocket opening, protein-ligand binding, and allosteric communications in a kinase. Metadynamics simulations using SFA as collective variables manage to sample 'deep' cryptic pocket opening within a few hundreds of nanoseconds which was beyond the reach of microsecond long unbiased molecular dynamics simulations. SFA augmented metadynamics also managed to capture conformational plasticity of protein upon ligand binding/unbinding and provided novel insights into allosteric communication in receptor-interacting protein kinase 2 (RIPK2) which dictates protein-protein interaction. Taken together, our results show how SFA acts as a dimensionality reduction tool which bridges the gap between AlphaFold, molecular dynamics simulation and metadynamics in context of capturing rare events in biomolecules, extending the scope of structure-based drug discovery in the era of AlphaFold.

AlphaFold and Protein Folding: Not Dead Yet! The Frontier Is Conformational Ensembles

Like the black knight in the classic Monty Python movie, grand scientific challenges such as protein folding are hard to finish off. Notably, AlphaFold is revolutionizing structural biology by bringing highly accurate structure prediction to the masses and opening up innumerable new avenues of research. Despite this enormous success, calling structure prediction, much less protein folding and related problems, "solved" is dangerous, as doing so could stymie further progress. Imagine what the world would be like if we had declared flight solved after the first commercial airlines opened and stopped investing in further research and development. Likewise, there are still important limitations to structure prediction that we would benefit from addressing. Moreover, we are limited in our understanding of the enormous diversity of different structures a single protein can adopt (called a conformational ensemble) and the dynamics by which a protein explores this space. What is clear is that conformational ensembles are critical to protein function, and understanding this aspect of protein dynamics will advance our ability to design new proteins and drugs.

Relevance of G protein-coupled receptor (GPCR) dynamics for receptor activation, signalling bias and allosteric modulation

G protein-coupled receptors (GPCRs) are one of the major drug targets. In recent years, computational drug design for GPCRs has mainly focused on static structures obtained through X-ray crystallography, cryogenic electron microscopy (cryo-EM) or in silico modelling as a starting point for virtual screening campaigns. However, GPCRs are highly flexible entities with the ability to adopt different conformational states that elicit different physiological responses. Including this knowledge in the drug discovery pipeline can help to tailor novel conformation-specific drugs with an improved therapeutic profile. In this review, we outline our current knowledge about GPCR dynamics that is relevant for receptor activation, signalling bias and allosteric modulation. Ultimately, we highlight new technological implementations such as time-resolved X-ray crystallography and cryo-EM as well as computational algorithms that can contribute to a more comprehensive understanding of receptor dynamics and its relevance for GPCR functionality.

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.

The Secret Ballet Inside Multivesicular Bodies

Lipid bilayers possess the capacity for self-assembly due to the amphipathic nature of lipid molecules, which have both hydrophobic and hydrophilic regions. When confined, lipid bilayers exhibit astonishing versatility in their forms, adopting diverse shapes that are challenging to observe through experimental means. Exploiting this adaptability, lipid structures motivate the development of bio-inspired mechanomaterials and integrated nanobio-interfaces that could seamlessly merge with biological entities, ultimately bridging the gap between synthetic and biological systems. In this work, we demonstrate how, in numerical simulations of multivesicular bodies, a fascinating evolution unfolds from an initial semblance of order toward states of higher entropy over time. We observe dynamic rearrangements in confined vesicles that reveal unexpected limit shapes of distinct geometric patterns. We identify five structures as the basic building blocks that systematically repeat under various conditions of size and composition. Moreover, we observe more complex and less frequent shapes that emerge in confined spaces. Our results provide insights into the dynamics of multivesicular systems, offering a richer understanding of how confined lipid bodies spontaneously self-organize.

Similarities and Differences in Ligand Binding to Protein and RNA Targets: The Case of Riboflavin

It is nowadays clear that RNA molecules can play active roles in several biological processes. As a result, an increasing number of RNAs are gradually being identified as potentially druggable targets. In particular, noncoding RNAs can adopt highly organized conformations that are suitable for drug binding. However, RNAs are still considered challenging targets due to their complex structural dynamics and high charge density. Thus, elucidating relevant features of drug-RNA binding is fundamental for advancing drug discovery. Here, by using Molecular Dynamics simulations, we compare key features of ligand binding to proteins with those observed in RNA. Specifically, we explore similarities and differences in terms of (i) conformational flexibility of the target, (ii) electrostatic contribution to binding free energy, and (iii) water and ligand dynamics. As a test case, we examine binding of the same ligand, namely riboflavin, to protein and RNA targets, specifically the riboflavin (RF) kinase and flavin mononucleotide (FMN) riboswitch. The FMN riboswitch exhibited enhanced fluctuations and explored a wider conformational space, compared to the protein target, underscoring the importance of RNA flexibility in ligand binding. Conversely, a similar electrostatic contribution to the binding free energy of riboflavin was found. Finally, greater stability of water molecules was observed in the FMN riboswitch compared to the RF kinase, possibly due to the different shape and polarity of the pockets.

Probing allosteric communication with combined molecular dynamics simulations and network analysis

Understanding the allosteric mechanisms within biomolecules involved in diseases is of paramount importance for drug discovery. Indeed, characterizing communication pathways and critical hotspots in signal transduction can guide a rational approach to leverage allosteric modulation for therapeutic purposes. While the atomistic signatures of allosteric processes are difficult to determine experimentally, computational methods can be a remarkable resource. Network analysis built on Molecular Dynamics simulation data is particularly suited in this respect and is gradually becoming of routine use. Herein, we collect the recent literature in the field, discussing different aspects and available options for network construction and analysis. We further highlight interesting refinements and extensions, eventually providing our perspective on this topic.

Ligandability assessment of the C-terminal Rel-homology domain of NFAT1

Transcription factors are generally considered challenging, if not "undruggable", targets but they promise new therapeutic options due to their fundamental involvement in many diseases. In this study, we aim to assess the ligandability of the C-terminal Rel-homology domain of nuclear factor of activated T cells 1 (NFAT1), a TF implicated in T-cell regulation. Using a combination of experimental and computational approaches, we demonstrate that small molecule fragments can indeed bind to this protein domain. The newly identified binder is the first small molecule binder to NFAT1 validated with biophysical methods and an elucidated binding mode by X-ray crystallography. The reported eutomer/distomer pair provides a strong basis for potential exploration of higher potency binders on the path toward degrader or glue modalities.

Leveraging Bidirectional Nature of Allostery To Inhibit Protein-Protein Interactions (PPIs): A Case Study of PCSK9-LDLR Interaction

The protein PCSK9 (proprotein convertase subtilisin/Kexin type 9) negatively regulates the recycling of LDLR (low-density lipoprotein receptor), leading to an elevated plasma level of LDL. Inhibition of PCSK9-LDLR interaction has emerged as a promising therapeutic strategy to manage hypercholesterolemia. However, the large interaction surface area between PCSK9 and LDLR makes it challenging to identify a small molecule competitive inhibitor. An alternative strategy would be to identify distal cryptic sites as targets for allosteric inhibitors that can remotely modulate PCSK9-LDLR interaction. Using several microseconds long molecular dynamics (MD) simulations, we demonstrate that on binding with LDLR, there is a significant conformational change (population shift) in a distal loop (residues 211-222) region of PCSK9. Consistent with the bidirectional nature of allostery, we establish a clear correlation between the loop conformation and the binding affinity with LDLR. Using a thermodynamic argument, we establish that the loop conformations predominantly present in the apo state of PCSK9 would have lower LDLR binding affinity, and they would be potential targets for designing allosteric inhibitors. We elucidate the molecular origin of the allosteric coupling between this loop and the LDLR binding interface in terms of the population shift in a set of salt bridges and hydrogen bonds. Overall, our work provides a general strategy toward identifying allosteric hotspots: compare the conformational ensemble of the receptor between the apo and bound states of the protein and identify distal conformational changes, if any. The inhibitors should be designed to bind and stabilize the apo-specific conformations.

Some mechanistic underpinnings of molecular adaptations of SARS-COV-2 spike protein by integrating candidate adaptive polymorphisms with protein dynamics

We integrate evolutionary predictions based on the neutral theory of molecular evolution with protein dynamics to generate mechanistic insight into the molecular adaptations of the SARS-COV-2 spike (S) protein. With this approach, we first identified candidate adaptive polymorphisms (CAPs) of the SARS-CoV-2 S protein and assessed the impact of these CAPs through dynamics analysis. Not only have we found that CAPs frequently overlap with well-known functional sites, but also, using several different dynamics-based metrics, we reveal the critical allosteric interplay between SARS-CoV-2 CAPs and the S protein binding sites with the human ACE2 (hACE2) protein. CAPs interact far differently with the hACE2 binding site residues in the open conformation of the S protein compared to the closed form. In particular, the CAP sites control the dynamics of binding residues in the open state, suggesting an allosteric control of hACE2 binding. We also explored the characteristic mutations of different SARS-CoV-2 strains to find dynamic hallmarks and potential effects of future mutations. Our analyses reveal that Delta strain-specific variants have non-additive (i.e., epistatic) interactions with CAP sites, whereas the less pathogenic Omicron strains have mostly additive mutations. Finally, our dynamics-based analysis suggests that the novel mutations observed in the Omicron strain epistatically interact with the CAP sites to help escape antibody binding.

Unveiling the Mechanisms of EGCG-p53 Interactions through Molecular Dynamics Simulations

Green tea consumption is associated with protective and preventive effects against various types of cancer. These effects are attributed to polyphenols, particularly epigallocatechin-3-gallate (EGCG). EGCG acts by directly inhibiting tumor suppressor protein p53. The binding mechanism by which EGCG inhibits p53 activity is associated with residues Trp23-Lys24 and Pro47-Thr55 within the p53 N-terminal domain (NTD). However, the structural and thermodynamic aspects of the interaction between EGCG and p53 are poorly understood. Therefore, based on crystallographic data, we combine docking, molecular dynamics (MD) simulations, and molecular mechanics generalized Born surface area approaches to explore the intricacies of the EGCG-p53 binding mechanism. A triplicate microsecond MD simulation for each system is initially performed to capture diverse p53 NTD conformations. From the start, the most populated cluster of the second run (R2-1) stands out due to a unique opening between Trp23 and Trp53. During MD simulations, this conformation allows EGCG to sustain a high level of stability and affinity while interacting with both regions of interest and deepening the binding pocket. Structural analysis emphasizes the significance of pyrogallol motifs in EGCG binding. Therefore, the conformational shift in this gap is pivotal, enabling EGCG to impede p53 interactions and manifest its anticancer properties. These findings enhance the present comprehension of the anticancer properties of green tea polyphenols and pave the way for future therapeutic developments.

Data-oriented protein kinase drug discovery

The continued growth of data from biological screening and medicinal chemistry provides opportunities for data-driven experimental design and decision making in early-phase drug discovery. Approaches adopted from data science help to integrate internal and public domain data and extract knowledge from historical in-house data. Protein kinase (PK) drug discovery is an exemplary area where large amounts of data are accumulating, providing a valuable knowledge base for discovery projects. Herein, the evolution of PK drug discovery and development of small molecular PK inhibitors (PKIs) is reviewed, highlighting milestone developments in the field and discussing exemplary studies providing a basis for increasing data orientation of PK discovery efforts.

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.

Graph Attention Site Prediction (GrASP): Identifying Druggable Binding Sites Using Graph Neural Networks with Attention

Identifying and discovering druggable protein binding sites is an important early step in computer-aided drug discovery, but it remains a difficult task where most campaigns rely on a priori knowledge of binding sites from experiments. Here, we present a binding site prediction method called Graph Attention Site Prediction (GrASP) and re-evaluate assumptions in nearly every step in the site prediction workflow from data set preparation to model evaluation. GrASP is able to achieve state-of-the-art performance at recovering binding sites in PDB structures while maintaining a high degree of precision which will minimize wasted computation in downstream tasks such as docking and free energy perturbation.

Drugging the entire human proteome: Are we there yet?

Each of the ∼20,000 proteins in the human proteome is a potential target for compounds that bind to it and modify its function. The 3D structures of most of these proteins are now available. Here, we discuss the prospects for using these structures to perform proteome-wide virtual HTS (VHTS). We compare physics-based (docking) and AI VHTS approaches, some of which are now being applied with large databases of compounds to thousands of targets. Although preliminary proteome-wide screens are now within our grasp, further methodological developments are expected to improve the accuracy of the results.

Exploring Lead-Like Molecules of Traditional Chinese Medicine for Treatment Quest against Aliarcobacter butzleri: In Silico Toxicity Assessment, Dynamics Simulation, and Pharmacokinetic Profiling

Background

Aliarcobacter butzleri is a Gram-negative, curved or spiral-shaped, microaerophilic bacterium and causes human infections, specifically diarrhea, fever, and sepsis. The research objective of this study was to employ computer-aided drug design techniques to identify potential natural product inhibitors of a vital enzyme in this bacterium. The pyrimidine biosynthesis pathway in its core genome fraction is crucial for its survival and presents a potential target for novel therapeutics. Hence, novel small molecule inhibitors were identified (from traditional Chinese medicinal (TCM) compound library) against it, which may be used for possible curbing of infection by A. butzleri. Methods. A comprehensive subtractive genomics approach was utilized to identify a key enzyme (orotidine-5'-phosphate decarboxylase) cluster conserved in the core genome fraction of A. butzleri. It was selected for inhibitor screening due to its vital role in pyrimidine biosynthesis. TCM library (n > 36,000 compounds) was screened against it using pharmacophore model based on orotidylic acid (control), and the obtained lead-like molecules were subjected to structural docking using AutoDock Vina. The top-scoring compounds, ZINC70454134, ZINC85632684, and ZINC85632721, underwent further scrutiny via a combination of physiological-based pharmacokinetics, toxicity assessment, and atomic-scale dynamics simulations (100 ns).

Results

Among the screened compounds, ZINC70454134 displayed the most favorable characteristics in terms of binding, stability, absorption, and safety parameters. Overall, traditional Chinese medicine (TCM) compounds exhibited high bioavailability, but in diseased states (cirrhosis, renal impairment, and steatosis), there was a significant decrease in absorption, Cmax, and AUC of the compounds compared to the healthy state. Furthermore, MD simulation demonstrated that the ODCase-ZINC70454134 complex had a superior overall binding affinity, supported by PCA proportion of variance and eigenvalue rank analysis. These favorable characteristics underscore its potential as a promising drug candidate.

Conclusion

The computer-aided drug design approach employed for this study helped expedite the discovery of antibacterial compounds against A. butzleri, offering a cost-effective and efficient approach to address infection by it. It is recommended that ZINC70454134 should be considered for further experimental analysis due to its indication as a potential therapeutic agent for combating A. butzleri infections. This study provides valuable insights into the molecular basis of biophysical inhibition of A. butzleri through TCM compounds.

Protein dynamics underlying allosteric regulation

Allostery is the mechanism by which information and control are propagated in biomolecules. It regulates ligand binding, chemical reactions, and conformational changes. An increasing level of experimental resolution and control over allosteric mechanisms promises a deeper understanding of the molecular basis for life and powerful new therapeutics. In this review, we survey the literature for an up-to-date biological and theoretical understanding of protein allostery. By delineating five ways in which the energy landscape or the kinetics of a system may change to give rise to allostery, we aim to help the reader grasp its physical origins. To illustrate this framework, we examine three systems that display these forms of allostery: allosteric inhibitors of beta-lactamases, thermosensation of TRP channels, and the role of kinetic allostery in the function of kinases. Finally, we summarize the growing power of computational tools available to investigate the different forms of allostery presented in this review.

Some mechanistic underpinnings of molecular adaptations of SARS-COV-2 spike protein by integrating candidate adaptive polymorphisms with protein dynamics

We integrate evolutionary predictions based on the neutral theory of molecular evolution with protein dynamics to generate mechanistic insight into the molecular adaptations of the SARS-COV-2 Spike (S) protein. With this approach, we first identified Candidate Adaptive Polymorphisms (CAPs) of the SARS-CoV-2 Spike protein and assessed the impact of these CAPs through dynamics analysis. Not only have we found that CAPs frequently overlap with well-known functional sites, but also, using several different dynamics-based metrics, we reveal the critical allosteric interplay between SARS-CoV-2 CAPs and the S protein binding sites with the human ACE2 (hACE2) protein. CAPs interact far differently with the hACE2 binding site residues in the open conformation of the S protein compared to the closed form. In particular, the CAP sites control the dynamics of binding residues in the open state, suggesting an allosteric control of hACE2 binding. We also explored the characteristic mutations of different SARS-CoV-2 strains to find dynamic hallmarks and potential effects of future mutations. Our analyses reveal that Delta strain-specific variants have non-additive (i.e., epistatic) interactions with CAP sites, whereas the less pathogenic Omicron strains have mostly additive mutations. Finally, our dynamics-based analysis suggests that the novel mutations observed in the Omicron strain epistatically interact with the CAP sites to help escape antibody binding.

Molecular Dynamics and Other HPC Simulations for Drug Discovery

High performance computing (HPC) is taking an increasingly important place in drug discovery. It makes possible the simulation of complex biochemical systems with high precision in a short time, thanks to the use of sophisticated algorithms. It promotes the advancement of knowledge in fields that are inaccessible or difficult to access through experimentation and it contributes to accelerating the discovery of drugs for unmet medical needs while reducing costs. Herein, we report how computational performance has evolved over the past years, and then we detail three domains where HPC is essential. Molecular dynamics (MD) is commonly used to explore the flexibility of proteins, thus generating a better understanding of different possible approaches to modulate their activity. Modeling and simulation of biopolymer complexes enables the study of protein-protein interactions (PPI) in healthy and disease states, thus helping the identification of targets of pharmacological interest. Virtual screening (VS) also benefits from HPC to predict in a short time, among millions or billions of virtual chemical compounds, the best potential ligands that will be tested in relevant assays to start a rational drug design process.

Computer-Aided Drug Discovery and Design: Recent Advances and Future Prospects

Computer-aided drug discovery and design involve the use of information technologies to identify and develop, on a rational ground, chemical compounds that align a set of desired physicochemical and biological properties. In its most common form, it involves the identification and/or modification of an active scaffold (or the combination of known active scaffolds), although de novo drug design from scratch is also possible. Traditionally, the drug discovery and design processes have focused on the molecular determinants of the interactions between drug candidates and their known or intended pharmacological target(s). Nevertheless, in modern times, drug discovery and design are conceived as a particularly complex multiparameter optimization task, due to the complicated, often conflicting, property requirements.This chapter provides an updated overview of in silico approaches for identifying active scaffolds and guiding the subsequent optimization process. Recent groundbreaking advances in the field have also analyzed the integration of state-of-the-art machine learning approaches in every step of the drug discovery process (from prediction of target structure to customized molecular docking scoring functions), integration of multilevel omics data, and the use of a diversity of computational approaches to assist target validation and assess plausible binding pockets.

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.

Unmasking an Allosteric Binding Site of the Papain-like Protease in SARS-CoV-2: Molecular Dynamics Simulations of Corticosteroids

To date, mechanistic insights into many clinical drugs against COVID-19 remain unexplored. Dexamethasone, a corticosteroid, is one of them. While treating the entire corticosteroid database, including vitamins D2 and D3, with cutting-edge computational techniques, several intriguing results are unfolded. From the top-notch candidates, dexamethasone is likely to inhibit the viral main protease (Mpro), with vitamin D3 exhibiting multitarget [Mpro, papain-like protease (PLpro), and nucleocapsid protein (N-pro)] roles and ciclesonide's dynamic flipping disinterring a cryptic allosteric site in the PLpro enzyme. The results rationalize why these drugs improve the health of COVID-19 patients. Understanding an enzyme's secret binding site is essential to understanding how the enzyme works and how to inhibit its function. Ciclesonide's allosteric inhibition could not only jeopardize PLpro's catalytic role in polyprotein processing but also make it less vulnerable to the host body's defense machinery. Hotspot residues in the identified allosteric site could be considered for effective therapeutic designs against PLpro.

Precision Oncology Comes of Age: Designing Best-in-Class Small Molecules by Integrating Two Decades of Advances in Chemistry, Target Biology, and Data Science

Small-molecule drugs have enabled the practice of precision oncology for genetically defined patient populations since the first approval of imatinib in 2001. Scientific and technology advances over this 20-year period have driven the evolution of cancer biology, medicinal chemistry, and data science. Collectively, these advances provide tools to more consistently design best-in-class small-molecule drugs against known, previously undruggable, and novel cancer targets. The integration of these tools and their customization in the hands of skilled drug hunters will be necessary to enable the discovery of transformational therapies for patients across a wider spectrum of cancers.

SIGNIFICANCE

Target-centric small-molecule drug discovery necessitates the consideration of multiple approaches to identify chemical matter that can be optimized into drug candidates. To do this successfully and consistently, drug hunters require a comprehensive toolbox to avoid following the "law of instrument" or Maslow's hammer concept where only one tool is applied regardless of the requirements of the task. Combining our ever-increasing understanding of cancer and cancer targets with the technological advances in drug discovery described below will accelerate the next generation of small-molecule drugs in oncology.

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.

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.

Graph Attention Site Prediction (GrASP): Identifying Druggable Binding Sites Using Graph Neural Networks with Attention

Identifying and discovering druggable protein binding sites is an important early step in computer-aided drug discovery but remains a difficult task where most campaigns rely on a priori knowledge of binding sites from experiments. Here we present a novel binding site prediction method called Graph Attention Site Prediction (GrASP) and re-evaluate assumptions in nearly every step in the site prediction workflow from dataset preparation to model evaluation. GrASP is able to achieve state-of-the-art performance at recovering binding sites in PDB structures while maintaining a high degree of precision which will minimize wasted computation in downstream tasks such as docking and free energy perturbation.

Protein dynamics provide mechanistic insights about epistasis among common missense polymorphisms

Sequencing of the protein coding genome has revealed many different missense mutations of human proteins and different population frequencies of corresponding haplotypes, which consist of different sets of those mutations. Here, we present evidence for pairwise intramolecular epistasis (i.e., nonadditive interactions) between many such mutations through an analysis of protein dynamics. We suggest that functional compensation for conserving protein dynamics is a likely evolutionary mechanism that maintains high-frequency mutations that are individually nonneutral but epistatically compensating within proteins. This analysis is the first of its type to look at human proteins with specific high population frequency mutations and examine the relationship between mutations that make up that observed high-frequency protein haplotype. Importantly, protein dynamics revealed a separation between high and low frequency haplotypes within a target protein cytochrome P450 2A7, with the high-frequency haplotypes showing behavior closer to the wild-type protein. Common protein haplotypes containing two mutations display dynamic compensation in which one mutation can correct for the dynamic effects of the other. We also utilize a dynamics-based metric, EpiScore, that evaluates the epistatic interactions and allows us to see dynamic compensation within many other proteins.

Spatial Layouts of Low-Entropy Hydration Shells Guide Protein Binding

Protein-protein binding enables orderly biological self-organization and is therefore considered a miracle of nature. Protein‒protein binding is driven by electrostatic forces, hydrogen bonding, van der Waals force, and hydrophobic interactions. Among these physical forces, only hydrophobic interactions can be considered long-range intermolecular attractions between proteins due to the electrostatic shielding of surrounding water molecules. Low-entropy hydration shells around proteins drive hydrophobic attraction among them that essentially coordinate protein‒protein binding. Here, an innovative method is developed for identifying low-entropy regions of hydration shells of proteins by screening off pseudohydrophilic groups on protein surfaces and revealing that large low-entropy regions of the hydration shells typically cover the binding sites of individual proteins. According to an analysis of determined protein complex structures, shape matching between a large low-entropy hydration shell region of a protein and that of its partner at the binding sites is revealed as a universal law. Protein‒protein binding is thus found to be mainly guided by hydrophobic collapse between the shape-matched low-entropy hydration shells that is verified by bioinformatics analyses of hundreds of structures of protein complexes, which cover four test systems. A simple algorithm is proposed to accurately predict protein binding sites.

Combining structural and coevolution information to unveil allosteric sites

Understanding allosteric regulation in biomolecules is of great interest to pharmaceutical research and computational methods emerged during the last decades to characterize allosteric coupling. However, the prediction of allosteric sites in a protein structure remains a challenging task. Here, we integrate local binding site information, coevolutionary information, and information on dynamic allostery into a structure-based three-parameter model to identify potentially hidden allosteric sites in ensembles of protein structures with orthosteric ligands. When tested on five allosteric proteins (LFA-1, p38-α, GR, MAT2A, and BCKDK), the model successfully ranked all known allosteric pockets in the top three positions. Finally, we identified a novel druggable site in MAT2A confirmed by X-ray crystallography and SPR and a hitherto unknown druggable allosteric site in BCKDK validated by biochemical and X-ray crystallography analyses. Our model can be applied in drug discovery to identify allosteric pockets.

Deconstructing allostery by computational assessment of the binding determinants of allosteric PTP1B modulators

Fragment-based drug discovery is an established methodology for finding hit molecules that can be elaborated into lead compounds. However it is currently challenging to predict whether fragment hits that do not bind to an orthosteric site could be elaborated into allosteric modulators, as in these cases binding does not necessarily translate into a functional effect. We propose a workflow using Markov State Models (MSMs) with steered molecular dynamics (sMD) to assess the allosteric potential of known binders. sMD simulations are employed to sample protein conformational space inaccessible to routine equilibrium MD timescales. Protein conformations sampled by sMD provide starting points for seeded MD simulations, which are combined into MSMs. The methodology is demonstrated on a dataset of protein tyrosine phosphatase 1B ligands. Experimentally confirmed allosteric inhibitors are correctly classified as inhibitors, whereas the deconstructed analogues show reduced inhibitory activity. Analysis of the MSMs provide insights into preferred protein-ligand arrangements that correlate with functional outcomes. The present methodology may find applications for progressing fragments towards lead molecules in FBDD campaigns.

Local Xenon-Protein Interaction Produces Global Conformational Change and Allosteric Inhibition in Lysozyme

Noble gases have well-established biological effects, yet their molecular mechanisms remain poorly understood. Here, we investigated, both experimentally and computationally, the molecular modes of xenon (Xe) action in bacteriophage T4 lysozyme (T4L). By combining indirect gassing methods with a colorimetric lysozyme activity assay, a reversible, Xe-specific (20 ± 3)% inhibition effect was observed. Accelerated molecular dynamic simulations revealed that Xe exerts allosteric inhibition on the protein by expanding a C-terminal hydrophobic cavity. Xe-induced cavity expansion results in global conformational changes, with long-range transduction distorting the active site where peptidoglycan binds. Interestingly, the peptide substrate binding site that enables lysozyme specificity does not change conformation. Two T4L mutants designed to reshape the C-terminal Xe cavity established a correlation between cavity expansion and enzyme inhibition. This work also highlights the use of Xe flooding simulations to identify new cryptic binding pockets. These results enrich our understanding of Xe-protein interactions at the molecular level and inspire further biochemical investigations with noble gases.

Binding pocket dynamics along the recovery stroke of human β-cardiac myosin

The druggability of small-molecule binding sites can be significantly affected by protein motions and conformational changes. Ligand binding, protein dynamics and protein function have been shown to be closely interconnected in myosins. The breakthrough discovery of omecamtiv mecarbil (OM) has led to an increased interest in small molecules that can target myosin and modulate its function for therapeutic purposes (myosin modulators). In this work, we use a combination of computational methods, including steered molecular dynamics, umbrella sampling and binding pocket tracking tools, to follow the evolution of the OM binding site during the recovery stroke transition of human β-cardiac myosin. We found that steering two internal coordinates of the motor domain can recapture the main features of the transition and in particular the rearrangements of the binding site, which shows significant changes in size, shape and composition. Possible intermediate conformations were also identified, in remarkable agreement with experimental findings. The differences in the binding site properties observed along the transition can be exploited for the future development of conformation-selective myosin modulators.

A pH-dependent cluster of charges in a conserved cryptic pocket on flaviviral envelopes

Flaviviruses are enveloped viruses which include human pathogens that are predominantly transmitted by mosquitoes and ticks. Some, such as dengue virus, exhibit the phenomenon of antibody-dependent enhancement (ADE) of disease, making vaccine-based routes of fighting infections problematic. The pH-dependent conformational change of the envelope (E) protein required for fusion between the viral and endosomal membranes is an attractive point of inhibition by antivirals as it has the potential to diminish the effects of ADE. We examined six flaviviruses by employing large-scale molecular dynamics (MD) simulations of raft systems that represent a substantial portion of the flaviviral envelope. We utilised a benzene-mapping approach that led to a discovery of shared hotspots and conserved cryptic sites. A cryptic pocket previously shown to bind a detergent molecule exhibited strain-specific characteristics. An alternative conserved cryptic site at the E protein domain interfaces showed a consistent dynamic behaviour across flaviviruses and contained a conserved cluster of ionisable residues. Constant-pH simulations revealed cluster and domain-interface disruption under low pH conditions. Based on this, we propose a cluster-dependent mechanism that addresses inconsistencies in the histidine-switch hypothesis and highlights the role of cluster protonation in orchestrating the domain dissociation pivotal for the formation of the fusogenic trimer.