key: cord-0747783-u9o8695z
authors: Zhang, Liqun; Ghosh, Santosh K.; Basavarajappa, Shrikanth C.; Chen, Yinghua; Shrestha, Pravesh; Penfield, Jackson; Brewer, Ann; Ramakrishnan, Parameswaran; Buck, Matthias; Weinberg, Aaron
title: HBD-2 binds SARS-CoV-2 RBD and blocks viral entry: Strategy to combat COVID-19
date: 2022-02-02
journal: iScience
DOI: 10.1016/j.isci.2022.103856
sha: 0539b76be8920bcd67d30e1791428c106511262c
doc_id: 747783
cord_uid: u9o8695z

New approaches to complement vaccination are needed to combat the spread of SARS-CoV-2 and stop COVID-19 related deaths and medical complications. Human beta defensin 2 (hBD-2) is a naturally occurring epithelial cell derived host defense peptide that has antiviral properties. Our comprehensive in-silico studies demonstrate that hBD-2 binds the site on the CoV-2-RBD that docks with the ACE2 receptor. Biophysical measurements confirm that hBD-2 indeed binds to the CoV-2-receptor binding domain (RBD) (KD ∼ 2μM by surface plasmon resonance), preventing it from binding to ACE2 expressing cells. Importantly, hBD-2 shows specificity by blocking CoV-2/spike pseudoviral infection, but not VSVG mediated infection, of ACE2 expressing human cells with an IC50 of 2.8 + 0.4 μM. These promising findings offer opportunities to develop hBD-2 and/or its derivatives and mimetics to safely and effectively use as agents to prevent SARS-CoV-2 infection.

The SARS-CoV-2 (CoV-2)-spike protein receptor binding domain (RBD) is the interaction site that binds the angiotensin converting enzyme II (ACE2) receptor for entry into vulnerable cells.

Zhang et al visualize the interaction of human beta defensin 2 (hBD-2) with the RBD. Biophysical experiments suggest and then a cell based assay demonstrate the ability of hBD-2 to block CoV-2-spike protein expressing pseudovirions from entering ACE2 expressing human cells.

J o u r n a l P r e -p r o o f INTRODUCTION SARS-CoV-2 (CoV-2) continues to impact worldwide health; having claimed over 4 million lives as of August 2021 (worldometers.info/coronavirus/). The discovery of new approaches to complement vaccinations and impede the virus' relentless spread remains an urgent challenge. CoV-2 expresses the Spike (S) protein (Siu et al., 2008; Yoshimoto, 2020) , which is responsible for binding to the receptor angiotensin converting enzyme 2 (ACE2), followed by fusion of the viral and cellular membranes. To engage ACE2, the receptor-binding domain (RBD) of the S protein undergoes hinge-like conformational movements that transiently hide or expose its determinants for receptor binding (Wrapp et al., 2020) .

Since this is a critical initial event in the infection cascade, the RBD is a key target for therapeutic strategies (Tai et al., 2020) . The high degree of dynamics of the RBD:ACE2 complex (Brielle et al., 2020; Ghorbani et al., 2020; Spinello et al., 2020) suggests that binding of small flexible peptides may inhibit spike protein:host cell receptor interactions, which can be interrogated by computational modeling and simulations most suitable for exploring these interactions (Amaro and Mulholland, 2020 ).

Nature's own antimicrobial peptides (AMPs) have been proposed as multifunctional defenses that participate in the elimination of pathogenic microorganisms, including bacteria, fungi, and viruses (Mookherjee et al., 2020) . Exhibiting antimicrobial and immunomodulatory properties, AMPs have been intensively studied as alternatives and/or adjuncts to antibiotics in bacterial infections and have also gained substantial attention as anti-viral agents (Mulder et al., 2013) . Human beta defensins (hBDs), the major AMP group expressed in mucosal epithelium, provide a first-line of defense against various infectious pathogens, including enveloped viruses (Leikina et al., 2005; Quiñones-Mateu et al., 2003; Ryan et al., 2011) . HBD-2 is expressed throughout the respiratory epithelium from the oral cavity to the lungs, and plays an important role in defense against respiratory infections (Mookherjee et al., 2020) . It is an amphipathic, beta-sheeted, cationic (+6 charge) molecule of 41 amino acids, and is stabilized by three intramolecular disulfide bonds that protect it from degradation by proteases (Sawai et al., 2001) .

In the present study, we examined the ability of hBD-2 to act as a blocking agent against CoV-2. Through extensive in silico docking and molecular dynamics simulation analyses we discovered that hBD-2 binds to the receptor binding motif (RBM) of the RBD of CoV-2 that associates with the ACE2 receptor. Biophysical and biochemical studies, respectively, confirm that hBD-2 binds the RBD and also prevents it from binding ACE2. By utilizing a CoV-2 spike-pseudotyped luciferase assay platform, we revealed that hBD-2 effectively blocks CoV-2 spike expressing pseudovirions from entering ACE2 expressing human cells. Harnessing the utility of naturally occurring AMPs, such as hBD-2, and their derived smaller peptides, could be a viable approach in developing novel CoV-2 therapeutics.

Interaction of SARS-CoV-2 RBD with ACE2 and hBD-2 using in silico docking and molecular dynamics simulations.

A. RBD:ACE2 complex: We ran, as a reference, a 50 ns all-atom molecular dynamics (MD) simulation of the ACE2:RBD complex. The final structure is compared with the initial experimental structure (Lan et al., 2020) in Figure 1A . Only moderate deviations are seen in some of the loop regions and at the N-and C-termini of both proteins ( Figure 1B ). The solvent accessible surface area, buried between the RBD and ACE2 proteins in the complex, fluctuates only modestly (750-1000 Å 2 around an average of ~900 Å 2 over the last 25 ns of the trajectory), further indicating the overall stability of the complex in the simulations.

B. RBD:hBD-2 (monomer): We carried out docking with CLUSPRO and HADDOCK (see Methods) to explore the initial possible bound structures between the RBD and hBD-2. In the case of CLUSPRO, the entire proteins were assessed for possible interaction surfaces. The best predicted models involved the receptor-binding motif (RBM) of the RBD, and these structures were used as starting structures for allatom MD and carried out for up to 500 ns. We also ran repeat simulations with different starting seeds J o u r n a l P r e -p r o o f (initial velocity assignments). All simulations are summarized in Table S1 . By comparing the initial and final structures, and presenting the most stable trajectory ( Figure 1C & D), we observed visually, and by RMSD, a slight rotation of hBD-2 relative to the initial structure at 75 ns. However, for the remaining 425 ns, hBD-2 stayed in the same location but made a range of interactions. This observation is remarkable since it is becoming clear that many protein complexes with the moderate binding affinity of ACE2:RBD and hBD-2:RBD are considerably dynamic and that the initial association involves many encounter states, with a majority not even involving the final interfaces (e.g. Li & Buck, 2019) . The persistency of the complex was also confirmed by the amount of surface area buried between the two proteins, fluctuating moderately around a value of 700 ± 150 Å 2 . The value is smaller than that seen on binding to ACE2 (900 Å 2 ), indicating that less area is covered, consistent with the smaller size of hBD-2. The three repeat trajectories exhibited comparable interactions of hBD-2 with RBD ( Figure S2 ).

A distance map, comparing residues which are on average closer than 5 Å in the RBD:ACE2 and RBD:hBD-2 complexes is shown in Figures 1E and 1F . For the RBD:ACE2, interaction residues 20 to 45, 75 to 85, as well as a short stretch of residues around 327, 355 and 387 on ACE2 bind with the RBD, whose binding interface ranges from residue 445 to 505. Some of the RBD residues are in loop regions; e.g., 404 and 417, which also come close to ACE2 over the course of the simulation. In total, 5 regions of the RBD (shaded in green) are involved. Remarkably, essentially all residues of the RBD, which contact ACE2 (either the same ones or their close neighbors), are also in contact with hBD-2 ( Figure 1F ; also shaded green). There are, however, some subtle shifts; i.e., RBD residues 475-478 make contact with ACE2 (shaded region 4; Figure. 1E) but not with hBD-2, where these interactions may have shifted to residue 473. Also, a regional area of RBD residues 438-444 contacts hBD-2; an area not seen interacting with ACE2, possibly because it is less dynamic. These alternative interactions may provide a mechanistic entry for hBD-2 in replacing/competing away ACE2 from the spike trimer.

C. Dynamics in the RBD:hBD-2 compared to the RBD:ACE2 complex: The rms fluctuation (RMSF Ca) of each residue in the RBD and ACE2, both when they are in the complex and as differences between free J o u r n a l P r e -p r o o f and bound proteins are shown in Figure S2A and S2B, respectively. As expected, most regions at the RBD:ACE2 interaction interface become less flexible, while other changes, including increases in fluctuations are seen further away from the interface, consistent with the recent description of allostery in the spike protein (Gross et al., 2020) . Upon complex formation, the RBD and ACE2 proteins form intermolecular hydrogen bonds ( Figure S2C ), which is one of the driving forces for their binding.

Importantly, these bonds are highly dynamic with occupancy between 20-40% ( Figure S1D ). In total, 7 of the 9 H-bonds of the RBD:ACE2 interface in the crystal structure (Lan et al., 2020) are populated with reasonable occupancy in the simulations. Similar behavior has been seen in other simulations (Ghorbani et al., 2020; Spinello et al., 2020) with the difference likely explained by solution vs. crystallization conditions. Water molecules were observed at the interface in other simulations and are likely bridging the interactions (Ali and Vijayan, 2020) , also underscoring the dynamic nature of the interactions (see below).

The comparison of main-chain fluctuations in RBD:hBD-2, again between the bound and free states of the proteins, is shown in Figure 2A . Overall, the binding region becomes less flexible on the RBD in similar key regions which exhibited dampened dynamics by ACE2 binding, while on the side of hBD-2 a number of main-chain sites also see their fluctuations decreased ( Figure 2B , showing results mapped to the final structure). Intermolecular hydrogen bonds ( Figure 2C ) are fewer, with an average 4 ± 1, compared to those bridging the RBD:ACE2 complex ( Figure S2C ). With exception of the Glu 484:

Arg 23 hydrogen bond with a greater than 50% persistency ( Figure 2D ), the occupancy of other hydrogen bonds is reduced compared to the reference complex, suggesting that they are somewhat more dynamic (see discussion below) and are likely accompanied by indirect H-bond interactions with water molecules near or at the interface bridging the interactions. D. RBD:hBD-2 (dimer): Although the dimerization of hBD-2 by itself is modest (Hoover et al., 2000) , it is possible that binding to the RBD stabilizes the dimeric form. We, therefore, also docked the hBD-2 dimer to the RBD and carried out simulations. Results are given in Figure S3A -D, suggesting that binding of an J o u r n a l P r e -p r o o f hBD-2 dimer is slightly stronger, H-bonds are more persistent, and that both units of the dimer make contact with the RBD.

E. RBD:hBD-2-interaction energy calculation: The average interaction energies, calculated with GBSA (see Methods) have similar magnitudes and all are slightly negative ( Figure S4 ). The average binding energy of the RBD:ACE2 is -45 ± 7 kcal/mol whereas average binding energy of RBD:hBD-2 (monomers)

is -34 ± 8 kcal/mol (and is similar to that of the RBD:hBD-2 dimer). These estimates for the free energy of binding have a greater negative value (suggesting stronger binding) than the experimental binding energy which is measured to be -10 kcal/mol to -12 kcal/mol. Nevertheless, our prediction is consistent with the reports of others (e.g. He et al, 2020) , using the same computational method, who calculated a free energy for the RBD:ACE2 interaction of -50.4 kcal/mol. However, it is likely that the entropy change upon binding of the RBD to hBD-2 is more favorable, at least from a protein dynamics and solvation perspective, compared of that of RBD binding to ACE2. This is because hBD-2 is more dynamic in the bound state, giving less of an entropy penalty upon binding. A caveat is that the surface buried is smaller between hBD-2 and the RBD than between ACE2 and the RBD and that in the latter complex there may be long range dynamic allosteric changes.

We used plasmon surface resonance (SPR) to measure binding of recombinant hBD-2 (rhBD-2) with recombinant RBD (rRBD) and determined an affinity of 1.8 ± 1.1 µM ( Figure 3A ). For measurements using microscale thermophoresis (MST) see Figure S5 . We then followed up using a functional ELISA assay, and found that rhBD-2 specifically bound to immobilized rRBD ( Figure 3B ). Incubating rhBD-2 (1.75 µM) with His-tagged rRBD (His-rRBD; 0.19 µM), followed by nickel bead immunoprecipitation and probing for hBD-2 in Western blots, resulted in significant binding of hBD-2 to His-RBD ( Figure 3C ).

We utilized HEK 293T cells that overexpress the human ACE2 receptor in the assays and incubated these cells with HIS-RBD containing culture supernatant with and without rhBD-2. We immunoprecipitated RBD through the HIS-tag and examined the co-precipitation of ACE2. HIS-RBD effectively precipitated ACE2 and the addition of rhBD-2 competitively decreased RBD-ACE2 binding ( Figure 3D ).

A luciferase reporter expressing CoV-2 spike-dependent lentiviral system was used to study the competitive inhibitory effects of rhBD-2 on CoV-2 spike-mediated infection. To demonstrate that the inhibitory effect of hBD-2 is specific to spike-mediated infection, we used a vesicular stomatitis virus glycoprotein (VSVG) pseudotyped virus as a control. Viruses pseudotyped with VSVG are pantropic; i.e., they can infect all cell types (Lever et al., 2004) , and do not depend on ACE2 for entry.

We infected ACE2 expressing HEK 293T cells using different dilutions of the luciferase expressing spike or VSVG pseudotyped virus as a control, to determine similar relative RLUs that correspond to a comparable multiplicity of infection (MOI). We found that the infection of undiluted spike pseudotyped virus was similar to the infection obtained using a 1:1000 dilution of the VSVG pseudotyped virus ( Figure   4A and 4B), which was used in the subsequent experiments. To study the hBD-2 mediated inhibition of viral entry, we co-infected cells with spike or VSVG pseudotyped virus with varying amounts of rhBD-2.

To better represent real life infection, rhBD-2, spike pseudovirions and cells were incubated simultaneously, thereby preventing any advantage preincubation would have offered hBD-2. Additionally, preincubation of rhBD-2 with pseudovirions could have confounded results due to possible perturbation of viral membranes by the AMP. We found that rhBD-2 dose-dependently decreased entry of CoV-2 spike mediated infection ( Figure 4C ; red bars) with the inhibitory concentration50 (IC50) of 2.8 ± 0.4 µM ( Figure   4D ; red line). rhBD-2 showed maximum inhibition of spike-mediated pseudoviral infection by 80% at a concentration of 12.8 µM and plateaued thereafter ( Figure 

Beta defensins are localized to the mucosa of the oral cavity, nares and upper airway Ghosh et al., 2007; Khurshid et al., 2017; Lee et al., 2002; Mathews et al., 1999; Singh et al., 1998) ; sites deemed vulnerable to CoV-2 entry and initial infection. Indeed, hBD-2 has been shown to interrupt viral infection of various viruses, including coronaviruses (Kim et al., 2018) . However, when a mucosal site becomes overwhelmed by a microbial threat, replenishing the AMP armamentarium locally after initial release takes multiple hours. Nasopharyngeal swab samples have revealed that CoV-2infected patients have lower mRNA levels of several defensins including hBD-2 (Idris et al., 2020) . Thus, hBD-2 or their mimetics, if administered exogenously, could be a sound therapeutic strategy to protect the host at vulnerable mucosal sites without eliciting an immunological response.

We chose to investigate hBD-2's ability to block CoV-2 because of its innate role in protecting the oral cavity and the upper airway, and because its mouse ortholog has been shown to inhibit other coronaviruses (Zhao et al., 2016) . The computer simulations showed remarkable stability of the RBD:hBD-2 complex even after 500 ns, with a clear overlap of binding sites when compared to the RBD:ACE2 complex. Complementary methods involving SPR, ELISA and immunoprecipitation/western blotting independently verified that hBD-2 binds to the RBD, thereby validating our in silico docking and dynamics simulations. Competitive inhibition assays were able to show that hBD-2 reduced RBD:ACE2 binding, possibly by directly binding to the RBD as shown in our biophysical and biochemical assays ( Fig   3A-3C) . Finally, by incorporating a luciferase reporter expressing CoV-2 spike-dependent lentiviral system , we demonstrated that hBD-2 inhibited viral entry into ACE2 expressing HEK 293T cells in a dose dependent manner, with an IC50 of ~2.8 µM. Additionally, hBD-2 begins to show hemolytic activity only at a concentration 25 times greater (70 µM) than our IC50 (Koeninger et al., 2020) , J o u r n a l P r e -p r o o f and shows no signs of cytotoxic effects against various other human cells at over twice our IC50 (Herrera et al., 2016; Mi et al., 2018; Sakamoto et al., 2005) ; suggesting a favorable therapeutic window for hBD-2 against CoV-2 at low micro molar concentrations as shown in this study.

It is important to note that a recent paper by Xu et al, that appeared while our manuscript was under revision, found no blocking activity of hBD-2 in their pseudovirion platform (Xu et al., 2021) . While we do not know exactly why our results and theirs differ, we can conjecture that this may have been due to differences in hBD-2 itself; i.e., their hBD-2 was synthesized while ours was recombinantly produced in E. coli. While both were presumably one isoform of hBD-2, they may not have been the same isoform.

Post-translational modification in folding of the molecule due to varying oxidative conditions that are not controlled in recombinantly generated products, can lead to a conformationally different molecule than one that is synthesized (Tokmakov et al., 2012; Wu et al., 2003) . Moreover, various isoforms have been reported for beta defensins, and interestingly, have been shown to vary in their biological activity (Wu et al., 2003) . Therefore, it would be very interesting to compare the two molecules and determine if indeed conformation plays a role in anti-CoV-2 blocking activity.

Of relevance to our study is the use of smaller mimetics of mouse beta defensin 4 (Zhao et al., 2016) , the ortholog of hBD-2, that when administered intra-nasally, rescued 100% of mice from the lethal challenge of human and avian influenza A, SARS-CoV and MERS-CoV (LeMessurier et al., 2016) . The fact that the peptide is endogenous to humans and would not elicit an immunogenic response, give it a high probability of being safe and a quicker route to human clinical trials. In fact, several AMPs, as well as AMP mimetics are currently undergoing clinical trials for multiple different diseases (Mookherjee et al., 2020) .

Recent biophysical studies using bilayer interferometry predicted a strong binding interaction between LL-37, another human AMP, and the RBD (Wang et al., 2021a) . This finding supports the idea that more than one AMP could be utilized, possibly in a "cocktail" to act as a potent blocking agent (Ghosh and Weinberg, 2021) . Recent findings also highlight that neuropilin-1 (NRP1), is being utilized by CoV-J o u r n a l P r e -p r o o f 2 to facilitate entry and infection (Cantuti-Castelvetri et al., 2020; Daly et al., 2020; Li and Buck, 2021 )and time will tell if also blocking entry via NRP1 will be required to reduce CoV-2 infection.

Ongoing studies have shown that CoV-2 has developed a number of mutations of which 89 have been associated with the RBD Wang et al., 2020b) . Furthermore, 52 out of 89 mutations are in the ACE2 receptor-binding motif (RBM), indicating that the virus may be accumulating mutations in that region to increase transmissibility (Li et al., 2020) . The flexibility of peptide based blocking strategies could permit customization to address mutations that could interfere with hBD-2's blocking efficacy (Agarwal and Gabrani, 2020) . Additionally, combining hBD-2, or its derivatives, that target the RBD with AMPs (e.g. human defensin 5) that target ACE2 (Wang et al., 2020a) may be efficacious in maintaining therapeutic potential.

The interaction of the RBD with ACE2, and especially with hBD-2, is highly dynamic at the protein interaction molecular level. Although this has not yet been measured in the RBD:ACE2 or RBD:hBD-2 platforms, the entropy of the interaction is likely to be not as unfavorable as seen in complexes where one or both partner proteins have to become significantly rigid. It is now becoming clear that many proteinprotein complexes are inherently dynamic Zhang and Buck, 2017) , thus minimizing the unfavorable entropy change that would otherwise occur on binding. In fact, this latter indication suggests that peptides, which are initially unstructured in the unbound state could also maintain considerable flexibility in the bound state and may thus be powerful antagonists of the RBD:ACE2

interaction. An alternative would be a mechanism involving entropy-enthalpy compensation, which is observed in many natural systems (Peccati et al, 2021) . Detailed thermodynamics analyses, both experimental and computational are needed to clarify this point. Nevertheless, our work and that of others suggests that design of AMP derived peptides and peptide mimetics could be a useful strategy in combating CoV-2 variants (Bakovic et al., 2021; Diamond et al., 2021; Zhang et al., 2021) .

We see an AMP strategy as complementary to the new anti-COV-2 vaccines. While they are highly effective, many people are refusing vaccination (Pogue et al., 2020; Schwarzinger et al., 2010 ) (Fisher et al., 2020; Wang et al., 2021b) . Additionally, a significant number either fail to mount effective J o u r n a l P r e -p r o o f neutralizing antibodies or high enough titers (Goodwin et al., 2006; Ndifon et al., 2009; Ovsyannikova et al., 2017) . Because of the multiple advantages of using small peptides like hBD-2 and their derived smaller mimetics, such as high specificity, low toxicity, lack of immunogenicity and ease of administration, they possess the potential for both safety and efficacy (Ghosh and Weinberg, 2021 ) Such molecules could be delivered, in the future, intra-orally and/or intra-nasally as prophylactic aerosols, in early stages of infection, when telltale symptoms appear and in combinatorial therapeutic approaches for more severe situations.

A limitation of the study is that the in vitro binding experiments were carried out on the RBD, rather than on the entire SARS-Cov-2 Spike protein. This was done because the difference in size between the entire S-protein and hBD-2 is so large that the SPR technique would have been rendered insensitive. Also docking and MD simulations would have been computationally prohibitive given the resources we have access to. Molecular docking and dynamics calculations are well known to still be limited by the accuracy of the potential function and by the extent of sampling; i.e., number of simulations and their length.

Although we cannot be certain that the calculations are highly converged, the results have key predictive value concerning the location and the dynamics of the binding surfaces. A second consideration is that the biophysical affinity measurements by SPR require the labeling of the His-tag with a surface attachment of the RBD via biotin . Although we do not believe this affects the results, as countless control studies on other proteins have shown, we cannot be completely certain. A third consideration is that we showed that hBD-2 inhibits CoV-2 infection by using an experimental system comprising of CoV-2 spike protein expressing pseudovirus particles and HEK293T cells overexpressing ACE2. CoV-2 spike expressing pseudovirus particles showed only modest infectivity in naturally ACE2 expressing cells such as the human lung cell line A549 and the monkey kidney epithelial cell line Vero, which was not sufficient to derive meaningful conclusions in the competitive inhibition assays. However, since the main objective values in (C); IC50 was calculated by plotting rhBD2 concentration (µM log) against % inhibition observed (n=4). Values are Mean ± SEM. ***p < 0.001, **p < 0.01, *p < 0.05, and ns (non-significant) against CoV-2 spike-pseudotyped virus alone treated group.

Further information and requests for resources and reagents should be directed to, and will be fulfilled by, the lead contact, Aaron Weinberg (axw47@case.edu).

This study did not generate new unique reagents.

The data reported in this paper will be shared upon request [Lead contact: Aaron Weinberg (axw47@case.edu)]. This paper does not report original code. Any additional information required to reanalyze the data reported in this paper is available from the lead contact upon request.

HEK 293T and HEK 293T cells stably expressing ACE2 receptor (ACE2 HEK293T) were cultured in DMEM media containing 10% heat inactivated FBS, 100 U/ml penicillin/streptomycin and 4 mM L-Glutamine.

J o u r n a l P r e -p r o o f pHAGE-CMV-Luc2-IRES-ZSgreen-W, HDM-HgPM2, HDM-tat1b, pRC-CMV-Rev1b, and SARS-CoV-2 Spike-ALAYT plasmids . FLAG and HIS tagged RBD were expressed from a pcDNA3 vector with leader sequence and leucine zipper (Ramakrishnan et al., 2004) .

The structure of human beta defensin 2 (hBD-2) in the monomer and dimer form is available in the PDB with ID 1FD3 (Hoover et al., 2000) . The hBD-2 sequence is 41 residues long:

GIGDPVTCLKSGAICHPVFCPRRYKQIGTCGLPGTKCCKKP. The five residues in bold were found to form hydrogen bonds with the RBD during the simulations (see below/main paper). The structure of the RBD domain of the Spike protein is also available in complex with ACE2 at 2.45 Å resolution in the PDB with ID 6M0J (Lan et al., 2020) and was used in reference simulations.

Two kinds of docking programs were applied; CLUSPRO (Kozakov et al., 2006; Kozakov et al., 2017; Porter et al., 2017; Vajda et al., 2017) and HADDOCK (Dominguez et al., 2003; van Zundert et al., 2016) . The x-ray structures of hBD-2 and of the SARS-CoV-2 S-protein RBD were uploaded to the CLUSPRO docking webserver without additional preparation. The hBD-2 : RBD docked structures were clustered, and the best of these clustered structures, in terms of the population of clusters, involved sites on the side of the RBD which are utilized in binding ACE2. The best structure was selected based on the docking programs' score. The CLUSPRO method is based on a Fast Fourier Transform correlation approach, which makes it feasible to generate and evaluate billions of docked conformations by simple scoring functions as shown in Equation (1). It is a multistage protocol: rigid body docking using PIPER (Kozakov Interaction energy E=0.40Erep-0.40Eatt+600Eelec+1.00EDARS (1) Here, Erep and Eatt are contributions of the van der Waals interaction energy, and Eelec is an electrostatic energy term. EDARS is a pairwise structure-based potential constructed by the Decoys as the Reference State (DARS) method (Chuang et al., 2008) . Since an entropic term was not included in CLUSPRO docking, the energy result should not be used to rank clusters and the population of clusters was used to rank them.

The population of clusters was applied to rank the clusters. In our simulations, the RBD:hBD-2 complex structure from the top cluster was taken and continued with all-atom molecular dynamics simulations.

By contrast to CLUSPRO, the docking HADDOCK had to be more selective and since the binding interface between the ACE2 receptor and RBD are known, residues from 400 to 520 on the RBD were selected as the target binding sites, while the entire hBD-2 peptide taken as a potential binding site.

Default values for all other parameters were applied. After that, the best 5 structures, by HADDOCK scoring, were selected.

Based on the best 6 (including above 5 from HADDOCK and 1 from CLUSPRO docking) structures predicted above, all-atom molecular dynamics simulations were set up using the CHARMM36m (Huang et al., 2017) forcefield and VMD program (Humphrey et al., 1996) . One of the deprotonated states of histidine was used (denoted HSD), and the native disulfide bonding in the hBD-2 was set up. After solvating the protein with an equilibrated box of TIP3P water molecules, the closest distance between atoms on the proteins and the edge of simulation box is 12 Å. The equivalent of 0.15 M in Na + and Cl -J o u r n a l P r e -p r o o f ions was added into the box plus several ions to neutralize the net charge of the system. The simulation temperature was kept at 310 K and pressure at 1 atm, using standard thermo-and barostats. After a brief energy minimization using the conjugate gradient and line search algorithm, 4 ps of dynamics was run at 50 K, and then the system was brought up to 310 K over an equilibration period of 1 ns using NAMD program version 2.12 (Phillips et al., 2005) . This was followed by trajectories that continued for up to 200 or 500 ns at 1 atm and 310 K using the NPT ensemble. The trajectory shown in Fig.1B and Fig. 2 was started from the best CLUSPRO docked structure. The final structure was restarted with three different random seeds (velocity assignments) for another 500 ns, denoted C1, C2 and C3 in Fig. S2 .

As a comparison, we also simulated the RBD bound with ACE2 using the structures from (Lan et al., 2020) and the same method as above. HBD-2 can also form a non-covalent dimer at high concentration in solution (Hoover et al., 2000) (with PDB ID of 1FD3). The initial bound structure of the hBD-2 dimer with the RBD was predicted using targeted HADDOCK docking. The best structure predicted was used in all-atom MD simulations as detailed above.

The simulation systems, set up, the number of atoms and box size information are shown in Table S1 .

To analyze the trajectories, the Root Mean Square Deviation (RMSD) and Fluctuations (RMSF) of the proteins were calculated using the VMD program and an in-house analysis script based on the coordinates of the backbone Ca atoms after aligning the trajectories respectively, to the original crystal structure of the RBD, hBD-2, and to the initial complex structure of the RBD and hBD-2 predicted from docking. The buried surface area (BSA) for the complex was calculated in two steps using the VMD program and a script using the Richards and Lee method with the water probe size of 1.4 Å (Lee and Richards, 1971) . First, the total solvent accessible surface area of the complex (ASAcomplex) was calculated based on the complex's trajectory. Second, the accessible surface area of each protein in the complex (ASArbd, ASAhbd2) was calculated for each protein individually. Then, the buried surface area, BSA is calculated using Equation (2) (2)

The number of hydrogen bonds between the RBD and ACE2 or the RBD and hBD-2 were calculated using the VMD program with the heavy atom distance cutoff of 3.0 Å and the angle cutoff of 20 degrees deviation from H-bond linearity. The time a particular H-bond is formed over the course of the simulation is monitored and is expressed as % occupancy. In order to find out the residues on the binding interface, the closest distance between every residue atom (including hydrogen) between the RBD and hBD-2 was calculated and averaged over the trajectory run. The average distances between each residue on RBD and on hBD-2 are shaded by proximity on a red to white color-scale and were used to build the distance maps. Due to the caveats associated with calculations of free energy estimations from trajectories such as the ones run for this study, we carried out the binding interaction energy calculation for RBD binding with ACE2 and hBD-2 monomer/dimer, respectively, using the popular MM-GBSA method (Genheden and Ryde, 2015) . This total pairwise interaction energy was calculated over the entire trajectory by NAMD energy plugin of the VMD and results are shown in Supplementary Table S1 , (Humphrey et al., 1996) .

This interaction energy ( E_binding ) is calculated using Equation (3):

E_complex is the potential energy of protein-ligand complex, E_protein is the potential energy of protein, and E_ligand is the potential energy of ligand. < > is the ensemble average over simulation time.

In the MM-GBSA method, the solvent effect was counted using the generalized Born implicit solvent model (GBIS) (Tanner et al., 2011) . Measuring RBD:ACE2 association in vitro J o u r n a l P r e -p r o o f Untagged recombinant hBD-2 and C-terminally His-tagged biotinylated recombinant RBD were purchased from Peprotech, Inc. and SinoBiological Inc., respectively. Earlier experiments with Histagged, but not biotinylated rRBDs indicated a slower second kinetic phase in SPR, which did not show saturation behavior (possibly aggregation of the rhBD-2) and could not be used for fitting. Using, surface plasmon resonance (SPR), the measurements were carried out at 25⁰C using a BIAcore T200 (Cytiva, formerly GE Healthcare) in PBSP buffer with 0.3mg/mL of BSA. Biotinylated RBD was captured to an S series SA sensor chip (Cytiva) at a density of 2600RU. Samples of the rhBD-2 at 90. 23, 180.47, 360.94, 721.88, 1443.75, 2887.5, 5775, 11550 , and 23100 nM were injected at 30 µL/min over the rRBD. Data were analyzed using the BIAevaluation software using two-state binding model and redrawn by Origin software. A second experiment was carried out a day later and the two Kd's were averaged.

100 l of rhBD-2 (Peprotech, Inc.) (concentration ranges 0.00024 to 2M) in assay diluent buffer 2 (R&D system) containing 1% BSA in PBS, were incubated in an RBD coated plate (Ray biotech, Inc.) at 4 0 C for 18 hrs. To test specificity of hBD-2 binding to RBD, we coated the microplate with LL37 (2ug/ml) as a negative control. Plates were then washed 4 times with 300 l of 1X Quantikine ELISA wash buffer (R&D Systems, Inc.) followed by incubation with 100 l of biotinylated anti-hBD-2 (Peprotech, Inc.) [ 0.1g/ml] for 1 hr. Plates were then washed again as stated above, incubated with 100l of Streptavidin-HRP (R&D system, Inc.) for 20 minutes. Signal was developed using TMB substrate and measured at 450nm using a microplate reader.

To study the interaction between recombinant hBD-2 (Peprotech, Inc.) and recombinant HIS-tagged-RBD (Sino Biologicals, Inc.), rhBD-2 (1.75 µM) and rRBD (0.19 µM) were preincubated in a buffer J o u r n a l P r e -p r o o f containing 20 mM HEPES, 5mM MgCl2 and 150 mM NaCl (solution A; 400ul) at room temperature for 1 hour. The binding reaction was then mixed with 1:1 volume of buffer containing 20 mM HEPES, 5mM

MgCl2, 150 mM NaCl, 0.5 mM DTT, 1% triton X 100 and 1mm EDTA (solution B) and incubated with Ni-NTA agarose resin beads (25 µl; Qiagen) overnight at 4°C. Standard IP protocols call for adding low amounts of DTT and triton X 100 to reduce possible non-specific cysteine reactions and non-specific sticking of reagents in the tube as well as to decrease protein degradation (possibly by decreasing protease activity). After overnight incubation, Ni-NTA beads were collected by centrifugation at 1000 rpm for 1 min and washed thrice with binding buffer. Beads were boiled with 30 µl of Laemmli sample buffer and were analyzed by Western blotting (WB). Briefly, samples were separated on 20% SDS-Polyacrylamide gels and proteins were then transferred to nitrocellulose membrane (0.2 µm pore size)

at 70V for 40 min in cold. Membranes were blocked with 5% milk in TBST and then probed with goat antihuman BD2 antibody (0.2 µg/ml; Peprotech), followed by secondary antibody (1:5000) at room temperature, and visualized by enhanced chemiluminescence. To study the ability of hBD2 to compete with RBD binding to ACE2, ACE2 HEK 293T cells were seeded in 6 cm plates. At 50% confluency, media was replaced with conditioned media from HEK 293T cells transfected with secreted HIS-RBD plasmid or control media in the presence or absence of rhBD-2 (0.7-2.1 µM) and incubated at 37°C for 30 min.

Cells were washed and collected in PBS-EDTA solution and then lysed in Triton lysis buffer. Lysates were centrifuged at 12000 g for 10 min at 4°C, and immunoprecipitated using Ni-NTA beads for 2 hours at 4°C. Beads were collected, washed, and boiled with Laemmli sample buffer and analyzed by Western blotting.

Pseudotyped SARS-CoV-2 spike virus was generated and the luciferase assay was carried . HEK 293T cells were transfected with luciferase-IRES-Zs Green, HDM-HgPM2, HDM-tat1b, PRC-CMV-Rev1b, and SARS-CoV-2 Spike-ALAYT plasmids . Culture J o u r n a l P r e -p r o o f supernatants were harvested 48 hours after transfection and used to infect ACE2 HEK293T cells. To perform in-vitro infections, cells were incubated for 48h with different dilutions of spike or VSVGpseudotyped virus in culture media to arrive at similar RLUs. Cells were lysed and luminescence was measured using the luciferase assay system following the manufacturer's instructions (Promega, Inc.) in Spectramax i3 microplate detection platform (Molecular Devices, Inc.). Based on the preliminary experiments, 1:0 dilution for spike-pseudotyped virus and 1:1000 dilution for VSVG-pseudotyped virus, which exhibits similar RLUs, was selected for further experiments. To study the effect of hBD-2 on spike pseudotyped virus entry, rhBD-2 and pseudovirus with similar RLUs were added to cells simultaneously.

ACE2 expressing HEK293T cells were co-incubated with 100 µl of spike-pseudotyped virus media (1:0 diluted in culture media) or VSVG-pseudotyped virus media (1:1000 diluted in culture media) and varying concentration of rhBD-2 (0-25.6 µM) with the final assay volume of 200 µl for 48 hours. Cells were lysed and luminescence was measured.

Values are expressed as mean ± SEM and evaluated statistically using Instat software. Significance between the groups was analyzed by performing one-way ANOVA with Tukey post-t-test. P < 0.05 was considered statistically significant. ***p < 0.001, **p < 0.01, *p < 0.05, and ns (non-significant).

LZ, SKG, PR, MB and AW are co-inventors in using AMPs, and their derivatives, as anti-coronavirus agents, and have submitted an invention disclosure to Case Western Reserve University. SCB, YC, PS, JP and AB declare no competing interests.

J o u r n a l P r e -p r o o f

Antiviral Peptides: Identification and Validation

Dynamics of the ACE2-SARS-CoV-2/SARS-CoV spike protein interface reveal unique mechanisms

Biomolecular Simulations in the Time of COVID19, and After

Brilacidin Demonstrates Inhibition of SARS-CoV-2 in Cell Culture

The SARS-CoV-2 Exerts a Distinctive Strategy for Interacting with the ACE2 Human Receptor

Neuropilin-1 facilitates SARS-CoV-2 cell entry and infectivity

Mutations Strengthened SARS-CoV-2 Infectivity

DARS (Decoys As the Reference State) potentials for protein-protein docking

Protocol and Reagents for Pseudotyping Lentiviral Particles with SARS-CoV-2 Spike Protein for Neutralization Assays

Neuropilin-1 is a host factor for SARS-CoV-2 infection. Science, eabd3072

Beta-defensins: what are they really doing in the oral cavity?

HADDOCK: a protein-protein docking approach based on biochemical or biophysical information

Attitudes Toward a Potential SARS-CoV-2 Vaccine : A Survey of

The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities

Critical Sequence Hotspots for Binding of Novel Coronavirus to Angiotensin Converter Enzyme as Evaluated by Molecular Simulations

Quantification of human beta-defensin-2 and -3 in body fluids: application for studies of innate immunity

Ramping Up Antimicrobial Peptides Against Severe Acute Respiratory Syndrome Coronavirus-2. Front Mol Biosci 8, 620806

ACE2, the Receptor that Enables Infection by SARS-CoV-2: Biochemistry, Structure, Allostery and Evaluation of the Potential Development of ACE2 Modulators

Molecular Mechanism of Evolution and Human Infection with SARS-CoV-2

Human beta-defensins 2 and -3 cointernalize with human immunodeficiency virus via heparan sulfate proteoglycans and reduce infectivity of intracellular virions in tonsil epithelial cells

The structure of human beta-defensin-2 shows evidence of higher order oligomerization

CHARMM36m: an improved force field for folded and intrinsically disordered proteins

VMD: visual molecular dynamics

Downregulation of Defensin genes in SARS-CoV-2 infection. medRxiv

Significance and Diagnostic Role of Antimicrobial Cathelicidins (LL-37) Peptides in Oral Health

Human β-defensin 2 plays a regulatory role in innate antiviral immunity and is capable of potentiating the induction of antigen-specific immunity

Human β-Defensin 2 Mediated Immune Modulation as Treatment for Experimental Colitis

PIPER: an FFT-based protein docking program with pairwise potentials

The ClusPro web server for protein-protein docking

Structure of the SARS-CoV-2 spike receptor-binding domain bound to the ACE2 receptor

The interpretation of protein structures: estimation of static accessibility

Antimicrobial defensin peptides of the human nasal mucosa

Carbohydrate-binding molecules inhibit viral fusion and entry by crosslinking membrane glycoproteins

Antimicrobial peptides alter early immune response to influenza A virus infection in C57BL/6 mice

Lentiviral vectors

The Impact of Mutations in SARS-CoV-2 Spike on Viral Infectivity and Antigenicity

Modified Potential Functions Result in Enhanced Predictions of a Protein Complex by All-Atom Molecular Dynamics Simulations, Confirming a Stepwise Association Process for Native Protein-Protein Interactions

Neuropilin-1 assists SARS-CoV-2 infection by stimulating the separation of Spike protein S1 and S2

Production of beta-defensin antimicrobial peptides by the oral mucosa and salivary glands

The Designer Antimicrobial Peptide A-hBD-2 Facilitates Skin Wound Healing by Stimulating Keratinocyte Migration and Proliferation

Antimicrobial host defence peptides: functions and clinical potential

Current scenario of peptide-based drugs: the key roles of cationic antitumor and antiviral peptides

Differential neutralization efficiency of hemagglutinin epitopes, antibody interference, and the design of influenza vaccines

A large population-based association study between HLA and KIR genotypes and measles vaccine antibody responses

Scalable molecular dynamics with NAMD

Influences on Attitudes Regarding Potential COVID-19 Vaccination in the United States. Vaccines (Basel) 8

ClusPro PeptiDock: efficient global docking of peptide recognition motifs using FFT

Human epithelial beta-defensins 2 and 3 inhibit HIV-1 replication

Receptor-specific signaling for both the alternative and the canonical NF-kappaB activation pathways by NF-kappaB-inducing kinase

Modulation of human beta-defensin-1 (hBD-1) in plasmacytoid dendritic cells (PDC), monocytes, and epithelial cells by influenza virus, Herpes simplex virus, and Sendai virus and its possible role in innate immunity

Differential effects of alpha-and beta-defensin on cytokine production by cultured human bronchial epithelial cells

The NMR structure of human beta-defensin-2 reveals a novel alpha-helical segment

Low acceptability of A/H1N1 pandemic vaccination in French adult population: did public health policy fuel public dissonance?

Production of beta-defensins by human airway epithelia

structural proteins of the severe acute respiratory syndrome coronavirus are required for efficient assembly, trafficking, and release of virus-like particles

Is the Rigidity of SARS-CoV-2 Spike Receptor-Binding Motif the Hallmark for Its Enhanced Infectivity? Insights from All-Atom Simulations

Characterization of the receptorbinding domain (RBD) of 2019 novel coronavirus: implication for development of RBD protein as a viral attachment inhibitor and vaccine

Parallel Generalized Born Implicit Solvent Calculations with NAMD

Multiple post-translational modifications affect heterologous protein synthesis

New additions to the ClusPro server motivated by CAPRI

The HADDOCK2.2 Web Server: User-Friendly Integrative Modeling of Biomolecular Complexes

Human Cathelicidin Inhibits SARS-CoV-2 Infection: Killing Two Birds with One Stone

Human Intestinal Defensin 5 Inhibits SARS-CoV-2 Invasion by Cloaking ACE2

Willingness to receive SARS-CoV-2 vaccine among healthcare workers in public institutions of Zhejiang Province, China

Decoding SARS-CoV-2 Transmission and Evolution and Ramifications for COVID-19 Diagnosis, Vaccine, and Medicine

Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation

Engineering disulfide bridges to dissect antimicrobial and chemotactic activities of human beta-defensin 3

Human Defensins Inhibit SARS-CoV-2 Infection by Blocking Viral Entry

The Proteins of Severe Acute Respiratory Syndrome Coronavirus-2 (SARS CoV-2 or n-COV19), the Cause of COVID-19

Dissociation of a Dynamic Protein Complex Studied by All-Atom Molecular Simulations

Molecular Dynamics Simulations Reveal Isoform Specific Contact Dynamics between the Plexin Rho GTPase Binding Domain (RBD) and Small Rho GTPases Rac1 and Rnd1

Antimicrobial peptide DP7 with potential activity against SARS coronavirus infections

A novel peptide with potent and broad-spectrum antiviral activities against multiple respiratory viruses

 HBD-2 binds spike-RBD at the ACE2 interaction site in silico. Biophysical and biological assays confirm hBD-2 binding to spike-RBD. HBD-2 blocks spike-RBD:ACE2 binding. HBD-2 prevents CoV-2/spike pseudovirions from infecting ACE2 expressing human cells.J o u r n a l P r e -p r o o f