Simultaneous Inhibition of SARS-CoV-2 Entry Pathways by Cyclosporine

Cyclosporine may reduce the risk of symptomatic COVID-19 in patients with systemic lupus erythematosus: a retrospective cohort study

This study aimed to explore the effect of cyclosporine (CsA) on severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in systemic lupus erythematosus (SLE) patients to provide a valuable reference for clinical treatment strategies in the context of the long-term risk of SARS-CoV-2 infection. SLE patients who visited the Rheumatology Outpatient Department of Fujian Medical University Union Hospital between 1 May and 31 October 2022 were included. Data on SARS-CoV-2 infection in patients between 1 November 2022 and 31 July 2023 were obtained by telephone follow-up. Patients were divided into two groups according to whether CsA was used during the observation period: the glucocorticoid or hydroxychloroquine group and the CsA group. To assess the robustness of the results, Data sets 1 and 2 were established to be analyzed independently. Multivariate logistic regression was used to estimate the odds ratios (ORs) and 95% confidence intervals (CIs) for symptomatic coronavirus disease (COVID-19). A total of 184 patients were included, among whom 129 were definite symptomatic COVID-19 patients; 29 were presumptive symptomatic COVID-19 patients; and 4 had signs and symptoms of COVID-19, but tested negative for SARS-CoV-2 in a virological test. According to the multivariable-adjusted models, CsA was associated with lower odds of symptomatic COVID-19 (P = 0.042, OR = 0.316, 95% CI: 0.104-0.959 in Data set 1 and P = 0.021, OR = 0.257, 95% CI: 0.081-0.812 in Data set 2). CsA is associated with lower odds of contracting symptomatic COVID-19. The use of CsA may be considered an appropriate therapeutic option for disease management in patients with rheumatic diseases who have severe disease activity and persistent SARS-CoV-2 infection.

IMPORTANCE

Our study indicated that cyclosporine may reduce the risk of symptomatic COVID-19 in systemic lupus erythematosus patients in spite of its immunosuppressive effects. This study provides a reference for clinical treatment strategies for AIIRD patients in the context of the long-term risk of SARS-CoV-2 infection.

Nsp1 facilitates SARS-CoV-2 replication through calcineurin-NFAT signaling

SARS-CoV-2, the causative agent of COVID-19, has been intensely studied in search of effective antiviral treatments. The immunosuppressant cyclosporine A (CsA) has been suggested to be a pan-coronavirus inhibitor, yet its underlying mechanism remained largely unknown. Here, we found that non-structural protein 1 (Nsp1) of SARS-CoV-2 usurped CsA-suppressed nuclear factor of activated T cells (NFAT) signaling to drive the expression of cellular DEAD-box helicase 5 (DDX5), which facilitates viral replication. Nsp1 interacted with calcineurin A (CnA) to displace the regulatory protein regulator of calcineurin 3 (RCAN3) of CnA for NFAT activation. The influence of NFAT activation on SARS-CoV-2 replication was also validated by using the Nsp1-deficient mutant virus. Calcineurin inhibitors, such as CsA and VIVIT, inhibited SARS-CoV-2 replication and exhibited synergistic antiviral effects when used in combination with nirmatrelvir. Our study delineated the molecular mechanism of CsA-mediated inhibition of SARS-CoV-2 replication and the anti-SARS-CoV-2 action of calcineurin inhibitors.

IMPORTANCE

Cyclosporine A (CsA), commonly used to inhibit immune responses, is also known to have anti-SARS-CoV-2 activity, but its mode of action remains elusive. Here, we provide a model to explain how CsA antagonizes SARS-CoV-2 through three critical proteins: DDX5, NFAT1, and Nsp1. DDX5 is a cellular facilitator of SARS-CoV-2 replication, and NFAT1 controls the production of DDX5. Nsp1 is a viral protein absent from the mature viral particle and capable of activating the function of NFAT1 and DDX5. CsA and similar agents suppress Nsp1, NFAT1, and DDX5 to exert their anti-SARS-CoV-2 activity either alone or in combination with Paxlovid.

Calcitriol ameliorates motor deficits and prolongs survival of Chrne-deficient mouse, a model for congenital myasthenic syndrome, by inducing Rspo2

Signal transduction at the neuromuscular junction (NMJ) is compromised in a diverse array of diseases including congenital myasthenic syndromes (CMS). Germline mutations in CHRNE encoding the acetylcholine receptor (AChR) ε subunit are the most common cause of CMS. An active form of vitamin D, calcitriol, binds to vitamin D receptor (VDR) and regulates gene expressions. We found that calcitriol enhanced MuSK phosphorylation, AChR clustering, and myotube twitching in co-cultured C2C12 myotubes and NSC34 motor neurons. RNA-seq analysis of co-cultured cells showed that calcitriol increased the expressions of Rspo2, Rapsn, and Dusp6. ChIP-seq of VDR revealed that VDR binds to a region approximately 15 ​kbp upstream to Rspo2. Biallelic deletion of the VDR-binding site of Rspo2 by CRISPR/Cas9 in C2C12 myoblasts/myotubes nullified the calcitriol-mediated induction of Rspo2 expression and MuSK phosphorylation. We generated Chrne knockout (Chrne KO) mouse by CRISPR/Cas9. Intraperitoneal administration of calcitriol markedly increased the number of AChR clusters, as well as the area, the intensity, and the number of synaptophysin-positive synaptic vesicles, in Chrne KO mice. In addition, calcitriol ameliorated motor deficits and prolonged survival of Chrne KO mice. In the skeletal muscle, calcitriol increased the gene expressions of Rspo2, Rapsn, and Dusp6. We propose that calcitriol is a potential therapeutic agent for CMS and other diseases with defective neuromuscular signal transmission.

Comprehensive mapping of mutations in TDP-43 and α-Synuclein that affect stability and binding

Abnormal aggregation and amyloid inclusions of TAR DNA-binding protein 43 (TDP-43) and α-Synuclein (α-Syn) are frequently co-observed in amyotrophic lateral sclerosis, Parkinson's disease, and Alzheimer's disease. Several reports showed TDP-43 C-terminal domain (CTD) and α-Syn interact with each other and the aggregates of these two proteins colocalized together in different cellular and animal models. Molecular dynamics simulation was conducted to elucidate the stability of the TDP-43 and Syn complex structure. The interfacial mutations in protein complexes changes the stability and binding affinity of the protein that may cause diseases. Here, we have utilized the computational saturation mutagenesis approach including structure-based stability and binding energy calculations to compute the systemic effects of missense mutations of TDP-43 CTD and α-Syn on protein stability and binding affinity. Most of the interfacial mutations of CTD and α-Syn were found to destabilize the protein and reduced the protein binding affinity. The results thus shed light on the functional consequences of missense mutations observed in TDP-43 associated proteinopathies and may provide the mechanisms of co-morbidities involving these two proteins.Communicated by Ramaswamy H. Sarma.

In silico identification of microRNAs targeting the PPARα/γ: promising therapeutics for SARS-CoV‑2 infection

The SARS-CoV-2 lifecycle is dependent on the host metabolism machinery. It upregulates the PPARα and PPARγ genes in lipid metabolism, which supports the essential viral replication complex including lipid rafts and palmitoylation of viral protein. The use of PPAR ligands in SARS-CoV-2 infection may have positive effects by preventing cytokine storm and the ensuing inflammatory cascade. The inhibition of PPARα and PPARγ genes may alter the metabolism and may disrupt the lifecycle of SARS-CoV-2 and COVID-19 progression. In the present work, we have identified possible miRNAs targeting PPARα and PPARγ in search of modulators of PPARα and PPARγ genes expression. The identified miRNAs could possibly be viewed as new therapeutic targets against COVID-19 infection.

Harnessing and bioprospecting botanical-based herbal medicines against potential drug targets for COVID-19: a review coupled molecular docking studies

Since the end of February 2020, the world has come to a standstill due to the virus SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2). Since then, the global scientific community has explored various remedies and treatments against this virus, including natural products that have always been a choice because of their many benefits. Various known phytochemicals are well documented for their antiviral properties. Research is being carried out to discover new natural plant products or existing ones as a treatment measure for this disease. The three important targets in this regard are-papain like protease (PLpro), spike protein, and 3 chymotrypsin like proteases (3CLpro). Various docking studies are also being elucidated to identify the phytochemicals that modulate crucial proteins of the virus. The paper is simultaneously a comprehensive review that covers recent advances in the domain of the effect of various botanically derived natural products as an alternative treatment approach against Coronavirus Disease 2019 (COVID-19). Furthermore, the docking analyses revealed that rutin (inhibitor of the major protease of SARS-CoV-2), gallocatechin (e.g., interacting with 03 hydrogen bonds with a spike-like protein), lycorine (showing the best binding affinity with amino acids GLN498, THR500 and GLY446 of the spike-like protein), and quercetrin (inhabiting at its residues ASP216, PHE219, and ILE259) are promising inhibitors of SARS‑CoV‑2.Communicated by Ramaswamy H. Sarma.

An Enhanced Dissolving Cyclosporin-A Inhalable Powder Efficiently Reduces SARS-CoV-2 Infection In Vitro

This work illustrates the development of a dry inhalation powder of cyclosporine-A for the prevention of rejection after lung transplantation and for the treatment of COVID-19. The influence of excipients on the spray-dried powder's critical quality attributes was explored. The best-performing powder in terms of dissolution time and respirability was obtained starting from a concentration of ethanol of 45% (v/v) in the feedstock solution and 20% (w/w) of mannitol. This powder showed a faster dissolution profile (Weibull dissolution time of 59.5 min) than the poorly soluble raw material (169.0 min). The powder exhibited a fine particle fraction of 66.5% and an MMAD of 2.97 µm. The inhalable powder, when tested on A549 and THP-1, did not show cytotoxic effects up to a concentration of 10 µg/mL. Furthermore, the CsA inhalation powder showed efficiency in reducing IL-6 when tested on A549/THP-1 co-culture. A reduction in the replication of SARS-CoV-2 on Vero E6 cells was observed when the CsA powder was tested adopting the post-infection or simultaneous treatment. This formulation could represent a therapeutic strategy for the prevention of lung rejection, but is also a viable approach for the inhibition of SARS-CoV-2 replication and the COVID-19 pulmonary inflammatory process.

SARS-CoV-2 Receptors and Their Involvement in Cell Infection

The new coronavirus infection (COVID-19) pandemic caused by SARS-CoV-2 has many times surpassed the epidemics caused by SARS-CoV and MERS-CoV. The reason for this was the presence of sites in the protein sequence of SARS-CoV-2 that provide interaction with a broader range of receptor proteins on the host cell surface. In this review, we consider both already known receptors common to SARS-CoV and SARS-CoV-2 and new receptors specific to SARS-CoV-2.

Network-based approach for targeting human kinases commonly associated with amyotrophic lateral sclerosis and cancer

Background

Amyotrophic Lateral Sclerosis (ALS) is a rare progressive and chronic motor neuron degenerative disease for which at present no cure is available. In recent years, multiple genes encode kinases and other causative agents for ALS have been identified. Kinases are enzymes that show pleiotropic nature and regulate different signal transduction processes and pathways. The dysregulation of kinase activity results in dramatic changes in processes and causes many other human diseases including cancers.

Methods

In this study, we have adopted a network-based system biology approach to investigate the kinase-based molecular interplay between ALS and other human disorders. A list of 62 ALS-associated-kinases was first identified and then we identified the disease associated with them by scanning multiple disease-gene interaction databases to understand the link between the ALS-associated kinases and other disorders.

Results

An interaction network with 36 kinases and 381 different disorders associated with them was prepared, which represents the complexity and the comorbidity associated with the kinases. Further, we have identified 5 miRNAs targeting the majority of the kinases in the disease-causing network. The gene ontology and pathways enrichment analysis of those miRNAs were performed to understand their biological and molecular functions along with to identify the important pathways. We also identified 3 drug molecules that can perturb the disease-causing network by drug repurposing.

Conclusion

This network-based study presented hereby contributes to a better knowledge of the molecular underpinning of comorbidities associated with the kinases associated with the ALS disease and provides the potential therapeutic targets to disrupt the highly complex disease-causing network.

Isolation and In Silico Inhibitory Potential against SARS-CoV-2 RNA Polymerase of the Rare Kaempferol 3-O-(6″-O-acetyl)-Glucoside from Calligonum tetrapterum

The phytochemical constituents of Calligonum tetrapterum Jaub. & Spach (Family Polygonaceae) were studied for the first time. The study resulted in the isolation of the rare flavonol glycoside, kaempferol 3-O-(6″-O-acetyl)-glucoside,(K3G-A). The potential inhibitive activity of K3G-A toward SARS-CoV-2 was investigated utilizing several in silico approaches. First, molecular fingerprints and structural similarity experiments were carried out for K3G-A against nine co-crystallized ligands of nine proteins of SARS-CoV-2 to reveal if there is a structural similarity with any of them. The conducted studies showed the high similarity of K3G-A and remdesivir, the co-crystallized ligand of SARS-CoV-2 RNA-dependent RNA polymerase (PDB ID: 7BV2), RdRp. To validate these findings, a DFT study was conducted and confirmed the proposed similarity on the electronic and orbital levels. The binding of K3G-A against RdRp was confirmed through molecular docking studies exhibiting a binding energy of −27.43 kcal/mol, which was higher than that of remdesivir. Moreover, the RdRp-K3G-A complex was subjected to several MD studies at 100 ns that authenticated the accurate mode of binding and the correct dynamic behavior. Finally, in silico ADMET and toxicity evaluation of K3G-A was conducted and denoted the safety and the drug-likeness of K3G-A. In addition to K3G-A, two other metabolites were isolated and identified to be kaempferol (K) and β-sitosterol (β-S).

Methodology-Centered Review of Molecular Modeling, Simulation, and Prediction of SARS-CoV-2

Despite tremendous efforts in the past two years, our understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), virus-host interactions, immune response, virulence, transmission, and evolution is still very limited. This limitation calls for further in-depth investigation. Computational studies have become an indispensable component in combating coronavirus disease 2019 (COVID-19) due to their low cost, their efficiency, and the fact that they are free from safety and ethical constraints. Additionally, the mechanism that governs the global evolution and transmission of SARS-CoV-2 cannot be revealed from individual experiments and was discovered by integrating genotyping of massive viral sequences, biophysical modeling of protein-protein interactions, deep mutational data, deep learning, and advanced mathematics. There exists a tsunami of literature on the molecular modeling, simulations, and predictions of SARS-CoV-2 and related developments of drugs, vaccines, antibodies, and diagnostics. To provide readers with a quick update about this literature, we present a comprehensive and systematic methodology-centered review. Aspects such as molecular biophysics, bioinformatics, cheminformatics, machine learning, and mathematics are discussed. This review will be beneficial to researchers who are looking for ways to contribute to SARS-CoV-2 studies and those who are interested in the status of the field.

A Drug Discovery Approach for an Effective Pain Therapy through Selective Inhibition of Nav1.7

Chronic pain is a widespread disorder affecting millions of people and is insufficiently addressed by current classes of analgesics due to significant long-term or high dosage side effects. A promising approach that was recently proposed involves the systemic inhibition of the voltage-gated sodium channel Nav1.7, capable of cancelling pain perception completely. Notwithstanding numerous attempts, currently no drugs have been approved for the inhibition of Nav1.7. The task is complicated by the difficulty of creating a selective drug for Nav1.7, and avoiding binding to the many human paralogs performing fundamental physiological functions. In our work, we obtained a promising set of ligands with up to 5-40-fold selectivity and reaching 5.2 nanomolar binding affinity by employing a proper treatment of the problem and an innovative differential in silico screening procedure to discriminate for affinity and selectivity against the Nav paralogs. The absorption, distribution, metabolism, and excretion (ADME) properties of our top-scoring ligands were also evaluated, with good to excellent results. Additionally, our study revealed that the top-scoring ligand is a stereoisomer of an already-approved drug. These facts could reduce the time required to bring a new effective and selective Nav1.7 inhibitor to the market.

An inductive graph neural network model for compound-protein interaction prediction based on a homogeneous graph

Identifying the potential compound–protein interactions (CPIs) plays an essential role in drug development. The computational approaches for CPI prediction can reduce time and costs of experimental methods and have benefited from the continuously improved graph representation learning. However, most of the network-based methods use heterogeneous graphs, which is challenging due to their complex structures and heterogeneous attributes. Therefore, in this work, we transformed the compound–protein heterogeneous graph to a homogeneous graph by integrating the ligand-based protein representations and overall similarity associations. We then proposed an Inductive Graph AggrEgator-based framework, named CPI-IGAE, for CPI prediction. CPI-IGAE learns the low-dimensional representations of compounds and proteins from the homogeneous graph in an end-to-end manner. The results show that CPI-IGAE performs better than some state-of-the-art methods. Further ablation study and visualization of embeddings reveal the advantages of the model architecture and its role in feature extraction, and some of the top ranked CPIs by CPI-IGAE have been validated by a review of recent literature. The data and source codes are available at https://github.com/wanxiaozhe/CPI-IGAE.

Cyclosporine A Inhibits Viral Infection and Release as Well as Cytokine Production in Lung Cells by Three SARS-CoV-2 Variants

In December 2019, a new severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) started spreading worldwide causing the coronavirus disease 2019 (COVID-19) pandemic. The hyperactivation of the immune system has been proposed to account for disease severity and death in COVID-19 patients. Despite several approaches having been tested, no therapeutic protocol has been approved. Given that Cyclosporine A (CsA) is well-known to exert a strong antiviral activity on several viral strains and an anti-inflammatory role in different organs with relevant benefits in diverse pathological contexts, we tested its effects on SARS-CoV-2 infection of lung cells. We found that treatment with CsA either before or after infection of CaLu3 cells by three SARS-CoV-2 variants: (i) reduces the expression of both viral RNA and proteins in infected cells; (ii) decreases the number of progeny virions released by infected cells; (iii) dampens the virus-triggered synthesis of cytokines (including IL-6, IL-8, IL1α and TNF-α) that are involved in cytokine storm in patients. Altogether, these data provide a rationale for CsA repositioning for the treatment of severe COVID-19 patients.

IMPORTANCE SARS-CoV-2 is the most recently identified member of the betacoronavirus genus responsible for the COVID-19 pandemic. Repurposing of available drugs has been a “quick and dirty” approach to try to reduce mortality and severe symptoms in affected patients initially, and can still represent an undeniable and valuable approach to face COVID-19 as the continuous appearance and rapid diffusion of more “aggressive”/transmissible variants, capable of eluding antibody neutralization, challenges the effectiveness of some anti-SARS-CoV-2 vaccines. Here, we tested a known antiviral and anti-inflammatory drug, Cyclosporine A (CsA), and found that it dampens viral infection and cytokine release from lung cells upon exposure to three different SARS-CoV-2 variants. Knock down of the main intracellular target of CsA, Cyclophilin A, does not phenocopy the drug inhibition of viral infection. Altogether, these findings shed new light on the cellular mechanisms of SARS-CoV-2 infection and provide the rationale for CsA repositioning to treat severe COVID-19 patients.

Insights into the SARS-CoV-2-Mediated Alteration in the Stress Granule Protein Regulatory Networks in Humans

The rapidly and constantly evolving coronavirus, SARS-CoV-2, imposes a great threat to human health causing severe lung disease and significant mortality. Cytoplasmic stress granules (SGs) exert anti-viral activities due to their involvement in translation inhibition and innate immune signaling. SARS-CoV-2 sequesters important SG nucleator proteins and impairs SG formation, thus evading the host response for efficient viral replication. However, the significance of SGs in COVID-19 infection remains elusive. In this study, we utilize a protein-protein interaction network approach to systematically dissect the crosstalk of human post-translational regulatory networks governed by SG proteins due to SARS-CoV-2 infection. We uncovered that 116 human SG proteins directly interact with SARS-CoV-2 proteins and are involved in 430 different brain disorders including COVID-19. Further, we performed gene set enrichment analysis to identify the drugs against three important key SG proteins (DYNC1H1, DCTN1, and LMNA) and also looked for potential microRNAs (miRNAs) targeting these proteins. We identified bexarotene as a potential drug molecule and miRNAs, hsa-miR-615-3p, hsa-miR-221-3p, and hsa-miR-124-3p as potential candidates for the treatment of COVID-19 and associated manifestations.

Targeting cathepsins: A potential link between COVID-19 and associated neurological manifestations

Many studies have shown that the lysosomal cathepsins, especially cathepsins B/L (CTSB/L) are required for SARS-CoV-2 entry into host cells. Lysosomal proteases, cathepsins are indispensable for normal health and are involved in several brain disorders occurring at different development age periods. On the other hand, it has been well known that COVID-19 infection is largely associated with several neurological disorders. Taken together these findings and given the high levels of expression of CTSB/L in the brain, we here proposed a reasonable hypothesis about the involvement of CTSB/L in the neurological manifestations linked to COVID-19. Pharmacological inhibitions of the CTSB/L could be a potential therapeutic target to block the virus entry as well as to mitigate the brain disorders. To this end, we utilized the network-based drug repurposing analyses to identify the possible drugs that can target CTSB/L. This study identifies the molecules like cyclosporine, phenytoin, and paclitaxel as potential drugs with binding ability to the CTSB/L. Further, we have performed molecular docking and all-atom molecular dynamics (MD) simulations to investigate the stability of CTSL-drug complexes. The results showed strong and stable binding of drugs with CTSL.

Long-Term Clinical and Immunological Impact of Severe COVID-19 on a Living Kidney Transplant Recipient - A Case Report

The long-term impact of COVID-19 on transplant recipients remains unknown. We describe the case of a 30-year-old male kidney transplant recipient from Wuhan, China that was treated for severe COVID-19 in February 2020. He suffered an acute lung and renal injury and required systemic treatment including adjustment of his immunosuppressant regime. He was followed up to 1-year after discharge. No chronic lung fibrosis or deterioration of his pulmonary function was observed. Despite COVID-19 mediated damage to his renal tubular cells, no transplant rejection occurred. His immunological profile demonstrated both cellular anti-SARS-CoV-2 reactivity and specific humoral immunity, indicating that it is beneficial for the transplanted patients to be immunized with SARS-CoV-2 virus vaccine. This case will help guide clinical decision making for immunocompromised individuals that become infected with SARS-CoV-2.

Genomics-guided identification of potential modulators of SARS-CoV-2 entry proteases, TMPRSS2 and Cathepsins B/L

SARS-CoV-2 requires serine protease, transmembrane serine protease 2 (TMPRSS2), and cysteine proteases, cathepsins B, L (CTSB/L) for entry into host cells. These host proteases activate the spike protein and enable SARS-CoV-2 entry. We herein performed genomic-guided gene set enrichment analysis (GSEA) to identify upstream regulatory elements altering the expression of TMPRSS2 and CTSB/L. Further, medicinal compounds were identified based on their effects on gene expression signatures of the modulators of TMPRSS2 and CTSB/L genes. Using this strategy, estradiol and retinoic acid have been identified as putative SARS-CoV-2 alleviation agents. Next, we analyzed drug-gene and gene-gene interaction networks using 809 human targets of SARS-CoV-2 proteins. The network results indicate that estradiol interacts with 370 (45%) and retinoic acid interacts with 251 (31%) human proteins. Interestingly, a combination of estradiol and retinoic acid interacts with 461 (56%) of human proteins, indicating the therapeutic benefits of drug combination therapy. Finally, molecular docking analysis suggests that both the drugs bind to TMPRSS2 and CTSL with the nanomolar to low micromolar affinity. The results suggest that these drugs can simultaneously target both the entry pathways of SARS-CoV-2 and thus can be considered as a potential treatment option for COVID-19.

Atomic Resolution Homology Models and Molecular Dynamics Simulations of Plasmodium falciparum Tubulins

Microtubules are tubulin polymers present in the eukaryotic cytoskeleton essential for structural stability and cell division that are also roadways for intracellular transport of vesicles and organelles. In the human malaria parasite Plasmodium falciparum, apart from providing structural stability and cell division, microtubules also facilitate important biological activities crucial for parasite survival in hosts, such as egression and motility. Hence, parasite structures and processes involving microtubules are among the most important drug targets for discovering much-needed novel Plasmodium inhibitors. The current study aims to construct reliable and high-quality 3D models of α-, β-, and γ-tubulins using various modeling techniques. We identified a common binding pocket specific to Plasmodium α-, β-, and γ-tubulins. Molecular dynamics simulations confirmed the stability of the Plasmodium tubulin 3D structures. The models generated in the present study may be used for protein-protein and protein-drug interaction investigations targeted toward designing malaria parasite tubulin-specific inhibitors.

SARS-CoV-2 Infects Hamster Testes

The COVID-19 pandemic continues to affect millions of people worldwide. Although SARS-CoV-2 is a respiratory virus, there is growing concern that the disease could cause damage and pathology outside the lungs, including in the genital tract. Studies suggest that SARS-CoV-2 infection can damage the testes and reduce testosterone levels, but the underlying mechanisms are unknown and evidence of virus replication in testicular cells is lacking. We infected golden Syrian hamsters intranasally, a model for mild human COVID-19, and detected viral RNA in testes samples without histopathological changes up to one month post-infection. Using an ex vivo infection model, we detected SARS-CoV-2 replication in hamster testicular cells. Taken together, our data raise the possibility that testes damage observed in severe cases of COVID-19 could be partly explained by direct SARS-CoV-2 infection of the testicular cells.