Modeling the Opening SARS-CoV-2 Spike: an Investigation of its Dynamic Electro-Geometric Properties (preprint)

The recent COVID-19 pandemic has brought about a surge of crowd-sourced initiatives aimed at simulating the proteins of the SARS-CoV-2 virus. A bottleneck currently exists in translating these simulations into tangible predictions that can be leveraged for pharmacological studies. Here we report on extensive electrostatic calculations done on an exascale simulation of the opening of the SARS-CoV-2 spike protein, performed by the Folding@home initiative. We compute the electric potential as the solution of the non-linear Poisson-Boltzmann equation using a parallel sharp numerical solver. The inherent multiple length scales present in the geometry and solution are reproduced using highly adaptive Octree grids. We analyze our results focusing on the electro-geometric properties of the receptor-binding domain and its vicinity. This work paves the way for a new class of hybrid computational and data-enabled approaches, where molecular dynamics simulations are combined with continuum modeling to produce high-fidelity computational measurements serving as a basis for protein bio-mechanism investigations.

Large-scale single-cell analysis reveals critical immune characteristics of COVID-19 patients (preprint)

Dysfunctional immune response in the COVID-19 patients is a recurrent theme impacting symptoms and mortality, yet the detailed understanding of pertinent immune cells is not complete. We applied single-cell RNA sequencing to 284 samples from 205 COVID-19 patients and controls to create a comprehensive immune landscape. Lymphopenia and active T and B cell responses were found to coexist and associated with age, sex and their interactions with COVID-19. Diverse epithelial and immune cell types were observed to be virus-positive and showed dramatic transcriptomic changes. Elevation of ANXA1 and S100A9 in virus-positive squamous epithelial cells may enable the initiation of neutrophil and macrophage responses via the ANXA1-FPR1 and S100A8/9-TLR4 axes. Systemic up-regulation of S100A8/A9, mainly by megakaryocytes and monocytes in the peripheral blood, may contribute to the cytokine storms frequently observed in severe patients. Our data provide a rich resource for understanding the pathogenesis and designing effective therapeutic strategies for COVID-19.

SARS-CoV-2 Quasispecies provides insight into its genetic dynamics during infection (preprint)

A novel coronavirus disease (COVID-19) caused by SARS-CoV-2 has been pandemic worldwide. The genetic dynamics of quasispecies afford RNA viruses a great fitness on cell tropism and host range. However, no quasispecies data of SARS-CoV-2 have been reported yet. To explore quasispecies haplotypes and its transmission characteristics, we carried out single-molecule real-time (SMRT) sequencing of the full-length of SARS-CoV-2 spike gene within 14 RNA samples from 2 infection clusters, covering first-to third-generation infected-patients. We observed a special quasispecies structure of SARS-CoV-2 (modeled as One-King): one dominant haplotype (mean abundance ~70.15%) followed by numerous minor haplotypes (mean abundance < 0.10%). We not only discovered a novel dominant haplotype of F1040 but also realized that minor quasispecies were also worthy of attention. Notably, some minor haplotypes (like F1040 and currently pandemic one G614) could potentially reveal adaptive and converse into the dominant one. However, minor haplotypes exhibited a high transmission bottleneck (~6% could be stably transmitted), and the new adaptive/dominant haplotypes were likely originated from genetic variations within a host rather than transmission. The evolutionary rate was estimated as 2.68-3.86 x 10-3 per site per year, which was larger than the estimation at consensus genome level. The One-King model and conversion event expanded our understanding of the genetic dynamics of SARS-CoV-2, and explained the incomprehensible phenomenon at the consensus genome level, such as limited cumulative mutations and low evolutionary rate. Moreover, our findings suggested the epidemic strains may be multi-host origin and future traceability would face huge difficulties.

Spike protein up-regulates inflammatory axis of both thromboinflammation and leukotriene in severe COVID-19 (preprint)

Hypercytokinemia is a critically fatal factor in COVID-19. However, underlying pathogenic mechanisms are unknown. Here we show that fibrinogen and leukotriene-A4 hydrolase (LTA4H), two of the most potent inflammatory contributors, are elevated by 67.7 and astonishing 227.7% in the plasma of patients infected by SARS-CoV-2 and admitted to intensive care unit in comparison with healthy control, respectively. Conversely, transferrin identified as a fibrinogen immobilizer in our recent work and Spink6 are down-regulated by 40.3 and 25.9%, respectively. Furthermore, we identify Spink6 as the first endogenous inhibitor of LTA4H, a pro-inflammatory enzyme catalyzing final and rating limited step in biosynthesis of leukotriene-B4 that is an extremely inflammatory mediator and a target to design superior anti-inflammatory drugs. Additionally, virus Spike protein is found to evoke LTA4H and fibrinogen expression in vivo. Collectively, these findings identify the imbalance between inflammatory drivers and antagonists, which likely contributes to hypercytokinemia in COVID-19. Spink6 may have superior anti-inflammatory function because it specifically targets epoxide hydrolase of LTA4H to inhibit leukotriene-B4 biosynthesis without effecting LTA4H’s aminopeptidase activity. 

Viral Kinetics and Antibody Responses in Patients with COVID-19 (preprint)

Background A pandemic of coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been spreading over the world. However, the viral dynamics, host serologic responses, and their associations with clinical manifestations, have not been well described in prospective cohort. Methods We conducted a prospective cohort and enrolled 67 COVID-19 patients admitting between Jan 26 and Feb 5, 2020. Clinical specimens including nasopharyngeal swab, sputum, blood, urine and stool were tested periodically according to standardized case report form with final follow-up on February 27. The routes and duration of viral shedding, antibody response, and their associations with disease severity and clinical manifestations were systematically evaluated. Coronaviral particles in clinical specimens were observed by transmission electron microscopy (TEM). Results The median duration of SARS-CoV-2 RNA shedding were 12 (3-38), 19 (5-37), and 18 (7-26) days in nasopharyngeal swabs, sputum and stools, respectively. Only 13 urines (5.6%) and 12 plasmas (5.7%) were viral positive. Prolonged viral shedding was observed in severe patients than that of non-severe patients. Cough but not fever, aligned with viral shedding in clinical respiratory specimens, meanwhile the positive stool-RNA appeared to align with the proportion who concurrently had cough and sputum production, but not diarrhea. Typical coronaviral particles could be found directly in sputum by TEM. The anti-nucleocapsid-protein IgM started on day 7 and positive rate peaked on day 28, while that of IgG was on day 10 and day 49 after illness onset. IgM and IgG appear earlier, and their titers are significantly higher in severe patients than non-severe patients (p<0.05). The weak responders for IgG had a significantly higher viral clearance rate than that of strong responders (p= 0.011). Conclusions Nasopharyngeal, sputum and stools rather than blood and urine, were the major shedding routes for SARS-CoV-2, and meanwhile sputum had a prolonged viral shedding. Symptom cough seems to be aligned with viral shedding in clinical respiratory and fecal specimens. Stronger antibody response was associated with delayed viral clearance and disease severity.

Diagnosis of Acute Respiratory Syndrome Coronavirus 2 Infection by Detection of Nucleocapsid Protein (preprint)

BACKGROUND Nucleic acid test and antibody assay have been employed in the diagnosis for SARS-CoV-2 infection, but the use of viral antigen for diagnosis has not been successfully developed. Theoretically, viral antigen is the specific marker of the virus and precedes antibody appearance within the infected population. There is a clear need of detection of viral antigen for rapid and early diagnosis. METHODS We included a cohort of 239 participants with suspected SARS-CoV-2 infection from 7 centers for the study. We measured nucleocapsid protein in nasopharyngeal swab samples in parallel with the nucleic acid test. Nucleic acid test was taken as the reference standard, and statistical evaluation was taken in blind. We detected nucleocapsid protein in 20 urine samples in another center, employing nasopharyngeal swab nucleic acid test as reference standard. RESULTS We developed a fluorescence immunochromatographic assay for detecting nucleocapsid protein of SARS-CoV-2 in nasopharyngeal swab sample and urine within 10 minutes. 100% of nucleocapsid protein positive and negative participants accord with nucleic acid test for same samples. Further, earliest participant after 3 days of fever can be identified by the method. In an additional preliminary study, we detected nucleocapsid protein in urine in 73.6% of diagnosed COVID-19 patients. CONCLUSIONS Those findings indicate that nucleocapsid protein assay is an accurate, rapid, early and simple method for diagnosis of COVID-19. Appearance of nucleocapsid protein in urine coincides our finding of the SARS-CoV-2 invading kidney and might be of diagnostic value.