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.

Genome-wide analysis of copy-number variation in humans with cleft lip and/or cleft palate identifies COBLL1, RIC1, and ARHGEF38 as clefting genes.

Cleft lip with or without cleft palate (CL/P) is a common birth defect with a complex, heterogeneous etiology. It is well established that common and rare sequence variants contribute to the formation of CL/P, but the contribution of copy-number variants (CNVs) to cleft formation remains relatively understudied. To fill this knowledge gap, we conducted a large-scale comparative analysis of genome-wide CNV profiles of 869 individuals from the Philippines and 233 individuals of European ancestry with CL/P with three primary goals: first, to evaluate whether differences in CNV number, amount of genomic content, or amount of coding genomic content existed within clefting subtypes; second, to assess whether CNVs in our cohort overlapped with known Mendelian clefting loci; and third, to identify unestablished Mendelian clefting genes. Significant differences in CNVs across cleft types or in individuals with non-syndromic versus syndromic clefts were not observed; however, several CNVs in our cohort overlapped with known syndromic and non-syndromic Mendelian clefting loci. Moreover, employing a filtering strategy relying on population genetics data that rare variants are on the whole more deleterious than common variants, we identify several CNV-associated gene losses likely driving non-syndromic clefting phenotypes. By prioritizing genes deleted at a rare frequency across multiple individuals with clefts yet enriched in our cohort of individuals with clefts compared to control subjects, we identify COBLL1, RIC1, and ARHGEF38 as clefting genes. CRISPR-Cas9 mutagenesis of these genes in Xenopus laevis and Danio rerio yielded craniofacial dysmorphologies, including clefts analogous to those seen in human clefting disorders.

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.

Maternal Wnt11b regulates cortical rotation during Xenopus axis formation: analysis of maternal-effect wnt11b mutants.

Asymmetric signalling centres in the early embryo are essential for axis formation in vertebrates. These regions (e.g. amphibian dorsal morula, mammalian anterior visceral endoderm) require stabilised nuclear ß-catenin, but the role of localised Wnt ligand signalling activity in their establishment remains unclear. In Xenopus, dorsal ß-catenin is initiated by vegetal microtubule-mediated symmetry breaking in the fertilised egg, known as 'cortical rotation'. Localised wnt11b mRNA and ligand-independent activators of ß-catenin have been implicated in dorsal ß-catenin activation, but the extent to which each contributes to axis formation in this paradigm remains unclear. Here, we describe a CRISPR-mediated maternal-effect mutation in Xenopus laevis wnt11b.L. We find that wnt11b is maternally required for robust dorsal axis formation and for timely gastrulation, and zygotically for left-right asymmetry. Importantly, we show that vegetal microtubule assembly and cortical rotation are reduced in wnt11b mutant eggs. In addition, we show that activated Wnt coreceptor Lrp6 and Dishevelled lack behaviour consistent with roles in early ß-catenin stabilisation, and that neither is regulated by Wnt11b. This work thus implicates Wnt11b in the distribution of putative dorsal determinants rather than in comprising the determinants themselves. This article has an associated 'The people behind the papers' interview.

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.

Insights into the Effects of Cancer Associated Mutations at the UPF2 and ATP-Binding Sites of NMD Master Regulator: UPF1.

Nonsense-mediated mRNA decay (NMD) is a quality control mechanism that recognizes post-transcriptionally abnormal transcripts and mediates their degradation. The master regulator of NMD is UPF1, an enzyme with intrinsic ATPase and helicase activities. The cancer genomic sequencing data has identified frequently mutated residues in the CH-domain and ATP-binding site of UPF1. In silico screening of UPF1 stability change as a function over 41 cancer mutations has identified five variants with significant effects: K164R, R253W, T499M, E637K, and E833K. To explore the effects of these mutations on the associated energy landscape of UPF1, molecular dynamics simulations (MDS) were performed. MDS identified stable H-bonds between residues S152, S203, S205, Q230/R703, and UPF2/AMPPNP, and suggest that phosphorylation of Serine residues may control UPF1-UPF2 binding. Moreover, the alleles K164R and R253W in the CH-domain improved UPF1-UPF2 binding. In addition, E637K and E833K alleles exhibited improved UPF1-AMPPNP binding compared to the T499M variant; the lower binding is predicted from hindrance caused by the side-chain of T499M to the docking of the tri-phosphate moiety (AMPPNP) into the substrate site. The dynamics of wild-type/mutant systems highlights the flexible nature of the ATP-binding region in UPF1. These insights can facilitate the development of drug discovery strategies for manipulating NMD signaling in cell systems using chemical tools.

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.

Topologically correct central projections of tetrapod inner ear afferents require Fzd3.

Inner ear sensory afferent connections establish sensory maps between the inner ear hair cells and the vestibular and auditory nuclei to allow vestibular and sound information processing. While molecular guidance of sensory afferents to the periphery has been well studied, molecular guidance of central projections from the ear is only beginning to emerge. Disorganized central projections of spiral ganglion neurons in a Wnt/PCP pathway mutant, Prickle1, suggest the Wnt/PCP pathway plays a role in guiding cochlear afferents to the cochlear nuclei in the hindbrain, consistent with known expression of the Wnt receptor, Frizzled3 (Fzd3) in inner ear neurons. We therefore investigated the role of Wnt signaling in central pathfinding in Fzd3 mutant mice and Fzd3 morpholino treated frogs and found aberrant central projections of vestibular afferents in both cases. Ear transplantations from knockdown to control Xenopus showed that it is the Fzd3 expressed within the ear that mediates this guidance. Also, cochlear afferents of Fzd3 mutant mice lack the orderly topological organization observed in controls. Quantification of Fzd3 expression in spiral ganglion neurons show a gradient of expression with Fzd3 being higher in the apex than in the base. Together, these results suggest that a gradient of Fzd3 in inner ear afferents directs projections to the correct dorsoventral column within the hindbrain.

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.

A single KH domain in Bicaudal-C links mRNA binding and translational repression functions to maternal development.

Bicaudal-C (Bicc1) is a conserved RNA-binding protein that represses the translation of selected mRNAs to control development. In Xenopus embryos, Bicc1 binds and represses specific maternal mRNAs to control anterior-posterior cell fates. However, it is not known how Bicc1 binds its RNA targets or how binding affects Bicc1-dependent embryogenesis. Focusing on the KH domains, we analyzed Bicc1 mutants for their ability to bind RNA substrates in vivo and in vitro Analyses of these Bicc1 mutants demonstrated that a single KH domain, KH2, was crucial for RNA binding in vivo and in vitro, while the KH1 and KH3 domains contributed minimally. The Bicc1 mutants were also assayed for their ability to repress translation, and results mirrored the RNA-binding data, with KH2 being the only domain essential for repression. Finally, maternal knockdown and rescue experiments indicated that the KH domains were essential for the regulation of embryogenesis by Bicc1. These data advance our understanding of how Bicc1 selects target mRNAs and provide the first direct evidence that the RNA binding functions of Bicc1 are essential for both Bicc1-dependent translational repression and maternal vertebrate development.