Amphiphilic Heparinoids as Potent Antiviral Agents against SARS-CoV-2.

Herein, we report the synthesis and biological evaluation of a novel series of heparinoid amphiphiles as inhibitors of heparanase and SARS-CoV-2. By employing a tailor-made synthetic strategy, a library of highly sulfated homo-oligosaccharides bearing d-glucose or a C5-epimer (i.e., l-idose or l-iduronic acid) conjugated with various lipophilic groups was synthesized and investigated for antiviral activity. Sulfated higher oligosaccharides of d-glucose or l-idose with lipophilic aglycones displayed potent anti-SARS-CoV-2 and antiheparanse activity, similar to or better than pixatimod (PG545), and were more potent than their isosteric l-iduronic acid congeners. Lipophilic groups such as cholestanol and C18-aliphatic substitution are more advantageous than functional group appended lipophilic moieties. These findings confirm that fine-tuning of higher oligosaccharides, degree of sulfation, and lipophilic groups can yield compounds with potent anti-SARS-CoV-2 activity.

Sequence- and Structure-Dependent Cytotoxicity of Phosphorothioate and 2'-O-Methyl Modified Single-Stranded Oligonucleotides.

Single-stranded oligonucleotides (SSOs) are a rapidly expanding class of therapeutics that comprises antisense oligonucleotides, microRNAs, and aptamers, with ten clinically approved molecules. Chemical modifications such as the phosphorothioate backbone and the 2'-O-methyl ribose can improve the stability and pharmacokinetic properties of therapeutic SSOs, but they can also lead to toxicity in vitro and in vivo through nonspecific interactions with cellular proteins, gene expression changes, disturbed RNA processing, and changes in nuclear structures and protein distribution. In this study, we screened a mini library of 277 phosphorothioate and 2'-O-methyl-modified SSOs, with or without mRNA complementarity, for cytotoxic properties in two cancer cell lines. Using circular dichroism, nucleic magnetic resonance, and molecular dynamics simulations, we show that phosphorothioate- and 2'-O-methyl-modified SSOs that form stable hairpin structures through Watson-Crick base pairing are more likely to be cytotoxic than those that exist in an extended conformation. In addition, moderate and highly cytotoxic SSOs in our dataset have a higher mean purine composition than pyrimidine. Overall, our study demonstrates a structure-cytotoxicity relationship and indicates that the formation of stable hairpins should be a consideration when designing SSOs toward optimal therapeutic profiles.

Characterization of Posttranslationally Modified PHF-1 Tau Peptides Using Gaussian Accelerated Molecular Dynamics Simulation.

The microtubule-associated protein, Tau, is an intrinsically disordered protein that plays a crucial role in neurodegenerative diseases like Alzheimer's disease. The posttranslational modifications across the Tau protein domains are involved in regulating Tau protein's function and disease onset. Of the various posttranslational modifications at Ser, Thr, and Tyr sites, O-GlcNAcylation and phosphorylation are the most critical ones, playing a vital role in Tau aggregation and tauopathies. To understand the function, it is essential to characterize the structural changes associated with Tau modification. Previous experimental studies have focused on high-resolution nuclear magnetic resonance techniques to structurally characterize the effect of phosphorylation, O-GlcNAcylation, and combination of both PTMs on Tau conformation in small peptides centered on the PHF-1 epitope from amino acid 392 to 411. The structural characterization using atomistic molecular dynamics simulation of such disordered peptides requires long simulation time, proper sampling method, and utilization of appropriate force fields for accurate determination of conformational ensembles, resembling the experimental data. This chapter details the protocol for the structural characterization of modified Tau peptides using the CHARMM36m force field and enhanced sampling methods like Gaussian accelerated molecular dynamics (GaMD) simulation. We have focused on a detailed explanation of the GaMD method and analyses of molecular dynamics trajectories to explain the relationship between two modifications, phospho- and glyco-, at C-terminus of Tau protein and its stable conformation over the longer simulation timeframes. The analyses involve energetics reweighting, clustering of simulation trajectories, and characterization of secondary structure using circular dichroism data from the simulation. The reader can utilize this protocol to investigate the structures of complex proteins, especially the disordered ones.

Structural and mechanistic characterization of bifunctional heparan sulfate N-deacetylase-N-sulfotransferase 1.

Heparan sulfate (HS) polysaccharides are major constituents of the extracellular matrix, which are involved in myriad structural and signaling processes. Mature HS polysaccharides contain complex, non-templated patterns of sulfation and epimerization, which mediate interactions with diverse protein partners. Complex HS modifications form around initial clusters of glucosamine-N-sulfate (GlcNS) on nascent polysaccharide chains, but the mechanistic basis underpinning incorporation of GlcNS itself into HS remains unclear. Here, we determine cryo-electron microscopy structures of human N-deacetylase-N-sulfotransferase (NDST)1, the bifunctional enzyme primarily responsible for initial GlcNS modification of HS. Our structures reveal the architecture of both NDST1 deacetylase and sulfotransferase catalytic domains, alongside a non-catalytic N-terminal domain. The two catalytic domains of NDST1 adopt a distinct back-to-back topology that limits direct cooperativity. Binding analyses, aided by activity-modulating nanobodies, suggest that anchoring of the substrate at the sulfotransferase domain initiates the NDST1 catalytic cycle, providing a plausible mechanism for cooperativity despite spatial domain separation. Our data shed light on key determinants of NDST1 activity, and describe tools to probe NDST1 function in vitro and in vivo.

Targeting the hSSB1-INTS3 Interface: A Computational Screening Driven Approach to Identify Potential Modulators.

Human single-stranded DNA binding protein 1 (hSSB1) forms a heterotrimeric complex, known as a sensor of single-stranded DNA binding protein 1 (SOSS1), in conjunction with integrator complex subunit 3 (INTS3) and C9ORF80. This sensory protein plays an important role in homologous recombination repair of double-strand breaks in DNA to efficiently recruit other repair proteins at the damaged sites. Previous studies have identified elevated hSSB1-mediated DNA repair activities in various cancers, highlighting its potential as an anticancer target. While prior efforts have focused on inhibiting hSSB1 by targeting its DNA binding domain, this study seeks to explore the inhibition of the hSSB1 function by disrupting its interaction with the key partner protein INTS3 in the SOSS1 complex. The investigative strategy entails a molecular docking-based screening of a specific compound library against the three-dimensional structure of INTS3 at the hSSB1 binding interface. Subsequent assessments involve in vitro analyses of protein-protein interaction (PPI) disruption and cellular effects through co-immunoprecipitation and immunofluorescence assays, respectively. Moreover, the study includes an evaluation of the structural stability of ligands at the INTS3 hot-spot site using molecular dynamics simulations. The results indicate a potential in vitro disruption of the INTS3-hSSB1 interaction by three of the tested compounds obtained from the virtual screening with one impacting the recruitment of hSSB1 and INTS3 to chromatin following DNA damage. To our knowledge, our results identify the first set of drug-like compounds that functionally target INTS3-hSSB1 interaction, and this provides the basis for further biophysical investigations that should help to speed up PPI inhibitor discovery.

Dynamics and recognition of homeodomain containing protein-DNA complex of IRX4.

Iroquois Homeobox 4 (IRX4) belongs to a family of homeobox TFs having roles in embryogenesis, cell specification, and organ development. Recently, large scale genome-wide association studies and epigenetic studies have highlighted the role of IRX4 and its associated variants in prostate cancer. No studies have investigated and characterized the structural aspect of the IRX4 homeodomain and its potential to bind to DNA. The current study uses sequence analysis, homology modeling, and molecular dynamics simulations to explore IRX4 homeodomain-DNA recognition mechanisms and the role of somatic mutations affecting these interactions. Using publicly available databases, gene expression of IRX4 was found in different tissues, including prostate, heart, skin, vagina, and the protein expression was found in cancer cell lines (HCT166, HEK293), B cells, ascitic fluid, and brain. Sequence conservation of the homeodomain shed light on the importance of N- and C-terminal residues involved in DNA binding. The specificity of IRX4 homodimer bound to consensus human DNA sequence was confirmed by molecular dynamics simulations, representing the role of conserved amino acids including R145, A194, N195, S190, R198, and R199 in binding to DNA. Additional N-terminal residues like T144 and G143 were also found to have specific interactions highlighting the importance of N-terminus of the homeodomain in DNA recognition. Additionally, the effects of somatic mutations, including the conserved Arginine (R145, R198, and R199) residues on DNA binding elucidated the importance of these residues in stabilizing the protein-DNA complex. Secondary structure and hydrogen bonding analysis showed the roles of specific residues (R145, T191, A194, N195, R198, and R199) in maintaining the homogeneity of the structure and its interaction with DNA. The differences in relative binding free energies of all the mutants shed light on the structural modularity of this protein and the dynamics behind protein-DNA interaction. We also have predicted that the C-terminal sequence of the IRX4 homeodomain could act as a potential cell-penetrating peptide, emphasizing the role these small peptides could play in targeting homeobox TFs.

An Exploration of Small Molecules That Bind Human Single-Stranded DNA Binding Protein 1.

Human single-stranded DNA binding protein 1 (hSSB1) is critical to preserving genome stability, interacting with single-stranded DNA (ssDNA) through an oligonucleotide/oligosaccharide binding-fold. The depletion of hSSB1 in cell-line models leads to aberrant DNA repair and increased sensitivity to irradiation. hSSB1 is over-expressed in several types of cancers, suggesting that hSSB1 could be a novel therapeutic target in malignant disease. hSSB1 binding studies have focused on DNA; however, despite the availability of 3D structures, small molecules targeting hSSB1 have not been explored. Quinoline derivatives targeting hSSB1 were designed through a virtual fragment-based screening process, synthesizing them using AlphaLISA and EMSA to determine their affinity for hSSB1. In parallel, we further screened a structurally diverse compound library against hSSB1 using the same biochemical assays. Three compounds with nanomolar affinity for hSSB1 were identified, exhibiting cytotoxicity in an osteosarcoma cell line. To our knowledge, this is the first study to identify small molecules that modulate hSSB1 activity. Molecular dynamics simulations indicated that three of the compounds that were tested bound to the ssDNA-binding site of hSSB1, providing a framework for the further elucidation of inhibition mechanisms. These data suggest that small molecules can disrupt the interaction between hSSB1 and ssDNA, and may also affect the ability of cells to repair DNA damage. This test study of small molecules holds the potential to provide insights into fundamental biochemical questions regarding the OB-fold.

Targeting the COMMD4-H2B protein complex in lung cancer.

BACKGROUND: Lung cancer is the biggest cause of cancer-related deaths worldwide. Non-small cell lung cancer (NSCLC) accounts for 85-90% of all lung cancers. Identification of novel therapeutic targets are required as drug resistance impairs chemotherapy effectiveness. COMMD4 is a potential NSCLC therapeutic target. The aims of this study were to investigate the COMMD4-H2B binding pose and develop a short H2B peptide that disrupts the COMMD4-H2B interaction and mimics COMMD4 siRNA depletion. METHODS: Molecular modelling, in vitro binding and site-directed mutagenesis were used to identify the COMMD4-H2B binding pose and develop a H2B peptide to inhibit the COMMD4-H2B interaction. Cell viability, DNA repair and mitotic catastrophe assays were performed to determine whether this peptide can specially kill NSCLC cells. RESULTS: Based on the COMMD4-H2B binding pose, we have identified a H2B peptide that inhibits COMMD4-H2B by directly binding to COMMD4 on its H2B binding binding site, both in vitro and in vivo. Treatment of NSCLC cell lines with this peptide resulted in increased sensitivity to ionising radiation, increased DNA double-strand breaks and induction of mitotic catastrophe in NSCLC cell lines. CONCLUSIONS: Our data shows that COMMD4-H2B represents a novel potential NSCLC therapeutic target.

Integration of molecular modelling and in vitro studies to inhibit LexA proteolysis.

Introduction: As antibiotic resistance has become more prevalent, the social and economic impacts are increasingly pressing. Indeed, bacteria have developed the SOS response which facilitates the evolution of resistance under genotoxic stress. The transcriptional repressor, LexA, plays a key role in this response. Mutation of LexA to a non-cleavable form that prevents the induction of the SOS response sensitizes bacteria to antibiotics. Achieving the same inhibition of proteolysis with small molecules also increases antibiotic susceptibility and reduces drug resistance acquisition. The availability of multiple LexA crystal structures, and the unique Ser-119 and Lys-156 catalytic dyad in the protein enables the rational design of inhibitors. Methods: We pursued a binary approach to inhibit proteolysis; we first investigated ß-turn mimetics, and in the second approach we tested covalent warheads targeting the Ser-119 residue. We found that the cleavage site region (CSR) of the LexA protein is a classical Type II ß-turn, and that published 1,2,3-triazole compounds mimic the ß-turn. Generic covalent molecule libraries and a ß-turn mimetic library were docked to the LexA C-terminal domain using molecular modelling methods in FlexX and CovDock respectively. The 133 highest-scoring molecules were screened for their ability to inhibit LexA cleavage under alkaline conditions. The top molecules were then tested using a RecA-mediated cleavage assay. Results: The ß-turn library screen did not produce any hit compounds that inhibited RecA-mediated cleavage. The covalent screen discovered an electrophilic serine warhead that can inhibit LexA proteolysis, reacting with Ser-119 via a nitrile moiety. Discussion: This research presents a starting point for hit-to-lead optimisation, which could lead to inhibition of the SOS response and prevent the acquisition of antibiotic resistance.

Comparative Assessment of Docking Programs for Docking and Virtual Screening of Ribosomal Oxazolidinone Antibacterial Agents.

Oxazolidinones are a broad-spectrum class of synthetic antibiotics that bind to the 50S ribosomal subunit of Gram-positive and Gram-negative bacteria. Many crystal structures of the ribosomes with oxazolidinone ligands have been reported in the literature, facilitating structure-based design using methods such as molecular docking. It would be of great interest to know in advance how well docking methods can reproduce the correct ligand binding modes and rank these correctly. We examined the performance of five molecular docking programs (AutoDock 4, AutoDock Vina, DOCK 6, rDock, and RLDock) for their ability to model ribosomal-ligand interactions with oxazolidinones. Eleven ribosomal crystal structures with oxazolidinones as the ligands were docked. The accuracy was evaluated by calculating the docked complexes' root-mean-square deviation (RMSD) and the program's internal scoring function. The rankings for each program based on the median RMSD between the native and predicted were DOCK 6 > AD4 > Vina > RDOCK >> RLDOCK. Results demonstrate that the top-performing program, DOCK 6, could accurately replicate the ligand binding in only four of the eleven ribosomes due to the poor electron density of said ribosomal structures. In this study, we have further benchmarked the performance of the DOCK 6 docking algorithm and scoring in improving virtual screening (VS) enrichment using the dataset of 285 oxazolidinone derivatives against oxazolidinone binding sites in the S. aureus ribosome. However, there was no clear trend between the structure and activity of the oxazolidinones in VS. Overall, the docking performance indicates that the RNA pocket's high flexibility does not allow for accurate docking prediction, highlighting the need to validate VS. protocols for ligand-RNA before future use. Later, we developed a re-scoring method incorporating absolute docking scores and molecular descriptors, and the results indicate that the descriptors greatly improve the correlation of docking scores and pMIC values. Morgan fingerprint analysis was also used, suggesting that DOCK 6 underpredicted molecules with tail modifications with acetamide, n-methylacetamide, or n-ethylacetamide and over-predicted molecule derivatives with methylamino bits. Alternatively, a ligand-based approach similar to a field template was taken, indicating that each derivative's tail groups have strong positive and negative electrostatic potential contributing to microbial activity. These results indicate that one should perform VS. campaigns of ribosomal antibiotics with care and that more comprehensive strategies, including molecular dynamics simulations and relative free energy calculations, might be necessary in conjunction with VS. and docking.

Integrative competing endogenous RNA network analyses identify novel lncRNA and genes implicated in metastatic breast cancer.

Competing endogenous RNAs (ceRNAs) have gained attention in cancer research owing to their involvement in microRNA-mediated gene regulation. Previous studies have identified ceRNA networks of individual cancers. Nevertheless, none of these studies has investigated different cancer stages. We identify stage-specific ceRNAs in breast cancer using the cancer genome atlas data. Moreover, we investigate the molecular functions and prognostic ability of ceRNAs involved in stage I-IV networks. We identified differentially expressed candidate ceRNAs using edgeR and limma R packages. A three-step analysis was used to identify statistically significant ceRNAs of each stage. Survival analysis and functional enrichment analysis were conducted to identify molecular functions and prognostic ability. We found five genes and one long non-coding RNA unique to the stage IV ceRNA network. These genes have been described in previous breast cancer studies. Genes acted as ceRNAs are enriched in cancer-associated pathways. Two, three, and three microRNAs from stages I, II, and III were prognostic from the Kaplan-Meier survival analysis. Our results reveal a set of unique ceRNAs in metastatic breast cancer. Further experimental work is required to evaluate their role in metastasis. Moreover, identifying stage-specific ceRNAs will improve the understanding of personalised therapeutics in breast cancer.

Structural investigation of CDCA3-Cdh1 protein-protein interactions using in vitro studies and molecular dynamics simulation.

The anaphase-promoting complex/cyclosome (APC/C) ubiquitin ligase and its cofactor, Cdh1, regulate the expression of several cell-cycle proteins and their functions during mitosis. Levels of the protein cell division cycle-associated protein 3 (CDCA3), which is functionally required for mitotic entry, are regulated by APC/CCdh1 . CDCA3 is an intrinsically disordered protein and contains both C-terminal KEN box and D-box recognition motifs, enabling binding to Cdh1. Our previous findings demonstrate that CDCA3 has a phosphorylation-dependent non-canonical ABBA-like motif within the linker region bridging these two recognition motifs and is required for efficient binding to Cdh1. Here, we sought to identify and further characterize additional residues that participate within this ABBA-like motif using detailed in vitro experiments and in silico modeling studies. We identified the role of H-bonds, hydrophobic and ionic interactions across the CDCA3 ABBA-like motif in the linker region between KEN and D-box motifs. This linker region adopts a well-defined structure when bound to Cdh1 in the presence of phosphorylation. Upon alanine mutation, the structure of this region is lost, leading to higher flexibility, and alteration in affinities due to binding to alternate sites on Cdh1. Our findings identify roles for the anchoring residues in the non-canonical ABBA-like motif to promote binding to the APC/CCdh1 and regulation of CDCA3 protein levels.

Evidence for Multiple Binding Modes in the Initial Contact Between SARS-CoV-2 Spike S1 Protein and Cell Surface Glycans.

Infection of host cells by SARS-CoV-2 begins with recognition by the virus S (spike) protein of cell surface heparan sulfate (HS), tethering the virus to the extracellular matrix environment, and causing the subunit S1-RBD to undergo a conformational change into the 'open' conformation. These two events promote the binding of S1-RBD to the angiotensin converting enzyme 2 (ACE2) receptor, a preliminary step toward viral-cell membrane fusion. Combining ligand-based NMR spectroscopy with molecular dynamics, oligosaccharide analogues were used to explore the interactions between S1-RBD of SARS CoV-2 and HS, revealing several low-specificity binding modes and previously unidentified potential sites for the binding of extended HS polysaccharide chains. The evidence for multiple binding modes also suggest that highly specific inhibitors will not be optimal against protein S but, rather, diverse HS-based structures, characterized by high affinity and including multi-valent compounds, may be required.

In silico analysis of crustacean hyperglycemic hormone family G protein-coupled receptor candidates.

Ecdysteroid molting hormone synthesis is directed by a pair of molting glands or Y-organs (YOs), and this synthesis is inhibited by molt-inhibiting hormone (MIH). MIH is a member of the crustacean hyperglycemic hormone (CHH) neuropeptide superfamily, which includes CHH and insect ion transport peptide (ITP). It is hypothesized that the MIH receptor is a Class A (Rhodopsin-like) G protein-coupled receptor (GPCR). The YO of the blackback land crab, Gecarcinus lateralis, expresses 49 Class A GPCRs, three of which (Gl-CHHR-A9, -A10, and -A12) were provisionally assigned as CHH-like receptors. CrusTome, a transcriptome database assembled from 189 crustaceans and 12 ecdysozoan outgroups, was used to deorphanize candidate MIH/CHH GPCRs, relying on sequence homology to three functionally characterized ITP receptors (BNGR-A2, BNGR-A24, and BNGR-A34) in the silk moth, Bombyx mori. Phylogenetic analysis and multiple sequence alignments across major taxonomic groups revealed extensive expansion and diversification of crustacean A2, A24, and A34 receptors, designated CHH Family Receptor Candidates (CFRCs). The A2 clade was divided into three subclades; A24 clade was divided into five subclades; and A34 was divided into six subclades. The subclades were distinguished by conserved motifs in extracellular loop (ECL) 2 and ECL3 in the ligand-binding region. Eleven of the 14 subclades occurred in decapod crustaceans. In G. lateralis, seven CFRC sequences, designated Gl-CFRC-A2α1, -A24α, -A24ß1, -A24ß2, -A34α2, -A34ß1, and -A34ß2, were identified; the three A34 sequences corresponded to Gl-GPCR-A12, -A9, and A10, respectively. ECL2 in all the CFRC sequences had a two-stranded ß-sheet structure similar to human Class A GPCRs, whereas the ECL2 of decapod CFRC-A34ß1/ß2 had an additional two-stranded ß-sheet. We hypothesize that this second ß-sheet on ECL2 plays a role in MIH/CHH binding and activation, which will be investigated further with functional assays.

Dehydroepiandrosterone Sulfate and Colorectal Cancer Risk: A Mendelian Randomization Analysis.

Colorectal cancer is the third most common and second most deadly type of cancer worldwide, with approximately 1.9 million cases and 0.9 million deaths worldwide in 2020. Previous studies have shown that estrogen and testosterone hormones are associated with colorectal cancer risk and mortality. However, the potential effect of their precursor, dehydroepiandrosterone sulfate (DHEAS), on colorectal cancer risk has not been investigated. Therefore, evaluating DHEAS's effect on colorectal cancer will expand our understanding of the hormonal contribution to colorectal cancer risk. In this study, we conducted a two-sample Mendelian randomization (MR) analysis to investigate the causal effect of DHEAS on colorectal cancer. We obtained DHEAS and colorectal cancer genomewide association study (GWAS) summary statistics from the Leipzig Health Atlas and the GWAS catalog and conducted MR analyses using the TwoSampleMR R package. Our results suggest that higher DHEAS levels are causally associated with decreased colorectal cancer risk (odds ratio per unit increase in DHEAS levels z score = 0.70; 95% confidence interval [0.51, 0.96]), which is in line with previous observations in a case-control study of colon cancer. The outcome of this study will be beneficial in developing plasma DHEAS-based biomarkers in colorectal cancer. Further studies should be conducted to interpret the DHEAS-colorectal cancer association among different ancestries and populations.

Molecular basis of transport of surface functionalised gold nanoparticles to pulmonary surfactant.

Ligands like alkanethiol (e.g. dodecanethiol, hexadecanethiol, etc.) and polymers (e.g. poly(vinyl pyrrolidone), polyethylene glycol-thiol) capped to the gold nanoparticles (AuNPs) are widely used in biomedical field as drug carriers and as promising materials for probing and manipulating cellular processes. Ligand functionalised AuNPs are known to interact with the pulmonary surfactant (PS) monolayer once reaching the alveolar region. Therefore, it is crucial to understand the interaction between AuNPs and PS monolayers. Using coarse-grained molecular dynamics simulations, the effect of ligand density, and ligand length have been studied for two classes of ligands on a PS model monolayer consisting of DPPC, POPG, cholesterol and SP-B (mini-peptide). The ligands considered in this study are alkanethiol and polyethylene glycol (PEG) thiol as examples of hydrophobic and hydrophilic ligands, respectively. It was observed that the interaction between AuNPs and PS changes the biophysical properties of PS monolayer in compressed and expanded states. The AuNPs with hydrophilic ligand, can penetrate through the monolayer more easily, while the AuNPs with hydrophobic ligand are embedded in the monolayer and participated in deforming the monolayer structure particularly the monolayer in the compressed state. The bare AuNPs hinder to lower the monolayer surface tension value at the interface, however introducing ligand to the bare AuNPs or increasing the ligand length and density have an impact of lowering of monolayer surface tension to a minor extent. The simulation results guide the design of ligand protected NPs as drug carriers and can identify the nanoparticles' potential side effects on lung surfactant.

A supervised machine learning approach identifies gene-regulating factor-mediated competing endogenous RNA networks in hormone-dependent cancers.

Competing endogenous RNAs (ceRNAs) have become an emerging topic in cancer research due to their role in gene regulatory networks. To date, traditional ceRNA bioinformatic studies have investigated microRNAs as the only factor regulating gene expression. Growing evidence suggests that genomic (e.g., copy number alteration [CNA]), transcriptomic (e.g., transcription factors [TFs]), and epigenomic (e.g., DNA methylation [DM]) factors can influence ceRNA regulatory networks. Herein, we used the Least absolute shrinkage and selection operator regression, a machine learning approach, to integrate DM, CNA, and TFs data with RNA expression to infer ceRNA networks in cancer risk. The gene-regulating factors-mediated ceRNA networks were identified in four hormone-dependent (HD) cancer types: prostate, breast, colorectal, and endometrial. The shared ceRNAs across HD cancer types were further investigated using survival analysis, functional enrichment analysis, and protein-protein interaction network analysis. We found two (BUB1 and EXO1) and one (RRM2) survival-significant ceRNA(s) shared across breast-colorectal-endometrial and prostate-colorectal-endometrial combinations, respectively. Both BUB1 and BUB1B genes were identified as shared ceRNAs across more than two HD cancers of interest. These genes play a critical role in cell division, spindle-assembly checkpoint signalling, and correct chromosome alignment. Furthermore, shared ceRNAs across multiple HD cancers have been involved in essential cancer pathways such as cell cycle, p53 signalling, and chromosome segregation. Identifying ceRNAs' roles across multiple related cancers will improve our understanding of their shared disease biology. Moreover, it contributes to the knowledge of RNA-mediated cancer pathogenesis.

Synthetic Heparan Sulfate Mimetic Pixatimod (PG545) Potently Inhibits SARS-CoV-2 by Disrupting the Spike-ACE2 Interaction.

Heparan sulfate (HS) is a cell surface polysaccharide recently identified as a coreceptor with the ACE2 protein for the S1 spike protein on SARS-CoV-2 virus, providing a tractable new therapeutic target. Clinically used heparins demonstrate an inhibitory activity but have an anticoagulant activity and are supply-limited, necessitating alternative solutions. Here, we show that synthetic HS mimetic pixatimod (PG545), a cancer drug candidate, binds and destabilizes the SARS-CoV-2 spike protein receptor binding domain and directly inhibits its binding to ACE2, consistent with molecular modeling identification of multiple molecular contacts and overlapping pixatimod and ACE2 binding sites. Assays with multiple clinical isolates of SARS-CoV-2 virus show that pixatimod potently inhibits the infection of monkey Vero E6 cells and physiologically relevant human bronchial epithelial cells at safe therapeutic concentrations. Pixatimod also retained broad potency against variants of concern (VOC) including B.1.1.7 (Alpha), B.1.351 (Beta), B.1.617.2 (Delta), and B.1.1.529 (Omicron). Furthermore, in a K18-hACE2 mouse model, pixatimod significantly reduced SARS-CoV-2 viral titers in the upper respiratory tract and virus-induced weight loss. This demonstration of potent anti-SARS-CoV-2 activity tolerant to emerging mutations establishes proof-of-concept for targeting the HS-Spike protein-ACE2 axis with synthetic HS mimetics and provides a strong rationale for clinical investigation of pixatimod as a potential multimodal therapeutic for COVID-19.

Structural Insights into Pixatimod (PG545) Inhibition of Heparanase, a Key Enzyme in Cancer and Viral Infections.

Pixatimod (PG545), a heparan sulfate (HS) mimetic and anticancer agent currently in clinical trials, is a potent inhibitor of heparanase. Heparanase is an endo-ß-glucuronidase that degrades HS in the extracellular matrix and basement membranes and is implicated in numerous pathological processes such as cancer and viral infections, including SARS-CoV-2. To understand how PG545 interacts with heparanase, we firstly carried out a conformational analysis through a combination of NMR experiments and molecular modelling which showed that the reducing end ß-D-glucose residue of PG545 adopts a distorted conformation. This was followed by docking and molecular dynamics simulations to study the interactions of PG545 with heparanase, revealing that PG545 is able to block the active site by binding in different conformations, with the cholestanol side-chain making important hydrophobic interactions. While PG545 blocks its natural substrate HS from binding to the active site, small synthetic heparanase substrates are only partially excluded, and thus pentasaccharide or larger substrates are preferred for assaying this class of inhibitor. This study provides new insights for the design of next-generation heparanase inhibitors and substrates.

The Impact of Rare Human Variants on Barrier-To-Auto-Integration Factor 1 (Banf1) Structure and Function.

Barrier-to-Autointegration Factor 1 (Banf1/BAF) is a critical component of the nuclear envelope and is involved in the maintenance of chromatin structure and genome stability. Banf1 is a small DNA binding protein that is conserved amongst multicellular eukaryotes. Banf1 functions as a dimer, and binds non-specifically to the phosphate backbone of DNA, compacting the DNA in a looping process. The loss of Banf1 results in loss of nuclear envelope integrity and aberrant chromatin organisation. Significantly, mutations in Banf1 are associated with the severe premature ageing syndrome, Néstor-Guillermo Progeria Syndrome. Previously, rare human variants of Banf1 have been identified, however the impact of these variants on Banf1 function has not been explored. Here, using in silico modelling, biophysical and cell-based approaches, we investigate the effect of rare human variants on Banf1 structure and function. We show that these variants do not significantly alter the secondary structure of Banf1, but several single amino acid variants in the N- and C-terminus of Banf1 impact upon the DNA binding ability of Banf1, without altering Banf1 localisation or nuclear integrity. The functional characterisation of these variants provides further insight into Banf1 structure and function and may aid future studies examining the potential impact of Banf1 function on nuclear structure and human health.