AutoDock-SS: AutoDock for Multiconformational Ligand-Based Virtual Screening.

Ligand-based virtual screening (LBVS) can be pivotal for identifying potential drug leads, especially when the target protein's structure is unknown. However, current LBVS methods are limited in their ability to consider the ligand conformational flexibility. This study presents AutoDock-SS (Similarity Searching), which adapts protein-ligand docking for use in LBVS. AutoDock-SS integrates novel ligand-based grid maps and AutoDock-GPU into a novel three-dimensional LBVS workflow. Unlike other approaches based on pregenerated conformer libraries, AutoDock-SS's built-in conformational search optimizes conformations dynamically based on the reference ligand, thus providing a more accurate representation of relevant ligand conformations. AutoDock-SS supports two modes: single and multiple ligand queries, allowing for the seamless consideration of multiple reference ligands. When tested on the Directory of Useful Decoys─Enhanced (DUD-E) data set, AutoDock-SS surpassed alternative 3D LBVS methods, achieving a mean AUROC of 0.775 and an EF1% of 25.72 in single-reference mode. The multireference mode, evaluated on the augmented DUD-E+ data set, demonstrated superior accuracy with a mean AUROC of 0.843 and an EF1% of 34.59. This enhanced performance underscores AutoDock-SS's ability to treat compounds as conformationally flexible while considering the ligand's shape, pharmacophore, and electrostatic potential, expanding the potential of LBVS methods.

From intuition to AI: evolution of small molecule representations in drug discovery.

Within drug discovery, the goal of AI scientists and cheminformaticians is to help identify molecular starting points that will develop into safe and efficacious drugs while reducing costs, time and failure rates. To achieve this goal, it is crucial to represent molecules in a digital format that makes them machine-readable and facilitates the accurate prediction of properties that drive decision-making. Over the years, molecular representations have evolved from intuitive and human-readable formats to bespoke numerical descriptors and fingerprints, and now to learned representations that capture patterns and salient features across vast chemical spaces. Among these, sequence-based and graph-based representations of small molecules have become highly popular. However, each approach has strengths and weaknesses across dimensions such as generality, computational cost, inversibility for generative applications and interpretability, which can be critical in informing practitioners' decisions. As the drug discovery landscape evolves, opportunities for innovation continue to emerge. These include the creation of molecular representations for high-value, low-data regimes, the distillation of broader biological and chemical knowledge into novel learned representations and the modeling of up-and-coming therapeutic modalities.

Modelling of Hepatitis E virus RNA-dependent RNA polymerase genotype 3 from a chronic patient and in silico interaction analysis by molecular docking with Ribavirin.

Hepatitis E Virus (HEV) infection is an emergent zoonotic disease, where chronic hepatitis E associated to solid organ transplant (SOT) recipients, related to genotype 3, is the clinical manifestation of major concern. In this setting, ribavirin (RBV) treatment is the only available therapy, though drug-resistant variants could emerge leading to a therapeutic failure. Crystallographic structures have not been reported for most of the HEV proteins, including the RNA-polymerase (RdRp). Therefore, the mechanism of action of RBV against HEV and the molecular interactions between this drug and RdRp are largely unknown. In this work, we aimed to model in silico the 3 D structure of a novel HEV3 RdRp (HEV_C1_Uy) from a chronically HEV infected-SOT recipient treated with RBV and to perform a molecular docking simulation between RBV triphosphate (RBVT), 7-methyl-guanosine-5'-triphosphate and the modelled protein. The models were generated using I-TASSER server and validated with multiple bioinformatics tools. The docking analysis were carried out with AutoDock Vina and LeDock software. We obtained a suitable model for HEV_C1_Uy (C-Score=-1.33, RMSD = 10.4 ± 4.6 Å). RBVT displayed a binding affinity of -7.6 ± 0.2 Kcal/mol by molecular docking, mediated by 6 hydrogen-bonds (Q195-O14, S198-O11, E257-O13, S260-O2, O3, S311-O11) between the finger's-palm-domains and a free binding energy of 31.26 ± 16.81 kcal/mol by molecular dynamics simulations. We identified the possible HEV RdRp interacting region for incoming nucleotides or analogs and provide novel insights that will contribute to better understand the molecular interactions of RBV and the enzyme and the mechanism of action of this antiviral drug.Communicated by Ramaswamy H. Sarma.

SCORCH: Improving structure-based virtual screening with machine learning classifiers, data augmentation, and uncertainty estimation.

INTRODUCTION: The discovery of a new drug is a costly and lengthy endeavour. The computational prediction of which small molecules can bind to a protein target can accelerate this process if the predictions are fast and accurate enough. Recent machine-learning scoring functions re-evaluate the output of molecular docking to achieve more accurate predictions. However, previous scoring functions were trained on crystalised protein-ligand complexes and datasets of decoys. The limited availability of crystal structures and biases in the decoy datasets can lower the performance of scoring functions. OBJECTIVES: To address key limitations of previous scoring functions and thus improve the predictive performance of structure-based virtual screening. METHODS: A novel machine-learning scoring function was created, named SCORCH (Scoring COnsensus for RMSD-based Classification of Hits). To develop SCORCH, training data is augmented by considering multiple ligand poses and labelling poses based on their RMSD from the native pose. Decoy bias is addressed by generating property-matched decoys for each ligand and using the same methodology for preparing and docking decoys and ligands. A consensus of 3 different machine learning approaches is also used to improve performance. RESULTS: We find that multi-pose augmentation in SCORCH improves its docking power and screening power on independent benchmark datasets. SCORCH outperforms an equivalent scoring function trained on single poses, with a 1 % enrichment factor (EF) of 13.78 vs. 10.86 on 18 DEKOIS 2.0 targets and a mean native pose rank of 5.9 vs 30.4 on CSAR 2014. Additionally, SCORCH outperforms widely used scoring functions in virtual screening and pose prediction on independent benchmark datasets. CONCLUSION: By rationally addressing key limitations of previous scoring functions, SCORCH improves the performance of virtual screening. SCORCH also provides an estimate of its uncertainty, which can help reduce the cost and time required for drug discovery.

Comparison of ATP-binding pockets and discovery of homologous recombination inhibitors.

The ATP binding sites of many enzymes are structurally related, which complicates their development as therapeutic targets. In this work, we explore a diverse set of ATPases and compare their ATP binding pockets using different strategies, including direct and indirect structural methods, in search of pockets attractive for drug discovery. We pursue different direct and indirect structural strategies, as well as ligandability assessments to help guide target selection. The analyses indicate human RAD51, an enzyme crucial in homologous recombination, as a promising, tractable target. Inhibition of RAD51 has shown promise in the treatment of certain cancers but more potent inhibitors are needed. Thus, we design compounds computationally against the ATP binding pocket of RAD51 with consideration of multiple criteria, including predicted specificity, drug-likeness, and toxicity. The molecules designed are evaluated experimentally using molecular and cell-based assays. Our results provide two novel hit compounds against RAD51 and illustrate a computational pipeline to design new inhibitors against ATPases.

Evidence that nuclear receptors are related to terpene synthases.

Ligand-activated nuclear receptors (NRs) orchestrate development, growth, and reproduction across all animal lifeforms - the Metazoa - but how NRs evolved remains mysterious. Given the NR ligands including steroids and retinoids are predominantly terpenoids, we asked whether NRs might have evolved from enzymes that catalyze terpene synthesis and metabolism. We provide evidence suggesting that NRs may be related to the terpene synthase (TS) enzyme superfamily. Based on over 10,000 3D structural comparisons, we report that the NR ligand-binding domain and TS enzymes share a conserved core of seven α-helical segments. In addition, the 3D locations of the major ligand-contacting residues are also conserved between the two protein classes. Primary sequence comparisons reveal suggestive similarities specifically between NRs and the subfamily of cis-isoprene transferases, notably with dehydrodolichyl pyrophosphate synthase and its obligate partner, NUS1/NOGOB receptor. Pharmacological overlaps between NRs and TS enzymes add weight to the contention that they share a distant evolutionary origin, and the combined data raise the possibility that a ligand-gated receptor may have arisen from an enzyme antecedent. However, our findings do not formally exclude other interpretations such as convergent evolution, and further analysis will be necessary to confirm the inferred relationship between the two protein classes.

A 5-FU Precursor Designed to Evade Anabolic and Catabolic Drug Pathways and Activated by Pd Chemistry In Vitro and In Vivo.

5-Fluorouracil (5-FU) is an antineoplastic antimetabolite that is widely administered to cancer patients by bolus injection, especially to those suffering from colorectal and pancreatic cancer. Because of its suboptimal route of administration and dose-limiting toxicities, diverse 5-FU prodrugs have been developed to confer oral bioavailability and increase the safety profile of 5-FU chemotherapy regimens. Our contribution to this goal is presented herein with the development of a novel palladium-activated prodrug designed to evade the metabolic machinery responsible for 5-FU anabolic activation and catabolic processing. The new prodrug is completely innocuous to cells and highly resistant to metabolization by primary hepatocytes and liver S9 fractions (the main metabolic route for 5-FU degradation), whereas it is rapidly converted into 5-FU in the presence of a palladium (Pd) source. In vivo pharmokinetic analysis shows the prodrug is rapidly and completely absorbed after oral administration and exhibits a longer half-life than 5-FU. In vivo efficacy studies in a xenograft colon cancer model served to prove, for the first time, that orally administered prodrugs can be locally converted to active drugs by intratumorally inserted Pd implants.

Correction: Bypassing the requirement for aminoacyl-tRNA by a cyclodipeptide synthase enzyme.

[This corrects the article DOI: 10.1039/D0CB00142B.].

Bypassing the requirement for aminoacyl-tRNA by a cyclodipeptide synthase enzyme.

Cyclodipeptide synthases (CDPSs) produce a variety of cyclic dipeptide products by utilising two aminoacylated tRNA substrates. We sought to investigate the minimal requirements for substrate usage in this class of enzymes as the relationship between CDPSs and their substrates remains elusive. Here, we investigated the Bacillus thermoamylovorans enzyme, BtCDPS, which synthesises cyclo(l-Leu-l-Leu). We systematically tested where specificity arises and, in the process, uncovered small molecules (activated amino esters) that will suffice as substrates, although catalytically poor. We solved the structure of BtCDPS to 1.7 Å and combining crystallography, enzymatic assays and substrate docking experiments propose a model for how the minimal substrates interact with the enzyme. This work is the first report of a CDPS enzyme utilizing a molecule other than aa-tRNA as a substrate; providing insights into substrate requirements and setting the stage for the design of improved simpler substrates.

Design of drug-like hepsin inhibitors against prostate cancer and kidney stones.

Hepsin, a transmembrane serine protease abundant in renal endothelial cells, is a promising therapeutic target against several cancers, particularly prostate cancer. It is involved in the release and polymerization of uromodulin in the urine, which plays a role in kidney stone formation. In this work, we design new potential hepsin inhibitors for high activity, improved specificity towards hepsin, and promising ADMET properties. The ligands were developed in silico through a novel hierarchical pipeline. This pipeline explicitly accounts for off-target binding to the related serine proteases matriptase and HGFA (human hepatocyte growth factor activator). We completed the pipeline incorporating ADMET properties of the candidate inhibitors into custom multi-objective optimization functions. The ligands designed show excellent prospects for targeting hepsin via the blood stream and the urine and thus enable key experimental studies. The computational pipeline proposed is remarkably cost-efficient and can be easily adapted for designing inhibitors against new drug targets.

PRL3-DDX21 Transcriptional Control of Endolysosomal Genes Restricts Melanocyte Stem Cell Differentiation.

Melanocytes, replenished throughout life by melanocyte stem cells (MSCs), play a critical role in pigmentation and melanoma. Here, we reveal a function for the metastasis-associated phosphatase of regenerating liver 3 (PRL3) in MSC regeneration. We show that PRL3 binds to the RNA helicase DDX21, thereby restricting productive transcription by RNAPII at master transcription factor (MITF)-regulated endolysosomal vesicle genes. In zebrafish, this mechanism controls premature melanoblast expansion and differentiation from MSCs. In melanoma patients, restricted transcription of this endolysosomal vesicle pathway is a hallmark of PRL3-high melanomas. Our work presents the conceptual advance that PRL3-mediated control of transcriptional elongation is a differentiation checkpoint mechanism for activated MSCs and has clinical relevance for the activity of PRL3 in regenerating tissue and cancer.

Recognition Dynamics of Cancer Mutations on the ERp57-Tapasin Interface.

Down regulation of the major histocompatibility class (MHC) I pathway plays an important role in tumour development, and can be achieved by suppression of HLA expression or mutations in the MHC peptide-binding pocket. The peptide-loading complex (PLC) loads peptides on the MHC-I molecule in a dynamic multi-step assembly process. The effects of cancer variants on ERp57 and tapasin components from the MHC-I pathway is less known, and they could have an impact on antigen presentation. Applying computational approaches, we analysed whether the ERp57-tapasin binding might be altered by missense mutations. The variants H408R(ERp57) and P96L, D100A, G183R(tapasin) at the protein-protein interface improved protein stability (ΔΔG) during the initial screen of 14 different variants. The H408R(ERp57) and P96L(tapasin) variants, located close to disulphide bonds, were further studied by molecular dynamics (MD). Identifying intramolecular a-a' domain interactions, MD revealed open and closed conformations of ERp57 in the presence and absence of tapasin. In wild-type and mutant ERp57-tapasin complexes, residues Val97, Ser98, Tyr100, Trp405, Gly407(ERp57) and Asn94, Cys95, Arg97, Asp100(tapasin) formed common H-bond interactions. Moreover, comparing the H-bond networks for P96L and H408R with each other, suggests that P96L(tapasin) improved ERp57-tapasin binding more than the H408R(ERp57) mutant. During MD, the C-terminus domain (that binds MHC-I) in tapasin from the ERp57(H408R)-tapasin complex moved away from the PLC, whereas in the ERp57-tapasin(P96L) system was oppositely displaced. These findings can have implications for the function of PLC and, ultimately, for the presentation of MHC-I peptide complex on the tumour cell surface.

Lamin A molecular compression and sliding as mechanisms behind nucleoskeleton elasticity.

Lamin A is a nuclear intermediate filament protein critical for nuclear architecture and mechanics and mutated in a wide range of human diseases. Yet little is known about the molecular architecture of lamins and mechanisms of their assembly. Here we use SILAC cross-linking mass spectrometry to determine interactions within lamin dimers and between dimers in higher-order polymers. We find evidence for a compression mechanism where coiled coils in the lamin A rod can slide onto each other to contract rod length, likely driven by a wide range of electrostatic interactions with the flexible linkers between coiled coils. Similar interactions occur with unstructured regions flanking the rod domain during oligomeric assembly. Mutations linked to human disease block these interactions, suggesting that this spring-like contraction can explain in part the dynamic mechanical stretch and flexibility properties of the lamin polymer and other intermediate filament networks.

Intrinsic lipid binding activity of ATG16L1 supports efficient membrane anchoring and autophagy.

Membrane targeting of autophagy-related complexes is an important step that regulates their activities and prevents their aberrant engagement on non-autophagic membranes. ATG16L1 is a core autophagy protein implicated at distinct phases of autophagosome biogenesis. In this study, we dissected the recruitment of ATG16L1 to the pre-autophagosomal structure (PAS) and showed that it requires sequences within its coiled-coil domain (CCD) dispensable for homodimerisation. Structural and mutational analyses identified conserved residues within the CCD of ATG16L1 that mediate direct binding to phosphoinositides, including phosphatidylinositol 3-phosphate (PI3P). Mutating putative lipid binding residues abrogated the localisation of ATG16L1 to the PAS and inhibited LC3 lipidation. On the other hand, enhancing lipid binding of ATG16L1 by mutating negatively charged residues adjacent to the lipid binding motif also resulted in autophagy inhibition, suggesting that regulated recruitment of ATG16L1 to the PAS is required for its autophagic activity. Overall, our findings indicate that ATG16L1 harbours an intrinsic ability to bind lipids that plays an essential role during LC3 lipidation and autophagosome maturation.

An inter-subunit protein-peptide interface that stabilizes the specific activity and oligomerization of the AAA+ chaperone Reptin.

Reptin is a member of the AAA+ superfamily whose members can exist in equilibrium between monomeric apo forms and ligand bound hexamers. Inter-subunit protein-protein interfaces that stabilize Reptin in its oligomeric state are not well-defined. A self-peptide binding assay identified a protein-peptide interface mapping to an inter-subunit "rim" of the hexamer bridged by Tyrosine-340. A Y340A mutation reduced ADP-dependent oligomer formation using a gel filtration assay, suggesting that Y340 forms a dominant oligomer stabilizing side chain. The monomeric ReptinY340A mutant protein exhibited increased activity to its partner protein AGR2 in an ELISA assay, further suggesting that hexamer formation can preclude certain protein interactions. Hydrogen-deuterium exchange mass spectrometry (HDX-MS) demonstrated that the Y340A mutation attenuated deuterium suppression of Reptin in this motif in the presence of ligand. By contrast, the tyrosine motif of Reptin interacts with a shallower pocket in the hetero-oligomeric structure containing Pontin and HDX-MS revealed no obvious role of the Y340 side chain in stabilizing the Reptin-Pontin oligomer. Molecular dynamic simulations (MDS) rationalized how the Y340A mutation impacts upon a normally stabilizing inter-subunit amino acid contact. MDS also revealed how the D299N mutation can, by contrast, remove oligomer de-stabilizing contacts. These data suggest that the Reptin interactome can be regulated by a ligand dependent equilibrium between monomeric and hexameric forms through a hydrophobic inter-subunit protein-protein interaction motif bridged by Tyrosine-340. SIGNIFICANCE: Discovering dynamic protein-protein interactions is a fundamental aim of research in the life sciences. An emerging view of protein-protein interactions in higher eukaryotes is that they are driven by small linear polypeptide sequences; the linear motif. We report on the use of linear-peptide motif screens to discover a relatively high affinity peptide-protein interaction for the AAA+ and pro-oncogenic protein Reptin. This peptide interaction site was shown to form a 'hot-spot' protein-protein interaction site, and validated to be important for ligand-induced oligomerization of the Reptin protein. These biochemical data provide a foundation to understand how single point mutations in Reptin can impact on its oligomerization and protein-protein interaction landscape.

Oligopeptide Signaling through TbGPR89 Drives Trypanosome Quorum Sensing.

Trypanosome parasites control their virulence and spread by using quorum sensing (QS) to generate transmissible "stumpy forms" in their host bloodstream. However, the QS signal "stumpy induction factor" (SIF) and its reception mechanism are unknown. Although trypanosomes lack G protein-coupled receptor signaling, we have identified a surface GPR89-family protein that regulates stumpy formation. TbGPR89 is expressed on bloodstream "slender form" trypanosomes, which receive the SIF signal, and when ectopically expressed, TbGPR89 drives stumpy formation in a SIF-pathway-dependent process. Structural modeling of TbGPR89 predicts unexpected similarity to oligopeptide transporters (POT), and when expressed in bacteria, TbGPR89 transports oligopeptides. Conversely, expression of an E. coli POT in trypanosomes drives parasite differentiation, and oligopeptides promote stumpy formation in vitro. Furthermore, the expression of secreted trypanosome oligopeptidases generates a paracrine signal that accelerates stumpy formation in vivo. Peptidase-generated oligopeptide QS signals being received through TbGPR89 provides a mechanism for both trypanosome SIF production and reception.

ALDH1 Bio-activates Nifuroxazide to Eradicate ALDHHigh Melanoma-Initiating Cells.

5-Nitrofurans are antibiotic pro-drugs that have potential as cancer therapeutics. Here, we show that 5-nitrofurans can be bio-activated by aldehyde dehydrogenase (ALDH) 1A1/1A3 enzymes that are highly expressed in a subpopulation of cancer-initiating (stem) cells. We discover that the 5-nitrofuran, nifuroxazide, is selective for bio-activation by ALDH1 isoforms over ALDH2, whereby it both oxidizes ALDH1 and is converted to cytotoxic metabolites in a two-hit pro-drug mechanism. We show that ALDH1High melanoma cells are sensitive to nifuroxazide, while ALDH1A3 loss-of-function mutations confer drug resistance. In tumors, nifuroxazide targets ALDH1High melanoma subpopulations with the subsequent loss of melanoma-initiating cell potential. BRAF and MEK inhibitor therapy increases ALDH1 expression in patient melanomas, and effectively combines with nifuroxazide in melanoma cell models. The selective eradication of ALDH1High cells by nifuroxazide-ALDH1 activation goes beyond current strategies based on inhibiting ALDH1 and provides a rational basis for the nifuroxazide mechanism of action in cancer.

Fatty-acylation target sequence in the ligand-binding domain of vertebrate steroid receptors demarcates evolution from estrogen-related receptors.

Present-day nuclear receptors (NRs) responding to adrenal and sex steroids are key regulators of reproduction and growth in mammals, and are thought to have evolved from an ancestral NR most closely related to extant estrogen-related receptors (ERRs). The molecular events (and ligands) that distinguish steroid-activated NRs (SRs) from their inferred ancestor, that gave rise to both the ERRs and SRs, remain unknown. We report that target sequences for fatty-acylation (palmitoylation) at a key cysteine residue (corresponding to Cys447 in human estrogen receptor ERα) in helix 8 of the ligand-binding domain accurately demarcate SRs from ERRs. Docking studies are consistent with the hypothesis that palmitate embeds into a key groove in the receptor surface. The implications of lipidation, and of potential alternative ligands for the key cysteine residue, for receptor function and the evolution of SRs are discussed.

Development of Potent Inhibitors of Receptor Tyrosine Kinases by Ligand-Based Drug Design and Target-Biased Phenotypic Screening.

Pyrazolopyrimidines with potent antiproliferative properties were developed by an adaptive strategy that applies ligand-based design and phenotypic screening iteratively and is informed by biochemical assays. To drive development toward specific oncopathways, compounds were tested against cancer cells that overexpress, or not, AXL kinase. Identified phenotypic hits were found to inhibit oncotargets AXL, RET, and FLT3. Subsequent optimization generated antiproliferative lead compounds with unique selectivity profiles, including selective AXL inhibitors and a highly potent inhibitor of FLT3.

Mono-Substituted Hydrocarbon Diastereomer Combinations Reveal Stapled Peptides with High Structural Fidelity.

Modified peptides, such as stapled peptides, which replicate the structure of α-helical protein segments, represent a potential therapeutic advance. However, the 3D solution structure of these stapled peptides is rarely explored beyond the acquisition of circular dichroism (CD) data to quantify bulk peptide helicity; the detailed backbone structure, which underlies this, is typically undefined. Diastereomeric stapled peptides based on helical sections of three proteins (αSyn, Cks1 and CK1α) were generated; their overall helicity was quantified by CD; and the most helical peptide from each series was selected for structural analysis. Solution-phase models for the optimised peptides were generated using NMR-derived restraints and a modified CHARMM22 force field. Comparing these models with PDB structures allowed deviation between the stapled peptides and critical helical regions to be evaluated. These studies demonstrate that CD alone is not sufficient to assess the structural fidelity of a stapled peptide.