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.

AI is a viable alternative to high throughput screening: a 318-target study

High throughput screening (HTS) is routinely used to identify bioactive small molecules. This requires physical compounds, which limits coverage of accessible chemical space. Computational approaches combined with vast on-demand chemical libraries can access far greater chemical space, provided that the predictive accuracy is sufficient to identify useful molecules. Through the largest and most diverse virtual HTS campaign reported to date, comprising 318 individual projects, we demonstrate that our AtomNet® convolutional neural network successfully finds novel hits across every major therapeutic area and protein class. We address historical limitations of computational screening by demonstrating success for target proteins without known binders, high-quality X-ray crystal structures, or manual cherry-picking of compounds. We show that the molecules selected by the AtomNet® model are novel drug-like scaffolds rather than minor modifications to known bioactive compounds. Our empirical results suggest that computational methods can substantially replace HTS as the first step of small-molecule drug discovery.

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.

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.

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.

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.

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.

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.

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.

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/C^Cdh1 . 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/C^Cdh1 and regulation of CDCA3 protein levels.

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.

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.

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.

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.

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.

Molecular-level observations of the behavior of ligand functionalised gold nanoparticles with a lipid monolayers.

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.

Identifying Complex lncRNA/Pseudogene-miRNA-mRNA Crosstalk in Hormone-Dependent Cancers

Simple Summary

Competing endogenous RNAs (ceRNAs) have gained attention in cancer research owing to their involvement in microRNA-mediated gene regulation. Here, we identified a shared ceRNA network across five hormone-dependent (HD) cancers (prostate, breast, colon, rectal, and endometrial), that contain two long non-coding RNAs, nine mRNAs, and seventy-four microRNAs. Among them, two mRNAs and forty-one microRNAs were associated with at least one HD cancer survival. A similar analytical approach can be applied to identify shared ceRNAs across a group of related cancers, which will significantly contribute to understanding their shared disease biology.

Abstract

The discovery of microRNAs (miRNAs) has fundamentally transformed our understanding of gene regulation. The competing endogenous RNA (ceRNA) hypothesis postulates that messenger RNAs and other RNA transcripts, such as long non-coding RNAs and pseudogenes, can act as natural miRNA sponges. These RNAs influence each other’s expression levels by competing for the same pool of miRNAs through miRNA response elements on their target transcripts, thereby modulating gene expression and protein activity. In recent years, these ceRNA regulatory networks have gained considerable attention in cancer research. Several studies have identified cancer-specific ceRNA networks. Nevertheless, prior bioinformatic analyses have focused on long-non-coding RNA-associated ceRNA networks. Here, we identify an extended ceRNA network (including both long non-coding RNAs and pseudogenes) shared across a group of five hormone-dependent (HD) cancers, i.e., prostate, breast, colon, rectal, and endometrial cancers, using data from The Cancer Genome Atlas (TCGA). We performed a functional enrichment analysis for differentially expressed genes in the shared ceRNA network of HD cancers, followed by a survival analysis to determine their prognostic ability. We identified two long non-coding RNAs, nine genes, and seventy-four miRNAs in the shared ceRNA network across five HD cancers. Among them, two genes and forty-one miRNAs were associated with at least one HD cancer survival. This study is the first to investigate pseudogene-associated ceRNAs across a group of related cancers and highlights the value of this approach to understanding the shared molecular pathogenesis in a group of related diseases.

Integrative Transcriptome-Wide Analyses Uncover Novel Risk-Associated MicroRNAs in Hormone-Dependent Cancers

Background

Hormone-dependent cancers (HDC) are among the leading causes of death worldwide among both men and women. Some of the established risk factors of HDC include unhealthy lifestyles, environmental factors, and genetic influences. Numerous studies have been conducted to understand gene-cancer associations. Transcriptome-wide association studies (TWAS) integrate data from genome-wide association studies (GWAS) and gene expression (expression quantitative trait loci - eQTL) to yield meaningful information on biological pathways associated with complex traits/diseases. Recently, TWAS have enabled the identification of novel associations between HDC risk and protein-coding genes.

Methods

In the present study, we performed a TWAS analysis using the summary data-based Mendelian randomization (SMR)-heterogeneity in dependent instruments (HEIDI) method to identify microRNAs (miRNAs), a group of non-coding RNAs (ncRNAs) associated with HDC risk. We obtained eQTL and GWAS summary statistics from the ncRNA-eQTL database and the National Human Genome Research Institute-European Bioinformatics Institute (NHGRI-EBI) GWAS Catalog.

Results

We identified 13 TWAS-significant miRNAs at cis regions (±1 Mb) associated with HDC risk (two, five, one, two, and three miRNAs for prostate, breast, ovarian, colorectal, and endometrial cancers, respectively). Among them, eight novel miRNAs were recognized in HDC risk. Eight protein-coding genes targeted by TWAS-identified miRNAs (SIRT1, SOX4, RUNX2, FOXA1, ABL2, SUB1, HNRNPH1, and WAC) are associated with HDC functions and signaling pathways.

Conclusion

Overall, identifying risk-associated miRNAs across a group of related cancers may help to understand cancer biology and provide novel insights into cancer genetic mechanisms. This customized approach can be applied to identify significant miRNAs in any trait/disease of interest.

From Cancer to COVID-19: A Perspective on Targeting Heparan Sulfate-Protein Interactions

Heparan sulfate (HS) is a complex, polyanionic polysaccharide ubiquitously expressed on cell surfaces and in the extracellular matrix. HS interacts with numerous proteins to mediate a vast array of biological and pathological processes. Inhibition of HS‐protein interactions is thus an attractive approach for new therapeutic development for cancer and infectious diseases, including COVID‐19; however, synthesis of well‐defined native HS oligosaccharides remains challenging. This has aroused significant interest in the development of HS mimetics which are more synthetically tractable and have fewer side effects, such as undesired anticoagulant activity. This account provides a perspective on the design and synthesis of different classes of HS mimetics with useful properties, and the development of various assays and molecular modelling tools to progress our understanding of their interactions with HS‐binding proteins.

Phosphorylation and O-GlcNAcylation of the PHF-1 Epitope of Tau Protein Induce Local Conformational Changes of the C-Terminus and Modulate Tau Self-Assembly Into Fibrillar Aggregates

Phosphorylation of the neuronal microtubule-associated Tau protein plays a critical role in the aggregation process leading to the formation of insoluble intraneuronal fibrils within Alzheimer's disease (AD) brains. In recent years, other posttranslational modifications (PTMs) have been highlighted in the regulation of Tau (dys)functions. Among these PTMs, the O-β-linked N-acetylglucosaminylation (O-GlcNAcylation) modulates Tau phosphorylation and aggregation. We here focus on the role of the PHF-1 phospho-epitope of Tau C-terminal domain that is hyperphosphorylated in AD (at pS396/pS404) and encompasses S400 as the major O-GlcNAc site of Tau while two additional O-GlcNAc sites were found in the extreme C-terminus at S412 and S413. Using high resolution NMR spectroscopy, we showed that the O-GlcNAc glycosylation reduces phosphorylation of PHF-1 epitope by GSK3β alone or after priming by CDK2/cyclin A. Furthermore, investigations of the impact of PTMs on local conformation performed in small peptides highlight the role of S404 phosphorylation in inducing helical propensity in the region downstream pS404 that is exacerbated by other phosphorylations of PHF-1 epitope at S396 and S400, or O-GlcNAcylation of S400. Finally, the role of phosphorylation and O-GlcNAcylation of PHF-1 epitope was probed in in-vitro fibrillization assays in which O-GlcNAcylation slows down the rate of fibrillar assembly while GSK3β phosphorylation stimulates aggregation counteracting the effect of glycosylation.

Evidence of a putative glycosaminoglycan binding site on the glycosylated SARS-CoV-2 spike protein N-terminal domain

SARS-CoV-2 has rapidly spread throughout the world's population since its initial discovery in 2019. The virus infects cells via a glycosylated spike protein located on its surface. The protein primarily binds to the angiotensin-converting enzyme-2 (ACE2) receptor, using glycosaminoglycans (GAGs) as co-receptors. Here, we performed bioinformatics and molecular dynamics simulations of the spike protein to investigate the existence of additional GAG binding sites on the receptor-binding domain (RBD), separate from previously reported heparin-binding sites. A putative GAG binding site in the N-terminal domain (NTD) of the protein was identified, encompassing residues 245-246. We hypothesized that GAGs of a sufficient length might bridge the gap between this site and the PRRARS furin cleavage site, including the mutation S247R. Docking studies using GlycoTorch Vina and subsequent MD simulations of the spike trimer in the presence of dodecasaccharides of the GAGs heparin and heparan sulfate supported this possibility. The heparan sulfate chain bridged the gap, binding the furin cleavage site and S247R. In contrast, the heparin chain bound the furin cleavage site and surrounding glycosylation structures, but not S247R. These findings identify a site in the spike protein that favors heparan sulfate binding that may be particularly pertinent for a better understanding of the recent UK and South African strains. This will also assist in future targeted therapy programs that could include repurposing clinical heparan sulfate mimetics.

Identification of Proteins Deregulated by Platinum-Based Chemotherapy as Novel Biomarkers and Therapeutic Targets in Non-Small Cell Lung Cancer

Platinum-based chemotherapy remains the cornerstone of treatment for most people with non-small cell lung cancer (NSCLC), either as adjuvant therapy in combination with a second cytotoxic agent or in combination with immunotherapy. Resistance to therapy, either in the form of primary refractory disease or evolutionary resistance, remains a significant issue in the treatment of NSCLC. Hence, predictive biomarkers and novel combinational strategies are required to improve the effectiveness and durability of treatment response 6for people with NSCLC. The aim of this study was to identify novel biomarkers and/or druggable proteins from deregulated protein networks within non-oncogene driven disease that are involved in the cellular response to cisplatin. Following exposure of NSCLC cells to cisplatin, in vitro quantitative mass spectrometry was applied to identify altered protein response networks. A total of 65 proteins were significantly deregulated following cisplatin exposure. These proteins were assessed to determine if they are druggable targets using novel machine learning approaches and to identify whether these proteins might serve as prognosticators of platinum therapy. Our data demonstrate novel candidates and drug-like molecules warranting further investigation to improve response to platinum agents in NSCLC.