Jonathan M Labriola1, Shane Miersch1, Gang Chen1, Chao Chen1, Alevtina Pavlenco1, Reza Saberianfar1, Francesca Caccuri2, Alberto Zani2, Nitin Sharma3, Annie Feng3,4, Daisy W Leung4, Arnaldo Caruso2, Giuseppe Novelli5,6,7, Gaya K Amarasinghe3, Sachdev S Sidhu1. 1. Department of Molecular Genetics, The Donnelly Centre, University of Toronto, 160 College St., M5S 3E1 Toronto, Ontario, Canada. 2. Section of Microbiology, Department of Molecular and Translational Medicine, University of Brescia Medical School, 25123 Brescia, Italy. 3. Department of Pathology and Immunology, Washington University School of Medicine in St Louis, St Louis, Missouri 63110, United States. 4. Division of Infectious Diseases, John T. Milliken Department of Medicine, Washington University School of Medicine, St. Louis, Missouri 63110, United States. 5. Department of Biomedicine and Prevention, University of Rome Tor Vergata, Via Montpellier 1, 00133 Rome, Italy. 6. Italy 6 IRCCS Neuromed, Pozzilli (IS) 86077, Italy. 7. Department of Pharmacology, School of Medicine, University of Nevada, Reno, Nevada 89557, United States.
Abstract
The spread of COVID-19 has been exacerbated by the emergence of variants of concern (VoC). Many VoC contain mutations in the spike protein (S-protein) and are implicated in infection and response to therapeutics. Bivalent neutralizing antibodies (nAbs) targeting the S-protein receptor-binding domain (RBD) are promising therapeutics for COVID-19, but they are limited by low potency and vulnerability to RBD mutations in VoC. To address these issues, we used naïve phage-displayed peptide libraries to isolate and optimize 16-residue peptides that bind to the RBD or the N-terminal domain (NTD) of the S-protein. We fused these peptides to the N-terminus of a moderate-affinity nAb to generate tetravalent peptide-IgG fusions, and we showed that both classes of peptides were able to improve affinities for the S-protein trimer by >100-fold (apparent KD < 1 pM). Critically, cell-based infection assays with a panel of six SARS-CoV-2 variants demonstrated that an RBD-binding peptide was able to enhance the neutralization potency of a high-affinity nAb >100-fold. Moreover, this peptide-IgG was able to neutralize variants that were resistant to the same nAb in the bivalent IgG format, including the dominant B.1.1.529 (Omicron) variant that is resistant to most clinically approved therapeutic nAbs. To show that this approach is general, we fused the same peptide to a clinically approved nAb drug and showed that it enabled the neutralization of a resistant variant. Taken together, these results establish minimal peptide fusions as a modular means to greatly enhance affinities, potencies, and breadth of coverage of nAbs as therapeutics for SARS-CoV-2.
The spread of COVID-19 has been exacerbated by the emergence of variants of concern (VoC). Many VoC contain mutations in the spike protein (S-protein) and are implicated in infection and response to therapeutics. Bivalent neutralizing antibodies (nAbs) targeting the S-protein receptor-binding domain (RBD) are promising therapeutics for COVID-19, but they are limited by low potency and vulnerability to RBD mutations in VoC. To address these issues, we used naïve phage-displayed peptide libraries to isolate and optimize 16-residue peptides that bind to the RBD or the N-terminal domain (NTD) of the S-protein. We fused these peptides to the N-terminus of a moderate-affinity nAb to generate tetravalent peptide-IgG fusions, and we showed that both classes of peptides were able to improve affinities for the S-protein trimer by >100-fold (apparent KD < 1 pM). Critically, cell-based infection assays with a panel of six SARS-CoV-2 variants demonstrated that an RBD-binding peptide was able to enhance the neutralization potency of a high-affinity nAb >100-fold. Moreover, this peptide-IgG was able to neutralize variants that were resistant to the same nAb in the bivalent IgG format, including the dominant B.1.1.529 (Omicron) variant that is resistant to most clinically approved therapeutic nAbs. To show that this approach is general, we fused the same peptide to a clinically approved nAb drug and showed that it enabled the neutralization of a resistant variant. Taken together, these results establish minimal peptide fusions as a modular means to greatly enhance affinities, potencies, and breadth of coverage of nAbs as therapeutics for SARS-CoV-2.
SARS-CoV-2 has become endemic, necessitating development of various COVID-19 treatment
strategies beyond vaccines.[1] To a large extent, this is due to emergent
variants of concern (VoC) that have proven to be more infectious and partially resistant to
approved vaccines.[2−4] Consequently, there is an
urgent need for alternative therapeutic strategies to complement vaccine campaigns.SARS-CoV-2 uses its surface spike glycoprotein (S-protein) to interact with host surface
receptors and enter host cells. The virus surface displays 25–100 copies of the
S-protein homotrimer. Each S-protein contains two subunits: the N-terminal subunit (S1) that
mediates host cell recognition and the C-terminal subunit (S2) that mediates membrane
fusion.[5] The S1 subunit itself contains an N-terminal domain (NTD)
followed by a receptor-binding domain (RBD)[6] that interacts with the host
cell–surface protein angiotensin-converting enzyme 2 (ACE2) to initiate
infection.[7]Most natural neutralizing antibodies (nAbs) target the S1 subunit. Many of these bind to
the RBD and compete with ACE2,[8−14] but
a distinct subset has been shown to target a neutralizing epitope on the NTD.[9] Several natural nAbs have been produced recombinantly and engineered further
to develop therapeutics to treat SARS-CoV-2 infection in patients.[15]
However, the current approved antibody drugs must be administered at very high doses and
have proven to be ineffective against many VoC that have arisen since the original COVID-19
outbreak.[2,4] Indeed,
most VoC that resist the current therapeutic nAbs contain mutations within the RBD that
disrupt binding to the nAbs[2] but not to ACE2.[16]
Notably, the B.1.1.529 (Omicron) variant, which emerged in South Africa and contains
numerous RBD mutations, is particularly resistant to vaccination and clinical nAbs.To address these limitations on the potency and breadth of coverage of current bivalent IgG
therapies, several groups have developed higher valence protein-based inhibitors.[17] These include small modular Ab domains or non-Ab scaffolds that can be
assembled as multimers with enhanced potency due to simultaneous engagement with all three
RBDs on an S-protein trimer.[18,19] Alternatively, we have shown that the fusion of additional Fab arms to
either the N- or C-terminus of an IgG heavy chain results in tetravalent IgG-like molecules
with enhanced potency and effectiveness against VoC that resist bivalent IgGs.[20]Notably, small peptides that target the S-protein with submicromolar affinities have been
developed and have shown promise as diagnostic tools.[21] Here, we explored
whether synthetic peptides that bind to the S-protein could be used to augment the
neutralization potency of IgGs in the form of tetravalent peptide–IgG fusions that
combine the binding sites of potent nAbs with the small, modular binding sites of peptides.
We used naïve phage-displayed peptide libraries to derive synthetic peptides that
bind to neutralizing epitopes on the RBD or the NTD. We showed that these small peptides
could be fused to a moderate-affinity nAb to develop peptide–IgG fusions with
affinities enhanced by over 2 orders of magnitude. Most importantly, one such peptide fusion
was able to greatly enhance the neutralization potency against SARS-CoV-2 and
VoC—including the Omicron variant—for a high-affinity nAb we had engineered
earlier[20] and also for a clinically approved therapeutic nAb developed
by others.[10] Thus, these synthetic peptides hold great promise to enhance
the potency and breadth of coverage of therapeutic nAbs against SARS-CoV-2 and its VoC.
Results
Isolation and Characterization of S-Protein-Binding Peptides
To isolate novel peptides that bind to the S-protein of SARS-CoV-2, we used
phage-displayed libraries of 16-residue peptides. We pooled together phage representing a
panel of 10 peptide libraries in which diversified positions were encoded by an equimolar
mixture of 19 codons representing all genetically encoded amino acids except cysteine. An
“unconstrained” library (X16) contained 16 diversified positions
with no fixed positions,[22] whereas the remaining libraries contained 14
diversified positions and two fixed cysteine residues separated by 4–12 diversified
positions. These “constrained” libraries were designed to display peptides
containing disulfide-bonded loops, which promote tertiary structures that can enhance
binding to proteins.[23]Phages representing the pooled libraries were cycled through five rounds of binding
selections with immobilized S-protein ectodomain (ECD) trimer or RBD (unless otherwise
noted, virus proteins were from the B.1 variant), and several hundred clones were analyzed
for binding to the ECD. Clones that exhibited strong binding signals in phage ELISAs with
the ECD and negligible signals with bovine serum albumin (BSA) and neutravidin (NAV) were
subjected to DNA sequence analysis. This process yielded 160 and 128 unique peptide
sequences from the ECD and RBD selections, respectively. Alignment of the sequences
revealed that most (85%) of the ECD-selected peptides were derived from two libraries: 63%
from the X16 library and 22% from the C-X4-C library (Figure A). In contrast, most (69%) of the
RBD-selected peptides were from two different libraries: 42% from the C-X10-C
library and 27% from the C-X10-C library (Figure B). Inspection of the aligned sequences revealed significant
homology within each peptide family, which enabled us to derive consensus motifs (Figure ). Notably, the peptides derived from the
C-X10-C and C-X12-C libraries exhibited similar consensus motifs,
suggesting that they bind to the S-protein in a similar manner.
Figure 1
Sequence alignments of phage-derived S-protein-binding peptides. (A) Peptides
selected for binding to the S-protein ECD originated from either the X16
library (left) or the C-X4-C library
(right). (B) Peptides selected for binding to the S-protein RBD
originated from either the C-X10-C library (left) or the
C-X12-C library (right). The positions are numbered at
the top, followed by the consensus sequence for positions that exhibited high sequence
conservation (>25% for X16, >50% for others), and the unique selected sequences
are aligned below the numbering and consensus. Sequences that match the consensus are
shaded gray, and fixed cysteines are shaded black. Peptides that were characterized in
detail (N1, N2, R1, R2) are labeled and shown at the top of each alignment.
Sequence alignments of phage-derived S-protein-binding peptides. (A) Peptides
selected for binding to the S-protein ECD originated from either the X16
library (left) or the C-X4-C library
(right). (B) Peptides selected for binding to the S-protein RBD
originated from either the C-X10-C library (left) or the
C-X12-C library (right). The positions are numbered at
the top, followed by the consensus sequence for positions that exhibited high sequence
conservation (>25% for X16, >50% for others), and the unique selected sequences
are aligned below the numbering and consensus. Sequences that match the consensus are
shaded gray, and fixed cysteines are shaded black. Peptides that were characterized in
detail (N1, N2, R1, R2) are labeled and shown at the top of each alignment.For more detailed characterization, we chose four peptides that closely matched the
consensus motifs for their respective families. These included a peptide from the
X16 library (N1) and a peptide from the C-X4-C library (N2), which
were selected for binding to the ECD, and two peptides (R1 and R2) from the
C-X12-C library, which were selected for binding to the RBD (Figure ). As expected, phage ELISAs showed that all four
peptide phages bound to immobilized ECD, but only peptides R1 and R2 bound to immobilized
RBD (Figure A). To further define where each
peptide bound, we assessed whether the binding of peptide phage to the ECD could be
blocked by nAbs that recognized known epitopes on either the N-terminal domain (NTD) (IgGs
5–24 and 4–8)[9] or the RBD (IgGs 15033-7 and
REGN10933).[10,20]
Binding of peptide phages N1 and N2, but not R1 and R2, was blocked by the two antibodies
that bound to a neutralizing epitope on the NTD but not by those that bound to the
ACE2-binding site of the RBD. Conversely, the binding of peptide phages R1 and R2, but not
N1 and N2, was blocked by the two antibodies that bound to the RBD but not by those that
bound to the NTD (Figure B). Taken together,
these results showed that peptides N1 and N2 likely bind to the NTD, whereas peptides R1
and R2 bind to the RBD, and all peptides likely bind to sites that overlap with epitopes
of nAbs.
Figure 2
Characterization of epitopes for phage-displayed S-protein-binding peptides. (A)
Phage ELISAs for peptide-phage binding to immobilized S-protein ECD (blue bars) or RBD
(red bars). Binding was quantified by the optical density at 450 nm
(y-axis) for 0.5 nM phage, displaying the indicated peptide
(x-axis). (B) Phage ELISAs for peptide-phage
(x-axis) binding to the immobilized S-protein ECD in the presence of
saturating concentrations of antibodies that bind to the NTD (5–24, dark blue;
4–8, light blue) or RBD (15033-7, dark red; RGN10933, light red). The binding
signal (y-axis) was normalized to binding in the absence of an
antibody.
Characterization of epitopes for phage-displayed S-protein-binding peptides. (A)
Phage ELISAs for peptide-phage binding to immobilized S-protein ECD (blue bars) or RBD
(red bars). Binding was quantified by the optical density at 450 nm
(y-axis) for 0.5 nM phage, displaying the indicated peptide
(x-axis). (B) Phage ELISAs for peptide-phage
(x-axis) binding to the immobilized S-protein ECD in the presence of
saturating concentrations of antibodies that bind to the NTD (5–24, dark blue;
4–8, light blue) or RBD (15033-7, dark red; RGN10933, light red). The binding
signal (y-axis) was normalized to binding in the absence of an
antibody.
Optimization of S-Protein-Binding Peptides
We took advantage of the large families of related sequences to design biased
peptide-phage libraries to further optimize peptides N1 and R1 (Figure
). For each peptide, we used the sequence logo derived from
its family to identify highly conserved positions that were fixed, and variable positions,
which were randomized using degenerate codons to encode a mix of parent and other amino
acids that were prevalent in the sequence logo. Phage pools representing each biased
library were cycled through five rounds of selection for binding to the S-protein ECD
trimer, positive clones were identified by phage ELISA, and DNA sequencing revealed 86 and
25 unique sequences from the N1 and R1 libraries, respectively (Figure
). We used the unique sequences to derive a sequence logo for
each selection pool (Figure ). We chose unique
peptide sequences that closely matched the logo for further characterization, reasoning
that these likely represented high-affinity binders.
Figure 3
Optimization of S-protein-binding peptides. Biased peptide-phage libraries were
designed to optimize (A) peptide N1 or (B) peptide R1. For each peptide, the sequence
logo at the top was derived from the family of peptides selected from naïve
peptide-phage libraries (Figure ). Below the
logo, the parent peptide sequence is shown with positions that were fixed or
diversified shaded black or gray, respectively. For each diversified position,
additional amino acids present in the displayed library are shown below the parent
amino acid. Output clones obtained from selection of the biased library against the
S-protein ECD were used to construct the sequence logo at the bottom of each panel
(Figure ).
Figure 4
Sequence alignments of S-protein-binding peptides derived from biased phage-displayed
libraries. Peptides selected for binding to the S-protein ECD originated from the
biased library based on the sequence of (A) peptide N1 or (B) peptide R1 (see Figure ). Only a subset of peptides from the N1
selection is shown. The positions are numbered at the top, followed by the sequence of
the parent peptide. Amino acids that match the parental sequence at each position are
shaded black. Amino acids that do not match the parental but are most prevalent at
each position are shaded gray. Asterisks indicate positions at which sequences
diverged from the parental sequence, suggesting that sequence differences here may
enhance the affinity of the variants relative to the parent. The sequences of peptide
variants chosen for chemical synthesis and further characterization are labeled to the
left of each alignment.
Optimization of S-protein-binding peptides. Biased peptide-phage libraries were
designed to optimize (A) peptide N1 or (B) peptide R1. For each peptide, the sequence
logo at the top was derived from the family of peptides selected from naïve
peptide-phage libraries (Figure ). Below the
logo, the parent peptide sequence is shown with positions that were fixed or
diversified shaded black or gray, respectively. For each diversified position,
additional amino acids present in the displayed library are shown below the parent
amino acid. Output clones obtained from selection of the biased library against the
S-protein ECD were used to construct the sequence logo at the bottom of each panel
(Figure ).Sequence alignments of S-protein-binding peptides derived from biased phage-displayed
libraries. Peptides selected for binding to the S-protein ECD originated from the
biased library based on the sequence of (A) peptide N1 or (B) peptide R1 (see Figure ). Only a subset of peptides from the N1
selection is shown. The positions are numbered at the top, followed by the sequence of
the parent peptide. Amino acids that match the parental sequence at each position are
shaded black. Amino acids that do not match the parental but are most prevalent at
each position are shaded gray. Asterisks indicate positions at which sequences
diverged from the parental sequence, suggesting that sequence differences here may
enhance the affinity of the variants relative to the parent. The sequences of peptide
variants chosen for chemical synthesis and further characterization are labeled to the
left of each alignment.Based on this process, we identified four variants of peptide N1 (N1a, b, c, d) and two
variants of peptide R1 (R1a, b) for chemical synthesis (Figure ). These were synthesized as 22-residue peptides consisting of the
following: a Gly residue, followed by the selected 16-residue sequence, followed by a
5-residue hydrophilic linker (Gly–Gly–Lys–Gly–Lys), followed
by biotin. We used biolayer interferometry (BLI) to determine binding kinetics and
affinities. Biotinylated peptides were immobilized on sensor chips that were incubated
with solution-phase S-protein ECD (Figure S1). Peptide N1 did not exhibit detectable binding, presumably due to
low affinity, but its variants bound with moderate apparent affinities (apparent
KD = 26–68 nM). Peptide R1 and its variants exhibited
high apparent affinities in the single-digit nanomolar range for the S-proteins from both
the B.1 variant (apparent KD = 0.8–5.6 nM) and the
B.1.1.529 variant (apparent KD = 3–7 nM) (Table ). We note that polyvalent avidity effects
likely enhance these apparent affinities, and the intrinsic monovalent affinities are
likely lower. Nonetheless, this format does mimic the polyvalent nature of the
virus-neutralizing molecules we ultimately intended to produce, as described below.
Table 1
Kinetics of Synthetic Peptide Binding to S-Protein ECDa,b,c
Peptides were synthesized with a C-terminal biotin group. The 16-residue sequence
selected for binding to the S-protein ECD or RBD is shaded gray and numbered as in
Figures and 4. For
optimized variants of peptide N1 (N1a, b, c, d) and peptide R1 (R1a, b), dashes
indicate the identity with the parent amino acid. A biotin tag was added to the
C-terminus of each peptide for capture on streptavidin-coated sensors.
Kinetic binding constants (k and
k) and the equilibrium
dissociation constant (KD) were determined by the global
fitting of response traces obtained when dipping peptide-coated sensors into
solutions of S-protein. For R1, values in parentheses were determined from
measurements against B.1.1.529 S-protein.
N.D., not determined. Affinity for the CoV-2-spike was too low to be measured by
BLI.
Peptides were synthesized with a C-terminal biotin group. The 16-residue sequence
selected for binding to the S-protein ECD or RBD is shaded gray and numbered as in
Figures and 4. For
optimized variants of peptide N1 (N1a, b, c, d) and peptide R1 (R1a, b), dashes
indicate the identity with the parent amino acid. A biotin tag was added to the
C-terminus of each peptide for capture on streptavidin-coated sensors.Kinetic binding constants (k and
k) and the equilibrium
dissociation constant (KD) were determined by the global
fitting of response traces obtained when dipping peptide-coated sensors into
solutions of S-protein. For R1, values in parentheses were determined from
measurements against B.1.1.529 S-protein.N.D., not determined. Affinity for the CoV-2-spike was too low to be measured by
BLI.
Production and Characterization of Peptide–IgG Fusions
We next explored whether peptide fusions could enhance the affinity of a neutralizing IgG
targeting the ACE2-binding site of the RBD. For this purpose, we used the
moderate-affinity IgG 15033 that we had selected from a naïve phage-displayed
synthetic antibody library[20] to accurately discern affinity
differences. Peptide N1 or R1 was fused to the N-terminus of either the light chain (LC)
or heavy chain (HC) of IgG 15033 with an intervening 20-residue Gly/Ser linker (Figure S2). The resulting peptide–IgG fusion proteins and IgG 15033
were purified by transient transfection in mammalian Expi293F cells.[20]
All proteins could be purified to near-homogeneity by affinity chromatography with
protein-A resin, as evidenced by SDS-PAGE (Figure A). As expected, under nonreducing conditions, the single bands for the intact
peptide–IgG molecules migrated slightly slower than the band for IgG 15033. Under
reducing conditions, either the HC band or the LC band migrated slower for
peptide–IgG fusions compared with that for IgG 15033, as expected for HC or LC
peptide fusions, respectively. Like IgG 15033 and the highly specific IgG trastuzumab, the
four peptide–IgG fusions did not bind to seven immobilized, heterologous proteins
that are known to exhibit high nonspecific binding to some IgGs (Figure
B). Lack of binding in this assay is a predictor of good
pharmacokinetics in vivo.(24)
Figure 5
Characterization of peptide–IgG fusion proteins. Peptides were fused to the
N-terminus of the HC or LC of IgG 15033, and the resulting peptide–IgG fusions
were named as follows: N1 fused to HC, 33HN1; N1 fused to LC,
33LN1; R1 fused to HC, 33HR1; and R1 fused to LC,
33LR1. (A) SDS-PAGE analysis of peptide–IgG fusion proteins under
nonreducing (top) or reducing conditions (bottom). (B) Assessment of nonspecific
binding of peptide–IgG fusion proteins to immobilized antigens or a goat
antihuman Fc Ab (positive control). (C) BLI sensor traces for 20 nM IgG 15033 or
peptide–IgG fusions binding to the immobilized S-protein ECD. Association and
dissociation kinetics were monitored for 600 seconds each. The derived binding
constants are shown in Table .
Characterization of peptide–IgG fusion proteins. Peptides were fused to the
N-terminus of the HC or LC of IgG 15033, and the resulting peptide–IgG fusions
were named as follows: N1 fused to HC, 33HN1; N1 fused to LC,
33LN1; R1 fused to HC, 33HR1; and R1 fused to LC,
33LR1. (A) SDS-PAGE analysis of peptide–IgG fusion proteins under
nonreducing (top) or reducing conditions (bottom). (B) Assessment of nonspecific
binding of peptide–IgG fusion proteins to immobilized antigens or a goat
antihuman Fc Ab (positive control). (C) BLI sensor traces for 20 nM IgG 15033 or
peptide–IgG fusions binding to the immobilized S-protein ECD. Association and
dissociation kinetics were monitored for 600 seconds each. The derived binding
constants are shown in Table .
Table 2
Kinetics of Peptide–IgG Fusion Protein Binding to S-Protein ECD
IgGa
kon
(104M–1s–1)
koff
(10–6S–1)
KD(pM)b
15033
170 ± 5
500 ± 20
300 ± 20
33HN1
160 ± 5
38 ± 20
24 ± 10
33LN1
110 ± 3
<0.1
<1
33HR1
120 ± 4
<0.1
<1
33LR1
35 ± 1
2.2 ± 1
6.1 ± 3
Peptide–IgG fusion proteins consisted of phage-derived peptides fused to the
N-terminus of the light chain or heavy chain of IgG 15033. The following
nomenclature was used: L, light-chain fusion; H, heavy-chain fusion; N1, NTD-binding
peptide 1; and R1, RBD-binding peptide 1 (see Table for sequences).
Kinetic binding constants (kon and
koff) and the equilibrium dissociation constant
(KD) were determined by the global fitting of response
traces obtained when dipping S-protein-coated sensors into solutions of indicated
IgG.
Kinetic binding analysis by BLI showed that the tetravalent peptide–IgG fusions
exhibited greatly reduced off-rates compared with IgG 15033 (Figure S3), and consequently, affinities for the S-protein ECD were greatly
improved (Figure C and Table ). These results
were as expected since the addition of a peptide ligand to an RBD-binding IgG should
reduce the overall free energy of binding and thus increase the binding efficiency. In
particular, the peptide–IgG with peptide R1 fused to its HC (33HR1) and
the peptide–IgG with peptide N1 fused to its LC (33LN1) exhibited
extremely slow off-rates that were beyond the sensitivity of the instrument, and
consequently, the apparent dissociation constants were estimated to be subpicomolar
(apparent KD <1 pM), which was >100-fold improved
relative to IgG 15033 (apparent KD = 100 pM). We note that the
response signals of the N1 peptide–IgG fusions were much greater than those of the
R1 peptide–IgG fusions. This is likely due to peptide N1 being unconstrained with a
higher koff for the S-protein compared to R1 (Table ). This may allow the N1 peptide–IgG fusions to
adopt more extended conformations upon interaction with the S-protein, which would
manifest as greater signals in the BLI experiment.Peptide–IgG fusion proteins consisted of phage-derived peptides fused to the
N-terminus of the light chain or heavy chain of IgG 15033. The following
nomenclature was used: L, light-chain fusion; H, heavy-chain fusion; N1, NTD-binding
peptide 1; and R1, RBD-binding peptide 1 (see Table for sequences).Kinetic binding constants (kon and
koff) and the equilibrium dissociation constant
(KD) were determined by the global fitting of response
traces obtained when dipping S-protein-coated sensors into solutions of indicated
IgG.
Inhibition of Virus Infection in Cell-Based Assays
To assess the efficacy of the peptides and peptide–IgG fusion proteins for the
neutralization of SARS-CoV-2 VoC, we first used mammalian cell infection assays with
pseudoviruses consisting of HIV-gag-based, lentivirus-like particles pseudotyped with the
SARS-CoV-2 S-protein.[25] Neither peptide, nor their optimized variants,
exhibited any effect on pseudovirus infection at concentrations up to 10 μM (Figure A). However, peptide–IgG fusion
proteins with the N1 peptide or the R1 peptide fused to the N-terminus of the HC
(33HN1 or 33HR1, respectively) or LC of IgG 15033
(33LN1 or 33LR1, respectively) exhibited enhanced neutralization
potency compared with the IgG alone (Figure B).
Figure 6
Neutralization of SARS-CoV-2 pseudovirus. Neutralization assays were conducted with
HIV-gag-based, lentivirus-like particles pseudotyped with the S-protein of the
original B.1 strain from Wuhan. (A) Assays were conducted with synthetic peptides N1
(blue squares), R1 (red triangles), or the highest affinity optimized peptide for each
(clear square and triangle, respectively) and IgG 15033-7 (black circle) as a positive
control. (B) Assays were conducted with IgG 15033 (gray circle) or with
peptide–IgG 33HN1 (blue square), 33HR1 (red triangle,
upper), or 33LN1 (blue square), 33LR1 (red triangle, lower).
Neutralization of SARS-CoV-2 pseudovirus. Neutralization assays were conducted with
HIV-gag-based, lentivirus-like particles pseudotyped with the S-protein of the
original B.1 strain from Wuhan. (A) Assays were conducted with synthetic peptides N1
(blue squares), R1 (red triangles), or the highest affinity optimized peptide for each
(clear square and triangle, respectively) and IgG 15033-7 (black circle) as a positive
control. (B) Assays were conducted with IgG 15033 (gray circle) or with
peptide–IgG 33HN1 (blue square), 33HR1 (red triangle,
upper), or 33LN1 (blue square), 33LR1 (red triangle, lower).Finally, we explored the critical question of whether peptide fusions could enhance
potency of nAbs against authentic SARS-CoV-2 in ACE2 expressing Vero E6 cell-based assays.
For this purpose, we fused peptide R1 to the N-terminus of the HC of IgG 15033-7, a more
potent variant of IgG 15033 with an optimized LC.[20] Neutralization
efficiency of the resulting peptide–IgG fusion protein (33-7HR1) was
compared to IgG 15033-7 against a panel of seven SARS-CoV-2 variants: Wuhan (B.1), Italy
(B.1.1), the United Kingdom (B.1.1.7), South Africa (B.1.351 and B.1.529), Nigeria
(B.1.525), and Brazil (B.1.128). In every case, the potency of the peptide–IgG
fusion greatly exceeded that of the IgG (Figure A). Indeed, peptide–IgG 33-7HR1 neutralized three VoC in the
single-digit ng/mL range (IC50 = 6.5–7.6 ng/mL), and three other VOC in
the double-digit ng/mL range (IC50 = 12–49 ng/mL). Notably,
33-7HR1 was even able to neutralize the B.1.1.529 variant, which is resistant
to most clinically approved IgG drugs, albeit with reduced potency (IC50 = 260
ng/mL). In contrast, IgG 15033-7 was less effective than IgG 33-7HR1 against
three VOC (IC50 = 64–310 ng/mL), much less effective against three other
VoC (IC50 > 900 ng/mL), and completely ineffective against B.1.1.529.
Figure 7
Neutralization of SARS-CoV-2 variants. (A) Authentic virus neutralization assays to
test the efficacy of IgG 15033-7 and peptide–IgG 33-7HR1 (IgG
15033-7 with peptide R1 fused to the N-terminus of its HC) against isolated SARS-CoV-2
variants. The virus was pretreated with serial dilutions of IgG, followed by infection
of ACE2-expressing Vero E6 cells, and measured relative to the untreated control. At
the top, the graph shows the neutralization curves for seven SARS-CoV-2 variants
treated with IgG 15033-7 (black dashed lines) or peptide–IgG 33-7HR1
(red solid lines). At the bottom, the table shows the symbols used for the curves and
the country of origin for each virus variant. IC50 values calculated from
the curves are shown in ng/mL or nM units. (B) Focus reduction neutralization assays
to test the efficacy of IgG REGN10933 and peptide–IgG REGN10933HR1
(REGN10933 with peptide R1 fused to the N-terminus of its HC) against isogenic
SARS-CoV-2 variants D641G and B.1.351. The virus was pretreated with serial dilutions
of IgG followed by infection of ACE2-expressing Vero E6 cells, which was measured
relative to the untreated virus. Each value represents the mean of duplicate
measurements.
Neutralization of SARS-CoV-2 variants. (A) Authentic virus neutralization assays to
test the efficacy of IgG 15033-7 and peptide–IgG 33-7HR1 (IgG
15033-7 with peptide R1 fused to the N-terminus of its HC) against isolated SARS-CoV-2
variants. The virus was pretreated with serial dilutions of IgG, followed by infection
of ACE2-expressing Vero E6 cells, and measured relative to the untreated control. At
the top, the graph shows the neutralization curves for seven SARS-CoV-2 variants
treated with IgG 15033-7 (black dashed lines) or peptide–IgG 33-7HR1
(red solid lines). At the bottom, the table shows the symbols used for the curves and
the country of origin for each virus variant. IC50 values calculated from
the curves are shown in ng/mL or nM units. (B) Focus reduction neutralization assays
to test the efficacy of IgG REGN10933 and peptide–IgG REGN10933HR1
(REGN10933 with peptide R1 fused to the N-terminus of its HC) against isogenic
SARS-CoV-2 variants D641G and B.1.351. The virus was pretreated with serial dilutions
of IgG followed by infection of ACE2-expressing Vero E6 cells, which was measured
relative to the untreated virus. Each value represents the mean of duplicate
measurements.We also explored whether the modularity of the peptides could be exploited to enhance the
potency of a clinically approved therapeutic nAb (REGN10933; Figure B). We constructed a variant of REGN10933 by fusing peptide R1 to
the N-terminus of the HC. The resulting peptide–IgG fusion
(REGN10933HR1) proved to be more potent than REGN10933 in neutralization assays
against the SARS-CoV-2 variant D614G (IC50 < 10 or ∼10 ng/mL,
respectively). Most strikingly, peptide–IgG REGN10933HR1 was also
extremely potent against the South African variant B.1.351 (IC50 < 10
ng/mL), against which REGN10933 was completely ineffective (IC50 > 1000
ng/mL), consistent with previous reports.[2] Taken together, these
results showed that peptide fusions greatly enhanced the potency and breadth of coverage
of neutralizing IgGs against SARS-CoV-2 and its VoC.
Discussion
The trimeric structure of the SARS-CoV-2 S-protein can be exploited to engineer Ab-based
inhibitors with enhanced neutralization by engaging three neutralizing epitopes on a single
trimer. Structural studies from us and others have shown that potent neutralizing IgGs
engage two RBDs on an S-protein trimer,[26] and we have shown that the
addition of Fab arms to either end of an IgG can further enhance neutralization by enabling
engagement to the third RBD.[20] Here, rather than using large Fabs as the
additional binding domain, we used small peptides to greatly enhance the affinities and
potencies of tetravalent peptide–IgG fusions compared with our bivalent IgG 15033-7.
Most importantly, peptide–IgG 33-7HR1 acted as a potent inhibitor of virus
variants that resisted IgG 15033-7 (Figure A).
Indeed, we showed further that a peptide–IgG version of the approved drug RGN10933
was able to potently neutralize a VoC against which REGN10933 was completely ineffective
(Figure B).To gain insight into the structural basis for how a peptide within a peptide–IgG
fusion could enhance affinity, we examined our published model of two 15033-7 Fabs bound to
the S-protein trimer, reasoning that this likely provides an accurate view of how the two
Fab arms of a bivalent IgG would bind (Figure A).
In this model, the RBDs bound to Fabs are in an “up” conformation, whereas the
unbound RBD is in a “down” conformation. The C-termini of the two Fab HCs are
separated by 45 Å and are well positioned to be linked to an Fc in an IgG molecule,
whereas their N-termini are close to the unbound RBD. In particular, the N-terminus of one
HC is 45 Å from the center of the exposed region of the neutralizing epitope for Fab
15033-7 on the unbound RBD. Although peptide R1 alone does not appear to neutralize virus
infection (Figure A), it likely binds to a site
that overlaps with this epitope, given that the peptide and IgG compete for binding to the
S-protein (Figure B). Consequently, in
peptide–IgG 33-7HR1 with two Fabs bound to an S-protein timer, it is
reasonable to assume that one of the R1 peptides could bind simultaneously to a region close
to this third epitope, and the estimated length (∼70 Å) of the extended
20-residue linker that connects the peptide to the IgG is consistent with this model.
Figure 8
Structural models of multivalent ligands binding to the S-protein trimer. Structural
models are shown for (A) two 15033-7 Fabs or (B) two RGN10933 Fabs bound to the S1
subunit of the S-protein of SARS-CoV-2. The S1-protein is shown as a surface colored in
light gray, except for the following regions. The bound RBDs in the “up”
position are colored dark gray. The unbound RBD in the “down” position is
colored dark or light red for residues that are within or outside the epitope for Fab
15033-7, respectively. The Fabs are shown as ribbons, with the HC and LC colored forest
green or teal, respectively. The C-termini and N-termini of the Fab HCs are shown as
yellow or magenta spheres, respectively. Distances are shown as dashed black lines and
are labeled accordingly.
Structural models of multivalent ligands binding to the S-protein trimer. Structural
models are shown for (A) two 15033-7 Fabs or (B) two RGN10933 Fabs bound to the S1
subunit of the S-protein of SARS-CoV-2. The S1-protein is shown as a surface colored in
light gray, except for the following regions. The bound RBDs in the “up”
position are colored dark gray. The unbound RBD in the “down” position is
colored dark or light red for residues that are within or outside the epitope for Fab
15033-7, respectively. The Fabs are shown as ribbons, with the HC and LC colored forest
green or teal, respectively. The C-termini and N-termini of the Fab HCs are shown as
yellow or magenta spheres, respectively. Distances are shown as dashed black lines and
are labeled accordingly.We next examined an analogous structural model of two RGN10933 Fabs bound to an S-protein
trimer (Figure B). The epitopes for REGN10933 and
15033-7 Fabs share a significant overlap, but the REGN10933 Fabs are oriented such that the
N-termini of the two HCs are 34 or 66 Å from the 15033-7 epitope on the unbound RBD.
Consequently, the ∼70 Å linker that connects each peptide to the IgG would be
sufficient to enable either of the two fused peptides within
peptide–IgG10933HR1 to bind to a site close to this epitope.We also developed peptides that bound to the NTD, likely at a site that overlaps with
epitopes for other neutralizing nAbs (Figure B).
We showed that peptide N1 greatly enhanced the affinity of peptide–IgG
33LN1, and both peptide–IgG 33LN1 and 33HR1
exhibited enhanced potency in a pseudovirus neutralization assay (Figure
B). Notably, peptides R1 and N1 appear to work best as fusions
to the HC or LC, respectively, and given that they recognize completely distinct epitopes,
it will be interesting to explore whether dual fusions of the two peptides could enhance the
potencies of peptide–IgG fusions even further, and these studies are ongoing.The ultimate value of peptide–IgG fusions targeting SARS-CoV-2 and its variant
resides in the potential for the development of superior therapeutics for the treatment of
COVID-19 and other viral diseases. We are currently assessing and optimizing the biophysical
properties of the peptide–IgG fusions with the aim of developing drug-grade
biologics. However, there is already a strong precedent for peptide–IgG fusions as
clinical drugs in the form of biologics that have reached phase 2 clinical trials as cancer
therapeutics.[27] These anticancer biologics are very similar in format
to the antiviral peptide–IgG fusions we report here.[28] Critically,
anticancer peptide–IgG fusions have exhibited good pharmacokinetic profiles[29] and tissue penetration,[30] and they have been well
tolerated in the clinic.[31] Taken together, our results establish
peptide–IgG fusions as a powerful means for greatly enhancing the potency and
coverage of next-generation biologics for the treatment of diseases caused by SARS-CoV-2 and
its variants.
Materials and Methods
Protein Production
The SARS-CoV-2 S-protein ECD and RBD were produced and purified as described[20] and were a kind gift from Dr. James Rini. Purified proteins were
site-specifically biotinylated in a reaction with 200 μM biotin, 500 μM ATP,
500 μM MgCl2, 30 μg/mL BirA, 0.1% (v/v) protease inhibitor
cocktail, and not more than 100 μM of the protein–AviTag substrate. The
reactions were incubated at 30 °C for 2 h, and biotinylated proteins were purified by
size-exclusion chromatography.
Antibody Production
IgG and peptide–IgG fusion proteins were produced as described.[20] S-protein-binding peptides were fused to the N-terminus of the heavy or light chain
through a 20-residue linker (sequence: GGGGSGGGGSGGGGSGGGGS) using standard molecular
biology techniques.
Phage Display Selections
Naïve libraries were constructed as described.[22] For libraries
for peptide optimization, oligonucleotides were synthesized using degenerate codons
encoding for the amino acids at each position indicated in Figure . Oligonucleotide-directed mutagenesis was performed to introduce
randomized sequences fused to the gene-8 major coat protein of M-13 bacteriophage, as
described.[32] Each phage-displayed peptide library was selected for
binding to the immobilized S-protein ECD or RBD, as described,[32] but
with the following modifications. Biotinylated S-protein ECD or RBD was captured on wells
coated with NAV, followed by incubation with the phage-displayed peptide library. After
five rounds of binding selections, individual clones were picked and DNA was sequenced.
Clones that showed a significant binding signal to S-protein ECD and/or RBD, but not to
BSA or NAV, were selected for further analysis.
Peptide Synthesis
Peptides were synthesized by AbClonal. Peptide sequences selected by phage display (gray
shading in Table ) were synthesized with a Gly
residue added to the N-terminus and a Gly–Gly–Lys–Gly–Lys
linker added to the C-terminus. Biotin was conjugated via the C-terminal lysine
residue.
ELISAs
Phage ELISAs were performed, as described,[20] but with the following
modifications: 384-well Maxisorp plates (Sigma-Aldrich) were coated with NAV, or left
uncoated as a negative control, and blocked with PBS, 0.5% bovine serum albumin (BSA).
Biotinylated target protein was captured by incubation in NAV-coated and BSA-blocked
wells, or with buffer solution alone as a negative control, at room temperature. For
competition ELISAs, blocking IgG was incubated with coated and blocked wells for 1 h at
room temperature. Wells were incubated with peptide phage in PBS, 0.5% BSA for 1 h. Plates
were washed, incubated with anti-M-13-HRP antibody (Sino Biological, catalog number
11973-MM05T-H), and developed with TMB substrate (Mandel, catalog number KP-50-76-03).For assessing the specificity of antibodies (Figure B), the following proteins were used: KLH (Sigma-Aldrich, H8283), cardiolipin
(Sigma-Aldrich, C0563), BSA (Medstore, ALB001), insulin (Sigma-Aldrich, I0516), LPS
(Invivogen, tlrl-eblps), ssDNA (Sigma-Aldrich, D8899), dsDNA (Sigma-Aldrich, D4522), and
goat antihuman Fc (Jackson Labs—109-005-098). Proteins (5 μg/mL) were coated
for 2 h at room temperature on 96-well maxisorp plates (Sigma-Aldrich), followed by
blocking with 0.5% BSA in PBS for 1 h. Antibodies were applied to the blocked and coated
plate for 30 min followed by washing and detection with an anti-kappa-chain-HRP conjugate
(Southern Biotech, 2060-05).
Biolayer Interferometry
For IgGs and peptide–IgG fusions, binding kinetics for the S-protein were
determined by BLI with an Octet HTX instrument (ForteBio), as described.[20] For peptides, biotinylated peptides (Table ) were immobilized on streptavidin-coated sensors that were
subsequently blocked with excess biotin. Following equilibration with assay buffer, the
loaded biosensors were dipped for 600 s into wells containing 3-fold serial
dilutions of S-protein ECD and subsequently were transferred back into the assay buffer
for 600 s. Binding response data were corrected by the subtraction of response from
a reference and were fitted with a 1:1 binding model using the ForteBio Octet Systems
software 9.0. We determined KD values at various peptide
loading densities to ensure that high densities were not impacting kinetic measurements at
the sensor.
Sequence Alignment and Analysis
Sequences of binding peptides were imported into Geneious R9 software (Biomatters Ltd.).
Peptides were sorted based on the library from which they originated and aligned
separately using the MUSCLE algorithm. In the case of peptides lacking cysteines, a strict
penalty was imposed on the formation of gaps during alignment. Sequence logos were created
using the weblogo server (https://weblogo.berkeley.edu/logo.cgi).
Production of Pseudoviruses
HEK-293 cells (ATCC) were seeded in a 6-well plate at 3 × 105 cells/well
in DMEM (ThermoFisher, 11995-065) supplemented with 10% FBS and 1%
penicillin–streptomycin (Gibco, 15140122) and grown overnight at 37 °C with 5%
CO2. HEK-293 cells were cotransfected with 1 μg pNL4-3.luc.R-E-plasmid
(luciferase-expressing HIV-1 with a defective envelop protein) (NIH AIDS Reagent Program,
ARP2128) and a 0.06 μg CMV-promoter-driven plasmid encoding the S-protein using the
Lipofectamine 2000 transfection reagent (ThermoFisher, 11668027). Pseudovirus particles
were harvested by collecting the supernatant 48 h after transfection and were
filter-sterilized (0.44 μm, Millipore-Sigma, SLHA033SS).
Pseudovirus Infection Assays
HEK293T cells stably overexpressing the full-length human ACE2 protein were seeded in
96-well white polystyrene microplates (Corning, CLS3610) at 3 × 104
cells/well in DMEM (10% FBS and 1% penicillin–streptomycin) and were grown
overnight at 37 °C with 5% CO2. Pseudovirus particles were mixed with Ab,
incubated at room temperature for 10 min, and added to the cells. The cells were incubated
at 37 °C with 5% CO2 and the medium was replaced with fresh DMEM (10% FBS
and 1% penicillin–streptomycin) after 6 h and again every 24 h up to 72 h. To
measure the luciferase signal (pseudovirus entry), DMEM was removed and DPBS
(ThermoFisher) was added to cells before mixing with an equal volume of the ONE-Glo EX
Luciferase Assay System (Promega, E8130). Relative luciferase units were measured using a
BioTek Synergy Neo plate reader (BioTek Instruments Inc.). The data were analyzed by
GraphPad Prism Version 8.4.3 (GraphPad Software, LLC).
Authentic virus infection assays
To evaluate the neutralization of SARS-Cov-2 VoC, the authentic isolates of USA_WA1/2020,
B.1.1.7, B.1.351, B.1.617.2, B.1.617.2+, B.1.1.529 (BA.1), and the chimeric B.1.1.28 (that
was generated on the genetic background of WA1/2020) were evaluated for neutralization
using a focus reduction neutralization test and potency was estimated as
described.[20,33,34] All viruses were verified by next-generation sequencing.
Alternatively, authentic SARS-CoV-2 strains were used in a microneutralization assay, as
described.[35] In brief, serial fourfold dilutions of antibody,
starting from 10 μg/mL, were preincubated with 102 focus-forming units
per 100 μL at 37 °C for 1 h. The antibody–virus mixture was transferred
to 96-well tissue culture plates containing Vero-hACE2-TMPRSS2 cell monolayers in
duplicate and incubated at 37 °C and 5% CO2 for 1 h; subsequently, the
cells were overlaid with 1% (wt/vol) methylcellulose in MEM. Plates were harvested at 24 h
by the removal of overlays and fixation for 20 min with 4% paraformaldehyde in PBS at room
temperature. Plates were washed with permeabilization buffer (PBS supplemented with 0.1%
saponin and 0.1% BSA) and sequentially incubated overnight with the oligoclonal pool of
anti-SARS-Cov-2 Abs.[9] The plates with Omicron (B.1.1.529) were
additionally incubated with the cross-reactive pool of Abs for SARS-CoV-1 that bind to the
RBD and compete with Ab (VanBlargan et al.). The plates were incubated with HRP-conjugated
antimouse (Sigma-Aldrich, A5278), goat antihuman Abs (Sigma-Aldrich, A6029), and
visualized by the KPL TrueBlue substrate and quantitated by an Immunospot microanalyzer
(Cellular Technologies), as described.[36,37]
Structural Analysis
The model of two 15033-7 Fabs bound to the SARS-CoV-2 S-protein (PDB entry 7KMK) was
imported into PyMol (DeLano Scientific, LLC). Distances between the HC N-termini and the
15033-7 epitope on the unbound RBD were measured as the distance between the
Cα of the first residue of the HC and the Cβ of
Tyr489 of the S-protein. To build the model of two RGN10933 Fabs bound to the
S-protein, the model of RGN10933 and RGN10987 bound to the RBD (PDB entry 6XDG) was imported
into PyMol along with the data from PDB entry 7KMK. The data from PDB entry 6XDG were
duplicated, and the RBDs of the model were superposed with the two RBDs in the
“up” position in the model from PDB entry 7KMK. The RBDs in the model from PDB
entry 6XDG, RGN10987 Fabs, and 15033-7 Fabs were then eliminated from the model, leaving
only the two RNG10933 from PDB entry 6XDG bound to the S-protein from PDB entry 7KMK.
Distances were measured the same way as for 15033-7.
Authors: Raoul De Gasparo; Mattia Pedotti; Luca Simonelli; Petr Nickl; Frauke Muecksch; Irene Cassaniti; Elena Percivalle; Julio C C Lorenzi; Federica Mazzola; Davide Magrì; Tereza Michalcikova; Jan Haviernik; Vaclav Honig; Blanka Mrazkova; Natalie Polakova; Andrea Fortova; Jolana Tureckova; Veronika Iatsiuk; Salvatore Di Girolamo; Martin Palus; Dagmar Zudova; Petr Bednar; Ivana Bukova; Filippo Bianchini; Dora Mehn; Radim Nencka; Petra Strakova; Oto Pavlis; Jan Rozman; Sabrina Gioria; Josè Camilla Sammartino; Federica Giardina; Stefano Gaiarsa; Qiang Pan-Hammarström; Christopher O Barnes; Pamela J Bjorkman; Luigi Calzolai; Antonio Piralla; Fausto Baldanti; Michel C Nussenzweig; Paul D Bieniasz; Theodora Hatziioannou; Jan Prochazka; Radislav Sedlacek; Davide F Robbiani; Daniel Ruzek; Luca Varani Journal: Nature Date: 2021-03-25 Impact factor: 49.962
Authors: Luc R Desnoyers; Olga Vasiljeva; Jennifer H Richardson; Annie Yang; Elizabeth E M Menendez; Tony W Liang; Chihunt Wong; Paul H Bessette; Kathy Kamath; Stephen J Moore; Jason G Sagert; Daniel R Hostetter; Fei Han; Jason Gee; Jeanne Flandez; Kate Markham; Margaret Nguyen; Michael Krimm; Kenneth R Wong; Shouchun Liu; Patrick S Daugherty; James W West; Henry B Lowman Journal: Sci Transl Med Date: 2013-10-16 Impact factor: 17.956
Authors: James Brett Case; Paul W Rothlauf; Rita E Chen; Zhuoming Liu; Haiyan Zhao; Arthur S Kim; Louis-Marie Bloyet; Qiru Zeng; Stephen Tahan; Lindsay Droit; Ma Xenia G Ilagan; Michael A Tartell; Gaya Amarasinghe; Jeffrey P Henderson; Shane Miersch; Mart Ustav; Sachdev Sidhu; Herbert W Virgin; David Wang; Siyuan Ding; Davide Corti; Elitza S Theel; Daved H Fremont; Michael S Diamond; Sean P J Whelan Journal: Cell Host Microbe Date: 2020-07-03 Impact factor: 21.023
Authors: Johanna Hansen; Alina Baum; Kristen E Pascal; Vincenzo Russo; Stephanie Giordano; Elzbieta Wloga; Benjamin O Fulton; Ying Yan; Katrina Koon; Krunal Patel; Kyung Min Chung; Aynur Hermann; Erica Ullman; Jonathan Cruz; Ashique Rafique; Tammy Huang; Jeanette Fairhurst; Christen Libertiny; Marine Malbec; Wen-Yi Lee; Richard Welsh; Glen Farr; Seth Pennington; Dipali Deshpande; Jemmie Cheng; Anke Watty; Pascal Bouffard; Robert Babb; Natasha Levenkova; Calvin Chen; Bojie Zhang; Annabel Romero Hernandez; Kei Saotome; Yi Zhou; Matthew Franklin; Sumathi Sivapalasingam; David Chien Lye; Stuart Weston; James Logue; Robert Haupt; Matthew Frieman; Gang Chen; William Olson; Andrew J Murphy; Neil Stahl; George D Yancopoulos; Christos A Kyratsous Journal: Science Date: 2020-06-15 Impact factor: 47.728