Physics-driven structural docking and protein language models accelerate antibody screening and design for broad-spectrum antiviral therapy.

Therapeutic antibodies have become one of the most influential therapeutics in modern medicine to fight against infectious pathogens, cancer, and many other diseases. However, experimental screening for highly efficacious targeting antibodies is labor-intensive and of high cost, which is exacerbated by evolving antigen targets under selective pressure such as fast-mutating viral variants. As a proof-of-concept, we developed a machine learning-assisted antibody generation pipeline that greatly accelerates the screening and re-design of immunoglobulins G (IgGs) against a broad spectrum of SARS-CoV-2 coronavirus variant strains. These viruses infect human host cells via the viral spike protein binding to the host cell receptor angiotensin-converting enzyme 2 (ACE2). Using over 1300 IgG sequences derived from convalescent patient B cells that bind with spike's receptor binding domain (RBD), we first established protein structural docking models in assessing the RBD-IgG-ACE2 interaction interfaces and predicting the virus-neutralizing activity of each IgG with a confidence score. Additionally, employing Gaussian process regression (also known as Kriging) in a latent space of an antibody language model, we predicted the landscape of IgGs' activity profiles against individual coronaviral variants of concern. With functional analyses and experimental validations, we efficiently prioritized IgG candidates for neutralizing a broad spectrum of viral variants (wildtype, Delta, and Omicron) to prevent the infection of host cells in vitro and hACE2 transgenic mice in vivo. Furthermore, the computational analyses enabled rational redesigns of selective IgG clones with single amino acid substitutions at the RBD-binding interface to improve the IgG blockade efficacy for one of the severe, therapy-resistant strains - Delta (B.1.617). Our work expedites applications of artificial intelligence in antibody screening and re-design even in low-data regimes combining protein language models and Kriging for antibody sequence analysis, activity prediction, and efficacy improvement, in synergy with physics-driven protein docking models for antibody-antigen interface structure analyses and functional optimization.

Dynamic Glycoprotein Hyposialylation Promotes Chemotherapy Evasion and Metastatic Seeding of Quiescent Circulating Tumor Cell Clusters in Breast Cancer.

Most circulating tumor cells (CTC) are detected as single cells, whereas a small proportion of CTCs in multicellular clusters with stemness properties possess 20- to 100-times higher metastatic propensity than the single cells. Here we report that CTC dynamics in both singles and clusters in response to therapies predict overall survival for breast cancer. Chemotherapy-evasive CTC clusters are relatively quiescent with a specific loss of ST6GAL1-catalyzed α2,6-sialylation in glycoproteins. Dynamic hyposialylation in CTCs or deficiency of ST6GAL1 promotes cluster formation for metastatic seeding and enables cellular quiescence to evade paclitaxel treatment in breast cancer. Glycoproteomic analysis reveals newly identified protein substrates of ST6GAL1, such as adhesion or stemness markers PODXL, ICAM1, ECE1, ALCAM1, CD97, and CD44, contributing to CTC clustering (aggregation) and metastatic seeding. As a proof of concept, neutralizing antibodies against one newly identified contributor, PODXL, inhibit CTC cluster formation and lung metastasis associated with paclitaxel treatment for triple-negative breast cancer. SIGNIFICANCE: This study discovers that dynamic loss of terminal sialylation in glycoproteins of CTC clusters contributes to the fate of cellular dormancy, advantageous evasion to chemotherapy, and enhanced metastatic seeding. It identifies PODXL as a glycoprotein substrate of ST6GAL1 and a candidate target to counter chemoevasion-associated metastasis of quiescent tumor cells. This article is featured in Selected Articles from This Issue, p. 1949.

Unique molecular signatures sustained in circulating monocytes and regulatory T cells in convalescent COVID-19 patients.

Over two years into the COVID-19 pandemic, the human immune response to SARS-CoV-2 during the active disease phase has been extensively studied. However, the long-term impact after recovery, which is critical to advance our understanding SARS-CoV-2 and COVID-19-associated long-term complications, remains largely unknown. Herein, we characterized single-cell profiles of circulating immune cells in the peripheral blood of 100 patients, including convalescent COVID-19 and sero-negative controls. Flow cytometry analyses revealed reduced frequencies of both short-lived monocytes and long-lived regulatory T (Treg) cells within the patients who have recovered from severe COVID-19. sc-RNA seq analysis identifies seven heterogeneous clusters of monocytes and nine Treg clusters featuring distinct molecular signatures in association with COVID-19 severity. Asymptomatic patients contain the most abundant clusters of monocytes and Tregs expressing high CD74 or IFN-responsive genes. In contrast, the patients recovered from a severe disease have shown two dominant inflammatory monocyte clusters featuring S100 family genes: one monocyte cluster of S100A8 & A9 coupled with high HLA-I and another cluster of S100A4 & A6 with high HLA-II genes, a specific non-classical monocyte cluster with distinct IFITM family genes, as well as a unique TGF-ß high Treg Cluster. The outpatients and seronegative controls share most of the monocyte and Treg clusters patterns with high expression of HLA genes. Surprisingly, while presumably short-lived monocytes appear to have sustained alterations over 4 months, the decreased frequencies of long-lived Tregs (high HLA-DRA and S100A6) in the outpatients restore over the tested convalescent time (≥ 4 months). Collectively, our study identifies sustained and dynamically altered monocytes and Treg clusters with distinct molecular signatures after recovery, associated with COVID-19 severity.

Computational ranking-assisted identification of Plexin-B2 in homotypic and heterotypic clustering of circulating tumor cells in breast cancer metastasis.

Metastasis is the cause of over 90% of all deaths associated with breast cancer, yet the strategies to predict cancer spreading based on primary tumor profiles and therefore prevent metastasis are egregiously limited. As rare precursor cells to metastasis, circulating tumor cells (CTCs) in multicellular clusters in the blood are 20-50 times more likely to produce viable metastasis than single CTCs. However, the molecular mechanisms underlying various CTC clusters, such as homotypic tumor cell clusters and heterotypic tumor-immune cell clusters, are yet to be fully elucidated. Combining machine learning-assisted computational ranking with experimental demonstration to assess cell adhesion candidates, we identified a transmembrane protein Plexin- B2 (PB2) as a new therapeutic target that drives the formation of both homotypic and heterotypic CTC clusters. High PB2 expression in human primary tumors predicts an unfavorable distant metastasis-free survival and is enriched in CTC clusters compared to single CTCs in advanced breast cancers. Loss of PB2 reduces formation of homotypic tumor cell clusters as well as heterotypic tumor-myeloid cell clusters in triple-negative breast cancer. Interactions between PB2 and its ligand Sema4C on tumor cells promote homotypic cluster formation, and PB2 binding with Sema4A on myeloid cells (monocytes) drives heterotypic CTC cluster formation, suggesting that metastasizing tumor cells hijack the PB2/Sema family axis to promote lung metastasis in breast cancer. Additionally, using a global proteomic analysis, we identified novel downstream effectors of the PB2 pathway associated with cancer stemness, cell cycling, and tumor cell clustering in breast cancer. Thus, PB2 is a novel therapeutic target for preventing new metastasis.

Machine learning-assisted elucidation of CD81-CD44 interactions in promoting cancer stemness and extracellular vesicle integrity.

Tumor-initiating cells with reprogramming plasticity or stem-progenitor cell properties (stemness) are thought to be essential for cancer development and metastatic regeneration in many cancers; however, elucidation of the underlying molecular network and pathways remains demanding. Combining machine learning and experimental investigation, here we report CD81, a tetraspanin transmembrane protein known to be enriched in extracellular vesicles (EVs), as a newly identified driver of breast cancer stemness and metastasis. Using protein structure modeling and interface prediction-guided mutagenesis, we demonstrate that membrane CD81 interacts with CD44 through their extracellular regions in promoting tumor cell cluster formation and lung metastasis of triple negative breast cancer (TNBC) in human and mouse models. In-depth global and phosphoproteomic analyses of tumor cells deficient with CD81 or CD44 unveils endocytosis-related pathway alterations, leading to further identification of a quality-keeping role of CD44 and CD81 in EV secretion as well as in EV-associated stemness-promoting function. CD81 is coexpressed along with CD44 in human circulating tumor cells (CTCs) and enriched in clustered CTCs that promote cancer stemness and metastasis, supporting the clinical significance of CD81 in association with patient outcomes. Our study highlights machine learning as a powerful tool in facilitating the molecular understanding of new molecular targets in regulating stemness and metastasis of TNBC.

Unique molecular signatures sustained in circulating monocytes and regulatory T cells in Convalescent COVID-19 patients.

Over two years into the COVID-19 pandemic, the human immune response to SARS-CoV-2 during the active disease phase has been extensively studied. However, the long-term impact after recovery, which is critical to advance our understanding SARS-CoV-2 and COVID-19-associated long-term complications, remains largely unknown. Herein, we characterized multi-omic single-cell profiles of circulating immune cells in the peripheral blood of 100 patients, including covenlesent COVID-19 and sero-negative controls. The reduced frequencies of both short-lived monocytes and long-lived regulatory T (Treg) cells are significantly associated with the patients recovered from severe COVID-19. Consistently, sc-RNA seq analysis reveals seven heterogeneous clusters of monocytes (M0-M6) and ten Treg clusters (T0-T9) featuring distinct molecular signatures and associated with COVID-19 severity. Asymptomatic patients contain the most abundant clusters of monocyte and Treg expressing high CD74 or IFN-responsive genes. In contrast, the patients recovered from a severe disease have shown two dominant inflammatory monocyte clusters with S100 family genes: S100A8 & A9 with high HLA-I whereas S100A4 & A6 with high HLA-II genes, a specific non-classical monocyte cluster with distinct IFITM family genes, and a unique TGF-ß high Treg Cluster. The outpatients and seronegative controls share most of the monocyte and Treg clusters patterns with high expression of HLA genes. Surprisingly, while presumably short-ived monocytes appear to have sustained alterations over 4 months, the decreased frequencies of long-lived Tregs (high HLA-DRA and S100A6) in the outpatients restore over the tested convalescent time (>= 4 months). Collectively, our study identifies sustained and dynamically altered monocytes and Treg clusters with distinct molecular signatures after recovery, associated with COVID-19 severity.

Circulating ACE2-expressing extracellular vesicles block broad strains of SARS-CoV-2.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused the pandemic of the coronavirus induced disease 2019 (COVID-19) with evolving variants of concern. It remains urgent to identify novel approaches against broad strains of SARS-CoV-2, which infect host cells via the entry receptor angiotensin-converting enzyme 2 (ACE2). Herein, we report an increase in circulating extracellular vesicles (EVs) that express ACE2 (evACE2) in plasma of COVID-19 patients, which levels are associated with severe pathogenesis. Importantly, evACE2 isolated from human plasma or cells neutralizes SARS-CoV-2 infection by competing with cellular ACE2. Compared to vesicle-free recombinant human ACE2 (rhACE2), evACE2 shows a 135-fold higher potency in blocking the binding of the viral spike protein RBD, and a 60- to 80-fold higher efficacy in preventing infections by both pseudotyped and authentic SARS-CoV-2. Consistently, evACE2 protects the hACE2 transgenic mice from SARS-CoV-2-induced lung injury and mortality. Furthermore, evACE2 inhibits the infection of SARS-CoV-2 variants (α, ß, and δ) with equal or higher potency than for the wildtype strain, supporting a broad-spectrum antiviral mechanism of evACE2 for therapeutic development to block the infection of existing and future coronaviruses that use the ACE2 receptor.

ICAM1 initiates CTC cluster formation and trans-endothelial migration in lung metastasis of breast cancer.

Circulating tumor cell (CTC) clusters mediate metastasis at a higher efficiency and are associated with lower overall survival in breast cancer compared to single cells. Combining single-cell RNA sequencing and protein analyses, here we report the profiles of primary tumor cells and lung metastases of triple-negative breast cancer (TNBC). ICAM1 expression increases by 200-fold in the lung metastases of three TNBC patient-derived xenografts (PDXs). Depletion of ICAM1 abrogates lung colonization of TNBC cells by inhibiting homotypic tumor cell-tumor cell cluster formation. Machine learning-based algorithms and mutagenesis analyses identify ICAM1 regions responsible for homophilic ICAM1-ICAM1 interactions, thereby directing homotypic tumor cell clustering, as well as heterotypic tumor-endothelial adhesion for trans-endothelial migration. Moreover, ICAM1 promotes metastasis by activating cellular pathways related to cell cycle and stemness. Finally, blocking ICAM1 interactions significantly inhibits CTC cluster formation, tumor cell transendothelial migration, and lung metastasis. Therefore, ICAM1 can serve as a novel therapeutic target for metastasis initiation of TNBC.

ITGA2 promotes expression of ACLY and CCND1 in enhancing breast cancer stemness and metastasis.

Cancer metastasis is largely incurable and accounts for 90% of breast cancer deaths, especially for the aggressive basal-like or triple negative breast cancer (TNBC). Combining patient database analyses and functional studies, we examined the association of integrin family members with clinical outcomes as well as their connection with previously identified microRNA regulators of metastasis, such as miR-206 that inhibits stemness and metastasis of TNBC. Here we report that the integrin receptor CD49b-encoding ITGA2, a direct target of miR-206, promotes breast cancer stemness and metastasis. ITGA2 knockdown suppressed self-renewal related mammosphere formation and pluripotency marker expression, inhibited cell cycling, compromised migration and invasion, and therefore decreased lung metastasis of breast cancer. ITGA2 overexpression reversed miR-206-caused cell cycle arrest in G1. RNA sequencing analyses revealed that ITGA2 knockdown inhibits genes related to cell cycle regulation and lipid metabolism, including CCND1 and ACLY as representative targets, respectively. Knockdown of CCND1 or ACLY inhibits mammosphere formation of breast cancer cells. Overexpression of CCND1 rescues the phenotype of ITGA2 knockdown-induced cell cycle arrest. ACLY-encoded ATP citrate lyase is essential to maintain cellular acetyl-CoA levels. CCND1 knockdown further mimics ITGA2 knockdown in abolishing lung colonization of breast cancer cells. We identified that the low levels of miR-206 as well as high expression levels of ITGA2, ACLY and CCND1 are associated with an unfavorable relapse-free survival of the patients with estrogen receptor-negative or high grade breast cancer, especially basal-like or TNBC, possibly serving as potential biomarkers of cancer stemness and therapeutic targets of breast cancer metastasis.

EGFR inhibition blocks cancer stem cell clustering and lung metastasis of triple negative breast cancer.

Triple-negative breast cancer (TNBC) is one of the most aggressive and metastatic breast cancer subtypes lacking targeted therapy. Our recent work demonstrated that circulating tumor cell (CTC) clusters and polyclonal metastasis of TNBC are driven by aggregation of CD44+ cancer stem cells (CSC) and associated with an unfavorable prognosis, such as low overall survival. However, there is no existing therapeutic that can specifically block CTC or CSC cluster formation. Methods: Using patient-derived xenograft (PDX) models, we established an ex vivo tumor cell clustering assay for a pilot screening of blockade antibodies. After identifying EGFR as a target candidate, we modulated the gene expression and inhibited its kinase activity to determine its functional importance in tumor cell clustering and therapeutic inhibition of lung metastasis. We also examined the molecular regulation network of EGFR and a potential connection to CSC marker CD44 and microRNAs, which regulate CTC clustering. Results: We report here that EGFR inhibition successfully blocks circulating CSC (cCSC) clustering and lung metastasis of TNBC. EGFR enhances CD44-mediated tumor cell aggregation and CD44 stabilizes EGFR. Importantly, blocking EGFR by a novel anti-EGFR monoclonal antibody (clone LA1) effectively blocked cell aggregation in vitro and reduced lung metastasis in vivo. Furthermore, our data demonstrated that the tumor suppressor microRNA-30c serves as another negative regulator of cCSC clustering and lung metastasis by targeting CD44 as well as its downstream effector EGFR. Conclusion: Our studies identify a novel anti-EGFR therapeutic strategy to inhibit cCSC aggregation and therefore abolish cCSC cluster-mediated metastasis of TNBC.

Extracellular Domains I and II of cell-surface glycoprotein CD44 mediate its trans-homophilic dimerization and tumor cluster aggregation.

CD44 molecule (CD44) is a well-known surface glycoprotein on tumor-initiating cells or cancer stem cells. However, its utility as a therapeutic target for managing metastases remains to be fully evaluated. We previously demonstrated that CD44 mediates homophilic interactions for circulating tumor cell (CTC) cluster formation, which enhances cancer stemness and metastatic potential in association with an unfavorable prognosis. Furthermore, CD44 self-interactions activate the P21-activated kinase 2 (PAK2) signaling pathway. Here, we further examined the biochemical properties of CD44 in homotypic tumor cell aggregation. The standard CD44 form (CD44s) mainly assembled as intercellular homodimers (trans-dimers) in tumor clusters rather than intracellular dimers (cis-dimers) present in single cells. Machine learning-based computational modeling combined with experimental mutagenesis tests revealed that the extracellular Domains I and II of CD44 are essential for its trans-dimerization and predicted high-score residues to be required for dimerization. Substitutions of 10 these residues in Domain I (Ser-45, Glu-48, Phe-74, Cys-77, Arg-78, Tyr-79, Ile-88, Arg-90, Asn-94, and Cys-97) or 5 residues in Domain II (Ile-106, Tyr-155, Val-156, Gln-157, and Lys-158) abolished CD44 dimerization and reduced tumor cell aggregation in vitro Importantly, the substitutions in Domain II dramatically inhibited lung colonization in mice. The CD44 dimer-disrupting substitutions decreased downstream PAK2 activation without affecting the interaction between CD44 and PAK2, suggesting that PAK2 activation in tumor cell clusters is CD44 trans-dimer-dependent. These results shed critical light on the biochemical mechanisms of CD44-mediated tumor cell cluster formation and may help inform the development of therapeutic strategies to prevent tumor cluster formation and block cluster-mediated metastases.

New Advances and Challenges of Targeting Cancer Stem Cells.

The second International Cancer Stem Cell Conference in Cleveland, Ohio, on September 20-23, 2016, convened 330 attendees from academic, industrial, and clinical organizations. It featured a debate on the concepts and challenges of the cancer stem cells (CSC) as well as CSC-centered scientific sessions on clinical trials, genetics and epigenetics, tumor microenvironment, immune suppression, metastasis, therapeutic resistance, and emerging novel concepts. The conference hosted 35 renowned speakers, 100 posters, 20 short talks, and a preconference workshop. The reported advances of CSC research and therapies fostered new collaborations across national and international borders, and inspired the next generation's young scientists. Cancer Res; 77(19); 5222-7. ©2017 AACR.

Crosstalk between tumor suppressors p53 and PKCδ: Execution of the intrinsic apoptotic pathways.

p53 and PKCδ are tumor suppressors that execute apoptotic mechanisms in response to various cellular stresses. p53 is a transcription factor that is frequently mutated in human cancers; it regulates apoptosis in transcription-dependent and -independent ways in response to genotoxic stresses. PKCδ is a serine/threonine protein kinase and mutated in human cancers. Available evidence shows that PKCδ activates p53 by direct and/or indirect mechanisms. Moreover, PKCδ is also implicated in the transcriptional regulation of p53 in response to DNA damage. Recent findings demonstrated that p53, in turn, binds onto the PKCδ promoter and induces its expression upon DNA damage to facilitate apoptosis. Both p53 and PKCδ are associated with the apoptotic mechanisms in the mitochondria by regulating Bcl-2 family proteins to provide mitochondrial outer membrane permeabilization. This review discusses the crosstalk between p53 and PKCδ in the context of apoptotic cell death and cancer therapy.

Transcriptional induction of protein kinase C delta by p53 tumor suppressor in the apoptotic response to DNA damage.

Genetic alterations and aberrant gene expression trigger malignant tumors. Tumor suppressor p53 is the most altered gene in human cancers. p53 induces apoptosis by promoting pro-apoptotic genes in response to DNA damage. Protein kinase C delta (PKCδ) also induces apoptosis via various mechanisms including modification of p53. The PKCδ-dependent apoptotic mechanism has been extensively studied; however, the transcriptional regulation of PKCδ remains obscure. The current study demonstrates the transcriptional regulation of PKCδ by p53 upon genotoxic stress. The p53-binding site in the promoter region of PKCδ was detected by the ChIP-sequencing assay. Notably, the expression of PKCδ was increased upon DNA damage, which is required for the stabilization of p53. More importantly, targeting single guide RNA-driven dead Cas9 to the p53-binding site of PKCδ disturbed p53-promoted PKCδ expression and suppressed apoptosis following DNA damage. Thus, our findings suggest that the transcriptional regulation of PKCδ is controlled by p53 in a positive feedback mechanism to induce apoptosis in response to DNA damage.

Cell death decision by p53 via control of the mitochondrial membrane.

The tumor suppressor p53 is mutated in more than half of human cancers. Recent evidence has revealed that p53 not only regulates apoptosis but also regulates necrotic/necroptotic cell death via the mitochondria. The regulation of apoptosis by p53 is tightly connected to the mitochondrial outer membrane permeabilization and the induction of and interaction with Bcl-2 family members. Interestingly, p53-mediated regulation of necrosis/necroptosis is correlated with mitochondrial permeabilization pore opening via interactions with CypD and Drp1. This review discusses the p53-regulating molecules that induce apoptosis or necrosis/necroptosis via the mitochondria.

Tumor suppressor p53 induces miR-1915 processing to inhibit Bcl-2 in the apoptotic response to DNA damage.

The antiapoptotic protein Bcl-2 is overexpressed in human cancers, and confers resistance to antitumor agents in cancer cells. Bcl-2 is negatively regulated by the tumor suppressor p53 in response to DNA damage during apoptotic cell death. However, this molecular mechanism remains unclear. The available evidence indicates that miR-1915 represses Bcl-2 expression at the post-transcriptional level in human colorectal carcinoma cells, which is correlated with drug resistance. Here, we show that p53 controls miR-1915 expression in response to DNA damage. Induction of p53 affects the expression of precursor and mature, but not primary, miR-1915. Inhibition of miR-1915 abrogates downregulation of Bcl-2 expression following treatment with genotoxin. These findings demonstrate that p53 negatively regulates Bcl-2 expression by targeting miR-1915 processing from primary into precursor miRNA. Taken together, the findings of the current study reveal a novel mechanism whereby p53 negatively modulates Bcl-2 by controlling miR-1915.

Programmed cell death 6, a novel p53-responsive gene, targets to the nucleus in the apoptotic response to DNA damage.

The cellular response to genotoxic stress is multifaceted in nature. Following DNA damage, the tumor suppressor gene p53 activates and plays critical roles in cell cycle arrest, activation of DNA repair and in the event of irreparable damage, induction of apoptosis. The breakdown of apoptosis causes the accumulation of mutant cells. The elucidation of the mechanism for the p53-dependent apoptosis will be crucial in applying the strategy for cancer patients. However, the mechanism of p53-dependent apoptosis remains largely unclear. Here, we carried out ChIP followed by massively parallel DNA sequencing assay (ChIP-seq) to uncover mechanisms of apoptosis. Using ChIP-seq, we identified PDCD6 as a novel p53-responsive gene. We determined putative p53-binding sites that are important for p53 regulation in response to DNA damage in the promoter region of PDCD6. Knockdown of PDCD6 suppressed p53-dependent apoptosis. We also observed that cytochrome c release and the cleavage of PARP by caspase-3 were suppressed by depletion of PDCD6. We further observed that PDCD6 localizes in the nucleus in response to DNA damage. We identified the nuclear localization signal of PDCD6 and, importantly, the nuclear accumulation of PDCD6 significantly induced apoptosis after genotoxic stress. Therefore, we conclude that a novel p53-responsive gene PDCD6 is accumulated in the nucleus and induces apoptosis in response to DNA damage.