The ongoing outbreak of COVID-19 has been a serious threat to human health worldwide. The virus SARS-CoV-2 initiates its infection to the human body via the interaction of its spike (S) protein with the human Angiotensin-Converting Enzyme 2 (ACE2) of the host cells. Therefore, understanding the fundamental mechanisms of how SARS-CoV-2 S protein receptor binding domain (RBD) binds to ACE2 is highly demanded for developing treatments for COVID-19. Here we implemented multi-scale computational approaches to study the binding mechanisms of human ACE2 and S proteins of both SARS-CoV and SARS-CoV-2. Electrostatic features, including electrostatic potential, electric field lines, and electrostatic forces of SARS-CoV and SARS-CoV-2 were calculated and compared in detail. The results demonstrate that SARS-CoV and SARS-CoV-2 S proteins are both attractive to ACE2 by electrostatic forces even at different distances. However, the residues contributing to the electrostatic features are quite different due to the mutations between SARS-CoV S protein and SARS-CoV-2 S protein. Such differences are analyzed comprehensively. Compared to SARS-CoV, the SARS-CoV-2 binds with ACE2 using a more robust strategy: The electric field line related residues are distributed quite differently, which results in a more robust binding strategy of SARS-CoV-2. Also, SARS-CoV-2 has a higher electric field line density than that of SARS-CoV, which indicates stronger interaction between SARS-CoV-2 and ACE2, compared to that of SARS-CoV. Key residues involved in salt bridges and hydrogen bonds are identified in this study, which may help the future drug design against COVID-19.

The Mycobacterium tuberculosis virulence factor EsxA and its chaperone EsxB are secreted as a heterodimer (EsxA:B) and are crucial for mycobacterial escape from phagosomes and cytosolic translocation. Current findings support the idea that for EsxA to interact with host membranes, EsxA must dissociate from EsxB at low pH. However, the molecular mechanism by which the EsxA:B heterodimer separates is not clear. In the present study, using liposome-leakage and cytotoxicity assays, LC-MS/MS–based proteomics, and CCF-4 FRET analysis, we obtained evidence that the Nα-acetylation of the Thr-2 residue on EsxA, a post-translational modification that is present in mycobacteria but absent in Escherichia coli, is required for the EsxA:B separation. Substitutions at Thr-2 that precluded Nα-acetylation inhibited the heterodimer separation and hence prevented EsxA from interacting with the host membrane, resulting in attenuated mycobacterial cytosolic translocation and virulence. Molecular dynamics simulations revealed that at low pH, the Nα-acetylated Thr-2 makes direct and frequent “bind-and-release” contacts with EsxB, which generates a force that pulls EsxB away from EsxA. In summary, our findings provide evidence that the Nα-acetylation at Thr-2 of EsxA facilitates dissociation of the EsxA:B heterodimer required for EsxA membrane permeabilization and mycobacterial cytosolic translocation and virulence.

In the last three decades, many giant DNA viruses have been discovered. Giant viruses present a unique and essential research frontier for studies of self-assembly and regulation of supramolecular assemblies. The question on how these giant DNA viruses assemble thousands of proteins so accurately to form their protein shells, the capsids, remains largely unanswered. Revealing the mechanisms of giant virus assembly will help to discover the mysteries of many self-assembly biology problems. Paramecium bursaria Chlorella virus-1 (PBCV-1) is one of the most intensively studied giant viruses. Here, we implemented a multi-scale approach to investigate the interactions among PBCV-1 capsid building units called capsomers. Three binding modes with different strengths are found between capsomers around the relatively flat area of the virion surface at the icosahedral 2-fold axis. Furthermore, a capsomer structure manipulation package is developed to simulate the capsid assembly process. Using these tools, binding forces among capsomers were investigated and binding funnels were observed that were consistent with the final assembled capsid. In addition, total binding free energies of each binding mode were calculated. The results helped to explain previous experimental observations. Results and tools generated in this work established an initial computational approach to answer current unresolved questions regarding giant virus assembly mechanisms. Results will pave the way for studying more complicated process in other biomolecular structures.

SARS-CoV-2 that caused COVID-19 has spread since the end of 2019. Its major effects resulted in over four million deaths around the whole world by August 2021. Therefore, understanding virulence mechanisms is important to prevent future outbreaks and for COVID-19 drug development. The envelope (E) protein is an important structural protein, affecting virus assembly and budding. The E protein pentamer is a viroporin, serving as an ion transferring channel in cells. In this work, we applied molecular dynamic simulations and topological and electrostatic analyses to study the effects of palmitoylation on the E protein pentamer. The results indicate that the cation transferring direction is more from the lumen to the cytosol. The structure of the palmitoylated E protein pentamer is more stable while the loss of palmitoylation caused the pore radius to reduce and even collapse. The electrostatic forces on the two sides of the palmitoylated E protein pentamer are more beneficial to attract cations in the lumen and to release cations into the cytosol. The results indicate the importance of palmitoylation, which can help the drug design for the treatment of COVID-19.

A large population in the world has been infected by COVID-19. Understanding the mechanisms of Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2) is important for the management and treatment of COVID-19. When it comes to the infection process, one of the most important proteins in SARS-CoV-2 is the spike (S) protein, which is able to bind to human Angiotensin-Converting Enzyme 2 (ACE2) and initializes the entry of the host cell. In this study, we implemented multiscale computational approaches to study the electrostatic features of the interfaces of the SARS-CoV-2 S protein receptor binding domain and ACE2. The simulations and analyses were performed on highperformance computing resources in the Texas Advanced Computing Center. Our study identified key residues on SARS-CoV-2, which can be used as targets for future drug design. The results shed light on future drug design and therapeutic targets for COVID-19. & THE NUMBER OF confirmed cases of Coronavirus Disease 2019 (COVID-19) is increasing dramatically 1 due to the fast spread of SARS-CoV-2. The large coronavirus family includes hundreds of viruses that usually do not pose a threat to human health. SARS-CoV-2 is the seventh member of those coronaviruses that infect the human body. Of these, four (HCoV-229E, HCoV-OC43, HCoV-NL63, HKU1) 2 cause mild to moderate symptoms, while the other three can cause serious, even fatal diseases. SARS coronavirus (SARS-CoV) broke out in 2002 and caused Severe Acute Respiratory Syndrome (SARS). MERS coronavirus (MERS-CoV)