key: cord-1044719-5ftpxlfe authors: Romano, Maria; Ruggiero, Alessia; Squeglia, Flavia; Berisio, Rita title: An engineered stable mini-protein to plug SARS-Cov-2 Spikes date: 2020-04-29 journal: bioRxiv DOI: 10.1101/2020.04.29.067728 sha: 856c001e04bdb4eaaae851843a85930818fdfddf doc_id: 1044719 cord_uid: 5ftpxlfe The novel betacoronavirus SARS-CoV-2 is the etiological agent of the current pandemic COVID-19. Like other coronaviruses, this novel virus relies on the surface Spike glycoprotein to access the host cells, mainly through the interaction of its Receptor Binding Domain (RBD) with the human angiotensin-converting enzyme 2 (ACE2). Therefore, molecular entities able to interfere with binding of the SARS-CoV-2 Spike protein to ACE2 have a great potential to inhibit viral entry. Starting from the available structural data on the interaction between SARS-CoV-2 Spike protein and the host ACE2 receptor, we here engineered a mini-protein with the aim of creating a soluble and stable Spike interactor. This mini-protein, which was recombinantly produced in high yields, possesses a stable α helical conformation and is able to interact with the RBD of glycosylated Spike protein from SARS-CoV-2 with nanomolar affinity, as measured by microscale thermophoresis. By plugging the Spike protein, our mini-protein constitutes a valid tool for the development of treatments against different types of coronavirus. The novel coronavirus SARS-CoV-2 has spread widely and rapidly since it was first identified in Wuhan, China in December 2019 [1] [2] [3] . Its associated disease, COVID-19, causes severe respiratory difficulties, with aged patients at higher risk of mortality [1] . At the time of writing, over 3.000.000 confirmed cases and over 200.000 deaths have been registered worldwide. Given the dramatic public health emergency, there is a strong and urgent need for new antiviral agents to block human-to-human transmission and to treat infected patients. Like other coronaviruses, SARS-CoV-2 makes use of a densely glycosylated Spike (S) protein to gain access into host cells [4, 5] . The S protein forms homotrimers protruding from the viral envelope and binds with high affinity to the host receptor ACE2 (angiotensin-converting enzyme 2), mainly expressed by epithelial cells of the respiratory tract [6] . Recently, another human receptor, CD147, has been identified as a possible route of viral entrance, again mediated by the S protein [7] . CD147 is also known as Basigin or EMMPRIN, and is a transmembrane glycoprotein that belongs to the immunoglobulin superfamily involved in many processes, including tumor development, plasmodium invasion and virus infection [7] . After attachment, the human transmembrane protease serine 2 (TMPRSS2) cleaves and activates the S protein, thus allowing SARS-CoV-2 to enter the host cells [8] . Compared to SARS-CoV, an additional protease, possibly furin, is likely involved in priming of the SARS-Cov-2, since the S protein of SARS-CoV-2 contains four redundant furin cut Pro-Arg-Arg-Ala motifs, which are absent in SARS-CoV [6] . The S protein contains two subunits, S1 and S2. Of these, S1 comprises the receptor binding domain (RBD), which is responsible for recognizing and binding the cell surface receptor ACE2 [5] . Being essential for infection, the S protein is a promising target for antibodies and vaccines [9] . Indeed, preventing attachment of the S protein to either the ACE2 or the CD147 receptors would hamper infection at the early viral entry step. Based on the structural information available on the complex between the S protein and the ACE2 receptor [10] , we engineered a miniaturised mimic of ACE2, here named Spikeplug. This mini-protein, which we produced in high yields, is highly soluble, conformationally stable and displays nanomolar binding affinity with the S protein of SARS-CoV-2. Given its properties, Spikeplug is a promising molecule for the development of treatments based on inhibition of viral entry for different types of coronavirus. We used all available structural information on the interactions between S protein and the human ACE2 receptor to design Spike interactors. An analysis of the structure of SARS-CoV-2 complex with ACE2, reveals that most hydrogen bonds between protein S and ACE2 involve the helix H1 (residues 20-52) and the C-terminal part of helix H2 (residues 56-82) of ACE2 [10] . Outside this region, an important interaction between the S protein and ACE2 is mediated by Lys353 of ACE2. Indeed, the presence of histidine in the place of lysine at position 353 of rat ACE2 disfavours viral binding [11] . Starting from the region H1-H2 of ACE2 as a basic scaffold for the design of a Spike interactor, we included a number of mutations to increase its stability and solubility and with the aim to compensate for the missing interactions mediated by residues outside helices H1 and H2. An important issue with peptide or small protein design is their conformational stability, since peptides are usually not folded in solution. Therefore, to stabilise this scaffold, we included an extra helix, H3 (residues 91-101) which naturally caps helices H1 and H2 in the ACE2 structure, through a cluster of hydrophobic interactions involving Val92, Leu97 and Leu100 of H3, Phe28, Leu29, Phe32 of H1 and Ala80 of H2, therefore holding together the two helices H1 and H2 ( Figure 1 ). In addition, we noticed that Glu37 is important to stabilise the structure of the ACE2 receptor by forming a salt bridge interaction with Arg393 and a hydrogen bond with the backbone N of Lys353 ( Figure 2A ). However, this residue forms a weak hydrogen bond with Tyr505, which is considered outside the Receptor Binding Motif SARS-CoV-2 [10] , while it is not engaged in any contact in the complex of ACE2 with RBD of the S protein of SARS-CoV (74% homologous to SARS-CoV-2 and 83% on the RBD domain) [11] . Therefore, we mutated Glu37 to an arginine residue and energy minimised using GROMACS [12] . As shown in Figure 2A , Arg37 side chain can form hydrogen bonds with the main chain of Gly496 of the S protein ( Figure 2A ). As such, the conformation of Arg37, also stabilised by a salt bridge with Asp38, mimics the interactions with the S protein mediated by Lys353 of ACE2 ( Figure 2A ) [10, 11] . We further engineered this variant by including two mutations, Leu 91 to Gln and Leu 95 to Gln, to improve protein solubility by replacing hydrophobic with hydrophilic residues, and hamper aggregation in solution. In ACE2, these two residues are located on the opposite side of the S-binding region and form hydrophobic interactions with the core of ACE2 to stabilise it ( Figure 2B ,C). Their mutation to glutamine residues was also motivated by the characteristics of glutamine to further stabilize α helix structures. This recombinant mini-protein, here named Spikeplug (Table 1, Figure 2C ), was successfully over-expressed in E.coli, resulting in a high yield, of 70 mg of pure protein from one litre of bacterial culture. Using far-UV CD spectroscopy, we observed that the spectrum of Spikeplug is typical of a wellstructured fold with a high α-helical content ( Figure 3A) , with typical minima at 208 and 222 nm ( Figure 3 ). Spectra also showed that folding is fully reversible, with the CD spectrum after refolding fully superimposable to that recorded at 4°C ( Figure 3A ). To investigate the heat-induced changes in the protein secondary structure, thermal unfolding curves were recorded by following the CD signal at 208 nm as a function of temperature, thus providing a melting temperature Tm of 37°C ( Figure 3B ). Microscale thermophoresis (MST) analysis was performed to measure the binding affinity of Blocking the very first step of SARS-CoV-2 entry into host cells by hampering the interaction of the S protein with the cellular receptor, represents a highly promising therapeutic strategy against COVID-19. We focused on generating a Spike interactor that resembles the human ACE2 receptor as we argued that such type of molecule would act against all known coronaviruses using this receptor, such as SARS-CoV and HCoV-NL63, as well as other coronavirues possibily emerging in the future. Several structures of the complex between the S protein and the ACE2 receptor have been recently published. Based on this structural information, we designed a mini-protein which embeds all important Spike-interacting residues of the ACE2 receptor. This miniaturised mimic of ACE2, was properly engineered to be stable in solution and to enhance its affinity with the S protein. Indeed, an important issue with peptide-or protein-based drugs is its conformational stability. This mini-protein, Spikeplug, presents a highly stable α helical conformation in solution and a nanomolar binding affinity to the RBD domain of glycosylated S protein from SARS-CoV-2. The dissociation constant KD determined here using MST, 40 nM, is in the same range of the affinity measured between the S protein and the ACE2 receptor of 14.7 nM [5] . We are currently designing additional variants to further improve binding affinity to the S protein. Our basic idea is that mimics of the host receptor have a great advantage over other drugs for several reasons. First, if a novel coronavirus will emerge, able to infect humans through the ACE2 receptor, then its S proteins will be recognised and plugged by our receptor mimics ( Figure 5) , blocking viral entry. Importantly, Spikeplug is large enough to cover most important interactions with the S protein, but small enough to bind simultaneously to the three chains of the S protein, thus increasing its potential to hamper interactions with the ACE2 receptor ( Figure 5 ). Another point we considered in our research was drug resistance, as pathogens can easily mutate to escape therapeutic treatments. However, if the virus mutates to escape binding to Spikeplug, then it will also likely drop its affinity for the natural ACE2 receptor, resulting in suicide action. The validity of such an approach has been demonstrated by the development of the antiretroviral drug Enfuvirtide, which blocks the action of the HIV-1 gp41 fusion protein, thus preventing viral entry. As for all antiviral drugs, also Enfuvirtide can select for resistance mutations. However, it has been shown that resistance to the drug comes at a serious fitness cost for the virus since it involves structural alterations in an essential component of the virus-host cell fusion complex [13] . Hence, targeting the receptor binding domain of SARS-CoV-2 might also result in a high genetic barrier towards resistance. Last, we expect Spikeplug to be well tolerated because of its almost complete identity with the correspondent domains of the endogenous human protein. In conclusion, we believe that Spikeplug is a first-in class promising lead candidate for the development of molecular targeted therapies against SARS-CoV-2 and other coronaviruses. Figure 5 . Cartoon representation of a plugged Spike trimer. The structure of the S protein (pdb code 6vxx) [9] is reported in magenta, prune and blue. The position of Spikeplug molecules (green surface representation) on the S protein was determined upon superposition of RBD domain of our minimised Spikeplug-RBD complex on the RBD domain of the entire S protein (pdb code 6vxx). A and B panels report side and top views, respectively. Molecular modelling was initially performed using the cryo EM structure of the SARS-CoV-2 S protein (PDB code 6vsb) [5] and the crystal structure of the complex between ACE2 and the RBD domain of SARS-CoV S protein (PDB code 2ajf) [11] . To complete missing regions in the cryo EM structure, we computed the homology model of SARS-CoV-2 S protein using MODELLER [14] , and the structure 2ajf as a template. While this work was in progress, the crystal structure of the complex between ACE2 and the RBD domain of the S protein from SARS-CoV-2 (PDB code 6m0j) was released [10] . In this structure, most of the interactions between the S protein and ACE2 are conserved. Given the high sequence identity between RBD of SARS-CoV and SARS-CoV-2 (84%), we relied on both structures to take also possible crystal packing bias into account. Mutations in the basic scaffold were generated using the software Coot [15] . Models were energy minimised using the GROMACS package [12] . Figures were generated with Pymol [16] . The gene encoding Spikeplug was synthesised and codon-optimised for E. coli expression by GeneArt® Gene Synthesis company (Invitrogen) and subsequently sub-cloned into pETM-13 vector (EMBL, Germany) between the NcoI and XhoI restriction sites. This plasmid endows the placement of a histidine-tag at the C-terminal of the mini-protein. The expression level was optimised after carrying out small scale expression screening using different E. coli available strains. From the screenings, the recombinant protein was successfully over-expressed in E.coli BL21(DE3) cells, resulting in a yield of 70 mg of pure protein from one litre of bacterial culture. Briefly, an overnight starting culture of 10 mL was prepared for growth in 1L of Luria Bertani (LB) medium containing 50µg L −1 kanamycin, which was then induced with 0.8 mM of IPTG at 16°C for approximately 16 h. The protein was purified by sonicating bacterial cells resuspended in binding buffer (300mM NaCl, 50mM Tris-HCl, 10 mM imidazole, 2.5% (v/v) glycerol, pH 7.8) containing a protease-inhibitor cocktail (Roche Diagnostics). The lysate was cleared by centrifugation at 18000 rpm at 4°C and the supernatant was loaded onto 5 mL Ni-NTA resin (Qiagen, Milan, Italy) equilibrated with binding buffer. After washing with ten volumes of binding buffer, the protein was eluted by adding 300 mM imidazole to the binding buffer. The fractions containing the eluted protein were pooled, concentrated and then loaded onto a Superdex 75 HR 10/30 gel-filtration column (GE Healthcare) equilibrated with 150mM NaCl, 50mM Tris-HCl, 2.5% (v/v) glycerol pH 7.8 for a further purification step. The sample eluted in a single peak and was homogeneous, as judged by 18% SDS-PAGE analysis. The protein was concentrated using a centrifugal filter (Merck Millipore) and the concentration was determined by UV absorbance using the corresponding ɛ values (M -1 cm -1 ). CD measurements were carried out using a JASCO J-815 CD spectropolarimeter equipped with a The thermophoretic measurements were performed using Monolith NT.115 device with red detection channel (NanoTemper Technologies, Munich, Germany). Recombinant Spikeplug was produced in our laboratory, while SARS-CoV-2 Spike RBD was purchased from Sino Biological Inc. For MST recording, the SARS-CoV-2 Spike RBD was labelled with the fluorescent dye NT647 using the protocol suggested from the NanoTemper. Thermophoretic experiments were conducted using Monolith NT.115 (NanoTemper Technologies, Munich, Germany). Briefly, 90 µl of a 10 µM solution of SARS-CoV-2 Spike protein (RBD) in labelling buffer (130 mM NaHCO3, 50 mM NaCl, pH 8.2) was mixed with 10 µl of 300 µM NT647-N-hydroxysuccinimide fluorophore (NanoTemper Technologies) and incubated in the dark for 30 min at room temperature. Spike RBD concentration after labelling (NT647-Spike RBD) was measured using a UV-Vis spectrophotometer and the labelling efficiency was determined to be about 25%. The MST experiment was performed in a buffer containing 200 mM NaCl, 50mM Tris-HCl, 0.05% A pneumonia outbreak associated with a new coronavirus of probable bat origin Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding Spike protein recognition of mammalian ACE2 predicts the host range and an optimized ACE2 for SARS-CoV-2 infection Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation SARS-CoV-2 receptor ACE2 and TMPRSS2 are primarily expressed in bronchial transient secretory cells SARS-CoV-2 invades host cells via a novel route: CD147-spike protein SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor Structure, Function, and Antigenicity of the SARS-CoV-2 Spike Glycoprotein Structure of the SARS-CoV-2 spike receptor-binding domain bound to the ACE2 receptor Structure of SARS coronavirus spike receptorbinding domain complexed with receptor GROMACS: fast, flexible, and free Enfuvirtide resistance mutations: impact on human immunodeficiency virus envelope function, entry inhibitor sensitivity, and virus neutralization Protein Structure Modeling with MODELLER Coot: model-building tools for molecular graphics PyMOL and Inkscape Bridge the Data and the Data Visualization We thank Prof. Giovanni Maga, Institute of Molecular Genetics IGM-CNR, for careful reading of this manuscript.