Next Article in Journal
TGFβ-Treated Placenta-Derived Mesenchymal Stem Cells Selectively Promote Anti-Adipogenesis in Thyroid-Associated Ophthalmopathy
Next Article in Special Issue
Determination of IgG1 and IgG3 SARS-CoV-2 Spike Protein and Nucleocapsid Binding—Who Is Binding Who and Why?
Previous Article in Journal
Antiretroviral Therapy-Induced Dysregulation of Gene Expression and Lipid Metabolism in HIV+ Patients: Beneficial Role of Antioxidant Phytochemicals
Previous Article in Special Issue
Expression of Oxidative Stress and Inflammation-Related Genes in Nasal Mucosa and Nasal Polyps from Patients with Chronic Rhinosinusitis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Structure-Based Development of SARS-CoV-2 Spike Interactors

1
Institute of Biostructures and Bioimaging, IBB, CNR, 80131 Napoli, Italy
2
Department of Pharmacy, University of Salerno, 84084 Fisciano, Italy
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(10), 5601; https://doi.org/10.3390/ijms23105601
Submission received: 21 April 2022 / Revised: 13 May 2022 / Accepted: 14 May 2022 / Published: 17 May 2022
(This article belongs to the Collection Feature Papers in “Molecular Biology”)

Abstract

:
Coronaviruses, including SARS-CoV-2 (the etiological agent of the current COVID-19 pandemic), rely on the surface spike glycoprotein to access the host cells, mainly through the interaction of their receptor-binding domain (RBD) with the human angiotensin-converting enzyme 2 (ACE2). Therefore, molecular entities able to interfere with the binding of the SARS-CoV-2 spike protein to ACE2 have 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 engineered a set of soluble and stable spike interactors, here denoted as S-plugs. Starting from the prototype S-plug, we adopted a computational approach by combining stability prediction, associated to single-point mutations, with molecular dynamics to enhance both S-plug thermostability and binding affinity to the spike protein. The best developed molecule, S-plug3, possesses a highly stable α-helical con-formation (with melting temperature Tm of 54 °C) and can interact with the spike RBD and S1 domains with similar low nanomolar affinities. Importantly, S-plug3 exposes the spike RBD to almost the same interface as the human ACE2 receptor, aimed at the recognition of all ACE2-accessing coronaviruses. Consistently, S-plug3 preserves a low nanomolar dissociation constant with the delta B.1.617.2 variant of SARS-CoV-2 spike protein (KD = 29.2 ± 0.6 nM). Taken together, we provide valid starting data for the development of therapeutical and diagnostic tools against coronaviruses accessing through ACE2.

1. Introduction

The 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]. Given the dramatic public health emergency and the growing number of variants that have emerged so far, 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 protein to gain access into host cells [4,5]. The spike protein forms homo-trimers 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]. After attachment, the human transmembrane protease serine 2 (TMPRSS2) cleaves and activates the spike protein, thus allowing SARS-CoV-2 to enter the host cells [7]. Compared to SARS-CoV, an additional protease, possibly furin, is likely involved in priming the SARS-CoV-2, because the spike protein of SARS-CoV-2 contains four redundant furin cut Pro-Arg-Arg-Ala motifs that are absent in SARS-CoV [6]. Once inside the cell, the infecting RNA acts as a messenger RNA, which is translated by host ribosomes to produce components necessary to assemble the new viral particles [8,9].
The spike protein contains two subunits, S1 and S2. Of these, S1 comprises the receptor-binding domain (RBD), which is responsible for recognising and binding the cell surface receptor ACE2 [5]. As it is essential for infection, the spike protein has been exploited as a target for antibodies and vaccines [10]. It is the principal target of the neutralising antibodies generated following infection by SARS-CoV-2 [11,12]. Consistently, it is the antigen of both mRNA and adenovirus-based vaccines licensed for use and others awaiting regulatory approval [13]. As SARS-CoV-2 continues to spread and cause diseases, emerging variants are being identified around the globe. As a consequence, there is emerging evidence of reduced neutralisation of some SARS-CoV-2 variants [13]. In addition, regarding the effectiveness of neutralising antibodies, their therapeutic and diagnostic use has several limitations associated with their large molecular mass, which often causes poor tissue penetration [14,15]. These features prompted researchers to develop smaller and less expensive protein scaffolds as alternatives to antibodies [15]. In addition to this, the large number of spike mutations so far observed [16], and the potential impact of mutations on the antigenicity of the spike protein, still prompt a strong interest in identifying spike-interacting molecules with high resemblance to the ACE2 receptor as a tool to target all variants that are able access the host through the ACE2 receptor, both for therapeutic and diagnostic use.
Based on the structural information available on the complex between the spike protein and the ACE2 receptor [17], we developed high-affinity binders to the spike with the property to share the same interacting interface with ACE2. The first developed molecule, S-plug, was highly soluble and displayed nanomolar binding affinity with the spike protein of SARS-CoV-2. By combining mutational stability prediction analysis with molecular dynamics (MD) simulations, we identified possible mutations to enhance S-plug thermostability. The resulting proteins, S-plug2 and S-plug3, displayed significantly enhanced melting temperatures and low nanomolar dissociation constant (KD) to the spike protein.

2. Results

2.1. Design of S-plug and Thermostabilised S-Plugs

An analysis of the structure of SARS-CoV-2 complex with ACE2, reveals that most hydrogen bonds between the spike protein and ACE2 involve the helix H1 (residues 20–52) and the C-terminal part of helix H2 (residues 56–82) of ACE2 [17]. Starting from the region H1-H2 of ACE2 as a basic scaffold for the design of a spike interactor, we included several mutations to increase its stability and solubility, with the aim to compensate for the missing interactions that are mediated by residues outside helices H1 and H2. 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, and Phe32 of H1; and Ala80 of H2, holding together the two helices, H1 and H2 (Figure 1). We further engineered this variant by including two mutations, Leu91 to Gln and Leu95 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 1B,C). Their mutation to glutamine residues was also motivated by the characteristics of glutamine to further stabilise α-helix structures.
An important issue in peptide or small protein design is thermal stability, because peptides are usually not folded in solution. Therefore, our further goal was to identify thermo-stabilising mutations, located in the region that does not interact with the spike protein, to preserve its binding. Using Rosetta mutation analysis [19], we identified a series of predicted thermo-stabilising mutations, acting either by increasing the buried hydrophobic surface area, or by increasing the quality of side-chain packing. Of the identified mutations, 10 were located at the interface between S-plug and RBD and were not applied to prevent destabilisation of the S-plug–RBD complex. To validate the effect of the mutations identified as stabilising (but not disruptive to the complex) on the structural integrity and flexibility of the proteins, we combined this strategy with MD simulations of S-plug and its mutated form, S-plug2 (Table 1).

2.2. MD Simulations

Two independent 2 µs MD simulations computations were performed, starting from complexes of the RBD of the spike protein (residues 336-516) with either S-plug or S-plug2. The evaluation of root-mean-square deviations (RMSD) (calculated on the protein Cα atoms) between the starting models and the trajectory structures shows limited conformational variations of the entire system, with RMSD values close to four Å in both simulations (Figure S1A,B). However, when focusing on the S-plug chains, the analysis of secondary structure (SS) shows a clear loss of α-helical conformation for the C-terminal part of H1 (residues 40–51) and the N-terminal of H2 (residues 57–65) of S-plug in the MD simulation of RBD–S-plug (Figure 2A). By contrast, stable SS evolution characterises the MD simulation of the RBD–S-plug2 complex (Figure 2B). Consistently, the number of hydrogen bonds stabilising α-helices (n + 4) of S-plugs strongly decreases throughout the MD simulation trajectory of RBD–S-plug, whereas it remains constant in the RBD–S-plug2 MD simulation (Figure 2C). Additionally, the RMSD between the starting models and the trajectory structures, computed on the 40-65 region, shows a remarkable increase in the RBD–S-plug MD trajectory (from one to five Å), whereas it remains unaltered (close to one Å) in the RBD–S-plug2 MD trajectory (Figure S1C,D). Importantly, the region 40-65 is not predicted to affect the binding affinity of spike interactors because the corresponding region of ACE2 is not involved in spike binding. Therefore, better organisation of this region in an α-helical arrangement in S-plug2 is likely to strongly contribute to protein stability without decreasing its binding affinity to the spike protein.

2.3. S-Plugs Folding and Thermostability

Designed miniproteins were successfully over-expressed in E. coli, resulting in high yields of pure proteins. Far-UV CD spectroscopy spectra of S-plug and S-plug2 are typical of well-structured folds with high α-helical content, with typical minima at 208 and 222 nm (Figure 3). Spectra also show that folding is fully reversible, with the CD spectrum after refolding fully superimposable to that recorded at 20 °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 222 nm as a function of temperature, using a 1 °C/min heating rate. Thermal unfolding curves of S-plug and S-plug2 exhibit the characteristic sigmoidal profile expected of two-state systems (Figure 3A,B). Consistent with MD computations, S-plug2 displays a thermo-stabilised structure, with an increase in the melting temperature (Tm) from 37 °C of S-plug to 48 °C (Figure 3A,B).
Analysis of the oligomeric states of S-plug and S-plug2 shows that S-plug exists in solution in equilibrium between a monomeric and dimeric state, whereas S-plug2 is mostly dimeric (Figure 3C). This result indicates that part of the observed protein stabilisation achieved is due to the stabilisation of the dimeric state in S-plug2. Using the docking procedure implemented in Rosetta software, [19] we computed a three-dimensional model of the S-plug2 dimer (Figure 4). As shown in Figure 4, the S-plug2 dimer leaves the spike-interacting regions on the two opposite sides of the molecule, as two molecules are placed back-to-back. In this arrangement, the C-termini of the two monomers are in close proximity. Therefore, we further modified S-plug2 sequence to add a C-terminal arm ending with a cysteine residue (with sequence PPGC, Table 1), to induce a covalent bond formation (Figure 4). The miniprotein with the resulting sequence, S-plug3, was over-expressed in E. coli and checked for its dimeric state by mass spectrometry. As reported in Figure 3D, CD spectroscopy data show that the cysteine-lock induced in S-plug3 further stabilised the molecule, as its denaturation temperature further increased from 48 to 54 °C. All developed miniproteins were assayed for their ability to interact with the spike protein.

2.4. Binding Analyses

Real-time binding assays were performed to assess binding kinetics and affinity between developed miniproteins and either the RBD or the S1 domain of the spike protein from SARS-CoV-2, using Biacore T200 (Cytiva, Uppsala, Sweden). The spike RBD-His was stably captured at the surface of the CM5 sensor chip using an anti-histidine antibody that had been covalently bound to the surface, as recommended by the manufacturer. The SARS-CoV-2 spike protein (S1) was immobilised on a CM5 sensorchip using standard amine-coupling protocols (see methods). Following immobilisations, protein solutions (S-plug, S-plug2, S-plug3) were injected at various concentrations (from 0.005 to 5µM), using a flow rate of 30 μL/min for 120 s (association phase), and then the buffer alone for 600 s (dissociation phase). Equilibrium dissociation constant (KD) values were derived from the ratio between kinetic dissociation (kd) and association (ka) constants, obtained by fitting data from all injections at different concentrations of each compound, using the simple 1:1 Langmuir binding fit model of the BIAevaluation software (version 2.0.2).
Results showed strong binding of S-plug to either RBD or S1, in the range 50–200 nM. However, low reproducibility characterised these experiments, depending on the preparation and time of storage. This behaviour is likely due to the conformational variability we observed for this construct, both in terms of secondary structure stability and of oligomerisation state (Figure 2 and Figure 3). Highly reproducible binding was observed for all other miniproteins. In particular, SPR analysis shows that S-plug2 binds to either the spike RBD or the S1 domain with elevated affinity, with KD values of 61.7 ± 2.3 nM and 71.7 ± 1.8 nM, respectively (Table 2, Figure S3). The strongest binding characterises S-plug3, with KD to RBD of 39.5 ± 1.1 nM and to S1 of 31.7 ± 1.8 nM (Table 2, Figure 5). Given this result, we also tested the binding of S-plug3 to the delta variant of spike RBD (B.1.617.2, carrying the mutations L452R and T478K), one of the deadliest variants so far reported [20]. Results show strong binding of S-plug3, with KD of 29.2 ± 0.6 nM (Figure 5). We also measured the affinity of spike S1 domain to the ACE2 receptor as a control experiment, which provided a dissociation constant KD of 51.9 ±1.4 nM (Table 2 and Figure S2).

3. Discussion

Several structures of the spike protein in different functional states and of complexes between the spike protein and the ACE2 receptor have recently become available [5,17,21,22,23,24]. In the prefusion state, the RBDs of the spike protein alternate between open ‘up’ and closed ‘down’ conformations. Consequently, the receptor-binding site, which can bind to human angiotensin-converting enzyme 2 ACE2, is transiently exposed in the ‘up’ conformation [25]. Many neutralising antibodies are directed against the RBD, with the highest neutralising activity associated with antibodies binding the RBD in the down state [11,12,26]. However, severe limitations are associated with the use of antibodies as therapeutics due to their typical suboptimal pharmacokinetics [27]. Their use is also limited for diagnostic purposes, due to the elevated costs of production. Miniaturisation represents a major direction in the development of rationally designed, less expensive protein scaffolds as alternatives to antibodies [15,27].
Based on the available structural information of the spike–ACE2 complex, we developed miniproteins which embed key spike-interacting residues in the ACE2 receptor. After the development of our prototype miniprotein spike binder, S-plug [28], one study reported the development of other miniproteins binding to the spike with high affinity; the best ranking showed low nanomolar binding affinity and was based on large-scale de novo design [29]. We started from our prototype S-plug [28], containing H1-H2 and H3 helices of ACE2, with mutations on the H3 helix (L91Q and L95Q) to improve its solubility (Figure 1). Then, our protein-engineering goal was to identify thermo-stabilising mutations located in regions that do not interact with the spike protein, to preserve its binding mode. The rationale of this study is to generate spike interactors that resemble the human ACE2 receptor, as we expect that such molecules would recognise all known coronaviruses accessing the host through this receptor, such as SARS-CoV, HCoV-NL63, as well as other coronaviruses that may emerge in the future.
Flexible regions are typically considered “hot spots” for engineering protein thermostability. Using MD simulations, we predicted that the introduced mutations induce a higher stability of α-helical and loop regions (residues 40–65), which are not involved in spike binding. Consistently, we recombinantly produced the resulting miniprotein, S-plug2, and showed that it presents a highly stable α-helical conformation in solution, with a Tm of 48 °C, 11 °C higher than the original S-plug (Tm 37 °C). Additionally, S-plug2 displays a low nanomolar dissociation constant KD, determined here using surface plasmon resonance, to both the RBD of glycosylated S protein and to the entire S1 region of SARS-CoV-2 (Figure 5 and Figure 6). Based on docking analysis of S-plug2, which presents a dimeric organisation in solution, we further improved it by adding a C-terminal arm ending with a cysteine residue to induce a cysteine lock. This protein, which we named S-plug3, is more stable than S-plug2, with a Tm of 54 °C, and displays enhanced binding affinities with both the spike S1 domain and RBD (Figure 6). This affinity is comparable to those hitherto measured by SPR with the ACE2 receptor (Table S2) [5,17,24]. Additionally, molecular dimensions of S-plug3 are compatible to full saturation of the RBDs in the spike trimer (Figure 7).
In conclusion, we designed and developed highly stable spike interactors, sharing with the ACE2 receptor a highly resemblant spike interface, aimed at the recognition of all coronaviruses accessing the host through the ACE2 receptor. Consistently, we show that our best performing miniprotein, S-plug3, binds to the deadliest delta variant of spike RBD with a dissociation constant KD of 29.2 (±0.6) nM. The most immediate application of these molecules is for diagnostic purposes, as the validation of their therapeutic properties still awaits a long experimental assessment. Taken together, data reported here provide a protocol for the design and development of stable miniproteins that mimic larger host receptors and valid starting tools for the development of therapeutic and/or diagnostic measures against ACE2-accessing coronaviruses.

4. Experimental Methods

4.1. Computational Design of S-Plug and the Thermo-Stabilised S-Plugs

The crystal structure of the complex between ACE2 and the RBD of the S protein from SARS-CoV-2 (PDB code 6m0j) was used as a starting model to design the spike interactors [17]. Interactions between the spike protein and ACE2 were computed to identify the smallest region of the ACE2 receptor able to bind the S-protein. This region, including residues 19-103 of the ACE2 receptor, was suitably mutated to enhance protein solubility (Table 1). Mutations in the basic scaffold were generated using the software Coot [30] and models were energy minimised using the GROMACS package to check the structural effect of mutation [31].
To design spike interactors with increased thermal stability, we used a combined computational approach including (i) structure-guided rational design methods, using Rosetta [32] with (ii) molecular dynamics approaches, using GROMACS [31]. Starting from the model of S-plug, generated as explained above, Rosetta software was used to predict the relative change in folding free energy of residues in flexible regions due to point mutations [32]. Key residues for interactions with RBD were left unaltered. Inserted mutations were selected to either increase buried hydrophobic surface area or increase the quality of side-chain packing. To predict structural stability changes of the spike interactors, we adopted MD dynamics, starting from either S-plug or the modified variant (S-plug2) in complex with spike RBD.

4.2. Molecular Dynamics

Simulations were performed by the GROMACS package [31], using protocols already tested for other protein systems [33,34,35]. Starting models were obtained as described above. For each starting model, 2 µs MD simulations were carried out using the Amber99sb all-atom force field [36,37] with tip4pew water model. All simulations were run in the isobaric-isothermal (NPT) ensemble at 300 K, using periodic boundary condition.
The starting models were solvated in a dodecahedron box with minimal distance of the model to the box wall of 11 Å. Sodium or chloride ions were added by replacing water molecules to neutralise the overall charge. Details are listed in Table S1. Before starting MD simulations, the systems were subjected to energy minimisation, then the solvent was equilibrated first in NVT and then in NPT ensemble for 100 ps, during which the protein atoms were restrained to the energy-minimised initial coordinates. Subsequently, the thermalisation process consisted of a sequence of 100-ps runs where the temperature was increased from 50 K to 300 K in 50 K steps. At this stage, the Berendsen method for pressure control was used with a coupling constant τp = 1.0 ps, whereas the temperature was kept constant at 300 K by the V-rescale method with a τt = 0.1 ps. For the production run, the pressure coupling was employed using the Parrinello–Rahman barostat with τp = 2.0 ps. The integration time step was 0.002 ps and LINCS was the constraint algorithm used [38]. The electrostatic interactions were treated using the particle-mesh Ewald (PME) method with a real-space cutoff of 10 Å [39]. Lennard-Jones interactions were truncated at 10 Å [40]. All the trajectories were analysed using GROMACS routines. Structural properties, such as root mean-square deviation (RMSD), root mean-square fluctuation (RMSF), and secondary structure (SS), were calculated with GROMACS or Visual Molecular Dynamics (VMD) [41] standard analysis tools.

4.3. Protein Expression and Purification

The genes encoding the S-plug and S-plug2 proteins were synthesised by GeneArt® (ThermoFisher), and subsequently subcloned in pETM-11 expression vector (EMBL, Germany) between the NcoI and XhoI restriction sites. The S-plug3 variant was produced by polymerase chain reaction (PCR), by introducing a spacer consisting of three amino acids and terminating a cysteine residue (-PPGC) at the protein C-terminal end, using the following oligonucleotides: Forward 5′-CATGCCATGGCGAGCGCGATGG-3′; Reverse 5′-CCGCTCGAGTTAGCAGCCCGGCGGGCTGCTCGCGGTTTTTTCC-3′. The PCR product was cloned in pETM-11 expression vector. All proteins were expressed with a TEV-cleavable N-terminal histidine tag. The expression levels were optimised after carrying out small-scale expression screening using different strains of E. coli. The recombinant proteins were successfully over-expressed in E. coli BL21(DE3) cells. Briefly, an overnight starting culture of 10 mL was prepared for the growth with 1L of Luria Bertani (LB) medium containing 50 µg/mL kanamycin, which was then induced by adding 0.8 mM of IPTG at 18 °C for approximately 16 h. For protein purification, a similar procedure was adopted. The harvested bacterial cells were resuspended in the binding buffer (50 mM Tris-HCl, 300 mM NaCl, 5% (v/v) glycerol and 10 mM imidazole, pH 7.8) containing a complete protease inhibitor cocktail (Roche Diagnostics, Mannheim, Germany), and then the cell suspension was lysed by sonication. The lysate was cleared by centrifugation at 18,000× g at 4 °C and the supernatant was loaded onto 5 mL of Ni–NTA resin (Qiagen, Milan, Italy) equilibrated with the binding buffer. After washing with ten volumes of binding buffer, the protein was eluted by increasing the imidazole concentration. The fractions containing the eluted protein were pooled and concentrated for TEV cleavage, which was performed overnight at 4 °C. After the tag cleavage, proteins were further purified by gel filtration using a Superdex 75 increase 10/30 column (GE Healthcare) equilibrated with 50 mM Tris-HCl, 150 mM NaCl, and 2.5% (v/v) glycerol, pH 7.8. The purity of proteins was assessed by 18% SDS–PAGE analysis. The proteins were concentrated using a centrifugal filter (Merck Millipore) and the concentration was determined by UV absorbance using the corresponding ɛ values (M−1 cm−1).

4.4. CD Spectroscopy

CD measurements were obtained using a JASCO J-715 CD spectro-polarimeter equipped with a Peltier temperature controller (Model PTC-423S) at 4 °C in a 0.1 cm optical path length cell in the 200–260 nm wavelength range. Protein concentration was 0.2 mg/mL in 20 mM sodium phosphate buffer (pH 7.4). The data were recorded with a scanning speed of 20 nm/min and a band width of 1 nm. All spectra were averaged from three scans and baseline-corrected using a blank consisting of the protein buffer. The molar ellipticity per mean residue, [Ѳ] in deg·cm2·dmol1, was calculated from the following equation: [Ѳ] = [Ѳ]obs × mrw × (10 × l × C)1 where [Ѳ]obs is the ellipticity measured in degrees, mrw is the mean residue molecular mass, l is the optical path length of the cell in cm, and C is the protein concentration in g/L. During the melting experiments, CD spectra were collected with increasing temperatures, from 4 °C to 80 °C, with an average rate of 1 °C/min. Thermal denaturation was investigated by recording the CD signal at 222 nm.

4.5. Light Scattering Measurements

Light scattering measurements were performed using a size exclusion chromatography (Superdex 75 increase 10/30, Cytiva, Italy) coupled with a light-scattering detector (miniDAWN TREOS, Wyatt Technology Corporation, Goleta, CA, USA) and a differential refractive index detector (Shodex RI-101, Wyatt Technology Corporation, Goleta, CA, USA). The measurements were obtained by loading ~1.0 mg of protein sample on the column equilibrated in 50 mM Tris–HCl, 150 mM NaCl, 2.5 (v/v)% glycerol (pH 7.8) as running buffer, at a flow rate of 0.5 mL/min. Data were recorded and analysed using the Astra software (version 5.3.4, Wyatt Technology Corporation, Goleta, CA, USA).

4.6. Surface Plasmon Resonance

Assessments of binding kinetics and affinity between S-plugs with either the RBD or the S1 domain of the spike protein from SARS-CoV-2 were performed at 25 °C using Biacore T200 (Cytiva, Uppsala, Sweden). The spike RBD-His (SinoBiological, Cat: 40592-V08H) was stably captured at the surface of the CM5 sensor chip by means of an anti-histidine antibody (His Capture Kit, Cytiva) that had been covalently bound to the surface, as recommended by the manufacturer. In particular, the anti-histidine antibody provided in the His Capture Kit was diluted to 50 µg/mL in the immobilisation buffer included in the kit and covalently coupled to Sensor Chip CM5 by standard amine coupling to a level of approximately 12,000 RU. Then, the spike RBD-His was injected (25 μg·mL−1) over the anti-histidine antibody surface for 1 min. No protein was injected over the reference surface. The dissociation was monitored by injecting running buffer for 120 s. Surface regeneration was performed by injecting glycine buffer at a flow rate of 30 µL/min (10 mM, pH 1.5, 1 min).
The SARS-CoV-2 spike protein (S1) (Genscript, cat. Z03501-1) was immobilised on a CM5 sensorchip (98 µg/mL in 10 mM sodium acetate pH 4.59) at a flow rate of 10 µL/min using standard amine-coupling protocols to obtain densities of approximately 11000 RU. HBS-P+ buffer was used as a running buffer. All samples and buffers were degassed prior to use. For a positive control, ACE2-His protein (SinoBiological, Cat: 10108-H08H) was used.
Protein solutions (S-plug, S-plug2, S-plug3, ACE2) in HBS-P+ buffer at various concentrations were injected at 25 °C with a flow rate of 30 μL/min for 120 s (association phase), and then the buffer alone was injected for 600 s (dissociation phase).
All mathematical manipulations and fitting operations were performed using BIAevaluation software (v2.02) provided with the Biacore T200 instrument (Cytiva) and assuming a 1:1 Langmuir binding model.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijms23105601/s1.

Author Contributions

Conceptualisation, all authors; methodology, A.R., F.S., G.B. and M.S.; investigation, all authors; data curation, all authors; writing—original draft preparation, R.B.; writing—review and editing, all authors; supervision, R.B. and A.R. All authors have read and agreed to the published version of the manuscript.

Funding

Research was supported by the project “DELETE_COVID” (DEveLopmEnT of a rapid biosEnsor against COVID19). FISR2020IP_00014. 2021-2022 and “RECOVER-COVID19” (RicErCa e sviluppO VERsus COVID19 in Campania), POR FESR CAMPANIA 2014-2020- Asse III Obiettivo Specifico 1.3 – 2020-2021. Authors also acknowledge CINECA Supercomputing (framework ISCRA@CINECA—project code HP10CSBSQC) for computational support.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Plasmids are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

COVID-19Coronavirus disease 19
PDBprotein databank
RBDreceptor-binding domain
RBMreceptor-binding motif
ACE2angiotensin-converting enzyme 2
SPRsurface plasmon resonance

References

  1. Zhou, P.; Yang, X.-L.; Wang, X.-G.; Hu, B.; Zhang, L.; Zhang, W.; Si, H.-R.; Zhu, Y.; Li, B.; Huang, C.-L.; et al. A Pneumonia Outbreak Associated with a New Coronavirus of Probable Bat Origin. Nature 2020, 579, 270–273. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Wu, J.T.; Leung, K.; Leung, G.M. Nowcasting and Forecasting the Potential Domestic and International Spread of the 2019-NCoV Outbreak Originating in Wuhan, China: A Modelling Study. Lancet 2020, 395, 689–697. [Google Scholar] [CrossRef] [Green Version]
  3. Lu, R.; Zhao, X.; Li, J.; Niu, P.; Yang, B.; Wu, H.; Wang, W.; Song, H.; Huang, B.; Zhu, N.; et al. Genomic Characterisation and Epidemiology of 2019 Novel Coronavirus: Implications for Virus Origins and Receptor Binding. Lancet 2020, 395, 565–574. [Google Scholar] [CrossRef] [Green Version]
  4. Luan, J.; Lu, Y.; Jin, X.; Zhang, L. Spike Protein Recognition of Mammalian ACE2 Predicts the Host Range and an Optimized ACE2 for SARS-CoV-2 Infection. Biochem. Biophys. Res. Commun. 2020, 526, 165–169. [Google Scholar] [CrossRef] [PubMed]
  5. Wrapp, D.; Wang, N.; Corbett, K.S.; Goldsmith, J.A.; Hsieh, C.-L.; Abiona, O.; Graham, B.S.; McLellan, J.S. Cryo-EM Structure of the 2019-NCoV Spike in the Prefusion Conformation. Science 2020, 367, 1260–1263. [Google Scholar] [CrossRef] [Green Version]
  6. Lukassen, S.; Chua, R.L.; Trefzer, T.; Kahn, N.C.; Schneider, M.A.; Muley, T.; Eils, R. SARS-CoV-2 Receptor ACE2 and TMPRSS2 Are Primarily Expressed in Bronchial Transient Secretory Cells. EMBO J. 2020, 39, e105114. [Google Scholar] [CrossRef]
  7. Hoffmann, M.; Kleine-Weber, H.; Schroeder, S.; Krüger, N.; Herrler, T.; Erichsen, S.; Schiergens, T.S.; Herrler, G.; Wu, N.-H.; Nitsche, A.; et al. SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor. Cell 2020, 181, 271–280.e8. [Google Scholar] [CrossRef]
  8. Squeglia, F.; Romano, M.; Ruggiero, A.; Maga, G.; Berisio, R. Host DDX Helicases as Possible SARS-CoV-2 Proviral Factors: A Structural Overview of Their Hijacking Through Multiple Viral Proteins. Front. Chem. 2020, 8, 602162. [Google Scholar] [CrossRef]
  9. Romano, M.; Ruggiero, A.; Squeglia, F.; Maga, G.; Berisio, R. A Structural View of SARS-CoV-2 RNA Replication Machinery: RNA Synthesis, Proofreading and Final Capping. Cells 2020, 9, 1267. [Google Scholar] [CrossRef]
  10. Walls, A.C.; Park, Y.-J.; Tortorici, M.A.; Wall, A.; McGuire, A.T.; Veesler, D. Structure, Function, and Antigenicity of the SARS-CoV-2 Spike Glycoprotein. Cell 2020, 181, 281–292.e6. [Google Scholar] [CrossRef]
  11. Piccoli, L.; Park, Y.-J.; Tortorici, M.A.; Czudnochowski, N.; Walls, A.C.; Beltramello, M.; Silacci-Fregni, C.; Pinto, D.; Rosen, L.E.; Bowen, J.E.; et al. Mapping Neutralizing and Immunodominant Sites on the SARS-CoV-2 Spike Receptor-Binding Domain by Structure-Guided High-Resolution Serology. Cell 2020, 183, 1024–1042.e21. [Google Scholar] [CrossRef]
  12. Liu, L.; Wang, P.; Nair, M.S.; Yu, J.; Rapp, M.; Wang, Q.; Luo, Y.; Chan, J.F.-W.; Sahi, V.; Figueroa, A.; et al. Potent Neutralizing Antibodies against Multiple Epitopes on SARS-CoV-2 Spike. Nature 2020, 584, 450–456. [Google Scholar] [CrossRef]
  13. Harvey, W.T.; Carabelli, A.M.; Jackson, B.; Gupta, R.K.; Thomson, E.C.; Harrison, E.M.; Ludden, C.; Reeve, R.; Rambaut, A.; COVID-19 Genomics UK (COG-UK) Consortium; et al. SARS-CoV-2 Variants, Spike Mutations and Immune Escape. Nat. Rev. Microbiol. 2021, 19, 409–424. [Google Scholar] [CrossRef]
  14. Chames, P.; van Regenmortel, M.; Weiss, E.; Baty, D. Therapeutic Antibodies: Successes, Limitations and Hopes for the Future. Br. J. Pharm. 2009, 157, 220–233. [Google Scholar] [CrossRef]
  15. Moreira, M.; Ruggiero, A.; Iaccarino, E.; Barra, G.; Sandomenico, A.; Ruvo, M.; Berisio, R. A Structure-Based Approach for the Development of a Bicyclic Peptide Acting as a Miniaturized Anti-CD55 Antibody. Int. J. Biol. Macromol. 2021, 182, 1455–1462. [Google Scholar] [CrossRef]
  16. Balasco, N.; Damaggio, G.; Esposito, L.; Villani, F.; Berisio, R.; Colonna, V.; Vitagliano, L. A Global Analysis of Conservative and Non-Conservative Mutations in SARS-CoV-2 Detected in the First Year of the COVID-19 World-Wide Diffusion. Sci. Rep. 2021, 11, 24495. [Google Scholar] [CrossRef] [PubMed]
  17. Lan, J.; Ge, J.; Yu, J.; Shan, S.; Zhou, H.; Fan, S.; Zhang, Q.; Shi, X.; Wang, Q.; Zhang, L.; et al. Structure of the SARS-CoV-2 Spike Receptor-Binding Domain Bound to the ACE2 Receptor. Nature 2020, 581, 215–220. [Google Scholar] [CrossRef] [Green Version]
  18. Yuan, S.; Chan, H.C.S.; Filipek, S.; Vogel, H. PyMOL and Inkscape Bridge the Data and the Data Visualization. Structure 2016, 24, 2041–2042. [Google Scholar] [CrossRef] [Green Version]
  19. Leman, J.K.; Weitzner, B.D.; Lewis, S.M.; Adolf-Bryfogle, J.; Alam, N.; Alford, R.F.; Aprahamian, M.; Baker, D.; Barlow, K.A.; Barth, P.; et al. Macromolecular Modeling and Design in Rosetta: Recent Methods and Frameworks. Nat. Methods 2020, 17, 665–680. [Google Scholar] [CrossRef] [PubMed]
  20. Mohamed, K.; Yazdanpanah, N.; Saghazadeh, A.; Rezaei, N. Cognitive Biases Affecting the Maintenance of COVID-19 Pandemic. Acta Biomed. 2021, 92, e2021102. [Google Scholar] [CrossRef]
  21. Zhou, T.; Tsybovsky, Y.; Gorman, J.; Rapp, M.; Cerutti, G.; Chuang, G.-Y.; Katsamba, P.S.; Sampson, J.M.; Schön, A.; Bimela, J.; et al. Cryo-EM Structures of SARS-CoV-2 Spike without and with ACE2 Reveal a PH-Dependent Switch to Mediate Endosomal Positioning of Receptor-Binding Domains. Cell Host Microbe 2020, 28, 867–879.e5. [Google Scholar] [CrossRef] [PubMed]
  22. Zhou, T.; Tsybovsky, Y.; Olia, A.S.; Gorman, J.; Rapp, M.A.; Cerutti, G.; Chuang, G.-Y.; Katsamba, P.S.; Nazzari, A.; Sampson, J.M.; et al. Cryo-EM Structures Delineate a PH-Dependent Switch That Mediates Endosomal Positioning of SARS-CoV-2 Spike Receptor-Binding Domains. BioRxiv 2020. [Google Scholar] [CrossRef]
  23. Cai, Y.; Zhang, J.; Xiao, T.; Peng, H.; Sterling, S.M.; Walsh, R.M.; Rawson, S.; Rits-Volloch, S.; Chen, B. Distinct Conformational States of SARS-CoV-2 Spike Protein. Science 2020, 369, 1586–1592. [Google Scholar] [CrossRef]
  24. Shang, J.; Ye, G.; Shi, K.; Wan, Y.; Luo, C.; Aihara, H.; Geng, Q.; Auerbach, A.; Li, F. Structural Basis of Receptor Recognition by SARS-CoV-2. Nature 2020, 581, 221–224. [Google Scholar] [CrossRef] [Green Version]
  25. Juraszek, J.; Rutten, L.; Blokland, S.; Bouchier, P.; Voorzaat, R.; Ritschel, T.; Bakkers, M.J.G.; Renault, L.L.R.; Langedijk, J.P.M. Stabilizing the Closed SARS-CoV-2 Spike Trimer. Nat. Commun. 2021, 12, 244. [Google Scholar] [CrossRef]
  26. Robbiani, D.F.; Gaebler, C.; Muecksch, F.; Lorenzi, J.C.C.; Wang, Z.; Cho, A.; Agudelo, M.; Barnes, C.O.; Gazumyan, A.; Finkin, S.; et al. Convergent Antibody Responses to SARS-CoV-2 in Convalescent Individuals. Nature 2020, 584, 437–442. [Google Scholar] [CrossRef]
  27. Jin, S.; Sun, Y.; Liang, X.; Gu, X.; Ning, J.; Xu, Y.; Chen, S.; Pan, L. Emerging New Therapeutic Antibody Derivatives for Cancer Treatment. Signal. Transduct. Target. Ther. 2022, 7, 39. [Google Scholar] [CrossRef]
  28. Romano, M.; Ruggiero, A.; Squeglia, F.; Berisio, R. An Engineered Stable Mini-Protein to Plug SARS-Cov-2 Spikes. BioRxiv 2020. [Google Scholar] [CrossRef]
  29. Cao, L.; Goreshnik, I.; Coventry, B.; Case, J.B.; Miller, L.; Kozodoy, L.; Chen, R.E.; Carter, L.; Walls, A.C.; Park, Y.-J.; et al. De Novo Design of Picomolar SARS-CoV-2 Miniprotein Inhibitors. Science 2020, 370, 426–431. [Google Scholar] [CrossRef]
  30. Emsley, P.; Cowtan, K. Coot: Model-Building Tools for Molecular Graphics. Acta Cryst. D Biol. Cryst. 2004, 60, 2126–2132. [Google Scholar] [CrossRef] [Green Version]
  31. Van der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A.E.; Berendsen, H.J.C. GROMACS: Fast, Flexible, and Free. J. Comput. Chem. 2005, 26, 1701–1718. [Google Scholar] [CrossRef] [PubMed]
  32. Kellogg, E.H.; Leaver-Fay, A.; Baker, D. Role of Conformational Sampling in Computing Mutation-Induced Changes in Protein Structure and Stability. Proteins 2011, 79, 830–838. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Squeglia, F.; Romano, M.; Ruggiero, A.; Vitagliano, L.; De Simone, A.; Berisio, R. Carbohydrate Recognition by RpfB from Mycobacterium Tuberculosis Unveiled by Crystallographic and Molecular Dynamics Analyses. Biophys. J. 2013, 104, 2530–2539. [Google Scholar] [CrossRef] [Green Version]
  34. Ruggiero, A.; Marchant, J.; Squeglia, F.; Makarov, V.; De Simone, A.; Berisio, R. Molecular Determinants of Inactivation of the Resuscitation Promoting Factor B from Mycobacterium Tuberculosis. J. Biomol. Struct. Dyn. 2013, 31, 195–205. [Google Scholar] [CrossRef]
  35. Ruggiero, A.; Squeglia, F.; Romano, M.; Vitagliano, L.; De Simone, A.; Berisio, R. Structure and Dynamics of the Multi-Domain Resuscitation Promoting Factor RpfB from Mycobacterium Tuberculosis. J. Biomol. Struct. Dyn. 2017, 35, 1322–1330. [Google Scholar] [CrossRef] [Green Version]
  36. Jorgensen, W.L.; Maxwell, D.S.; Tirado-Rives, J. Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids. J. Am. Chem. Soc. 1996, 118, 11225–11236. [Google Scholar] [CrossRef]
  37. Kaminski, G.A.; Friesner, R.A.; Tirado-Rives, J.; Jorgensen, W.L. Evaluation and Reparametrization of the OPLS-AA Force Field for Proteins via Comparison with Accurate Quantum Chemical Calculations on Peptides. J. Phys. Chem. B 2001, 105, 6474–6487. [Google Scholar] [CrossRef]
  38. Hess, B.; Bekker, H.; Berendsen, H.J.C.; Fraaije, J.G.E.M. LINCS: A Linear Constraint Solver for Molecular Simulations. J. Comput. Chem. 1997, 18, 1463–1472. [Google Scholar] [CrossRef]
  39. Petersen, H.G. Accuracy and Efficiency of the Particle Mesh Ewald Method. J. Chem. Phys. 1995, 103, 3668–3679. [Google Scholar] [CrossRef]
  40. Darden, T.; York, D.; Pedersen, L. Particle Mesh Ewald: An N log(N) Method for Ewald Sums in Large Systems. J. Chem. Phys. 1993, 98, 10089–10092. [Google Scholar] [CrossRef] [Green Version]
  41. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual Molecular Dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
Figure 1. Design of S-plug based on spike RBD-ACE2 complex structure. (A) Cartoon representation of the complex between the RBD of SARS-CoV-2 spike protein (blue/cyan) and the human ACE2 receptor (grey/green). The receptor-binding motif (RBM) is drawn in cyan. The green portion of the ACE2 domain, including helices H1, H2 and H3, is drawn in green. (B) Cartoon representations showing the superposition of S-plug (green) on ACE2 (grey). Leu91 and Leu95, involved in hydrophobic interactions with V209 and V212 in the ACE2 structure, were mutated to glutamine residues in S-plug. (C) The structure of S-plug, mutated residues are shown in stick. Figures were generated with Pymol [18].
Figure 1. Design of S-plug based on spike RBD-ACE2 complex structure. (A) Cartoon representation of the complex between the RBD of SARS-CoV-2 spike protein (blue/cyan) and the human ACE2 receptor (grey/green). The receptor-binding motif (RBM) is drawn in cyan. The green portion of the ACE2 domain, including helices H1, H2 and H3, is drawn in green. (B) Cartoon representations showing the superposition of S-plug (green) on ACE2 (grey). Leu91 and Leu95, involved in hydrophobic interactions with V209 and V212 in the ACE2 structure, were mutated to glutamine residues in S-plug. (C) The structure of S-plug, mutated residues are shown in stick. Figures were generated with Pymol [18].
Ijms 23 05601 g001
Figure 2. Molecular dynamics (MD) studies. (A,B) Secondary structure time evolution of S-plug and S-plug2, respectively, in the MD simulations complexes; best superimposition of the starting and average structures of the trajectories are reported on the left panels; RBD is shown in cyan/blue, S-plugs in yellow/green (starting/average); average structures (on Cα atoms) are computed over the entire trajectories. (C) Evolution of the number of hydrogen bonds in helical conformation in S-plug (left) and S-plug2 (right).
Figure 2. Molecular dynamics (MD) studies. (A,B) Secondary structure time evolution of S-plug and S-plug2, respectively, in the MD simulations complexes; best superimposition of the starting and average structures of the trajectories are reported on the left panels; RBD is shown in cyan/blue, S-plugs in yellow/green (starting/average); average structures (on Cα atoms) are computed over the entire trajectories. (C) Evolution of the number of hydrogen bonds in helical conformation in S-plug (left) and S-plug2 (right).
Ijms 23 05601 g002
Figure 3. Thermostability and light scattering studies. (A,B,D) Far-UV CD spectra (left) and thermal denaturation curves monitored at 222 nm (right) for S-plugs. All spectra were measured at 0.2 mg mL-1 in in 20 mM sodium phosphate buffer (pH 7.4). The panel (C) reports the analytical SEC-LS of S-plug (grey) and S-plug2 (black). Relative Rayleigh ratios (left scale) and derived molar masses (right scale, blue) are plotted versus elution time.
Figure 3. Thermostability and light scattering studies. (A,B,D) Far-UV CD spectra (left) and thermal denaturation curves monitored at 222 nm (right) for S-plugs. All spectra were measured at 0.2 mg mL-1 in in 20 mM sodium phosphate buffer (pH 7.4). The panel (C) reports the analytical SEC-LS of S-plug (grey) and S-plug2 (black). Relative Rayleigh ratios (left scale) and derived molar masses (right scale, blue) are plotted versus elution time.
Ijms 23 05601 g003
Figure 4. Cartoon representations of S-plug2 (A) and S-plug3 (B) models obtained using Rosetta docking procedures [19]. The C-terminal arm (PGGC) inserted in S-plug3 to induce a disulphide bond (B) is drawn in red ribbon representation.
Figure 4. Cartoon representations of S-plug2 (A) and S-plug3 (B) models obtained using Rosetta docking procedures [19]. The C-terminal arm (PGGC) inserted in S-plug3 to induce a disulphide bond (B) is drawn in red ribbon representation.
Ijms 23 05601 g004
Figure 5. Surface plasmon resonance (SPR) analysis. Sensorgrams to evaluate the binding affinity between S-plug3 and (A) spike S1 domain, (B) spike RBD, (C) spike RBD B.1.617.2 (delta) variant.
Figure 5. Surface plasmon resonance (SPR) analysis. Sensorgrams to evaluate the binding affinity between S-plug3 and (A) spike S1 domain, (B) spike RBD, (C) spike RBD B.1.617.2 (delta) variant.
Ijms 23 05601 g005
Figure 6. A sketch of the S-plugs’ design strategy with a summary of their thermal stabilities and binding affinities to the spike RBD and S1 domain.
Figure 6. A sketch of the S-plugs’ design strategy with a summary of their thermal stabilities and binding affinities to the spike RBD and S1 domain.
Ijms 23 05601 g006
Figure 7. Cartoon representation of the spike trimer (closed conformation) in complex with S-plug3. The positions of S-plug3 molecules (orange) on the spike protein were determined upon superposition of RBD of our minimised S-plug3–RBD complex on the RBD of the entire S protein (pdb code 6m0j). (A,B) panels report side and top views, respectively.
Figure 7. Cartoon representation of the spike trimer (closed conformation) in complex with S-plug3. The positions of S-plug3 molecules (orange) on the spike protein were determined upon superposition of RBD of our minimised S-plug3–RBD complex on the RBD of the entire S protein (pdb code 6m0j). (A,B) panels report side and top views, respectively.
Ijms 23 05601 g007
Table 1. Sequences of developed S-plugs. Most relevant residues for ACE2 interactions with spike RBD are reported in red.
Table 1. Sequences of developed S-plugs. Most relevant residues for ACE2 interactions with spike RBD are reported in red.
SequenceMelting Temperature (Tm)
S-plug STIEE Q AK TF L DK FN HE A ED LF YQ SS L ASWNYNTNITEENVQNMNNAGDKWSAFLKEQSTLAQMYPLQEIQNQTVKQQLQALQQN 37 °C
S-plug2 SAMEE Q AK FF L DK LQ HE L ED TQ YQ LL Y KALALNKSQQEELQKQYKKTYDEYVTLAKKWAKAASSIPLNSITDEKLYKAMEKTASS 48 °C
S-plug3 SAMEE Q AK FF L DK LQ HE L ED TQ YQ LL Y KALALNKSQQEELQKQYKKTYDEYVTLAKKWAKAASSIPLNSITDEKLYKAMEKTASSPPGC 54 °C
Table 2. Binding parameters of SARS-CoV-2 spike to S-plugs and ACE2 receptor. Dissociation constants were measured by SPR.
Table 2. Binding parameters of SARS-CoV-2 spike to S-plugs and ACE2 receptor. Dissociation constants were measured by SPR.
LigandAnalyteKD (nM)
SARS-CoV-2 S1S-plug271.7 ± 1.8
RBDS-plug261.7 ± 2.3
SARS-CoV-2 S1S-plug331.7 ± 1.7
RBDS-plug339.5 ± 1.1
RBD delta (B.1.617.2)S-plug329.2 ± 0.6
SARS-CoV-2 S1ACE251.9 ±1.4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Squeglia, F.; Romano, M.; Esposito, L.; Barra, G.; Campiglia, P.; Sala, M.; Scala, M.C.; Ruggiero, A.; Berisio, R. Structure-Based Development of SARS-CoV-2 Spike Interactors. Int. J. Mol. Sci. 2022, 23, 5601. https://doi.org/10.3390/ijms23105601

AMA Style

Squeglia F, Romano M, Esposito L, Barra G, Campiglia P, Sala M, Scala MC, Ruggiero A, Berisio R. Structure-Based Development of SARS-CoV-2 Spike Interactors. International Journal of Molecular Sciences. 2022; 23(10):5601. https://doi.org/10.3390/ijms23105601

Chicago/Turabian Style

Squeglia, Flavia, Maria Romano, Luciana Esposito, Giovanni Barra, Pietro Campiglia, Marina Sala, Maria Carmina Scala, Alessia Ruggiero, and Rita Berisio. 2022. "Structure-Based Development of SARS-CoV-2 Spike Interactors" International Journal of Molecular Sciences 23, no. 10: 5601. https://doi.org/10.3390/ijms23105601

APA Style

Squeglia, F., Romano, M., Esposito, L., Barra, G., Campiglia, P., Sala, M., Scala, M. C., Ruggiero, A., & Berisio, R. (2022). Structure-Based Development of SARS-CoV-2 Spike Interactors. International Journal of Molecular Sciences, 23(10), 5601. https://doi.org/10.3390/ijms23105601

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop