extfieldsofstudy: Chemistry Medicine Biology

text:
Maraviroc inhibits SARS-CoV-2 multiplication and s-protein mediated cell fusion in cell culture

In an effort to identify therapeutic intervention strategies for the treatment of COVID-19, we have investigated a selection of FDA-approved small molecules and biologics that are commonly used to treat other human diseases. A investigation into 18 small molecules and 3 biologics was conducted in cell culture and the impact of treatment on viral titer was quantified by plaque assay. The investigation identified 4 FDA-approved small molecules, Maraviroc, FTY720 (Fingolimod), Atorvastatin and Nitazoxanide that were able to inhibit SARS-CoV-2 infection. Confocal microscopy with over expressed S-protein demonstrated that Maraviroc reduced the extent of S-protein mediated cell fusion as observed by fewer multinucleate cells in the context of drug-treatment. Mathematical modeling of drug-dependent viral multiplication dynamics revealed that prolonged drug treatment will exert an exponential decrease in viral load in a multicellular/tissue environment. Taken together, the data demonstrate that Maraviroc, Fingolimod, Atorvastatin and Nitazoxanide inhibit SARS-CoV-2 in cell culture.

Introduction
Severe Acute Respiratory Syndrome coronavirus 2 (SARS-CoV-2) is a novel coronavirus which is the causative agent of the coronavirus disease 2019 (COVID-19) pandemic (Zhou et al., 2020). As of November 2021, there are 254,036,628 confirmed cases of COVID-19 resulting in 5,111,176 deaths worldside (Johns Hopkins Coronavirus Resource Center, ). Coronaviruses (CoV) are generally the cause of common colds resulting in mild respiratory infections (Weiss and Leibowitz, 2011).The emergence of pathogenic coronavirus (CoV) strains with high mortality and morbidity has occurred multiple times over the past two decades, (De Wit et al., 2016). Severe acute respiratory virus (SARS) emerged in 2002 (De Wit et al., 2016;Stadler et al., 2003;Zhong and Wong, 2004). Several SARS-like coronaviruses have been identified in bat species (Banerjee et al., 2019;Ge et al., 2013;Menachery et al., 2015;Wang et al.,2006). Middle Eastern respiratory virus was detected in 2012 (De Wit et al., 2016;Mackay and Arden, 2015). In December of 2019, SARS-CoV-2 began circulating in Wuhan, China resulting in the current COVID-19 pandemic. SARS-CoV-2 is known to be transmitted from human-to-human with inhalational acquisition being the primary mode of transfer (Liu and Zhang, 2020;Nishiura et al., 2020;Song et al., 2020).
Coronaviruses belong to the superfamily Nidovirales and subfamily Orthocoronavirinae (Weiss and Leibowitz, 2011). CoVs are enveloped, nonsegmented, positive-sense single-stranded RNA viruses (Weiss and Leibowitz, 2011). They are classified in four sub-groups: alpha, beta, gamma and delta. The seven viruses that are known to infect humans belong to the alpha and beta subgroups. HCoV-229E and HCoV-NL62 are classified as alpha while HCoV-OC43, HCoV-HKU1, MERS-CoV, SARS-CoV, and SARS-CoV-2 are beta CoVs. The virus measures approximately 65-125nm in diameter, and the viral genome measures approximately 29.9 Kb (Astuti and Ysrafil, 2020). The viral genome encompasses 14 open reading frames (ORFs) that encode both structural and nonstructural viral proteins. Among the structural proteins, S-protein has achieved significant attention due to the critical role it plays in the interaction of the virus with the ACE2 receptor on host cells (Ahmed et al., 2020). In addition to ACE2, the type II transmembrane serine protease (TMPRSS2) is also required for SARS-CoV-2 entry into cells, thus making these two membrane-associated proteins as the primary determinants for viral entry. ACE2 receptor expression can be detected in various organs besides in the lungs, including heart, kidneys and the gastrointestinal tract. Notably, COVID-19 is characterized by disease manifestations that impact all of the ACE2 positive organs and tissues. The critical aspects of interactions between the viral spike protein and the host membrane proteins ACE2 and TMPRSS2 have led to several therapeutic candidates that interfere with this virus-host protein interactions.
Infection by coronaviruses such as infectious bronchitis virus (IBV) is known to result in cellcell fusion and formation of large, multinucleated cells referred to as syncytia (Fehr and Perlman, 2015;Sisk et al., 2018). Newly synthesized S-protein either in the context of coronavirus infected cells or cells that over express S-protein, is said to accumulate on the plasma membrane (Lontok et al., 2004).Such S-protein enriched sections of plasma membranes can fuse resulting in cell-cell fusion. Inhibition of IBV infected cells with Abl kinase inhibitors (Imatinib) resulted in reduced syncytia in addition to decreasing viral load (Sisk et al., 2018).This inhibition of S-protein mediated cell fusion by Imatinib could be achieved even in the absence of other viral proteins suggesting that the cell-cell fusion event in coronavirus infected cells is likely to be dependent only on S-protein function. The cell-cell fusion event mediated by coronavirus S-protein is suggested to be controlled by different host enzymatic components than those that influence virus-host membrane fusion. S-protein dependent cell-cell fusion was shown to be independent of cathepsin L which was essential for virus-cell fusion. A novel leupeptin-sensitive host cell protease activated S-protein dependent cell-cell fusion in target cells expressing high levels of ACE2 in the context of SARS-CoV-1 (Simmons et al., 2011).This mechanism of S-protein mediated cell-cell fusion was implicated in viral spread in the context of SARS-CoV-1 infection and the ability of the virus to evade host humoral immune responses, thus posing an important challenge for antibody-mediated viral control . Increased incidence of S-protein mediated syncytia were observed in the context of SARS-CoV-2 S-protein, as compared to SARS-CoV-1, thus highlighting the need to address S-protein mediated cell-cell fusion in SARS-CoV-2 to control viral spread in the infected host (Xia et al., 2020). A significant contributor towards the morbidity and mortality associated with COVID-19 is the host inflammatory response, with several pro-inflammatory cytokines known to be involved in the tissue damage sustained due to infection. Cytokine storm is suggested to be a significant cause of organ failure and Acute Respiratory Distress Syndrome (ARDS) (Ye et al., 2020). Several host signaling events including the NFkB pathway, JAK-STAT pathway and IFN pathway cumulatively contribute to the proinflammatory environment in target tissues such as the lungs. Onset of ARDS during the later stages of COVID-19 is associated with poor prognosis and is a direct outcome of inflammatory lung damage (Fanelli et al., 2013). Additional complications that are unique to COVID-19 include microclots in blood vessels that supply blood to the heart, lungs and the brain (Ciceri et al., 2020). This disseminated coagulation contributes to death resulting from strokes and heart complications in addition to ARDS (Bansal, 2020;Trejo-Gabriel-Galán, 2020).
It is becoming increasingly apparent that COVID-19 is a combined outcome of viral multiplication in multiple target sites in the infected individual and various host-based complications resulting from dysregulated thrombotic and inflammatory events. An effective therapeutic regimen will need to address the viral and host-based aspects of the disease in order to be able to decrease the health burden and increase survival of affected individuals. Furthermore, such an approach may be essential to address any potential long-term complications that may arise in survivors, especially those who are considered to be high-risk populations. High-risk individuals include those with preexisting health conditions such as coronary heart disease, and inflammatory states including asthma and diabetes. Seeking a therapeutic regimen that may be easily adapted to such high-risk individuals, we focused our attention on a subset of FDA-approved small molecules and biologics that are commonly used to treat heart disease, inflammation and infectious disease. We screened a total of 22 such molecules and have identified four small molecules (Maraviroc, FTY 720, Atorvastatin and Nitazoxanide) and three antibodies (Bevacizumab, Eculizumab and Anti-CD126/IL6R), all of which were able to decrease infectious virus titers in cell culture. Additional studies performed with the prioritized small molecule inhibitors demonstrated that these compounds moderately impacted viral entry and interfered with the interaction between the spike protein and ACE2 to varying degrees. Cell biology studies demonstrate that Maraviroc was able to decrease the extent of S-protein-mediated cell fusion. Our efforts to develop simulation models of the antiviral outcomes of Maraviroc on a tissue level demonstrated that sustained availability of drug will have a significant impact on decreasing viral load and cellular fusion and increasing cell survival. Simulations also show that drugs targeting viral fusion may increase efficacy of anti-replication drugs when used in combination.

Viral infection and inhibitor treatment.
Vero cells were seeded at a density of 3E4 cells/well. The inhibitors were resuspended in dimethyl sulfoxide (DMSO) prior to use and stored as recommended by the manufacturer as stock solutions (5 mM -100 mM, varied depending dry weight and molecular weight of inhibitor received). Inhibitors were diluted to stated concentrations in cell culture media. Cells were pretreated with the target inhibitors before infection for 2h. Virus was diluted in media to appropriate MOI and placed on cells. The plate was incubated for 1 hour at 37°C, 5% CO2 to allow the uptake of virus. Viral inoculum was removed, cells were gently washed three times with Dulbecco's phosphate-buffered saline (DPBS) without calcium and magnesium (Gibco 14190144) and conditioned media with the inhibitors was added back onto the cells. Cells were incubated at 37°C and supernatants were collected at indicated times post infection and stored at -80°C until required for further analyses.

Plaque assay.
Plaque assays were performed using Vero cells grown to a cell number of 2E5 cells/well in a 12well plate. Supernatants from infected cells were serially diluted in media and used to infect Vero cells for one hour. After infection a 1:1 overlay consisting of 2X Eagle's Minimum Essential Medium without phenol red (Quality Biological, 115-073-101), supplemented with 10% fetal bovine serum (FBS) (Gibco, 10,437,028), non-essential amino acids (Gibco, 11140-050), 1 mM sodium pyruvate (Corning, 25-000-Cl), 2 mM L-glutamine, 20 U/ml of penicillin and 20 μg/ml of streptomycin) and 0.6% agarose was added to each well. Plates were incubated at 37°C for 48 hours. Cells were fixed with 10% formaldehyde for 1 hour at room temperature. Formaldehyde was aspirated and the agarose overlay was removed. Cells were stained with crystal violet (1% CV w/v in a 20% ethanol solution). Viral titer (PFU/mL) of SARS-CoV-2 was determined by counting the number of plaques.

Cell viability assay.
Cell viability was determined using the CellTiterGlo assay (Promega, G7570). Inhibitors were resuspended in DMSO and stored appropriately. Inhibitors were diluted to required concentrations in cell culture media. Vero cells were seeded at a density of 3E4 in white-walled 96-well plates. Inhibitors were overlaid on cells and incubated at 37°C, 5% CO2 for 24 hours. CellTiter-Glo substrate was used according to the manufacturer's instructions and luminescence was determined using a DTX 880 multimode detector (Beckman Coulter) with an integration time of 100ms/well. Cell viability was determined as percent cell viability normalized to untreated control.
Small molecules inhibitors were screened for ability to disrupt SARS-CoV-2 spike proteins binding to ACE2 receptors. ACE2: Spike S1-Biotin (SARS-CoV-2) Inhibitor Screening Assay Kit (BPS Bioscience, Cat# 79945) was utilized per manufacturer's instructions. Plate was washed and blocked in between each addition step with provided wash buffer and blocking buffer. A nickel-coated plate was treated with his-tagged ACE2. Drugs were screened at a low concentration (25 μM) and a high concentration (100 μM). Inhibitor or DMSO was added to the plate to bring final dilution to indicated concentration. Plate was incubated for 1 hour at room temperature with shaking. SARS-CoV-2 spike S1-biotin was added to the plate. Plate was incubated for 1 hour at room temperature with shaking. Streptavidin-HRP is added to the plate. Plate was shaken for one hour at room temperature. HRP substrate was added to the plate and the Chemiluminescence was measured immediately using a DTX 880 multimode detector (Beckman Coulter) with an integration time of 1 second, and delay of 100ms.

rVSV Pseudotyped SARS-CoV2 Spike Neutralization Assay.
The vesicular stomatitis virus (rVSV) with deleted glycoprotein and replacement SARS-CoV-2 spike was acquired from IBT Bioservices. Vero cells were seed at a density of 6E4 in a black, 96-well tissue culture plate. Cells were maintained in EMEM supplemented with 10% fetal bovine serum (FBS) and 1x penicillin/streptomycin (P/S). Cells were incubated overnight 37°C and 5% CO2. Virus was diluted per manufacture's protocol and incubated with drug at a final concentration of 5 μM or DMSO control in unsupplemented media for one hour at room temperature. Media was removed from plate, replaced with treated or untreated virus. Plate was incubated for one hour at 37°C. EMEM supplemented with 2% FBS and 1x P/S was added to the wells and the plate was incubated overnight at 37°C and 5% CO2. Cells were lysed with 1X Passive Lysis Buffer for 30 minutes while shaking. Firefly Luciferase Assay substrate was added to the plate. Plate was read immediately using a DTX 880 multimode detector (Beckman Coulter) with an integration time of 1 second.

Transfection Methods.
Vero cells were grown on 6-well plates and were transiently transfected with Lipofectamine 2000 (Invitrogen, 11668019) based on the manufacturer's protocols. In brief, prior to treatment, Lipofectamine 2000/siRNA complexes were prepared in reduced serum medium, OptiMEM (Gibco, 31985088), at the recommended ratio. Cells were treated for 6 hours before being replaced with EMEM complete medium.

Microscopy analysis of cell-cell fusion.
African green monkey kidney (Vero) cell lines (ATCC Cat# CCL-81) were grown in Eagle's Minimum Essential Medium (EMEM) (Corning) supplemented with 10% heat-inactivated fetal bovine serum (Gibco, USA). Cells were cultured at 37°C in a humidified incubator supplied with 5% CO2. pTwist-EF1alpha-SARS-CoV-2-S-2xStrep (a gift from Nevan Krogan, Addgene plasmid # 141382) was used to express S-protein from the original SARS-CoV-2 strain (WT). To express D614G S-protein this S-protein construct was mutated to introduce the D614G substitution. Additionally, we also expressed WT and D614 S-protein-FLAG (Addgene plasmids 156420 and 156421, gift of Hyeryun Choe) for improved immunostaining. Lipofectamine . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2020.08.12.246389 doi: bioRxiv preprint 2000 (ThermoFisher) was used to transfect cells with the S-protein DNA (per manufacture instructions). Maraviroc (Millipore Sigma, PZ0002) was dissolved in DMSO at a 1 mM (stock). Before or following transfection, 5 µM of Maraviroc solution was added to each well with complete EMEM medium. Vero cells were transferred to 35 mm Glass bottom dishes (Cellvis, D35-20-1-N) 24 hours after transfection. Then cells were fixed with 4% para-formaldehyde (Electron Microscopy Sciences, 15710) in PBS (Gibco, 14040182) 24 hours after attachment and then incubated with PBS. Thereafter, cells were mounted using Prolong Diamond Antifade Kit (Invitrogen, P36970) according to manufacturer's instructions. Hoechst 33342 (ThermoFisher, H3570) were used to identify nuclei. 2xStrep tagged S-protein was stained using StrepMAB-Classic DY-488 conjugate (IBA Lifesciences, 2-1563-050) and FLAG tagged S-protein was stained using anti-FLAG (clone M2, Sigma-Aldrich, F1804). SARS-CoV-2 Spike Antibody (ProSci, 3525) was also used to detect expressed S-protein. Actin-stain 555 phalloidin (Cytoskeleton, PHDH1-A) was used to label the actin cytoskeleton. Cells were imaged using a Zeiss LSM 710 Laser Scanning Microscope. Images were analyzed by Fiji (U.S. National Institutes of Health, Bethesda, MD). Total nuclei and S-protein transfected nuclei were counted using Cell Counter plugin.

Simulation modeling.
To simulate Maraviroc treatments, we adapted an open source, community-developed simulator (Version 3.2) for SARS-CoV-2 infection dynamics in tissue that currently includes diffusive viral spread through the tissue, ACE2 receptor dynamics, intracellular viral replication, single-cell viral and antiviral responses, and immune response (Getz et al., 2020). The simulator uses an agent-based modeling approach (PhysiCell (Ghaffarizadeh et al., 2018)) to simulate individual epithelial and immune cells. After adapting to the in vitro experimental conditions (no immune components), we extended the model to simulate cell fusion under the assumption that infected cells express spike protein on the surface (modeled to be correlate with the infected cell's assembled virus) that can bind to to nearby cells, leading to cellular fusion (Braga et al., 2021;Stoyanov et al., 2020). While inhibition experiments suggest that fusion is not necessarily driven by binding to ACE2 receptors on nearby cells (see Results), we use unoccupied ACE2 receptors in our model to more broadly model possible interaction targets (Hoffmann et al., 2020) to mediate fusion. Thus, we model the probability !"#$%& of two cells fusing as where h1 and h2 are Hill powers, [AV] is assembled virus in the infected cell (correlated with surface expression of spike protein), [ACE2] is unoccupied ACE2 receptors on the cell to be fused (proxy model of all possible binding targets), [AV] half and [ACE2] half are respective halfmaximum values, and P maxfusion is the maximum fusion probability when the cells have plentiful surface expression of spike protein and binding targets. In this model, fused cells continue viral replication with the same mass-action kinetics as individually-fused cells.
. CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2020.08.12.246389 doi: bioRxiv preprint To prepare for treatment simualtions, we added a generalized pharmacodynamic response model to each cell agent, where local drug availability could drive a combination of responses (including dose-dependent inhibition of cellular fusion, viral endo-and exocytosis, or viral replication) in each simulated cell. This allows us to model both the specific targeted effects of anti-fusion drugs as well as off-target, secondary, or mixed effects that are typically observed for drugs (e.g., Braga et al. (2021)). See Getz et al. (2020) for full details on the ACE2 receptor binding and trafficking, viral replication and exocytosis, and infected cell death response. See Figure S1A for the model overview.
To model the individual-cell pharmacodynamics, we first modeled the rate of drug internalization as diffusion through the cell membrane, allowing us to track the accumulated total drug (n) in each cell (Equation 2). Similarly to Aljayyoussi et al. (2017) and Aljayyoussi et al. (2019), we model the total effect E for each cell based on the internal drug concentration $&/0( using a Hill function (Equation 3). The effect is then used to modulate a behavioral parameter p (endocytosis rate, exocytosis rate, replication rate) from its untreated value (p0) to its value pmax at maximum drug effect when E = 1. (Equation 4). Note that '() = 1 for any parameter not affected by the simulated drug. In addition, we can model the impact of partial drug efficacy by multiplying the effect by a sensitivity S, where S = 0 denotes a drug with no efficacy (e.g., it fails to bind its target, or binding its target does not affect cell phenotype), and S = 1 is an ideal drug that fully binds and inhibits its target. In addition, we used a nonliear interpolation for calculating fusion probability under drug effect (Equation 5). Taken together, we model: Here, 5 is drug internalization rate, 6)/0( and $&/0( are the extracellular and intracellular drug concentrations, 76--is cell volume, EC >1 is the drug dose that gives 50% of maximum effect, and ℎ is Hill power (Goutelle et al., 2008). <1 is fusion probability without drug effect (see Equation 1), m is the inhibition coefficient.
Full simulation source code is available as open source under the BSD license; see Appendix A. We used xml2jupyter (Heiland et al., 2019) to create a cloud-hosted version of this model; interested readers can interactively run the model in a web browser at https://nanohub.org/tools/vi-rion2pbpd ). See Figure S1B-C.

Simulation model fitting
We began with the default parameter values from Getz et al. (2020) and adjusted to match viral titer measurements at 16 hours under control conditions and a 5 µM dose of Maraviroc. First, we adjusted the receptor dynamics parameters (ACE2 binding rate, ACE2 endocytosis rate, ACE2 cargo release rate and ACE2 recycling rate), viral replication parameters (virion uncoating rate, RNA preparation rate, viral protein synthesis rate, and virion assembly rate) and export rate, and the . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2020.08.12.246389 doi: bioRxiv preprint fusion probability to match the viral titer in untreated Vero cells at 16 hours, assuming that 3000 cells are initially infected.
To calibrate the drug uptake rate, we looked at prior experiments which show that peak cytoplasmic drug concentrations are reached in about 3.1 hours (300mg dose) (Carter and Keating, 2007). To match this drug uptake time scale, we set the internalization rI at 0.0116 M • µ ? /min (1 hour to reach half maximum). We initially approximated EC50 as 2.5 M (half of experimental dose), the Hill power n as 1, and m as 1 as an intial parameter estimate.
We then sought a parsimonious pharmacodynamic model that could fit the full response of the 5 µM Maraviroc dose, with both a reduction in viral titer at 16 hours and cell fusion at 48 hours. First, we attempted to fit the simplest possible model where the drug inhibits cell fusion only but has no effect on viral endocytosis and exocytosis, and no effect on replication. While we were able to fit the model to match the reduction of cell fusion, this model could not match experimental measurements of reduced viral titer. Hypothesizing that a drug that impacts viral protein localization on the cell surface may also impact endocytosis and exocytosis (similarly to other recent observations (Hoffmann et al., 2020)), we next modeled the drug as inhibiting fusion, endocytosis, and exocytosis (with equal impact on endocytosis and exocytosis, for simplicity). To match with both the reduction of viral titer at 16 hours and the inhibtion of cellular fusion at 48 hour with a dose of 5 µM, we varied the Hill power n and inhibition coefficient m until the full dynamical model yielded a viral titer and cell fusion that matched experimental measurements. The resulting parameters were EC50 = 2.5 M, n = 2.5, and m=2.5. See Goutelle et al. (2008) for more details on common pharmacodynamic forms.
Because this was the simplest (most parsimonious) set of drug effect assumptions that could fit the data, we used it for the remainder of the study. We note that simulating Maraviroc as partially inhibiting replication as well as endo/exocytosis and fusion could also potentially fit the data, but as a less parsimonious model.

Computing.
We used a total of 13 parameters sets. Because the simulation model is stochastic, we ran 10 simulation replicates for each parameter set. All simulation images are chosen from one representative replicate, and all aggregate dynamical curves (e.g., viral titer versus time) are reported as the mean and confidence interval as 95% over all 10 simulation replicates. The simulations were performed on the Big Red 3 supercomputer at Indiana University. Each simulation was run on a single node using 24 threads. All jobs (130 runs) were submitted as a batch. Each simulation was completed in approximately 25 minutes, and with a total wall time of approximately of 9 hours.

Data and Code Availability
. CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2020.08.12.246389 doi: bioRxiv preprint The simulation source code is available at:

Results
Identification of anti-SARS-CoV-2 activity of FDA approved small molecules.
The investigation was conducted with 19 FDA-approved small molecules that are commonly prescribed to the general population which included antiviral and anti-parasitic compounds. The investigation also included compounds that were documented to alter host signaling cascade events and exhibit anti-inflammatory properties. The list of small molecules that were tested is included in Figure 1A. The investigation was conducted using Vero cells and the Washington strain of SARS-CoV-2 (2019-nCoV/USA-WA1/2020). The toxicity of the small molecules was first quantified by treating the cells with the inhibitor at a low, medium and high concentration (1 μM, 5μM and 10 μM) for 24 hours after which cell viability was quantified using Cell Titer Glo assay ( Figure 1B). Almost all of the small molecules demonstrated less than 10% cytotoxicity in this assay with the exception of FTY720 and Atorvastatin that showed 20% and 10% reduction in cell viability (compared to the untreated control) at the highest concentration. Anti-SARS-CoV-2 activities of these small molecules were quantified by plaque assay and the outcomes are included in Figure 1C. Vero cells were pre-incubated with nontoxic concentrations of the compounds for 1 hour prior to infection. Infection with SARS-CoV-2, Isolate USA-WA1/2020 was carried out by incubating the virus with cells (multiplicity of infection [MOI]: 0.1) in the presence of the designated inhibitor for 1 hour. After allowing time for the infection, the virus overlay was removed and media containing the appropriate inhibitor was added back to the cells. DMSO treated cells were maintained as controls alongside. The amount of infectious virus in the supernatants of the drug-treated cells was quantified by plaque assay at 16 hours post infection and compared to the DMSO-treated control.
Among the antiviral compounds that were tested, Maraviroc demonstrated a 1.7-fold decrease (p<0.0001) in viral load as compared to the DMSO control ( Figure 1C). Among the anti-parasitic compounds that were tested, Nitazoxanide demonstrated a 7.1-fold decrease (p<0.0001) when compared to the DMSO control. Interestingly, Nitazoxanide was suggested as a potential anti-SARS-CoV-2 therapeutic based on prior activity observed against MERS-CoV (Martins-Filho et al., 2020). Nitazoxanide was also suggested as a combinatorial therapeutic along with Azithromycin (Kelleni, 2020). Nitazoxanide is currently undergoing clinical trials for safety and efficacy to treat COVID-19 patients (ClinicalTrials.gov Identifier: NCT04435314). In the set of compounds that influenced host signaling and inflammatory events, FTY-720 showed a 38.4fold decrease (p<0.0001) of SARS-CoV-2 viral titer as compared to the DMSO control. FTY-. CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2020.08.12.246389 doi: bioRxiv preprint 720 is FDA-approved for the treatment of multiple sclerosis and is documented to exhibit antiinflammatory activities by modulating the NFkB signaling cascade (Pul et al., 2016). In this class of compounds, Atorvastatin, a commonly used statin for the treatment of high cholesterol, was noted to decrease viral load 1.6-fold (p<0.01) as compared to the DMSO control. The relative percent inhibitions of these prioritized compounds is indicated in Figure 1D, in comparison with the DMSO control. Maraviroc and FTY-720 were screened against SARS-CoV-2, Isolate Italy-INMI1 ( Figure 1E) in addition to Isolate USA-WA1/2020.Maraviroc decreased viral titers 45.2% while FTY-720 showed a 97.4% decrease.

Identification of anti-SARS-CoV-2 activity of FDA approved biologics.
Three FDA-approved biologics namely, Bevacizumab, Eculizumab and CD126/IL6R mAb were evaluated for anti-SARS-CoV-2 activity in a similar manner as described above (Figure 2A). Bevacizumab (Avastin) is an inhibitor of VEGF-R2 by binding to VEGF-A and is used for the treatment of several types of cancers (Mahfouz et al., 2017;Papachristos et al., 2019). Eculizumab (Soliris) is undergoing clinical trials for use as a therapeutic for the treatment of COVID-19 patients, potentially for the anti-inflammatory activities and the ability to moderate inflammation related complications (ClinicalTrials.gov Identifier: NCT04288713). Briefly, toxicity of the biologics was evaluated by incubating Vero cells with a low (1 μg) and high (2.5 μg) concentration of the target biologic for 24 hours after which cell viability was quantified by Cell Titer Glo ( Figure 2B). Less than 10% cytotoxicity was detected for these molecules at the concentrations tested, with the exception of CD126/IL-6R antibody. A concentration of 2.5 µg was selected for efficacy evaluation for the less toxic biologics and 1 µg was selected for the CD126/IL-6R antibody. Impact of inhibitor treatment on SARS-CoV-2 infectivity was assessed by performing infection of Vero cells (MOI: 0.1) in the presence of the biologic at nontoxic concentration for 1 hour. The infected cells were incubated with the biologic for an additional 16 hours after which infectious virus load in the supernatant was quantified in comparison with the DMSO control ( Figure 2C). Bevacizumab and Eculizumab treatment resulted in a 4.2-fold decrease (p<0.0001) and a 3.7-fold decrease of virus titers, respectively, which coorilate to 76.2% and 73.1% total reduction in virus respectively ( Figure 2D). CD126/IL6R antibody resulted in 1.9-fold decrease of virus load.

Interference of anti-SARS-CoV-2 small molecules with ACE2:Spike protein interaction.
A VSV-G pseudovirus with a firefly luciferase reporter which retains the SARS-CoV-2 spike protein on the surface was used to query if the prioritized compounds affected early steps of the infectious process. The pseudovirus is capable of ACE2-based viral entry which can be quantified by measuring intracellular luciferase expression. As the pseudovirus is not capable of viral RNA synthesis, the inhibitory effect is exclusive to the early entry and post-entry steps. The stock pseudovirus was diluted 1:5 in media as specified in the protocol and was used to infect Vero cells in the presence of the inhibitor and luciferase expression was measured at 16 hours post-infection ( Figure S2A). While Nitazoxanide exerted a 74.1% inhibition of luciferase expression, Maraviroc and FTY720 produced 28.6% and a 41.1% inhibition, respectively, of reporter expression suggesting that Maraviroc or FTY720 were likely to inhibit additional steps besides entry to achieve viral inhibition.
. CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10. 1101 To further ascertain if this observed inhibition to viral entry exerted by Nitazoxanide involved interference in the molecular interaction between the viral spike protein and the cellular ACE2 receptor, a biotinylated receptor binding assay was employed. In this assay, the potential for the inhibitor to interfere with the binding of spike protein and ACE2 receptor is quantified as a measure of biotin expression. The assay involves using a nickel-coated plate to which his-tagged ACE2 is added. The drug is then added to the plate and incubated for 1 hour at room temperature after which the SARS-CoV-2 spike protein bound to biotin is added and the plate is incubated for 1 hour after which Streptavidin-HRP is added. Chemiluminescence is measured following the addition of HRP substrate which is a reflection of the interaction between the ACE2 protein and S-protein. Quantification of luminescence in drug-treated wells as compared to the DMSOtreated control demonstrated that Maraviroc, FTY 720 or Nitazoxanide did not interfere with the interaction at lower concentration of the inhibitor (Figure S2B). Atorvastatin elicited a modest interference in the interaction as demonstrated by a reduction of 31.7% (p<0.01) of biotin expression ( Figure S2B). When the experiment was conducted at higher concentrations of the drug (100 µM), FTY720 was able to exert a modest interference (27% relative to control, p<0.01), while Maraviroc and Nitazoxanide treatment did not produce any impact on the interaction. The data suggest that the mechanism of inhibition for Maraviroc and Nitazoxanide is unlikely to involve disruption of receptor binding by the viral spike protein and is likely that early post-entry steps are impacted by the drug treatment ( Figure S2B).

Impact of Maraviroc on intracellular events in SARS-CoV-2 infected cells.
The observations that Maraviroc exhibited strong inhibition of SARS-CoV-2 when infection was performed in the presence of the inhibitor, without necessarily impacting the interaction of the viral surface spike protein with the ACE2 receptor suggested that alternate mechanisms that may be critical for the establishment of a productive viral infection may be impacted by Maraviroc. Evidence in the literature indicates that the coronavirus spike protein drives membrane fusion resulting in the formation of multinucleate cells (Fehr and Perlman, 2015;Lontok et al., 2004;Sisk et al., 2018). Notably, this phenomenon was observed to be elicited by over expressed S-protein (spike protein) as well, suggesting that it is independent of functionality from other proteins. As a next step, the impact of Maraviroc on S-protein based cellular fusion was assessed by confocal microscopy ( Figure 3A). Vero cells were transfected with an over expression construct that expressed the SARS-CoV-2 spike protein. Transfected cells were maintained in the presence of Maraviroc for up to 48 hours to permit expression of S-protein. DMSO-treated cells were maintained as controls to ascertain baseline numbers of multinucleate cells. The cells were fixed and stained at 48 hours post transfection and the extent of formation of multinucleate cells were quantified by confocal microscopy (Figure 3A-B). The data demonstrate that increased expression of S-protein in Vero cells results in the formation of large multinucleate cells that typically contain 4-8 nuclei per cell. Vero cells transfected with a GFP control plasmid did not form multinucleate cells (data not shown). Interestingly, maintaining the transfected cells in the presence of Maraviroc reduced the frequency of multinucleate cell formation by approximately 2-fold as compared to the DMSO-treated controls (Figure 3C), as well as decreasing the average number of nuclei in these clusters ( Figure 3D) and % fusion ( Figure 3E) for S-protein from wild-type and D614G strains . This outcome suggests the possibility that Maraviroc treatment may result in . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2020.08.12.246389 doi: bioRxiv preprint a decrease in viral load in the cell population by preventing membrane fusion events that impact intracellular viral multiplication and spread. The mechanism of action for S-protein to prevent cell-cell fusion is not yet identified. Staining of S-protein in Maraviroc-treated cells showed increased cytosolic staining ( Figure 3B), which may indicate that Maraviroc inhibits trafficking of S-protein to the extracellular surface of the cell.

Determination of synergy between the anti-SARS-CoV-2 small molecules.
The ability of Maraviroc to exert synergistic inhibitory outcomes when combined with FTY720, Atorvastatin or Nitazoxanide was assessed. As a first step, toxicity of the combination of drugs was assessed when the drugs were combined at 20%, 50% and 100% of their effective concentrations as observed in Figure 1D. To that end, the four compounds were assessed when they were combined at 1 μM each, 2.5 μM each and 5 μM each. The combinations were used to treat Vero cells and cell viability was quantified at 24 hours post treatment by Cell Titer Glo assay ( Figure  4A). The combinations that resulted in >90% cell survival (<10% toxicity) were considered for efficacy assessment. To assess potential synergy between these compounds, infection with SARS-CoV-2 were carried out in the presence of the drug combinations and the infected cells were maintained for 16 hours in the presence of the combinations. The viral load in the supernatants were quantified at 16 hours post infection by plaque assay and relative inhibition compared to that of the DMSO control ( Figure 4B). A modest increase in inhibitory potential was noted when Maraviroc was combined with any of the three compounds ( Figure 4C). Inhibitory potential was increased 2.6% for FTY720, 10.6% for Atorvastatin and 11.9% for Nitazoxanide. Remdesivir was also screened for antiviral activity and possible synergy with Maraviroc. Toxicity of Remdesivir alone was below 10% for the concentrations that were screened. Combining Remdesivir with Maraviroc did not substantially increase cytotoxicity, and cell viability remained above 90%. Remdesivir alone decreased viral titer by 91.1% (p<0.0001). When Maraviroc was added to the treatment viral titers were decreased by 94% (p<0.0001).

Modeling of Maraviroc mediated inhibition of SARS-CoV-2 in a multicellular environment.
To study the dynamics of SARS-CoV-2 spread under control and Maraviroc-inhibited conditions, we adapted a computational model of SARS-CoV-2 (Getz et al., 2020) to include cellular fusion driven by surface-bound spike protein, and integrated a mixed-effect pharmacodynamic response model that can simulate any combination of inhibiting cell fusion, viral entry and exit, and replication. After verifying that single-effect pharmacodynamics models that solely inhibit cell fusion or viral entry/exit cannot simultaneously match our experimental measurements and observations that Maraviroc reduces both viral load at 16 hours, and cell fusion at 48 hours, we fitted a mixed-effect model that inhibits both cellular fusion and viral entry and exit. (See Fig. S3 and Simulation model fitting above.) In particular, we tuned the viral binding, endocytosis, replication, exocytosis, and pharmacodynamics parameters to match the viral titer at 16 hours, and the number of nuclei per cell at 48 hours under no treatment and at 5 µM of Maraviroc.
. CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ;https://doi.org/10.1101https://doi.org/10. /2020 The fitted model suggests that exocytosis is not fully inhibited at 5 µM: because the original 3000 virions are quickly bound and endocytosed, replication and exocytosis are required to explain the presence of unbound virus at 16 hours. Moreover, because Maraviroc does not inhibit SARS-CoV-2 binding to ACE2 receptors on uninfected cells, sustained replication and exocytosis are required to outpace binding to ACE2 receptors. In spite of incomplete inhibition, fewer cells are infected by 16 hours, and fewer infected cells have a high viral load of assembled virions. (Figure S3) When we extended the simulations to 2 days (Figure 5A), we found that the 5 µM dose ( Figure  5B) delayed but did not prevent the spread of the infection. By 2 days, all viable cells were infected with high viral load (non-blue cells), but cell survival was increased over the untreated control ( Figure 5A). Increasing the simulated dose to 7.5 µM ( Figure 5C) or 10 µM (Figiure 5D) further increased cell survival, and some cells remained uninfected at 2 days.
We also examined the potential to combine a moderate (5 µM) dose of Maraviroc with an antireplication drug such as Remdesivir. First, we calibrated an anti-replication pharmacodynamic response to match the Remdesivir-only experiments, and simulated to 2 days ( Figure 5E). Next, we simulated a combination of Remdesivir with 5 µM Maraviroc without any additional calibration. Just as in experiments, the simulation model predicted that the combination would provide a modest improvement (Figure 5F).
After verifying that the computational model recapitulated these results, we further investigated the role of cellular fusion in viral spread in untreated conditions. When we completely disabled cellular fusion in untreated conditions, we found virtually no reduction in viral load ( Figure  S4A), suggesting that extracellular spread of virions may ordinarily be more efficient at transmission than fusion. When we reduced extracellular virion spread (e.g., as might be expected in a tissue where virions are quickly removed by flow or neutralizing antibodies), we found that fusion substantially increased viral transmission (Figure S4B-C). Thus, cellular fusion could play a role as a "backup transmission" method when extracellular transmission is impaired by treatment or an adaptive immune response.

Discussion
SARS-CoV-2, the causative agent of the ongoing COVID-19 pandemic, continues to have a devastating effect globally, increasing the morbidity and mortality burden on a daily basis and giving rise to significant long-term health impacts even in cases of mild infections (Hayes et al., 2021;Michelen et al., 2021;Nalbandian et al., 2021). From the perspective of bringing quick treatment options to the clinic, there continues to be heightened focus on repurposing FDA-approved drugs. When a therapeutic intervention approach to treat COVID-19 is to be evaluated, it cannot be ignored that the population that presents itself at the clinic is likely to fall in the more severe disease scale, with a high probability of co-morbidities and ongoing medication regimens. This possibility has to be considered from two angles: 1. the aggressive therapeutic strategy that may be used to treat SARS-CoV-2 infection should not interfere with ongoing drug regimens or . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2020.08.12.246389 doi: bioRxiv preprint exert unforeseen catastrophic outcomes; 2. the possibility of treatment strategies that are commonly used to treat disease conditions such as cardiovascular disease, asthma or immune system disorders may be used to also treat SARS-CoV-2 infection. The latter aspect offers not only the possibility of being able to safely use some of these therapeutics in the treatment of these high risk groups, but also paves the way for determining if any of these ongoing therapeutic regimens may offer a level of protection to the target group by preventing escalation of disease severity.
In this study, a total of 22 small molecules and biologics that are commonly used to treat disease states including viral, bacterial and parasitic infections, heart disease and inflammatory conditions such as multiple sclerosis were analyzed for potential antiviral activity against SARS-CoV-2 in a cell culture model. The investigation identified four small molecules that demonstrated statistically significant anti-SARS-CoV-2 activity, namely, Maraviroc, FTY720, Nitazoxanide and Atorvastatin. Maraviroc elicited a statistically significant and consistent inhibitory effect against SARS-CoV-2 ( Figure 1C). Maraviroc is an inhibitor of the cellular CCR5 receptor and is used to treat HIV infection. Interestingly, Maraviroc is under consideration for a clinical trial as a COVID-19 therapeutic as it is known to be well tolerated in human populations, even in the context of a compromised immune system (ClinicalTrials.gov Identifier: NCT04435522). The rationale behind the use of Maraviroc is that it would be beneficial in preventing immune cell depletion (or exhaustion) and increasing routing of inflammatory cells away from sites that sustain inflammatory damage during COVID-19 disease, thus additively contributing to lowering viral load and inflammatory damage in organs such as the lungs. FTY720 (Fingolimod) is a drug that is used to treat multiple sclerosis which demonstrated robust inhibitory activity against SARS-CoV-2 ( Figure 1C). It functions as an immune modulator that is currently undergoing clinical trials for use as a COVID-19 therapeutic (ClinicalTrials.gov Identifier: NCT04280588). The rationale for consideration of FTY720 is that it may be able to reduce the inflammatory manifestations of ARDS in the lung by decreasing pulmonary edema and hyaline membrane formation, possibly if used in conjunction with additional methods including ventilator support. Furthermore, the mechanisms of FTY720-mediated reduction of egress of autoreactive lymphocytes from lymphoid tissues may play additional roles in decreasing systemic inflammatory damage. FTY720 treatment may produce short term fluctuations in cardiovascular parameters and blood pressure due to agonistic and antagonistic impacts on Sphingosine-1-phosphate receptor (S1PR) which is differentially expressed in tissues including cardiac cells (Camm et al., 2014). FTY720 has also been implicated in cases of Fingolimod associated macular edema, especially when administered in the same regimen as for multiple sclerosis (Pul et al., 2016). These factors will require consideration if FTY720 is to be evaluated for COVID-19 treatment.
Nitazoxanide, an FDA-approved drug used to treat diarrhea, exerted antiviral activity in the context of SARS-CoV-2 infection ( Figure 1C). There have been recent published reports that refer to the use of Nitazoxanide in combination with other drugs including Azithromycin and Ivermectin, both of which were demonstrated to show anti-SARS-CoV-2 activity in cell culture (Kelleni, 2020). Nitazoxanide has been demonstrated to have antiviral activity against MERS-CoV and potentiate interferon responses (Rossignol, 2016). In the context of MERS-CoV infection model, . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10. 1101 Nitazoxanide was also shown to decrease the expression of the proinflammatory cytokine IL6, which is also implicated in the inflammatory damage elicited due to SARS-CoV-2 infection. Nitazoxanide is under consideration for a clinical trial as a COVID-19 treatment strategy (Clini-calTrials.gov Identifier: NCT04435314).
Atorvastatin, a commonly used small molecule for the treatment of cardiovascular disease demonstrated modest antiviral activity against SARS-CoV-2 when compared to the other effective compounds (Figure 1C). Recently published outcomes suggested that statin use in patients is likely to be linked to decreased mortality, ventilator use and hospitalizations in the case of COVID-19 patients . Atorvastatin is currently undergoing clinical trials as adjunctive therapy for potentially addressing the cardiovascular complications associated with COVID-19 (ClinicalTrials.gov Identifier: NCT04380402).
Studies performed by varying the infection condition, specifically the presence or absence of the candidate drug at the time of infection revealed that including the drug during infection produced a higher rate of inhibition (data not shown). This observation led to the hypothesis that the candidate compounds, Maraviroc, FTY720, Nitazoxanide and Atorvastatin may interfere with entry and/or early post-entry steps in the SARS-CoV-2 infectious process. To investigate this possibility further, a VSV-G pseudovirus which expresses the SARS-CoV-2 spike protein on the surface was used in the context of these inhibitors. The extent of inhibition observed with the pseudovirus was significantly lower than what was noted in the context of SARS-CoV-2 ( Figure S2A) suggesting that additional steps beyond viral entry are likely to be influenced by these small molecules. It may also be possible that the VSV glycoproteins that are also expressed on the surface of the pseudovirion may interfere with the inhibition that these compounds may exert on the SARS-CoV-2 spike protein. Nevertheless, the observation that the early entry steps are impacted by these compounds led to the question whether they interfered with the interaction between the viral spike protein and the host ACE2 receptor. Analysis of binding interaction by a reporter assay in the presence and absence of the inhibitor revealed that while Atorvastatin may elicit a modest level of interference in the interaction between the spike protein and ACE2, Maraviroc and FTY720 did not, even at higher concentrations ( Figure S2B). Maraviroc and FTY720 were the most potent of the inhibitors that were tested in this effort and the data supported the idea that the mechanism of inhibition is not likely to involve interference with receptor binding.
Coronavirus spike protein has been shown to be involved in membrane fusion, both in the context of fusion between viral and host cell membranes and in the context of host cell membranes (Tang et al., 2020). The membrane fusion event is believed to be mediated by the S2 domain of the extracellular stalk portion of the spike protein. Within the S2 domain, a short, loosely defined fusion peptide, a stretch of 15-25 amino acids is important to initiate the fusion event and is considered to be very sensitive to mutations (Madu et al., 2009). The SARS-CoV-2 fusion peptide region in S2 is considered to be 93% similar to that observed in SARS-CoV-1, thus indicating a high degree of conservation in the fusion peptide (Tang et al., 2020). SARS-CoV-1 S-protein was shown to initiate cell fusion between human cells, which was also observed in the case of SARS-CoV-2 (Ou et al., 2020) and the protease TRPRSS11a is thought to aid this fusion event because of S-protein processing in the context of SARS-CoV-1 (Kam et al., 2009). Studies that . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ;https://doi.org/10.1101https://doi.org/10. /2020 were conducted with an over expressed SARS-CoV-2 S-protein expression construct demonstrated an increased incidence of multinucleate cells with an average of 4-7 nucleic per multinucleate cell (Figure 3). When the transfected cells were maintained in the presence of Maraviroc, a statistically significant reduction in the number of multinucleate cells were observed ( Figure  3B). This observation raises the interesting possibility of viral spread in a multicellular environment mediated by cell fusion. Such a possibility holds significant relevance to the efficacy of convalescent serum use as an intervention strategy because the intracellular viral population will be shielded from the antibodies and will continue to persist. Synthetic antibodies or nanobodies that are membrane permeable may prove to be effective in such cases. If such a possibility is true, it will explain the observed cases of "relapse" or "prolonged persistence" of viral RNA even in individuals who do not exhibit any disease symptoms. Such an event will also have a significant impact on vaccine design, particularly those that rely on humoral responses to extracellular virus, drawing attention to CD8 mediated responses as being an important requirement for an efficacious vaccine.
An interdisciplinary coalition is actively developing an open source, multiscale simulation model for SARS-CoV-2 infection dynamics in tissue (Getz et al., 2020), which can model the heterogeneous dynamics of individual cells and the spread of viral particles through a tissue (or cell culture). This study adapted and extended the simulator by adding celllar fusion and a Maraviroc pharmacodynamics model to each cell and calibrating to the untreated control and 5 µM experiments at 16 hours. As employed here, the simulator emphasizes the early entry and post entry events in the context of Maraviroc treatment as a test case with a goal of implementing a full in vivo simulation model to assess the clinical applicability of the drug. The simulations revealed that the inhibition is likely to extend protective effects in a multicellular setting and decrease cell death in relation to a lowered viral load if the effective drug concentration can be maintained. Because the model explicitly simulates viral binding, endocytosis, replication, release and spread as well as cellular fusion, it can help "fill in the gaps" between limited experimental observations. In this effort, the model suggests that antiviral drugs such as Maraviroc can reduce viral titer by over two orders of magnitude even with partial inhibition of its target at moderate doses. Moreover, it also agrees with experimental findings that combining Maraviroc could complement existing drugs that target replication (e.g., Remdesivir) for further efficacy, and that targeting cellular fusion may be beneficial when contact-based viral transmission otherwise evades treatment. Future simulation studies will assess whether Maraviroc can act synergistically with innate immune responses to improve control of the virus and reduce tissue damage. While simulations were employed in the context of Maraviroc as an initial demonstration of utility, the additional inhibitors identified in this study will be queried independently for control of viral and inflammatory load, relative pharmacodynamics properties as related to antiviral and anti-inflammatory activities and synergistic outcomes. The simulation model could also afford the opportunity to assess the relative impact of cell-to-cell viral transmission (e.g., via fusion) compared to diffusive dissemination of endocytosed virus particles, and whether treatments like Maraviroc could ultimately halt these alternative forms of viral spread in tissues.
In conclusion, our data demonstrate that Maraviroc, FTY720, Nitazoxanide and Atorvastatin are viable candidates that may be used to treat individuals infected with SARS-CoV-2. These drugs . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2020.08.12.246389 doi: bioRxiv preprint are commonly used to treat a wide demographic of the human population and are known to be well tolerated in humans. The data presented in this effort will also provide support to the ongoing clinical trials that are exploring the utility of these compounds as COVID-19 treatment strategies. Xia, S., Liu, M., Wang, C., Xu, W., Lan, Q., Feng, S., Qi, F., Bao, L., Du, L., Liu, S., et al. (2020). Inhibition of SARS-CoV-2 (previously 2019-nCoV) infection by a highly potent pancoronavirus fusion inhibitor targeting its spike protein that harbors a high capacity to mediate membrane fusion. Cell Research 30, 343-355.   (B.) Drug toxicity in Vero cells at 24 hours. Cells were seeded at a density of 3E4 and overlaid with media conditioned with 0.1% DMSO or the indicated concentration of inhibitor. Cells were incubated for 24 hours at 37ºC at 5% CO2. Viability was determined using CellTiter Glo assay.

Trejo
. CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2020.08.12.246389 doi: bioRxiv preprint (C.) Efficacy of small molecule inhibitors in Vero cells 16 hpi. Cells were seeded at a density of 3E4 and overlaid with conditioned media containing 0.1% DMSO or 5 μM of inhibitor. Cells were incubated for 2 hours. Conditioned media was removed. Viral inoculum (MOI: 0.1) containing 5 μM of inhibitor was added to cells and incubated for 1 hour at 37ºC and 5% CO2. Viral inoculum was removed, cells were washed with DPBS, and conditioned media was returned to the cells. Culture was incubated for 16 hpi. Supernatants were collected and a viral plaque assay was used to quantify the viral titers.
The quantitative data are depicted as the mean of three biologically independent experiments ± SD. * p < 0.05; ** p < 0.01; *** p < 0.001; ****p < 0.0001. (B.) Antibody toxicity in Vero cells at 24 hours. Cells were seeded at a density of 3E4 and overlaid with media conditioned with 0.1% DMSO or the indicated concentration of antibody. Cells were incubated for 24 hours at 37ºC at 5% CO2. Viability was determined using CellTiter Glo assay.
(C.) Efficacy of antibodies in Vero cells 16 hpi. Cells were seeded at a density of 3E4 and overlaid with conditioned media containing 0.1% DMSO or maximum permissible dose of antibody. Cells were incubated for 2 hours. Conditioned media was removed. Viral inoculum (MOI: 0.1) was added to cells and incubated for 1 hour at 37ºC and 5% CO2. Viral inoculum was removed, cells were washed with DPBS, and conditioned media was returned to the cells. Culture was incubated for 16 hpi. Supernatants were collected and a viral plaque assay was used to quantify the viral titers.
(D.) Inhibition that was observed at 16 hpi, represented as relative percentages.

Figure 3 -Maraviroc inhibits Wild Type (WT) and D614G strain S-protein induced cellcell fusion.
(A.) Vero cells were transfected with S protein tagged with Strep-tag II and treated with DMSO (labeled as condition ND for no drug) or 5 µM Maraviroc for 48 hours (treated immediately after transfection; labeled as condition A), or 5 µM Maraviroc for 72 hours (treated 24 hours before . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2020.08.12.246389 doi: bioRxiv preprint transfection, labeled as condition B). Immunohistochemistry showed that S protein expression resulted in the formation of multinucleated cells, which was inhibited with Maraviroc treatment.
(B.) The frequency of multinucleated cells, normalized to total cell number, was reduced with Maraviroc treatment.
(C.) The average number of nuclei in multinucleated cells was also reduced with Maraviroc treatment.
The quantitative data are depicted as the mean of three biologically independent experiments ± Standard Error ** p < 0.01; *** p < 0.001, ***p < 0.0001. (A.) Combination drug toxicity in Vero cells at 24 hours. Cells were seeded at a density of 3E4 and overlaid with media conditioned with 0.1% DMSO or the indicated concentration of inhibitor. Cells were incubated for 24 hours at 37ºC at 5% CO2. Viability was determined using CellTiter Glo assay.
(B.) Efficacy of inhibitors combinations in Vero cells at 16 hpi. Cells were seeded at a density of 3E4 and overlaid with conditioned media containing 0.1% DMSO or 5:5 μM of Maraviroc plus indicated inhibitor. Cells were incubated for 2 hours and conditioned media was removed. Viral inoculum (MOI: 0.1) containing 5:5 μM of inhibitor combination was added to cells and incubated for 1 hour at 37ºC and 5% CO2. Viral inoculum was removed, cells were washed with DPBS, and conditioned media was returned to the cells. Culture was incubated for 16 hpi after which supernatants were collected and a viral plaque assay was used to quantify the viral titers.
(C.) Inhibition that was observed at 16 hpi, represented as relative percentages.
The quantitative data are depicted as the mean of three biologically independent experiments ± SD. * p < 0.05; ** p < 0.01; *** p < 0.001; ****p < 0.0001. We simulated a 2-day treatment with an anti-replication drug (e.g., Remdesivir) calibrated to approximately match experimental responses at 2 days (E). Without additional fitting, we then simulated a combination treatment with 5 μM Maraviroc and the same anti-replication treatment (F). Just as in the experiments, the computational model predicts an additional synergistic improvement in cell survival. We plot the impact on viral load for these virtual experiments in (G) . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2020.08.12.246389 doi: bioRxiv preprint


s2fieldsofstudy: Biology Medicine

text:
Knockdown of lncRNA HAGLROS inhibits metastasis and promotes apoptosis in nephroblastoma cells by inhibition of autophagy

ABSTRACT Nephroblastoma, or Wilms tumor, is a primary renal malignant tumor that easily occurs in children. Previous studies have revealed the regulatory functions of LncRNA in nephroblastoma. LncRNA HAGLROS functions as a tumor promotor in various cancers including hepatocellular carcinoma, ovarian cancer and colorectal cancer. In this study, the HAGLROS expression in nephroblastoma cells was assayed through qRT-PCR. Cell proliferation assessment employed CCK-8. Moreover, the migration and invasion of cells were examined separately through wound healing and transwell assay. Moreover, flow cytometric analysis and Western blot assay were applied to evaluate cell apoptosis. Rapamycin and 3-methyladenine were used to serve as autophagy activator or inhibitor, respectively. In addition, autophagy was identified by immunofluorescence staining and Western blot analysis. Experiment results showed that HAGLROS expressed highly in nephroblastoma cell lines. HAGLROS knockdown prevented cells from proliferating, and also showed suppressive impact on migration and invasion in HFWT cells. In addition, knockdown of HAGLROS showed a facilitative effect on apoptosis and an inhibitory effect on autophagy. Stimulation of autophagy alleviated HAGLROS silencing-induced apoptosis, while inhibition of autophagy reversed the effect in nephroblastoma cells. In summary, our results revealed that HAGLROS executed an oncogenic role in the progress of nephroblastoma, offering a new perspective on the strategy for nephroblastoma therapy.

Introduction
Nephroblastoma, or Wilms tumor, is a primary malignant tumor that occurs in children [1][2][3]. It accounts for 6-7% of malignant tumors in children younger than 15, about 75% of which are distributed in children under 5 years old [4]. Treatment for Wilms' tumor is administered according to its diagnosis and stage, usually surgical resection, chemotherapy, or targeted therapy are implemented [5]. Among these approach, comprehensive treatment based on the combination of surgery, chemotherapy and radiotherapy is usually a more practical treatment for all stages of nephroblastoma [6]. Although the total survival rate is over 90% at 5 years, plenty of side effects and sequelae of drug treatment such as malnutrition, cardiac dysfunction and secondary malignancies, still trouble the patients [7]. In addition, there are still a few of children died of relapse, metastasis and resistance to chemotherapy and other reasons [8]. Thus, great efforts on exploring the pathogenesis of nephroblastoma may throw light on the identification of molecular diagnostic markers or therapeutic targets for nephroblastoma.
Long non-coding RNAs (lncRNAs) are a class of RNAs that have no protein coding capacity with a length of over 200 nucleotides [9]. LncRNAs are capable of regulating the transcription, translation, and protein modification by interacting with DNAs, RNAs or proteins [10]. It is widely accepted that LncRNAs are intimately involved in the formation, development and metastasis of numerous cancers [11,12]. Earlier researches indicated that several LncRNAs are critical in cell proliferation, migration, apoptosis and cell cycle of nephroblastoma [13,14]. Based on the study of Ge et al., lncRNA deleted in lymphocytic leukemia 1 (DLEU1) expression was upregulated in nephroblastoma tissues and cell lines and knockdown of DLEU1 suppressed the proliferation, migration and invasion of GHINK-1 cells via targeting the miR-300/HOXC8 axis [15]. HAGLROS, one of the lncRNAs, has been reported to be abnormally expressed in many types of malignancies, presenting a carcinogenic role and contributing to the malignant processes of a number of cancers [16,17]. For instance, Wei et al. reported that HAGLROS promoted cell proliferation and cell autophagy but inhibited apoptosis through mediating miR-5095/ATG12 axis in hepatocellular carcinoma cells [18]. In addition, Wang et al. demonstrated that HAGLROS expression was correlated with tumor stage and lymphatic metastasis of lung cancer. Furthermore, HAGLROS accelerated proliferation and migration of Lung cancer cells by regulating microRNA-152 [19]. However, the biological role of HAGLROS in the occurrence and progression of nephroblastoma remains incompletely discovered.
This study was designed to look into the expression and biological role of lncRNA HAGLROS in nephroblastoma cells. We demonstrated that HAGLROS expression was higher in nephroblastoma cell lines and knockdown of HAGLROS yielded suppressive effects on the development of Wilms tumor, which provides a reliable new target for the therapy of nephroblastoma. We hypothesized that HAGLROS silencing inhibits nephroblastoma development by regulating cell proliferation, migration, invasion and apoptosis. In addition, we intend to explore whether HAGLROS silencing promotes nephroblastoma cell apoptosis by the inhibition of autophagy.

Cell culture and treatment
Normal human embryonic kidney cell line HEK-293 T and nephroblastoma cell lines HFWT, 17-94 and WiT49 were supplied by American Type Culture Collection (ATCC; Manassas, VA, USA). For the culture of these cells, they were usually conducted in Dulbecco's Modified Eagle Medium (DMEM; Invitrogen, Carlsbad, CA, USA) 10% fetal bovine serum (FBS) and 1% penicillin/streptomycin. The temperature of the incubator was controlled at 37°C and the CO 2 concentration was 5%. For the induction and inhibition of autophagy, cells were subjected to treatment separately with 0.2 μg/ml rapamycin (an autophagy activator; Sigma, St Louis, MO, USA) or 5 mM 3-methyladenine (3 M, an autophagy inhibitor; Sigma) for 24 h [20,21].

Cell transfection
The short-hairpin RNA against HAGLROS (shRNA-HAGLROS-1/2) and negative control (shRNA-NC) were constructed by GenePharma (Shanghai, China). HFWT cells were seeded in sixwell plates and transfected with these shRNAs with the help of Lipofectamine 2000 Reagent (Life Technologies Corporation, Grand Island, NY, USA) in line with the standard procedures of the vendor. 48 h later, cells were collected to make subsequent experiment studies [22].

Cell proliferation
Cell proliferation assessment was undertaken with the application of Cell counting kit-8 (CCK-8) assay [23]. After transfection with shRNA-HAGLROS-1 or shRNA-NC, HFWT cell line was cultured in 96-well plates at a density of 5 × 10 3 cells/well for various time periods (24 h, 48 h and 72 h). Then CCK-8 reagent was incubated with the plates seeded with cells at 37°C in 5% CO 2 for 2 h. The measurement of optical density values employed a microplate reader (RT-3001; Thermo Fisher Scientific, Waltham, MA, USA) under 450 nm wavelength.

Quantitative real-time polymerase chain reaction (qRT-PCR)
Extraction of total RNA from nephroblastoma cell lines and transfected cells employed TRIzol reagent (Invitrogen) and DNase I was used to remove genomic contamination. A reverse transcription for cDNA synthesis was performed by PrimeScript RT Master Mix (Takara Biomedical Technology, Beijing, China). The PCR reaction was undertaken making use of the SYBR Premix Ex Taq™ II kit (Takara, Shiga, Japan) at an ABI PRISM 7000 Sequence Detection System (Applied Biosystems, Foster City, CA, USA). The primer sequences for PCR are presented as below: HAGLROS: 5′-TGTCACCCTTAAATACCGCTCT-3′ (forward) and 5′-CTTCCTCCCACACAAATACTCC-3′ (reverse); GAPDH: 5′-GGGAAACTGTGGCGTGAT-3′ (forward) and 5′-GAGTGGGTGTCGCTGTTGA-3′ (reverse). GAPDH was worked as the control of HAGLROS and the results were measured by using the 2 −ΔΔCt method [24].

Wound healing assay
After transfection, HFWT cells were grown in 6-well plates at a density of 4 × 10 5 cells /well until 90% confluence. The cell monolayer was scratched in a straight line utilizing a white pipette tip. Washed three times using serum-free medium, the monolayer was cultured for 24 h and the area occupied by migrated cells in the scratch was evaluated by an inverted microscopy [25].

Transwell assay
The experiment of cell invasive ability of transfected HFWT cells was undertaken with the application of transwell chamber (Corning Costar, Cambridge, MA) [26]. Firstly, 0.1 mL of Matrigel was used to coat the chambers at 37°C for 1 h. After transfection, HFWT cells were suspended in DMEM without serum and placed in the upper wells. Meanwhile, the lower chamber was loaded with DMEM with 10% FBS. 24 h later, a swab was adopted to wipe off the non-migrated cells on the transwell membrane gently. Then the cells crossing the membrane were fixed, stained with 0.05% crystal violet and counted by a microscope under in five randomly chosen fields.

Cell apoptosis analysis
The treatment of HFWT cells with 0.2 μg/ml rapamycin or 5 mM 3 M was conducted for 24 h after transfected with shRNA-HAGLROS-1 or shRNA-NC. Treated differently, the cultured cells were harvested, rinsed twice with precooled PBS and resuspended in binding buffer. Then, the cells were double-labeled with Annexin V-FITC and PI apoptosis detection kits (Ribobio, Guangzhou, China) in keeping with the supplier's guidance. Cell apoptosis was examined by means of BD LSRII Flow Cytometer System (BD Biosciences) and analyzed by Flowjo Software (Tree Star, Ashland, OR, USA) [27].

Immunofluorescence staining
The fixation and permeabilization of cell sections were performed with paraformaldehyde (4%) and Triton X-100 (0.1%), respectively. Then, 5% skimmed milk was added in 0.01 M PBS for blocking the cells which were incubated with primary antibodies targeting LC3B (1 µg/ml; Abcam) overnight at 4°C. Twice-washed with PBS, the cells were subsequently cultured with Goat Anti-Rabbit IgG H&L secondary antibody (1:500; Abcam). The coverslips were then rinsed three times with PBS and visualized with the aid of a fluorescence microscope (Olympus, Tokyo, Japan) [28].

Western blot analysis
Total protein was extracted from different treated cells using RIP assay lysis buffer (Beyotime Biotechnology, Shanghai, China) supplemented with protease inhibitor and the protein concentration was measured by BCA Protein Assay Kit (Pierce, Appleton, WI). Then, a total of 30 μg protein per well were separated by 10% SDS-PAGE and transferring onto a PVDF membrane (Millipore, Billerica, MA), and the separated membrane was blocked by the use of 5% nonfat milk solution for 1 h followed by an overnight incubation with primary antibodies against Bcl-2, Bax, cleaved caspase3, LC3B, ATG7, Beclin1 and GAPDH (1:1000, Santa Cruz, USA) at 4°C. At the end of three washes, the membrane was subjected to incubation for another 1 h using horseradish peroxidase-conjugated secondary antibody. Visualization of the bands was carried out using the enhanced chemiluminescence (ECL) system (Beyotime, Shanghai, China) and quantification using Image J (NIH, USA) [29].

Statistical analysis
All data analysis was implemented through SPSS 19.0 (Chicago, IL, USA). An unpaired two-tailed Student's t-test was applied for the comparison of differences between two groups. Differences among multiple groups were analyzed using oneway ANOVA with a post hoc Bonferroni multiple comparison test. The experiment results were recorded in the way of mean ± SD. A statistical significance was identified when the P value was <0.05.

Results
In this study, we investigated the biological roles of HAGLROS and the potential mechanism in nephroblastoma HFWT cells. The data showed that the silencing of HAGLROS suppressed HFWT cellproliferation, migration and invasion. In addition, HAGLROS silencing repressed autophagy and induced cell apoptosis. Moreover, activation of autophagy inhibited HAGLROS silencing-mediated increased cell apoptosis, while the inhibition of autophagy promoted the apoptosis of HFWT cells.

HAGLROS expression is upregulated in nephroblastoma cell lines
To examine the function of lncRNA HAGLROS in nephroblastoma, the detection of HAGLROS expression in several nephroblastoma cell lines and normal human embryonic kidney cells employed qRT-PCR. As shown in Figure 1, HAGLROS mRNA expression in nephroblastoma cells was remarkably higher (vs Control). The results suggest that HAGLROS might be associated with progression of nephroblastoma.

Knockdown of HAGLROS suppressed the proliferation, migration and invasion of HFWT cells
For the exploration of the influence of HAGLROS on nephroblastoma cells, HAGLROS expression in HFWT cells after its knockdown was measured via qRT-PCR (Figure 2(a)). shRNA-HAGLROS-1 showed the better transfetion and was chosen to perform follow-up experiments. Next, we found that HAGLROS silencing inhibited HFWT cell proliferation in comparison with the control and negative control groups (Figure 2(b)). In addition, the migration and invasion capacity of HFWT cells were evidently repressed by transfection with shRNA-HAGLROS-1 (Figure 2 (c-f)). These data indicate that HAGLROS silencing exerts suppressive roles in proliferating, migrating and invading of nephroblastoma cells.

HAGLROS silencing promoted cell apoptosis in HFWT cells
To further study the effects of HAGLROS on nephroblastoma, cell apoptosis detection was undertaken. As indicated in Figure 3(a), flow cytometric analysis revealed that downregulation of HAGLROS greatly increased cell apoptosis rate by contrast to the control and negative control groups. Moreover, the results from Western blot analysis showed that Bcl-2 level was reduced but the levels of Bax and cleaved caspase 3 were significantly increased in shRNA-HAGLROS -1 group when compared with the controls (Figure 3 (b)). The data indicate that HAGLROS plays a regulatory role in cell apoptosis of nephroblastoma.

HAGLROS silencing inhibits the autophagy in nephroblastoma cells
The objective of this set of experiments was to identify the effect of HAGLROS on cell autophagy in HFWT cells. The results obtained from immunofluorescence staining showed that LC3II level was remarkably decreased in cells transfected with shRNA-HAGLROS-1 compared with transfection with shRNA-NC ( Figure  4(a)). Furthermore, HAGLROS silencing markedly enhanced the protein levels of LC3II/I, ATG7 and Beclin1 according to Western blot assay (Figure 4(b)). Thus, HAGLROS may be involved in nephroblastoma cell autophagy.

Effects of cell autophagy on the apoptosis of HFWT cells with HAGLROS silencing
Finally, we elucidated the relationship between autophagy and apoptosis in HFWT cells after transfection with shRNA-HAGLROS-1. As shown in Figure 5(a), rapamycin significantly increased the protein level of LC3II while 3 M reduced the LC3II level in HFWT cells transfected with shRNA-HAGLROS-1. Consistently, the protein levels of LC3II/I, ATG7 and Beclin1 were conspicuously elevated after treatment of rapamycin while 3 M reduced the protein levels ( Figure 5 (b)), indicating that autophagy can be regulated by rapamycin or 3 M in in HFWT cells. On the basis of flow cytometric analysis, activation of autophagy reduced cell apoptosis rate while inhibition of autophagy facilitated the apoptosis upon transfection with shRNA-HAGLROS-1 (Figure 6(a)). Moreover, rapamycin increased the Bcl-2 level and attenuated the activities of Bax and cleaved caspase 3. However, 3 M showed opposite effects on in HFWT cells after transfection with shRNA-HAGLROS-1 (Figure 6(b)). These results suggest that autophagy restrains apoptosis while autophagy inhibition promotes apoptosis after HAGLROS silencing.

Discussion
Nephroblastoma is one of the most common cancers with a significant direct impact on the survival and health for children [30]. The molecular mechanisms underlying nephroblastoma are usually found to be linked with noncoding RNAs, especially lncRNAs [31,32]. Nevertheless, the key lncRNAs that take part in the occurrence and development of nephroblastoma are still incomplete grabbed. It was found in our study that HAGLROS expressed higher in nephroblastoma cell lines. HAGLROS silencing suppressed the proliferation of nephroblastoma cell, repressed migration, invasion and autophagy, but promoted apoptosis. Moreover, a negative regulation between apoptosis and autophagy was elucidated in HFWT cells after HAGLROS was downregulated.
Multiple evidences have confirmed that many lncRNAs play crucial roles in cell biological activities including proliferation, migration and invasion in nephroblastoma [23,33]. Zhang et al. reported that SOX21-AS1 expression was markedly elevated in the tissues and cell lines of nephroblastoma. SOX21-AS1 knockdown diminished the proliferation and colony-forming ability of nephroblastoma cells and inhibited cell cycle [34]. Another study showed that SNHG6 was found to be overexpressed in nephroblastoma tissues and knockdown of SNHG6 suppressed cell proliferation, migration and incursion via elevating miR-15a expression [35]. Consistently, our study revealed that HAGLROS expression was increased in nephroblastoma cells, among which HFWT cells showed the highest HAGLROS level. In addition, we found that the proliferation, migration and invasion of HFWT cells were inhibited after downregulation of HAGLROS, suggesting the oncogenic role of HAGLROS in nephroblastoma. HAGLROS was also found to affect proliferation, migration and invasion in other types of cancers. Wei et al. revealed that HAGLROS level was significantly higher in hepatocellular carcinoma (HCC) tissues and was correlated with clinical characteristics, tumor stage or tumor differentiation of patients with HCC [18]. Chen et al. showed that HAGLROS expression was significantly elevated and related to outcomes of gastric cancer (GC) patients. Knockdown of HAGLROS considerably inhibited the cell proliferation, invasion and migration by regulation of mTOR signals [36].
Numerous studies conducted on apoptosis in recent years have demonstrated its role in the development of disease. Nevertheless, apoptosis is not the only factor that determines the destiny of a cell [37]. Autophagy has been shown to engage in complex interplay with apoptosis and proved to contribute to cell survival or death according to the different cell microenvironment [38]. The role of the double-edged sword autophagy modulation in cancers is researching continually and excessive autophagy can lead to cell death [39]. One study showed that inhibition of autophagy reduced cell activity and induced cell apoptosis in childhood nephroblastoma [40]. In our study, HAGLROS silencing induced cell apoptosis by reducing cell apoptosis rate, inhibited Bcl-2 level and enhanced the activities of Bax and cleaved caspase 3 in nephroblastoma cell line HFWT. Furthermore, we also found that downregulation of HAGLROS suppressed autophagy through inhibiting the expressions of autophagy-related proteins, including ATG7, Beclin1 and LC3II/I. For the purpose of studying the impacts of autophagy on cell apoptosis in nephroblastoma, we regulated autophagy level in HFWT cells by administration of rapamycin or 3 M. The results revealed that activating autophagy by rapamycin significantly repressed cell apoptosis while inhibition of autophagy by 3 M increased the rate of cell apoptosis. Similarly, Western blot assay showed an increased expression of Bcl-2 and reduced levels of Bax and cleaved caspase 3 after induction of autophagy whereas 3 M treatment caused the opposite phenomenon inversely, which was in line with previous results that lncRNA HAGLROS facilitated cell proliferation, suppressed apoptosis and accelerated autophagy in hepatocellular carcinoma cells through the regulation of miR-5095/ATG12 axis [18]. However, we did not perform animal experiments or examine HAGLROS expression in clinical samples in this study. In addition, we only focus on the biological functions of HAGLROS in one nephroblastoma cell line HFWT. Thus, we plan to explore roles of HAGLROS in other nephroblastoma cell lines to verify our results and perform animal and clinical trials about HAGLROS in the further experiments. In addition, we will explore the expression pattern of HAGLROS in nephroblastoma by bioinformatics analysis in further study.

Conclusion
In summary, our findings revealed that HAGLROS expression might be associated with pathogenesis of nephroblastoma. Knockdown of HAGLROS suppressed cell proliferation and metastasis, accelerated apoptosis and diminished autophagy. HAGLROS may regulate nephroblastoma cell apoptosis via mediating the activation of autophagy. Our study provides a new strategy for the targeted therapy of nephroblastoma.

Highlights
• HAGLROS is highly expressed in nephroblastoma cell lines.  and Shijie Tang performed the statistical analysis. Pugui Li and Weizhu Tang wrote the manuscript. All authors contributed to the article and approved the submitted version.

Disclosure statement
No potential conflict of interest was reported by the author(s).

Funding
The author(s) reported there is no funding associated with the work featured in this article.