Cosolvent Analysis Toolkit (CAT): a robust hotspot identification platform for cosolvent simulations of proteins to expand the druggable proteome

Enhancing SILCS-MC via GPU Acceleration and Ligand Conformational Optimization with Genetic and Parallel Tempering Algorithms

In the domain of computer-aided drug design, achieving precise and accurate estimates of ligand-protein binding is paramount in the context of screening extensive drug libraries and performing ligand optimization. A fundamental aspect of the SILCS (site identification by ligand competitive saturation) methodology lies in the generation of comprehensive 3D free-energy functional group affinity maps (FragMaps), encompassing the entirety of the target molecule structure. These FragMaps offer an intricate landscape of functional group affinities across the protein, bilayer, or RNA, acting as the basis for subsequent SILCS-Monte Carlo (MC) simulations wherein ligands are docked to the target molecule. To augment the efficiency and breadth of ligand sampling capabilities, we implemented an improved SILCS-MC methodology. By harnessing the parallel computing capability of GPUs, our approach facilitates concurrent calculations over multiple ligands and binding sites, markedly enhancing the computational efficiency. Moreover, the integration of a genetic algorithm (GA) with MC allows us to employ an evolutionary approach to perform ligand sampling, assuring enhanced convergence characteristics. In addition, the potential utility of parallel tempering (PT) to improve sampling was investigated. Implementation of SILCS-MC on GPU architecture is shown to accelerate the speed of SILCS-MC calculations by over 2-orders of magnitude. Use of GA and PT yield improvements over Markov-chain MC, increasing the precision of the resultant docked orientations and binding free energies, though the extent of improvements is relatively small. Accordingly, significant improvements in speed are obtained through the GPU implementation with minor improvements in the precision of the docking obtained via the tested GA and PT algorithms.

AAp-MSMD: Amino Acid Preference Mapping on Protein-Protein Interaction Surfaces Using Mixed-Solvent Molecular Dynamics

Peptides have attracted much attention recently owing to their well-balanced properties as drugs against protein-protein interaction (PPI) surfaces. Molecular simulation-based predictions of binding sites and amino acid residues with high affinity to PPI surfaces are expected to accelerate the design of peptide drugs. Mixed-solvent molecular dynamics (MSMD), which adds probe molecules or fragments of functional groups as solutes to the hydration model, detects the binding hotspots and cryptic sites induced by small molecules. The detection results vary depending on the type of probe molecule; thus, they provide important information for drug design. For rational peptide drug design using MSMD, we proposed MSMD with amino acid residue probes, named amino acid probe-based MSMD (AAp-MSMD), to detect hotspots and identify favorable amino acid types on protein surfaces to which peptide drugs bind. We assessed our method in terms of hotspot detection at the amino acid probe level and binding free energy prediction with amino acid probes at the PPI site for the complex structure that formed the PPI. In hotspot detection, the max-spatial probability distribution map (max-PMAP) obtained from AAp-MSMD detected the PPI site, to which each type of amino acid can bind favorably. In the binding free energy prediction using amino acid probes, ΔGFE obtained from AAp-MSMD roughly estimated the experimental binding affinities from the structure-activity relationship. AAp-MSMD, with amino acid probes, provides estimated binding sites and favorable amino acid types at the PPI site of a target protein.

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.

Mechanistic Investigation of the Androgen Receptor DNA-Binding Domain and Modulation via Direct Interactions with DNA Abasic Sites: Understanding the Mechanisms Involved in Castration-Resistant Prostate Cancer

The androgen receptor (AR) is an important drug target in prostate cancer and a driver of castration-resistant prostate cancer (CRPC). A significant challenge in designing effective drugs lies in targeting constitutively active AR variants and, most importantly, nearly all AR variants lacking the ligand-binding domain (LBD). Recent findings show that an AR's constitutive activity may occur in the presence of somatic DNA mutations within non-coding regions, but the role of these mutations remains elusive. The discovery of new drugs targeting CRPC is hampered by the limited molecular understanding of how AR binds mutated DNA sequences, frequently observed in prostate cancer, and how mutations within the protein and DNA regulate AR-DNA interactions. Using atomistic molecular dynamics (MD) simulations and quantum mechanical calculations, we focused our efforts on (i) rationalising the role of several activating DBD mutations linked to prostate cancer, and (ii) DBD interactions in the presence of abasic DNA lesions, which frequently occur in CRPC. Our results elucidate the role of mutations within DBD through their modulation of the intrinsic dynamics of the DBD-DNA ternary complex. Furthermore, our results indicate that the DNA apurinic lesions occurring in the androgen-responsive element (ARE) enhance direct AR-DNA interactions and stabilise the DBD homodimerisation interface. Moreover, our results strongly suggest that those abasic lesions may form reversible covalent crosslinks between DNA and lysine residues of an AR via a Schiff base. In addition to providing an atomistic model explaining how protein mutations within the AR DNA-binding domain affect AR dimerisation and AR-DNA interactions, our findings provide insight into how somatic mutations occurring in DNA non-coding regions may activate ARs. These mutations are frequently observed in prostate cancer and may contribute to disease progression by enhancing direct AR-DNA interactions.

The druggable genome: Twenty years later

The concept of the druggable genome has been with us for 20 years. During this time, researchers have developed several methods and resources to help assess a target’s druggability. In parallel, evidence for target-disease associations has been collated at scale by Open Targets. More recently, the Protein Data Bank in Europe (PDBe) have built a knowledge base matching per-residue annotations with available protein structure. While each resource is useful in isolation, we believe there is enormous potential in bringing all relevant data into a single knowledge graph, from gene-level to protein residue. Automation is vital for the processing and assessment of all available structures. We have developed scalable, automated workflows that provide hotspot-based druggability assessments for all available structures across large numbers of targets. Ultimately, we will run our method at a proteome scale, an ambition made more realistic by the arrival of AlphaFold 2. Bringing together annotations from the residue up to the gene level and building connections within the graph to represent pathways or protein-protein interactions will create complexity that mirrors the biological systems they represent. Such complexity is difficult for the human mind to utilise effectively, particularly at scale. We believe that graph-based AI methods will be able to expertly navigate such a knowledge graph, selecting the targets of the future.

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.

Inverse Mixed-Solvent Molecular Dynamics for Visualization of the Residue Interaction Profile of Molecular Probes

To ensure efficiency in discovery and development, the application of computational technology is essential. Although virtual screening techniques are widely applied in the early stages of drug discovery research, the computational methods used in lead optimization to improve activity and reduce the toxicity of compounds are still evolving. In this study, we propose a method to construct the residue interaction profile of the chemical structure used in the lead optimization by performing “inverse” mixed-solvent molecular dynamics (MSMD) simulation. Contrary to constructing a protein-based, atom interaction profile, we constructed a probe-based, protein residue interaction profile using MSMD trajectories. It provides us the profile of the preferred protein environments of probes without co-crystallized structures. We assessed the method using three probes: benzamidine, catechol, and benzene. As a result, the residue interaction profile of each probe obtained by MSMD was a reasonable physicochemical description of the general non-covalent interaction. Moreover, comparison with the X-ray structure containing each probe as a ligand shows that the map of the interaction profile matches the arrangement of amino acid residues in the X-ray structure.

Direct Intracellular Delivery of Benzene Diazonium Ions As Observed by Increased Tyrosine Phosphorylation

A challenge within the field of bioconjugation is developing probes to uncover novel information on proteins and other biomolecules. Intracellular delivery of these probes offers the promise of giving relevant context to this information, and these probes can serve as hypothesis-generating tools within complex systems. Leveraging the utility of triazabutadiene chemistry, herein, we discuss the development of a probe that undergoes reduction-mediated deprotection to rapidly deliver a benzene diazonium ion (BDI) into cells. The intracellular BDI resulted in an increase in global tyrosine phosphorylation levels. Seeing phosphatase dysregulation as a potential source of this increase, a tyrosine phosphatase (PTP1B) was tested and shown to be both inhibited and covalently modified by the BDI. In addition to the expected azobenzene formation at tyrosine side chains, key reactive histidine residues were also modified.

Ligand Binding Path Sampling Based on Parallel Cascade Selection Molecular Dynamics: LB-PaCS-MD

Parallel cascade selection molecular dynamics (PaCS-MD) is a rare-event sampling method that generates transition pathways between a reactant and product. To sample the transition pathways, PaCS-MD repeats short-time MD simulations from important configurations as conformational resampling cycles. In this study, PaCS-MD was extended to sample ligand binding pathways toward a target protein, which is referred to as LB-PaCS-MD. In a ligand-concentrated environment, where multiple ligand copies are randomly arranged around the target protein, LB-PaCS-MD allows for the frequent sampling of ligand binding pathways. To select the important configurations, we specified the center of mass (COM) distance between each ligand and the relevant binding site of the target protein, where snapshots generated by the short-time MD simulations were ranked by their COM distance values. From each cycle, snapshots with smaller COM distance values were selected as the important configurations to be resampled using the short-time MD simulations. By repeating conformational resampling cycles, the COM distance values gradually decreased and converged to constants, meaning that a set of ligand binding pathways toward the target protein was sampled by LB-PaCS-MD. To demonstrate relative efficiency, LB-PaCS-MD was applied to several proteins, and their ligand binding pathways were sampled more frequently than conventional MD simulations.

Predicting the Activities of Drug Excipients on Biological Targets using One-Shot Learning

Excipients are major components of drugs and are used to improve drug attributes such as stability and appearance. Excipients approved by the U.S. Food and Drug Administration (FDA) are regarded as safe for humans in allowed concentrations, but their potential interactions with drug targets have not been investigated systematically, which might influence a drug's efficacy. Deep learning models have been used for the identification of ligands that could bind to the drug targets. However, due to the limited available data, it is challenging to reliably estimate the likelihood of a ligand-protein interaction. One-shot learning techniques provide a potential approach to address this low data problem as these techniques require only one or a few examples to classify the new data. In this study, we apply one-shot learning models to data sets that include ligands binding to G-protein-coupled receptors (GPCRs) and kinases. The predicted results suggest that one-shot learning could be used for predicting ligand-protein interactions, and the models attain better performance when protein targets contain conserved binding pockets. The trained models are also used to predict interactions between excipients and drug targets, which provides a potential efficient strategy to explore the activities of drug excipients. We find that a large number of drug excipients could interact with biological targets and influence their function. The results demonstrate how one-shot learning can be used to make accurate predictions for excipient-protein interactions, and these methods could be used for selecting excipients with limited drug-protein interactions.

Application of Site-Identification by Ligand Competitive Saturation in Computer-Aided Drug Design

Site Identification by Ligand Competitive Saturation (SILCS) is a molecular simulation approach that uses diverse small solutes in aqueous solution to obtain functional group affinity patterns of a protein or other macromolecule. This involves employing a combined Grand Canonical Monte Carlo (GCMC)-molecular dynamics (MD) method to sample the full 3D space of the protein, including deep binding pockets and interior cavities from which functional group free energy maps (FragMaps) are obtained. The information content in the maps, which include contributions from protein flexibilty and both protein and functional group desolvation contributions, can be used in many aspects of the drug discovery process. These include identification of novel ligand binding pockets, including allosteric sites, pharmacophore modeling, prediction of relative protein-ligand binding affinities for database screening and lead optimization efforts, evaluation of protein-protein interactions as well as in the formulation of biologics-based drugs including monoclonal antibodies. The present article summarizes the various tools developed in the context of the SILCS methodology and their utility in computer-aided drug design (CADD) applications, showing how the SILCS toolset can improve the drug-development process on a number of fronts with respect to both accuracy and throughput representing a new avenue of CADD applications.

Divide and Conquer. Pocket-Opening Mixed-Solvent Simulations in the Perspective of Docking Virtual Screening Applications for Drug Discovery

The existence of a druggable binding pocket is a prerequisite for computational drug-target interaction studies including virtual screening. Retrospective studies have shown that extended sampling methods like Markov State Modeling and mixed-solvent simulations can identify cryptic pockets relevant for drug discovery. Here, we apply a combination of mixed-solvent molecular dynamics (MD) and time-structure independent component analysis (TICA) to four retrospective case studies: NPC2, the CECR2 bromodomain, TEM-1, and MCL-1. We compare previous experimental and computational findings to our results. It is shown that the successful identification of cryptic pockets depends on the system and the cosolvent probes. We used alternative TICA internal features such as the unbiased backbone coordinates or backbone dihedrals versus biased interatomic distances. We found that in the case of NPC2, TEM-1, and MCL-1, the use of unbiased features is able to identify cryptic pockets, although in the case of the CECR2 bromodomain, more specific features are required to properly capture a pocket opening. In the perspective of virtual screening applications, it is shown how docking studies with the parent ligands depend critically on the conformational state of the targets.

Carnitine is a pharmacological allosteric chaperone of the human lysosomal α-glucosidase

Pompe disease is an inherited metabolic disorder due to the deficiency of the lysosomal acid α-glucosidase (GAA). The only approved treatment is enzyme replacement therapy with the recombinant enzyme (rhGAA). Further approaches like pharmacological chaperone therapy, based on the stabilising effect induced by small molecules on the target enzyme, could be a promising strategy. However, most known chaperones could be limited by their potential inhibitory effects on patient's enzymes. Here we report on the discovery of novel chaperones for rhGAA, L- and D-carnitine, and the related compound acetyl-D-carnitine. These drugs stabilise the enzyme at pH and temperature without inhibiting the activity and acted synergistically with active-site directed pharmacological chaperones. Remarkably, they enhanced by 4-fold the acid α-glucosidase activity in fibroblasts from three Pompe patients with added rhGAA. This synergistic effect of L-carnitine and rhGAA has the potential to be translated into improved therapeutic efficacy of ERT in Pompe disease.

Cosolvent and Dynamic Effects in Binding Pocket Search by Docking Simulations

The lack of conformational sampling in virtual screening projects can lead to inefficient results because many of the potential drugs may not be able to bind to the target protein during the static docking simulations. Here, we performed ensemble docking for around 2000 United States Food and Drug Administration (FDA)-approved drugs with the RNA-dependent RNA polymerase (RdRp) protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) as a target. The representative protein structures were generated by clustering classical molecular dynamics trajectories, which were evolved using three solvent scenarios, namely, pure water, benzene/water and phenol/water mixtures. The introduction of dynamic effects in the theoretical model showed improvement in docking results in terms of the number of strong binders and binding sites in the protein. Some of the discovered pockets were found only for the cosolvent simulations, where the nonpolar probes induced local conformational changes in the protein that lead to the opening of transient pockets. In addition, the selection of the ligands based on a combination of the binding free energy and binding free energy gap between the best two poses for each ligand provided more suitable binders than the selection of ligands based solely on one of the criteria. The application of cosolvent molecular dynamics to enhance the sampling of the configurational space is expected to improve the efficacy of virtual screening campaigns of future drug discovery projects.

EXPRORER: Rational Cosolvent Set Construction Method for Cosolvent Molecular Dynamics Using Large-Scale Computation

Cosolvent molecular dynamics (CMD) simulations involve an MD simulation of a protein in the presence of explicit water molecules mixed with cosolvent molecules to perform hotspot detection, binding site identification, and binding energy estimation, while other existing methods (e.g., MixMD, SILCS, and MDmix) utilize small molecules that represent functional groups of compounds. However, the cosolvent selections employed in these methods differ and there are only a few cosolvents that are commonly used in these methods. In this study, we proposed a systematic method for constructing a set of cosolvents for drug discovery, termed the EXtended PRObes set construction by REpresentative Retrieval (EXPRORER). First, we extracted typical substructures from FDA-approved drugs, generated 138 cosolvent structures, and for each cosolvent molecule, we conducted CMD simulations to generate a spatial probability distribution map of cosolvent atoms (PMAP). Analyses of PMAP similarity revealed that a cosolvent pair with a PMAP similarity greater than 0.70-0.75 shared similar structural features. We present a method for the construction of a cosolvent subset that satisfies a similarity threshold for all cosolvents, and we tested the constructed sets for four proteins. To our knowledge, this is the first study to include a systematic proposal for cosolvent set construction, and thus, the EXPRORER cosolvents will provide deeper insights into ligand binding sites of various proteins.

The calcium sensing receptor (CaSR) promotes Rab27B expression and activity to control secretion in breast cancer cells

Chemotactic and angiogenic factors secreted within the tumor microenvironment eventually facilitate the metastatic dissemination of cancer cells. Calcium-sensing receptor (CaSR) activates secretory pathways in breast cancer cells via a mechanism driven by vesicular trafficking of this receptor. However, it remains to be elucidated how endosomal proteins in secretory vesicles are controlled by CaSR. In the present study, we demonstrate that CaSR promotes expression of Rab27B and activates this secretory small GTPase via PI3K, PKA, mTOR and MADD, a guanine nucleotide exchange factor, also known as DENN/Rab3GEP. Active Rab27B leads secretion of various cytokines and chemokines, including IL-6, IL-1β, IL-8, IP-10 and RANTES. Expression of Rab27B is stimulated by CaSR in MDA-MB-231 and MCF-7 breast epithelial cancer cells, but not in non-cancerous MCF-10A cells. This regulatory mechanism also occurs in HeLa and PC3 cells. Our findings provide insightful information regarding how CaSR activates a Rab27B-dependent mechanism to control secretion of factors known to intervene in paracrine communication circuits within the tumor microenvironment.

Inflammageing in the cardiovascular system: mechanisms, emerging targets, and novel therapeutic strategies

In the elderly population, pathological inflammation has been associated with ageing-associated diseases. The term 'inflammageing', which was used for the first time by Franceschi and co-workers in 2000, is associated with the chronic, low-grade, subclinical inflammatory processes coupled to biological ageing. The source of these inflammatory processes is debated. The senescence-associated secretory phenotype (SASP) has been proposed as the main origin of inflammageing. The SASP is characterised by the release of inflammatory cytokines, elevated activation of the NLRP3 inflammasome, altered regulation of acetylcholine (ACh) nicotinic receptors, and abnormal NAD+ metabolism. Therefore, SASP may be 'druggable' by small molecule therapeutics targeting those emerging molecular targets. It has been shown that inflammageing is a hallmark of various cardiovascular diseases, including atherosclerosis, hypertension, and adverse cardiac remodelling. Therefore, the pathomechanism involving SASP activation via the NLRP3 inflammasome; modulation of NLRP3 via α7 nicotinic ACh receptors; and modulation by senolytics targeting other proteins have gained a lot of interest within cardiovascular research and drug development communities. In this review, which offers a unique view from both clinical and preclinical target-based drug discovery perspectives, we have focused on cardiovascular inflammageing and its molecular mechanisms. We have outlined the mechanistic links between inflammageing, SASP, interleukin (IL)-1β, NLRP3 inflammasome, nicotinic ACh receptors, and molecular targets of senolytic drugs in the context of cardiovascular diseases. We have addressed the 'druggability' of NLRP3 and nicotinic α7 receptors by small molecules, as these proteins represent novel and exciting targets for therapeutic interventions targeting inflammageing in the cardiovascular system and beyond.

Advances in Molecular Dynamics Simulations and Enhanced Sampling Methods for the Study of Protein Systems

Molecular dynamics (MD) simulation is a rigorous theoretical tool that when used efficiently could provide reliable answers to questions pertaining to the structure-function relationship of proteins. Data collated from protein dynamics can be translated into useful statistics that can be exploited to sieve thermodynamics and kinetics crucial for the elucidation of mechanisms responsible for the modulation of biological processes such as protein-ligand binding and protein-protein association. Continuous modernization of simulation tools enables accurate prediction and characterization of the aforementioned mechanisms and these qualities are highly beneficial for the expedition of drug development when effectively applied to structure-based drug design (SBDD). In this review, current all-atom MD simulation methods, with focus on enhanced sampling techniques, utilized to examine protein structure, dynamics, and functions are discussed. This review will pivot around computer calculations of protein-ligand and protein-protein systems with applications to SBDD. In addition, we will also be highlighting limitations faced by current simulation tools as well as the improvements that have been made to ameliorate their efficiency.