Dromedary camel nanobodies broadly neutralize SARS-CoV-2 variants.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike is a trimer of S1/S2 heterodimers with three receptor-binding domains (RBDs) at the S1 subunit for human angiotensin-converting enzyme 2 (hACE2). Due to their small size, nanobodies can recognize protein cavities that are not accessible to conventional antibodies. To isolate high-affinity nanobodies, large libraries with great diversity are highly desirable. Dromedary camels (Camelus dromedarius) are natural reservoirs of coronaviruses like Middle East respiratory syndrome CoV (MERS-CoV) that are transmitted to humans. Here, we built large dromedary camel VHH phage libraries to isolate nanobodies that broadly neutralize SARS-CoV-2 variants. We isolated two VHH nanobodies, NCI-CoV-7A3 (7A3) and NCI-CoV-8A2 (8A2), which have a high affinity for the RBD via targeting nonoverlapping epitopes and show broad neutralization activity against SARS-CoV-2 and its emerging variants of concern. Cryoelectron microscopy (cryo-EM) complex structures revealed that 8A2 binds the RBD in its up mode with a long CDR3 loop directly involved in the ACE2 binding residues and that 7A3 targets a deeply buried region that uniquely extends from the S1 subunit to the apex of the S2 subunit regardless of the conformational state of the RBD. At a dose of ≥5 mg/kg, 7A3 efficiently protected transgenic mice expressing hACE2 from the lethal challenge of variants B.1.351 or B.1.617.2, suggesting its therapeutic use against COVID-19 variants. The dromedary camel VHH phage libraries could be helpful as a unique platform ready for quickly isolating potent nanobodies against future emerging viruses.

Conformational energy range of ligands in protein crystal structures: The difficult quest for accurate understanding.

In this review, we address a fundamental question: What is the range of conformational energies seen in ligands in protein-ligand crystal structures? This value is important biophysically, for better understanding the protein-ligand binding process; and practically, for providing a parameter to be used in many computational drug design methods such as docking and pharmacophore searches. We synthesize a selection of previously reported conflicting results from computational studies of this issue and conclude that high ligand conformational energies really are present in some crystal structures. The main source of disagreement between different analyses appears to be due to divergent treatments of electrostatics and solvation. At the same time, however, for many ligands, a high conformational energy is in error, due to either crystal structure inaccuracies or incorrect determination of the reference state. Aside from simple chemistry mistakes, we argue that crystal structure error may mainly be because of the heuristic weighting of ligand stereochemical restraints relative to the fit of the structure to the electron density. This problem cannot be fixed with improvements to electron density fitting or with simple ligand geometry checks, though better metrics are needed for evaluating ligand and binding site chemistry in addition to geometry during structure refinement. The ultimate solution for accurately determining ligand conformational energies lies in ultrahigh-resolution crystal structures that can be refined without restraints.

New insights into the enzymatic mechanism of human chitotriosidase (CHIT1) catalytic domain by atomic resolution X-ray diffraction and hybrid QM/MM.

Chitotriosidase (CHIT1) is a human chitinase belonging to the highly conserved glycosyl hydrolase family 18 (GH18). GH18 enzymes hydrolyze chitin, an N-acetylglucosamine polymer synthesized by lower organisms for structural purposes. Recently, CHIT1 has attracted attention owing to its upregulation in immune-system disorders and as a marker of Gaucher disease. The 39 kDa catalytic domain shows a conserved cluster of three acidic residues, Glu140, Asp138 and Asp136, involved in the hydrolysis reaction. Under an excess concentration of substrate, CHIT1 and other homologues perform an additional activity, transglycosylation. To understand the catalytic mechanism of GH18 chitinases and the dual enzymatic activity, the structure and mechanism of CHIT1 were analyzed in detail. The resolution of the crystals of the catalytic domain was improved from 1.65 Š(PDB entry 1waw) to 0.95-1.10 Šfor the apo and pseudo-apo forms and the complex with chitobiose, allowing the determination of the protonation states within the active site. This information was extended by hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. The results suggest a new mechanism involving changes in the conformation and protonation state of the catalytic triad, as well as a new role for Tyr27, providing new insights into the hydrolysis and transglycosylation activities.

Establishment and Characterization of Amitrole-Induced Mouse Thyroid Adenomatous Nodule-Derived Cell Lines.

Background: Thyroid cancer cell lines have been of great value for the study of thyroid cancer. However, the availability of benign thyroid adenoma cell lines is limited. Methods: Cell lines were established from thyroid adenomatous nodules that developed in mice treated with the goitrogen amitrole. Expression of epithelial, mesenchymal, and thyroid markers of these established cell lines was determined, and the effect of lentivirus-transduced overexpression of NKX2-1, a master regulator of thyroid development, on the thyroid marker expression was examined. Signal transduction and cell proliferation were evaluated after treatment with insulin-like growth factor-I (IGF-I) and the selective IGF-I receptor (IGF-IR) inhibitor NVP-ADW742. Xenograft studies were performed to examine tumorigenicity of the cells in mice. Whole-genome sequencing (WGS) was used to comprehensively determine the genetic mutations in the established two cell lines. Results: Five mouse thyroid adenomatous nodules-derived cell lines named CAT (cells from amitrole-treated thyroids) were established. Among these, two cell lines, CAT458/458s (CAT458s: a subline of CAT458) and CAT459, were found to be positive for epithelial markers and negative for a mesenchymal marker. NKX2-1-positive CAT459 cells showed higher messenger RNA (mRNA) expression of some thyroid differentiation markers than NKX2-1-negative CAT458s cells, and NKX2-1 overexpression increased and/or induced their expression. IGF-I signaling was transduced in thyrotropin receptor (Tshr)-negative CAT458s and 459 cells, and NVP-ADW742 suppressed their proliferation. No tumors developed in mice after subcutaneous injection of CAT458s or 459 cells. The WGS analysis revealed the presence of missense mutations in the tumor suppressor genes such as Polk (encoding DNA polymerase kappa) and Tgfb1 (encoding transforming growth factor beta 1), while no mutations were found in the prominent thyroid cancer-related genes Braf, Trp53 (encoding p53), and Tert (encoding telomerase reverse transcriptase). Conclusions: Two mouse thyroid adenomatous nodule-derived cell lines with different thyroid differentiation marker expression were established. NKX2-1 induced partial differentiation of these cell lines. They lacked tumorigenicity and prominent gene mutations involved in thyroid cancer development, while missense mutations were found in some tumor suppressors as revealed by WGS. The CAT458s and 459 provide a new tool to further clarify the process of thyroid multistep carcinogenesis and differentiation.

Generation of antitumor chimeric antigen receptors incorporating T cell signaling motifs.

Chimeric antigen receptor (CAR) T cells have been used to successfully treat various blood cancers, but adverse effects have limited their potential. Here, we developed chimeric adaptor proteins (CAPs) and CAR tyrosine kinases (CAR-TKs) in which the intracellular ζ T cell receptor (TCRζ) chain was replaced with intracellular protein domains to stimulate signaling downstream of the TCRζ chain. CAPs contain adaptor domains and the kinase domain of ZAP70, whereas CAR-TKs contain only ZAP70 domains. We hypothesized that CAPs and CAR-TKs would be more potent than CARs because they would bypass both the steps that define the signaling threshold of TCRζ and the inhibitory regulation of upstream molecules. CAPs were too potent and exhibited high tonic signaling in vitro. In contrast, CAR-TKs exhibited high antitumor efficacy and significantly enhanced long-term tumor clearance in leukemia-bearing NSG mice as compared with the conventional CD19-28ζ-CAR-T cells. CAR-TKs were activated in a manner independent of the kinase Lck and displayed slower phosphorylation kinetics and prolonged signaling compared with the 28ζ-CAR. Lck inhibition attenuated CAR-TK cell exhaustion and improved long-term function. The distinct signaling properties of CAR-TKs may therefore be harnessed to improve the in vivo efficacy of T cells engineered to express an antitumor chimeric receptor.

The in-silico evaluation of important GLUT9 residue for uric acid transport based on renal hypouricemia type 2.

Uric acid is the end product of purine metabolism. Uric acid transporters in the renal proximal tubule plays a key role in uric acid transport. Functional abnormalities in these transporters could lead to high or low levels of uric acid in the blood plasma, known as hyperuricemia and hypouricemia, respectively. GLUT9 has been reported as a key transporter for uric acid reuptake in renal proximal tubule. GLUT9 mutation is known as causal gene for renal hypouricemia due to defective uric acid uptake, with more severe cases resulting in urolithiasis and exercise induced acute kidney injury (EIAKI). However, the effect of mutation is not fully investigated and hard to predict the change of binding affinity. We comprehensively described the effect of GLUT9 mutation for uric acid transport using molecular dynamics and investigated the specific site for uric acid binding differences. R171C and R380W showed the significant disruption of the structure not affecting transport dynamics whereas L75R, G216R, N333S, and P412R showed the reduced affinity of the extracellular vestibular area towards urate. Interestingly, T125 M showed a significant increase in intracellular binding energy, associated with distorted geometries. We can use this classification to consider the effect mutations by comparing the transport profiles of mutants against those of chemical candidates for transport and providing new perspectives to urate lowering drug discovery using GLUT9.

The IgG4 hinge with CD28 transmembrane domain improves VHH-based CAR T cells targeting a membrane-distal epitope of GPC1 in pancreatic cancer.

Heterogeneous antigen expression is a key barrier influencing the activity of chimeric antigen receptor (CAR) T cells in solid tumors. Here, we develop CAR T cells targeting glypican-1 (GPC1), an oncofetal antigen expressed in pancreatic cancer. We report the generation of dromedary camel VHH nanobody (D4)-based CAR T cells targeting GPC1 and the optimization of the hinge (H) and transmembrane domain (TM) to improve activity. We find that a structurally rigid IgG4H and CD28TM domain brings the two D4 fragments in proximity, driving CAR dimerization and leading to enhanced T-cell signaling and tumor regression in pancreatic cancer models with low antigen density in female mice. Furthermore, single-cell-based proteomic and transcriptomic analysis of D4-IgG4H-CD28TM CAR T cells reveals specific genes (e.g., HMGB1) associated with high T-cell polyfunctionality. This study demonstrates the potential of VHH-based CAR T for pancreatic cancer therapy and provides an engineering strategy for developing potent CAR T cells targeting membrane-distal epitopes.

Anopheles salivary apyrase regulates blood meal hemostasis and drives malaria parasite transmission.

Mosquito salivary proteins play a crucial role in regulating hemostatic responses at the bite site during blood feeding. In this study, we investigate the function of Anopheles gambiae salivary apyrase (AgApyrase) in Plasmodium transmission. Our results demonstrate that salivary apyrase interacts with and activates tissue plasminogen activator, facilitating the conversion of plasminogen to plasmin, a human protein previously shown to be required for Plasmodium transmission. Microscopy imaging shows that mosquitoes ingest a substantial amount of apyrase during blood feeding which reduces coagulation in the blood meal by enhancing fibrin degradation and inhibiting platelet aggregation. Supplementation of Plasmodium infected blood with apyrase significantly enhanced Plasmodium infection in the mosquito midgut. In contrast, AgApyrase immunization inhibited Plasmodium mosquito infection and sporozoite transmission. This study highlights a pivotal role for mosquito salivary apyrase for regulation of hemostasis in the mosquito blood meal and for Plasmodium transmission to mosquitoes and to the mammal host, underscoring the potential for new strategies to prevent malaria transmission.

On the role of the SP1 domain in HIV-1 particle assembly: a molecular switch?

Expression of a retroviral protein, Gag, in mammalian cells is sufficient for assembly of immature virus-like particles (VLPs). VLP assembly is mediated largely by interactions between the capsid (CA) domains of Gag molecules but is facilitated by binding of the nucleocapsid (NC) domain to nucleic acid. We have investigated the role of SP1, a spacer between CA and NC in HIV-1 Gag, in VLP assembly. Mutational analysis showed that even subtle changes in the first 4 residues of SP1 destroy the ability of Gag to assemble correctly, frequently leading to formation of tubes or other misassembled structures rather than proper VLPs. We also studied the conformation of the CA-SP1 junction region in solution, using both molecular dynamics simulations and circular dichroism. Consonant with nuclear magnetic resonance (NMR) studies from other laboratories, we found that SP1 is nearly unstructured in aqueous solution but undergoes a concerted change to an α-helical conformation when the polarity of the environment is reduced by addition of dimethyl sulfoxide (DMSO), trifluoroethanol, or ethanol. Remarkably, such a coil-to-helix transition is also recapitulated in an aqueous medium at high peptide concentrations. The exquisite sensitivity of SP1 to mutational changes and its ability to undergo a concentration-dependent structural transition raise the possibility that SP1 could act as a molecular switch to prime HIV-1 Gag for VLP assembly. We suggest that changes in the local environment of SP1 when Gag oligomerizes on nucleic acid might trigger this switch.

PDB ligand conformational energies calculated quantum-mechanically.

We present here a greatly updated version of an earlier study on the conformational energies of protein-ligand complexes in the Protein Data Bank (PDB) [Nicklaus et al. Bioorg. Med. Chem. 1995, 3, 411-428], with the goal of improving on all possible aspects such as number and selection of ligand instances, energy calculations performed, and additional analyses conducted. Starting from about 357,000 ligand instances deposited in the 2008 version of the Ligand Expo database of the experimental 3D coordinates of all small-molecule instances in the PDB, we created a "high-quality" subset of ligand instances by various filtering steps including application of crystallographic quality criteria and structural unambiguousness. Submission of 640 Gaussian 03 jobs yielded a set of about 415 successfully concluded runs. We used a stepwise optimization of internal degrees of freedom at the DFT level of theory with the B3LYP/6-31G(d) basis set and a single-point energy calculation at B3LYP/6-311++G(3df,2p) after each round of (partial) optimization to separate energy changes due to bond length stretches vs bond angle changes vs torsion changes. Even for the most "conservative" choice of all the possible conformational energies-the energy difference between the conformation in which all internal degrees of freedom except torsions have been optimized and the fully optimized conformer-significant energy values were found. The range of 0 to ~25 kcal/mol was populated quite evenly and independently of the crystallographic resolution. A smaller number of "outliers" of yet higher energies were seen only at resolutions above 1.3 Å. The energies showed some correlation with molecular size and flexibility but not with crystallographic quality metrics such as the Cruickshank diffraction-component precision index (DPI) and R(free)-R, or with the ligand instance-specific metrics such as occupancy-weighted B-factor (OWAB), real-space R factor (RSR), and real-space correlation coefficient (RSCC). We repeated these calculations with the solvent model IEFPCM, which yielded energy differences that were generally somewhat lower than the corresponding vacuum results but did not produce a qualitatively different picture. Torsional sampling around the crystal conformation at the molecular mechanics level using the MMFF94s force field typically led to an increase in energy.

Identification of a dysfunctional exon-skipping splice variant in GLUT9/SLC2A9 causal for renal hypouricemia type 2.

Renal hypouricemia (RHUC) is a pathological condition characterized by extremely low serum urate and overexcretion of urate in the kidney; this inheritable disorder is classified into type 1 and type 2 based on causative genes encoding physiologically-important urate transporters, URAT1 and GLUT9, respectively; however, research on RHUC type 2 is still behind type 1. We herein describe a typical familial case of RHUC type 2 found in a Slovak family with severe hypouricemia and hyperuricosuria. Via clinico-genetic analyses including whole exome sequencing and in vitro functional assays, we identified an intronic GLUT9 variant, c.1419+1G>A, as the causal mutation that could lead the expression of p.Gly431GlufsTer28, a functionally-null variant resulting from exon 11 skipping. The causal relationship was also confirmed in another unrelated Macedonian family with mild hypouricemia. Accordingly, non-coding regions should be also kept in mind during genetic diagnosis for hypouricemia. Our findings provide a better pathogenic understanding of RHUC and pathophysiological importance of GLUT9.

Supramolecular complexes of quantum dots and a polyamidoamine (PAMAM)-folate derivative for molecular imaging of cancer cells.

Polyamidoamine (PAMAM) dendrimers and water-soluble 3-mercaptopropionic acid (MPA)-capped CdSe quantum dots (QDs) were combined to produce a new gel containing supramolecular complexes of QDs/PAMAM dendrimers. The formation of the QDs/PAMAM supramolecular complexes was confirmed by high resolution electron microscopy and Fourier transform infrared (FTIR) analyses. Molecular dynamics simulations corroborated the structure of the new QDs/PAMAM-based supramolecular compound. Finally, on the basis of the prominent fluorescent properties of the supramolecular complexes, PAMAM dendrimer was functionalized with folic acid to produce a new QDs/PAMAM-folate derivative that showed an efficient and selective performance as a marker for gastric cancer cells.

Nanoinformatics: an emerging area of information technology at the intersection of bioinformatics, computational chemistry and nanobiotechnology.

After the progress made during the genomics era, bioinformatics was tasked with supporting the flow of information generated by nanobiotechnology efforts. This challenge requires adapting classical bioinformatic and computational chemistry tools to store, standardize, analyze, and visualize nanobiotechnological information. Thus, old and new bioinformatic and computational chemistry tools have been merged into a new sub-discipline: nanoinformatics. This review takes a second look at the development of this new and exciting area as seen from the perspective of the evolution of nanobiotechnology applied to the life sciences. The knowledge obtained at the nano-scale level implies answers to new questions and the development of new concepts in different fields. The rapid convergence of technologies around nanobiotechnologies has spun off collaborative networks and web platforms created for sharing and discussing the knowledge generated in nanobiotechnology. The implementation of new database schemes suitable for storage, processing and integrating physical, chemical, and biological properties of nanoparticles will be a key element in achieving the promises in this convergent field. In this work, we will review some applications of nanobiotechnology to life sciences in generating new requirements for diverse scientific fields, such as bioinformatics and computational chemistry.

Quantum model of catalysis based on a mobile proton revealed by subatomic x-ray and neutron diffraction studies of h-aldose reductase.

We present results of combined studies of the enzyme human aldose reductase (h-AR, 36 kDa) using single-crystal x-ray data (0.66 A, 100K; 0.80 A, 15K; 1.75 A, 293K), neutron Laue data (2.2 A, 293K), and quantum mechanical modeling. These complementary techniques unveil the internal organization and mobility of the hydrogen bond network that defines the properties of the catalytic engine, explaining how this promiscuous enzyme overcomes the simultaneous requirements of efficiency and promiscuity offering a general mechanistic view for this class of enzymes.

Characterization of a Compound Heterozygous SLC2A9 Mutation That Causes Hypouricemia.

Renal hypouricemia is a rare genetic disorder. Hypouricemia can present as renal stones or exercise-induced acute renal failure, but most cases are asymptomatic. Our previous study showed that two recessive variants of SLC22A12 (p.Trp258*, pArg90His) were identified in 90% of the hypouricemia patients from two independent cohorts: the Korean genome and epidemiology study (KoGES) and the Korean Cancer Prevention Study (KCPS-II). In this work, we investigate the genetic causes of hypouricemia in the rest of the 10% of unsolved cases. We found a novel non-synonymous mutation of SLC2A9 (voltage-sensitive uric acid transporter) in the whole-exome sequencing (WES) results. Molecular dynamics prediction suggests that the novel mutation p.Met126Val in SLCA9b (p.Met155Val in SLC2A9a) hinders uric acid transport through a defect of the outward open geometry. Molecular analysis using Xenopus oocytes confirmed that the p.Met126Val mutation significantly reduced uric acid transport but does not affect the SLC2A9 protein expression level. Our results will shed light on a better understanding of SLC2A9-mediated uric acid transport and the development of a uric acid-lowering agent.

Camel nanobodies broadly neutralize SARS-CoV-2 variants.

With the emergence of SARS-CoV-2 variants, there is urgent need to develop broadly neutralizing antibodies. Here, we isolate two V H H nanobodies (7A3 and 8A2) from dromedary camels by phage display, which have high affinity for the receptor-binding domain (RBD) and broad neutralization activities against SARS-CoV-2 and its emerging variants. Cryo-EM complex structures reveal that 8A2 binds the RBD in its up mode and 7A3 inhibits receptor binding by uniquely targeting a highly conserved and deeply buried site in the spike regardless of the RBD conformational state. 7A3 at a dose of ≥5 mg/kg efficiently protects K18-hACE2 transgenic mice from the lethal challenge of B.1.351 or B.1.617.2, suggesting that the nanobody has promising therapeutic potentials to curb the COVID-19 surge with emerging SARS-CoV-2 variants. ONE-SENTENCE SUMMARY: Dromedary camel ( Camelus dromedarius ) V H H phage libraries were built for isolation of the nanobodies that broadly neutralize SARS-CoV-2 variants.

TNIK Is a Therapeutic Target in Lung Squamous Cell Carcinoma and Regulates FAK Activation through Merlin.

Lung squamous cell carcinoma (LSCC) is the second most prevalent type of lung cancer. Despite extensive genomic characterization, no targeted therapies are approved for the treatment of LSCC. Distal amplification of the 3q chromosome is the most frequent genomic alteration in LSCC, and there is an urgent need to identify efficacious druggable targets within this amplicon. We identify the protein kinase TNIK as a therapeutic target in LSCC. TNIK is amplified in approximately 50% of LSCC cases. TNIK genetic depletion or pharmacologic inhibition reduces the growth of LSCC cells in vitro and in vivo. In addition, TNIK inhibition showed antitumor activity and increased apoptosis in established LSCC patient-derived xenografts. Mechanistically, we identified the tumor suppressor Merlin/NF2 as a novel TNIK substrate and showed that TNIK and Merlin are required for the activation of focal adhesion kinase. In conclusion, our data identify targeting TNIK as a potential therapeutic strategy in LSCC. SIGNIFICANCE: Targeted therapies have not yet been approved for the treatment of LSCC, due to lack of identification of actionable cancer drivers. We define TNIK catalytic activity as essential for maintaining LSCC viability and validate the antitumor efficacy of TNIK inhibition in preclinical models of LSCC.This article is highlighted in the In This Issue feature, p. 1307.

CAR T cells targeting tumor-associated exons of glypican 2 regress neuroblastoma in mice.

Targeting solid tumors must overcome several major obstacles, in particular, the identification of elusive tumor-specific antigens. Here, we devise a strategy to help identify tumor-specific epitopes. Glypican 2 (GPC2) is overexpressed in neuroblastoma. Using RNA sequencing (RNA-seq) analysis, we show that exon 3 and exons 7-10 of GPC2 are expressed in cancer but are minimally expressed in normal tissues. Accordingly, we discover a monoclonal antibody (CT3) that binds exons 3 and 10 and visualize the complex structure of CT3 and GPC2 by electron microscopy. The potential of this approach is exemplified by designing CT3-derived chimeric antigen receptor (CAR) T cells that regress neuroblastoma in mice. Genomic sequencing of T cells recovered from mice reveals the CAR integration sites that may contribute to CAR T cell proliferation and persistence. These studies demonstrate how RNA-seq data can be exploited to help identify tumor-associated exons that can be targeted by CAR T cell therapies.

Using iCn3D and the World Wide Web for structure-based collaborative research: Analyzing molecular interactions at the root of COVID-19.

The COVID-19 pandemic took us ill-prepared and tackling the many challenges it poses in a timely manner requires world-wide collaboration. Our ability to study the SARS-COV-2 virus and its interactions with its human host in molecular terms efficiently and collaboratively becomes indispensable and mission-critical in the race to develop vaccines, drugs, and neutralizing antibodies. There is already a significant corpus of 3D structures related to SARS and MERS coronaviruses, and the rapid generation of new structures demands the use of efficient tools to expedite the sharing of structural analyses and molecular designs and convey them in their native 3D context in sync with sequence data and annotations. We developed iCn3D (pronounced "I see in 3D")1 to take full advantage of web technologies and allow scientists of different backgrounds to perform and share sequence-structure analyses over the Internet and engage in collaborations through a simple mechanism of exchanging "lifelong" web links (URLs). This approach solves the very old problem of "sharing of molecular scenes" in a reliable and convenient manner. iCn3D links are sharable over the Internet and make data and entire analyses findable, accessible, and reproducible, with various levels of interoperability. Links and underlying data are FAIR2 and can be embedded in preprints and papers, bringing a 3D live and interactive dimension to a world of text and static images used in current publications, eliminating at the same time the need for arcane supplemental materials. This paper exemplifies iCn3D capabilities in visualization, analysis, and sharing of COVID-19 related structures, sequence variability, and molecular interactions.

Analysis of SARS-CoV-2 ORF3a structure reveals chloride binding sites.

SARS-CoV-2 ORF3a is believed to form ion channels, which may be involved in the modulation of virus release, and has been implicated in various cellular processes like the up-regulation of fibrinogen expression in lung epithelial cells, downregulation of type 1 interferon receptor, caspase-dependent apoptosis, and increasing IFNAR1 ubiquitination. ORF3a assemblies as homotetramers, which are stabilized by residue C133. A recent cryoEM structure of a homodimeric complex of ORF3a has been released. A lower-resolution cryoEM map of the tetramer suggests two dimers form it, arranged side by side. The dimer's cryoEM structure revealed that each protomer contains three transmembrane helices arranged in a clockwise configuration forming a six helices transmembrane domain. This domain's potential permeation pathway has six constrictions narrowing to about 1 Å in radius, suggesting the structure solved is in a closed or inactivated state. At the cytosol end, the permeation pathway encounters a large and polar cavity formed by multiple beta strands from both protomers, which opens to the cytosolic milieu. We modeled the tetramer following the arrangement suggested by the low-resolution tetramer cryoEM map. Molecular dynamics simulations of the tetramer embedded in a membrane and solvated with 0.5 M of KCl were performed. Our simulations show the cytosolic cavity is quickly populated by both K+ and Cl-, yet with different dynamics. K+ ions moved relatively free inside the cavity without forming proper coordination sites. In contrast, Cl- ions enter the cavity, and three of them can become stably coordinated near the intracellular entrance of the potential permeation pathway by an inter-subunit network of positively charged amino acids. Consequently, the central cavity's electrostatic potential changed from being entirely positive at the beginning of the simulation to more electronegative at the end.