Protein decoys may battle COVID-19 and more.

Drugs designed to resemble viruses' cellular targets move into clinical trials.

Distinct Conformations of SARS-CoV-2 Omicron Spike Protein and Its Interaction with ACE2 and Antibody.

Since November 2021, Omicron has been the dominant severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant that causes the coronavirus disease 2019 (COVID-19) and has continuously impacted human health. Omicron sublineages are still increasing and cause increased transmission and infection rates. The additional 15 mutations on the receptor binding domain (RBD) of Omicron spike proteins change the protein conformation, enabling the Omicron variant to evade neutralizing antibodies. For this reason, many efforts have been made to design new antigenic variants to induce effective antibodies in SARS-CoV-2 vaccine development. However, understanding the different states of Omicron spike proteins with and without external molecules has not yet been addressed. In this review, we analyze the structures of the spike protein in the presence and absence of angiotensin-converting enzyme 2 (ACE2) and antibodies. Compared to previously determined structures for the wildtype spike protein and other variants such as alpha, beta, delta, and gamma, the Omicron spike protein adopts a partially open form. The open-form spike protein with one RBD up is dominant, followed by the open-form spike protein with two RBD up, and the closed-form spike protein with the RBD down. It is suggested that the competition between antibodies and ACE2 induces interactions between adjacent RBDs of the spike protein, which lead to a partially open form of the Omicron spike protein. The comprehensive structural information of Omicron spike proteins could be helpful for the efficient design of vaccines against the Omicron variant.

Host range and structural analysis of bat-origin RshSTT182/200 coronavirus binding to human ACE2 and its animal orthologs.

Bat-origin RshSTT182 and RshSTT200 coronaviruses (CoV) from Rhinolophus shameli in Southeast Asia (Cambodia) share 92.6% whole-genome identity with SARS-CoV-2 and show identical receptor-binding domains (RBDs). In this study, we determined the structure of the RshSTT182/200 receptor binding domain (RBD) in complex with human angiotensin-converting enzyme 2 (hACE2) and identified the key residues that influence receptor binding. The binding of the RshSTT182/200 RBD to ACE2 orthologs from 39 animal species, including 18 bat species, was used to evaluate its host range. The RshSTT182/200 RBD broadly recognized 21 of 39 ACE2 orthologs, although its binding affinities for the orthologs were weaker than those of the RBD of SARS-CoV-2. Furthermore, RshSTT182 pseudovirus could utilize human, fox, and Rhinolophus affinis ACE2 receptors for cell entry. Moreover, we found that SARS-CoV-2 induces cross-neutralizing antibodies against RshSTT182 pseudovirus. Taken together, these findings indicate that RshSTT182/200 can potentially infect susceptible animals, but requires further evolution to obtain strong interspecies transmission abilities like SARS-CoV-2.

Effect of acetic acid inactivation of SARS-CoV-2.

Effective measures are needed to prevent the spread and infectivity of SARS-CoV-2 that causes COVID-19. Chemical inactivation may help to prevent the spread and transmission of this and other viruses. Hence, we tested the SARS-CoV-2 antiviral activity of acetic acid, the main component of vinegar, in vitro. Inactivation and binding assays suggest that acetic acid is virucidal. We found that 6% acetic acid, a concentration typically found in white distilled vinegar, effectively inactivated SARS-CoV-2 after 15-min incubation with a complete loss of replication of competent virus as measured by TCID50. Transmission electron microscopy further demonstrated that 6% acetic acid disrupts SARS-CoV-2 virion structure. In addition, 6% acetic acid significantly inhibits and disrupts the binding of SARS-CoV-2 spike protein binding to ACE2, the primary SARS-CoV-2 cell receptor, after contact with spike protein for 5, 10, 30 and 60 minutes incubation. Taken together, our findings demonstrate that acetic acid possesses inactivating activity against SARS-CoV-2 and may represent a safe alternative to commonly used chemical disinfectants to effectively control the spread of SARS-CoV-2.

Protein decoys may battle COVID-19 and more.

Drugs designed to resemble viruses' cellular targets move into clinical trials.

Investigation on the interaction mechanism of different SARS-CoV-2 spike variants with hACE2: insights from molecular dynamics simulations.

Since the COVID-19 pandemic caused by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), SARS-CoV-2 has evolved by acquiring genomic mutations, resulting in the recent emergence of several SARS-CoV-2 variants with improved transmissibility and infectivity relative to the original strain. An underlying mechanism may be the increased ability of the mutants to bind the receptor proteins and infect the host cell. In this work, we implemented all-atom molecular dynamics (MD) simulations to study the binding and interaction of the receptor binding domain (RBD) of the SARS-CoV-2 spike protein singly (D614G), doubly (D614G + L452R and D614G + N501Y), triply (D614G + N501Y + E484K), and quadruply (D614G + N501Y + E484K + K417T) mutated variants with the human angiotensin-converting enzyme 2 (hACE2) receptor protein in the host cell. A combination of multiple analysis approaches elucidated the effects of mutations and the extent of molecular divergence from multiple perspectives, including the dynamic correlated motions, interaction patterns, dominant motions, free energy landscape, and charge distribution on the electrostatic potential surface between the hACE2 and all RBD variants. Moreover, free energy calculations using the MM/PBSA method evaluated the binding affinity between these RBD variants and hACE2. The results showed that the D614G + N501Y + E484K variant possessed the lowest free energy value (highest affinity) compared to the D614G + N501Y + E484K + K417T, D614G + L452R, D614G + N501Y, and D614G mutants. The residue-based energy decomposition also indicated that the energy contribution of residues at the mutation site to the total binding energy was highly variable. The interaction mechanisms between the different RBD variants and hACE2 elucidated in this study will provide some insights into the development of drugs targeting the new SARS-CoV-2 variants.

Evaluating the transmission feasibility of SARS-CoV-2 Omicron (B.1.1.529) variant to 143 mammalian hosts: insights from S protein RBD and host ACE2 interaction studies.

In comparison to previously known severe respiratory syndrome coronavirus 2 (SARS-CoV-2) variants, the newly emerged Omicron (B.1.1.529) variant shows higher infectivity in humans. Exceptionally high infectivity of this variant raises concern of its possible transmission via other intermediate hosts. The SARS-CoV-2 infectivity is established via the association of spike (S) protein receptor binding domain (RBD) with host angiotensin I converting enzyme 2 (hACE2) receptor. In the course of this study, we investigated the interaction between Omicron S protein RBD with the ACE2 receptor of 143 mammalian hosts including human by protein-protein interaction analysis. The goal of this study was to forecast the likelihood that the virus may infect other mammalian species that coexist with or are close to humans in the household, rural, agricultural, or zoological environments. The Omicron RBD was found to interact with higher binding affinity with the ACE2 receptor of 122 mammalian hosts via different amino acid residues from the human ACE2 (hACE2). The rat (Rattus rattus) ACE2 was found to show the strongest interaction with Omicron RBD with a binding affinity of -1393.6 kcal/mol. These distinct strong binding affinity of RBD of Omicron with host ACE2 indicates a greater potential of new host transmissibility and infection via intermediate hosts. Though expected but the phylogenetic position of the mammalian species may not dictate the Omicron RBD binding to the host ACE2 receptor suggesting an involvement of multiple factors in guiding host divergence of the variant.

Host range and structural analysis of bat-origin RshSTT182/200 coronavirus binding to human ACE2 and its animal orthologs.

Bat-origin RshSTT182 and RshSTT200 coronaviruses (CoV) from Rhinolophus shameli in Southeast Asia (Cambodia) share 92.6% whole-genome identity with SARS-CoV-2 and show identical receptor-binding domains (RBDs). In this study, we determined the structure of the RshSTT182/200 receptor binding domain (RBD) in complex with human angiotensin-converting enzyme 2 (hACE2) and identified the key residues that influence receptor binding. The binding of the RshSTT182/200 RBD to ACE2 orthologs from 39 animal species, including 18 bat species, was used to evaluate its host range. The RshSTT182/200 RBD broadly recognized 21 of 39 ACE2 orthologs, although its binding affinities for the orthologs were weaker than those of the RBD of SARS-CoV-2. Furthermore, RshSTT182 pseudovirus could utilize human, fox, and Rhinolophus affinis ACE2 receptors for cell entry. Moreover, we found that SARS-CoV-2 induces cross-neutralizing antibodies against RshSTT182 pseudovirus. Taken together, these findings indicate that RshSTT182/200 can potentially infect susceptible animals, but requires further evolution to obtain strong interspecies transmission abilities like SARS-CoV-2.

Potential use of the S-protein-Angiotensin converting enzyme 2 binding pathway in the treatment of coronavirus disease 2019.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the pathogen that causes coronavirus disease 2019 (COVID-19), infects humans through a strong interaction between the viral spike protein (S-protein) and angiotensin converting enzyme 2 (ACE2) receptors on the cell surface. The infection of host lung cells by SARS-CoV-2 leads to clinical symptoms in patients. However, ACE2 expression is not restricted to the lungs; altered receptors have been found in the nasal and oral mucosa, vessel, brain, pancreas, gastrointestinal tract, kidney, and heart. The future of COVID-19 is uncertain, however, new viral variants are likely to emerge. The SARS-CoV-2 Omicron variant has a total of 50 gene mutations compared with the original virus; 15 of which occur in the receptor binding domain (RBD). The RBD of the viral S-protein binds to the human ACE2 receptor for viral entry. Mutations of the ACE2-RBD interface enhance tight binding by increasing hydrogen bond interactions and expanding the accessible surface area. Extracorporeal membrane oxygenation, hyperbaric oxygen, and aggressive dialysis for the treatment of COVID-19 have shown various degrees of clinical success. The use of decoy receptors based on the ACE2 receptor as a broadly potent neutralizer of SARS-CoV-2 variants has potential as a therapeutic mechanism. Drugs such as 3E8 could block binding of the S1-subunit to ACE2 and restrict the infection of ACE2-expressing cells by a variety of coronaviruses. Here, we discuss the development of ACE2-targeted strategies for the treatment and prevention of COVID-19.

A broadly neutralizing monoclonal antibody overcomes the mutational landscape of emerging SARS-CoV-2 variants of concern.

The emergence of new variants of SARS-CoV-2 necessitates unremitting efforts to discover novel therapeutic monoclonal antibodies (mAbs). Here, we report an extremely potent mAb named P4A2 that can neutralize all the circulating variants of concern (VOCs) with high efficiency, including the highly transmissible Omicron. The crystal structure of the P4A2 Fab:RBD complex revealed that the residues of the RBD that interact with P4A2 are a part of the ACE2-receptor-binding motif and are not mutated in any of the VOCs. The pan coronavirus pseudotyped neutralization assay confirmed that the P4A2 mAb is specific for SARS-CoV-2 and its VOCs. Passive administration of P4A2 to K18-hACE2 transgenic mice conferred protection, both prophylactically and therapeutically, against challenge with VOCs. Overall, our data shows that, the P4A2 mAb has immense therapeutic potential to neutralize the current circulating VOCs. Due to the overlap between the P4A2 epitope and ACE2 binding site on spike-RBD, P4A2 may also be highly effective against a number of future variants.

Lucidenic acid A inhibits the binding of hACE2 receptor with spike protein to prevent SARS-CoV-2 invasion.

High infection caused by mutations of SARS-CoV-2 calls for new prevention strategy. Ganoderma lucidum known as a superior immunoenhancer exhibits various antiviral effects, whether it can resist SARS-CoV-2 remains unclear. Herein, virtual screening combined with in vitro hACE2 inhibition assays were used to investigate its anti SARS-CoV-2 effect. Potential 54 active components, 80 core targets and 20 crucial pathways were identified by the component-target-pathway network. The binding characters of these components to hACE2 and its complexes with spike protein including omicron variant was analyzed by molecular docking. Lucidenic acid A was selected as the top molecule with high affinity to all receptors by forming hydrogen bonds. Molecular dynamics simulation showed it had good binding stability with the receptor proteins. Finally, in vitro FRET test demonstrated it inhibited the hACE2 activity with IC50 2 µmol/mL. Therefore, lucidenic acid A can prevent the virus invasion by blocking hACE2 binding with SARS-CoV-2.

Omicron BA.1 and BA.2 variants increase the interactions of SARS-CoV-2 spike glycoprotein with ACE2.

SARS-CoV-2 infection is initiated by binding of the receptor-binding domain (RBD) of its spike glycoprotein to the peptidase domain (PD) of angiotensin-converting enzyme 2 (ACE2) receptors in host cells. Recently detected Omicron variant of SARS-CoV-2 (B.1.1.529) is heavily mutated on RBD. First the BA.1 and later the BA.2 variant became the most dominant strains of the Omicron variant. To investigate how the mutations of these strains affect RBD-PD interactions, we performed all-atom molecular dynamics simulations of the BA.1 and BA.2 RBD-PD in the presence of full-length glycans, explicit water, and ions. Simulations revealed that RBDs of BA.1 and BA.2 variants exhibit a more dispersed interaction network and make an increased number of salt bridges and hydrophobic interactions with PD compared to wild-type RBD. Although BA.1 and BA.2 differ in two residues at the RBD-ACE2 interface, no major difference in RBD-PD interactions and binding strengths were observed between these variants. Using the conformations sampled in each trajectory, the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) method estimated ∼34% and ∼51% stronger binding free energies to PD for BA.1 and BA.2 RBD, respectively, than wild-type RBD, which may result in higher binding efficiency of the Omicron variant to infect host cells.

Biosynthetic proteins targeting the SARS-CoV-2 spike as anti-virals.

The binding of the SARS-CoV-2 spike to angiotensin-converting enzyme 2 (ACE2) promotes virus entry into the cell. Targeting this interaction represents a promising strategy to generate antivirals. By screening a phage-display library of biosynthetic protein sequences build on a rigid alpha-helicoidal HEAT-like scaffold (named αReps), we selected candidates recognizing the spike receptor binding domain (RBD). Two of them (F9 and C2) bind the RBD with affinities in the nM range, displaying neutralisation activity in vitro and recognizing distinct sites, F9 overlapping the ACE2 binding motif. The F9-C2 fusion protein and a trivalent αRep form (C2-foldon) display 0.1 nM affinities and EC50 of 8-18 nM for neutralization of SARS-CoV-2. In hamsters, F9-C2 instillation in the nasal cavity before or during infections effectively reduced the replication of a SARS-CoV-2 strain harbouring the D614G mutation in the nasal epithelium. Furthermore, F9-C2 and/or C2-foldon effectively neutralized SARS-CoV-2 variants (including delta and omicron variants) with EC50 values ranging from 13 to 32 nM. With their high stability and their high potency against SARS-CoV-2 variants, αReps provide a promising tool for SARS-CoV-2 therapeutics to target the nasal cavity and mitigate virus dissemination in the proximal environment.

Missense mutations in spike protein of SARS-CoV-2 delta variant contribute to the alteration in viral structure and interaction with hACE2 receptor.

INTRODUCTION: Many of the global pandemics threaten human existence over the decades among which coronavirus disease (COVID-19) is the newest exposure circulating worldwide. The RNA encoded severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus is referred as the pivotal agent of this deadly disease that induces respiratory tract infection by interacting host ACE2 receptor with its spike glycoprotein. Rapidly evolving nature of this virus modified into new variants helps in perpetrating immune escape and protection against host defense mechanism. Consequently, a new isolate, delta variant originated from India is spreading perilously at a higher infection rate. METHODS: In this study, we focused to understand the conformational and functional significance of the missense mutations found in the spike glycoprotein of SARS-CoV-2 delta variant performing different computational analysis. RESULTS: From physiochemical analysis, we found that the acidic isoelectric point of the virus elevated to basic pH level due to the mutations. The targeted mutations were also found to change the interactive bonding pattern and conformational stability analyzed by the molecular dynamic's simulation. The molecular docking study also revealed that L452R and T478K mutations found in the RBD domain of delta variant spike protein contributed to alter interaction with the host ACE2 receptor. CONCLUSIONS: Overall, this study provided insightful evidence to understand the morphological and attributive impact of the mutations on SARS-CoV-2 delta variant.

Characterization of Entry Pathways, Species-Specific Angiotensin-Converting Enzyme 2 Residues Determining Entry, and Antibody Neutralization Evasion of Omicron BA.1, BA.1.1, BA.2, and BA.3 Variants.

The SARS-CoV-2 Omicron variants were first detected in November 2021, and several Omicron lineages (BA.1, BA.2, BA.3, BA.4, and BA.5) have since rapidly emerged. Studies characterizing the mechanisms of Omicron variant infection and sensitivity to neutralizing antibodies induced upon vaccination are ongoing by several groups. In the present study, we used pseudoviruses to show that the transmembrane serine protease 2 (TMPRSS2) enhances infection of BA.1, BA.1.1, BA.2, and BA.3 Omicron variants to a lesser extent than ancestral D614G. We further show that Omicron variants have higher sensitivity to inhibition by soluble angiotensin-converting enzyme 2 (ACE2) and the endosomal inhibitor chloroquine compared to D614G. The Omicron variants also more efficiently used ACE2 receptors from 9 out of 10 animal species tested, and unlike the D614G variant, used mouse ACE2 due to the Q493R and Q498R spike substitutions. Finally, neutralization of the Omicron variants by antibodies induced by three doses of Pfizer/BNT162b2 mRNA vaccine was 7- to 8-fold less potent than the D614G. These results provide insights into the transmissibility and immune evasion capacity of the emerging Omicron variants to curb their ongoing spread. IMPORTANCE The ongoing emergence of SARS-CoV-2 Omicron variants with an extensive number of spike mutations poses a significant public health and zoonotic concern due to enhanced transmission fitness and escape from neutralizing antibodies. We studied three Omicron lineage variants (BA.1, BA.2, and BA.3) and found that transmembrane serine protease 2 has less influence on Omicron entry into cells than on D614G, and Omicron exhibits greater sensitivity to endosomal entry inhibition compared to D614G. In addition, Omicron displays more efficient usage of diverse animal species ACE2 receptors than D614G. Furthermore, due to Q493R/Q498R substitutions in spike, Omicron, but not D614G, can use the mouse ACE2 receptor. Finally, three doses of Pfizer/BNT162b2 mRNA vaccination elicit high neutralization titers against Omicron variants, although the neutralization titers are still 7- to 8-fold lower those that against D614G. These results may give insights into the transmissibility and immune evasion capacity of the emerging Omicron variants to curb their ongoing spread.

Intermediates in SARS-CoV-2 spike-mediated cell entry.

SARS-CoV-2 cell entry is completed after viral spike (S) protein-mediated membrane fusion between viral and host cell membranes. Stable prefusion and postfusion S structures have been resolved by cryo-electron microscopy and cryo-electron tomography, but the refolding intermediates on the fusion pathway are transient and have not been examined. We used an antiviral lipopeptide entry inhibitor to arrest S protein refolding and thereby capture intermediates as S proteins interact with hACE2 and fusion-activating proteases on cell-derived target membranes. Cryo-electron tomography imaged both extended and partially folded intermediate states of S2, as well as a novel late-stage conformation on the pathway to membrane fusion. The intermediates now identified in this dynamic S protein-directed fusion provide mechanistic insights that may guide the design of CoV entry inhibitors.

Computational studies evidenced the potential of steroidal lactone to disrupt surface interaction of SARS-CoV-2 spike protein and hACE2.

The critical event in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pathogenesis is recognition of host cells by the virus, which is facilitated by protein-protein interaction (PPI) of viral Spike-Receptor Binding Domain (S-RBD) and Human Angiotensin Converting Enzyme 2-Receptor (hACE2-R). Thus, disrupting the interaction between S-RBD and hACE2-R is widely accepted as a primary strategy for managing COVID-19. The purpose of this study is to assess the ability of three steroidal lactones (SL) (4-Dehydrowithaferin A, Withaferin A, and Withalongolide A) derived from plants to disrupt the PPI of S-RBD and hACE2-R under two conditions (CON-I and CON-II) using in-silico methods. Under CON-I, 4-Dehydrowithaferin A destabilizing the interactions between S-RBD and hACE2-R, as indicated by an increase in binding energy (BE) from -1028.5 kJ/mol (control) to -896.12 kJ/mol 4-Dehydrowithaferin A exhibited a strong interaction with S-RBD GLY496 with a hydrogen bond occupancy (HBO) of 37.33%. Under CON-II, Withalongolide A was capable of disrupting all types of PPI, as evidenced by an increased BE from -913 kJ/mol (control) to -133.69 kJ/mol and an increased distance (>3.55 nm) between selected AAR combinations of S-RBD and hACE2-R. Withalongolide A formed a hydrogen bond with TYR453 (97%, HBO) of S-RBD, which is required for interaction with hACE2-R's HIS34. Our studies demonstrated that SL molecules have the potential to disrupt the S-RBD and hACE2-R interaction, thereby preventing SARS-CoV-2 from recognizing host cells. The SL molecules can be considered for additional in-vitro and in-vivo studies with this research evidence.

ACE2-binding exposes the SARS-CoV-2 fusion peptide to broadly neutralizing coronavirus antibodies.

The coronavirus spike glycoprotein attaches to host receptors and mediates viral fusion. Using a broad screening approach, we isolated seven monoclonal antibodies (mAbs) that bind to all human-infecting coronavirus spike proteins from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) immune donors. These mAbs recognize the fusion peptide and acquire affinity and breadth through somatic mutations. Despite targeting a conserved motif, only some mAbs show broad neutralizing activity in vitro against alpha- and betacoronaviruses, including animal coronaviruses WIV-1 and PDF-2180. Two selected mAbs also neutralize Omicron BA.1 and BA.2 authentic viruses and reduce viral burden and pathology in vivo. Structural and functional analyses showed that the fusion peptide-specific mAbs bound with different modalities to a cryptic epitope hidden in prefusion stabilized spike, which became exposed upon binding of angiotensin-converting enzyme 2 (ACE2) or ACE2-mimicking mAbs.

Interfacial Water in the SARS Spike Protein: Investigating the Interaction with Human ACE2 Receptor and In Vitro Uptake in A549 Cells.

The severity of global pandemic due to severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has engaged the researchers and clinicians to find the key features triggering the viral infection to lung cells. By utilizing such crucial information, researchers and scientists try to combat the spread of the virus. Here, in this work, we performed in silico analysis of the protein-protein interactions between the receptor-binding domain (RBD) of the viral spike protein and the human angiotensin-converting enzyme 2 (hACE2) receptor to highlight the key alteration that happened from SARS-CoV to SARS-CoV-2. We analyzed and compared the molecular differences between spike proteins of the two viruses using various computational approaches such as binding affinity calculations, computational alanine, and molecular dynamics simulations. The binding affinity calculations showed that SARS-CoV-2 binds a little more firmly to the hACE2 receptor than SARS-CoV. The major finding obtained from molecular dynamics simulations was that the RBD-ACE2 interface is populated with water molecules and interacts strongly with both RBD and ACE2 interfacial residues during the simulation periods. The water-mediated hydrogen bond by the bridge water molecules is crucial for stabilizing the RBD and ACE2 domains. Near-ambient pressure X-ray photoelectron spectroscopy (NAP-XPS) confirmed the presence of vapor and molecular water phases in the protein-protein interfacial domain, further validating the computationally predicted interfacial water molecules. In addition, we examined the role of interfacial water molecules in virus uptake by lung cell A549 by binding and maintaining the RBD/hACE2 complex at varying temperatures using nanourchins coated with spike proteins as pseudoviruses and fluorescence-activated cell sorting (FACS) as a quantitative approach. The structural and dynamical features presented here may serve as a guide for developing new drug molecules, vaccines, or antibodies to combat the COVID-19 pandemic.

Computational Analysis of the SARS-CoV-2 RBD-ACE2-Binding Process Based on MD and the 3D-RISM Theory.

The binding process of angiotensin-converting enzyme 2 (ACE2) to the receptor-binding domain (RBD) of the severe acute respiratory syndrome-like coronavirus 2 spike protein was investigated using molecular dynamics simulation and the three-dimensional reference interaction-site model theory. The results suggested that the protein-binding process consists of a protein-protein approaching step, followed by a local structural rearrangement step. In the approaching step, the interprotein interaction energy decreased as the proteins approached each other, whereas the solvation free energy increased. As the proteins approached, the glycan of ACE2 first established a hydrogen bond with the RBD. Thereafter, the number of interprotein hydrogen bonds increased rapidly. The solvation free energy increased because of the desolvation of the protein as it approached its partner. The spatial distribution function of the solvent revealed the presence of hydrogen bonds bridged by water molecules on the RBD-ACE2 interface. Finally, principal component analysis revealed that ACE2 showed a pronounced conformational change, whereas there was no significant change in RBD.