Pedro F N Souza1,2, Maurício F vanTilburg3, Felipe P Mesquita2, Jackson L Amaral1, Luina B Lima2, Raquel C Montenegro2, Francisco E S Lopes1, Rafael X Martins4, Leonardo Vieira4, Davi F Farias4, Ana C O Monteiro-Moreira5, Cleverson D T Freitas1, Arnaldo S Bezerra3, Maria I F Guedes3, Débora S C M Castelo-Branco6, Jose T A Oliveira1. 1. Department of Biochemistry and Molecular Biology, Federal University of Ceará, Av Mister Hull, S/n-Pici, P.O. Box 60440-593, Fortaleza, Ceará 60020-181, Brazil. 2. Drug Research and Development Center, Department of Physiology and Pharmacology, Federal University of Ceará, Rua Coronel Nunes de Melo 100, P.O. Box 60430-275, Fortaleza, Ceará 60020-181, Brazil. 3. Biotechnology and Molecular Biology Laboratory, Renorbio, State University of Ceará, Av. Dr. Silas Munguba, 1700-Itaperi, P.O. Box 60714-903, Fortaleza, Ceará 60020-181, Brazil. 4. Laboratory for Risk Assessment of Novel Technologies (LabRisk), Department of Molecular Biology, Federal University of Paraiba, Campus I Lot. Cidade Universitaria, P.O. Box 58051-900, João Pessoa, Paraíba 58051-900, Brazil. 5. School of Pharmacy, University of Fortaleza, Av. Washington Soares, 1321, Edson Queiroz, P.O. Box 60811-905, Fortaleza, Fortaleza, Ceará 60811-690, Brazil. 6. Department of Pathology and Legal Medicine, Federal University of Ceará, Rodolfo Teófilo, P.O. Box 60010-681, Fortaleza, Ceará 60020-181, Brazil.
Abstract
The outbreak caused by SARS-CoV-2 has taken many lives worldwide. Although vaccination has started, the development of drugs to either alleviate or abolish symptoms of COVID-19 is still necessary. Here, four synthetic peptides were assayed regarding their ability to protect Vero E6 cells from SARS-CoV-2 infection and their toxicity to human cells and zebrafish embryos. All peptides had some ability to protect cells from infection by SARS-CoV-2 with the D614G mutation. Molecular docking predicted the ability of all peptides to interact with and induce conformational alterations in the spike protein containing the D614G mutation. PepKAA was the most effective peptide, by having the highest docking score regarding the spike protein and reducing the SARS-CoV-2 plaque number by 50% (EC50) at a concentration of 0.15 mg mL-1. Additionally, all peptides had no toxicity to three lines of human cells as well as to zebrafish larvae and embryos. Thus, these peptides have potential activity against SARS-CoV-2, making them promising to develop new drugs to inhibit cell infection by SARS-CoV-2.
The outbreak caused by SARS-CoV-2 has taken many lives worldwide. Although vaccination has started, the development of drugs to either alleviate or abolish symptoms of COVID-19 is still necessary. Here, four synthetic peptides were assayed regarding their ability to protect Vero E6 cells from SARS-CoV-2 infection and their toxicity to human cells and zebrafish embryos. All peptides had some ability to protect cells from infection by SARS-CoV-2 with the D614G mutation. Molecular docking predicted the ability of all peptides to interact with and induce conformational alterations in the spike protein containing the D614G mutation. PepKAA was the most effective peptide, by having the highest docking score regarding the spike protein and reducing the SARS-CoV-2 plaque number by 50% (EC50) at a concentration of 0.15 mg mL-1. Additionally, all peptides had no toxicity to three lines of human cells as well as to zebrafish larvae and embryos. Thus, these peptides have potential activity against SARS-CoV-2, making them promising to develop new drugs to inhibit cell infection by SARS-CoV-2.
The
outbreak of the coronavirus disease in late 2019 (COVID-19)
is still ongoing and has already claimed more than 4 million lives
worldwide.[1,2] The causative agent of COVID-19 was later
identified as SARS-CoV-2 (Severe Acute Respiratory Syndrome Coronavirus
2), which is 96 and 79% similar, respectively, to BatCoV RatG13 (a
coronavirus from bats) and another pandemic coronavirus, SARS-CoV.[3] Coronaviruses are a group of viruses belonging
to the family Coronaviridae, a group
of enveloped positive-stranded RNA viruses with both medical and veterinary
importance. The RNA genome produces a larger polyprotein composed
of non-structural (NSPs) and structural proteins (SPs). The NSPs are
classified from 1 to 16 and the SPs are identified as spike (S), envelope
(E), membrane (M), and nucleocapsid (N) proteins.[4]Due to the huge effort of scientists worldwide, an
unprecedented
breakthrough was achieved. In less than a year, vaccination campaigns
began in virtually all countries.[5] However,
two problems have arisen: (1) the rate of vaccination is not fast
enough to reach herd immunity worldwide, and (2) countries with low
vaccination rates are a repository of SARS-CoV-2 mutants that threaten
the efficacy of the vaccines, causing vaccinated people to be jeopardized
by COVID-19.[4,6] Most of the mutations in the SARS-CoV-2
genome are concentrated in the RBD domain of S protein. S protein
is the main focus of most vaccines developed, so mutations on it could
reduce their efficacy.[6] Other points are
that poor countries that cannot afford to buy vaccines are still struggling
with COVID-19, and the development of specific drugs for clinical
treatment of the symptoms is still necessary. Therefore, the development
of drugs for COVID-19, despite the availability of vaccines, is still
urgent.[7] Drug repositioning has been employed
to accelerate the process of drug development. This entails testing
drugs already approved for the treatment of other diseases against
SARS-CoV-2. Despite reaching good results in computational simulations,
in vitro tests have so far revealed the inefficiency of the targeted
drugs.[7]In this context, in two previous
studies, our research group employed
computational simulations to drive antimicrobial peptides toward the
S protein[8] and Mpro[9] of SARS-CoV-2. In the first study, Souza et al.[8] reported that out of 8 peptides, two peptides, Mo-CBP3-PepII and PepKAA, strongly bonded to
the S protein, leading to changes in structural conformation and interaction
with the ACE2 receptor. In the second study, Amaral et al.[9] employed the same peptides against Mpro. Of these,
three peptides, RcAlb-PepI, PepGAT, and PepKAA, interacted
the best with Mpro, leading to conformational changes and reduction
in the catalytic site. The results of both studies suggested that
those peptides could have anti-SARS-CoV-2 action in vitro, making
them potential candidates for the development of a specific drug to
treat COVID-19 symptoms.The four peptides tested here were
designed from plant proteins.[10−12]Mo-CBP3-PepII and RcAlb-PepI were designed, respectively,
from a chitin-binding protein
of Moringa oleifera and a 2S albumin
of Ricinus communis.[10,11] PepGAT and PepKAA were designed from a chitinase from Arabidopsis thaliana.[12] All peptides are positively charged, have hydrophobic ratio from
40 to 65%, were predicted to have antiviral action, and are resistant
to intestinal enzymes, indicating the potential for oral administration.[10−12] Here, we report the in vitro activity of those peptides against
SARS-CoV-2 in addition to their level of toxicity.
Results
Anti-SARS-CoV-2
Potential of Synthetic Peptides
To
test whether the peptides can suppress virus-induced cytopathic effects,
a simple and fast MTT [3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium]
bioassay for screening drugs against SARS-CoV-2 was used. The MTT
assay is a simple colorimetric test to assess cell metabolic activity
in a microplate reader. It is based on the activity of the NAD(P)H-dependent
cellular oxidoreductase, which transforms the tetrazolium bromide
from MTT into formazan, which has a purple color, allowing the determination
of presence of the number of viable cells.[13] Based on that, we expected that the cells infected by SARS-CoV-2
would not convert the MTT into formazan because they had been killed
by SARS-CoV-2, and the control cells (without SARS-CoV-2 exposure)
would convert all MTT into formazan because they would be metabolically
active (Figure ).
Figure 1
Inhibitory
activity of synthetic peptides against SARS-CoV-2 infection.
(A) Mo-CBP3-PepII, (B) RcAlb-PepI,
(C) PepGAT, and (D) PepKAA inhibiting Vero E6 cell infection by SARS-CoV-2
as revealed by MTT assay for cell viability. Peptides at different
concentrations were incubated with SARS-CoV-2 for 30 min; then, SARS-CoV-2
was mixed with Vero E6 cells for 4 h. After that, cells were incubated
for 96 h, and cell viability was evaluated by MTT assay. Vero E6 cells
with and without SARS-CoV-2 infection were used as controls for cell
viability. Data are shown as mean ± standard deviation of three
independent experiments. ***P < 0.001. The neutralizing
effect of peptides against SARS-CoV-2.
Inhibitory
activity of synthetic peptides against SARS-CoV-2 infection.
(A) Mo-CBP3-PepII, (B) RcAlb-PepI,
(C) PepGAT, and (D) PepKAA inhibiting Vero E6 cell infection by SARS-CoV-2
as revealed by MTT assay for cell viability. Peptides at different
concentrations were incubated with SARS-CoV-2 for 30 min; then, SARS-CoV-2
was mixed with Vero E6 cells for 4 h. After that, cells were incubated
for 96 h, and cell viability was evaluated by MTT assay. Vero E6 cells
with and without SARS-CoV-2 infection were used as controls for cell
viability. Data are shown as mean ± standard deviation of three
independent experiments. ***P < 0.001. The neutralizing
effect of peptides against SARS-CoV-2.All peptides were able to inhibit virus-induced cytopathic effects,
with EC50 values in the microgram/milliliter range (Table ). Thus, we performed
a test with a fixed concentration to investigate the inhibition of
virus-induced cytopathic effects at 0.15 and 0.30 mg mL–1. The positive control cells [exposed to dimethyl sulfoxide (DMSO)]
presented 100% viability, while the cells infected with SARS-CoV-2
presented about 40% viability (Figure ). At 0.15 mg mL–1, the peptides Mo-CBP3-PepII, RcAlb-PepI, PepGAT, and
PepKAA reduced the SARS-CoV-2 cytopathic effects by 60, 75, 65, and
90%, respectively (Figure ). At 0.30 mg mL–1, the reductions of SARS-CoV-2
cytopathic effects were 75, 90, 85, and 80%, respectively, for Mo-CBP3-PepII, RcAlb-PepI, PepGAT, and
PepKAA (Figure ).
The best peptides and concentrations were PepKAA at 0.15 mg mL–1 and RcAlb-PepI at 0.30 mg mL–1, in both cases reaching 90% cell viability.
Table 1
EC50 Values for Peptides Treatment
in Vero E6 Cell Culture Infected with SARS-CoV-2a
peptides
EC50 (mg mL–1)
CI95%
MoCBP3-PepII
0.09
0.04–0.19
RcAlbPepI
0.05
0.03–0.07
PepGAT
0.08
0.07–0.09
PepKAA
0.04
0.02–0.07
EC50: concentration of a drug that
gives a half-maximal response. CI95%: 95% confidence interval.
EC50: concentration of a drug that
gives a half-maximal response. CI95%: 95% confidence interval.To confirm the neutralizing effect
of peptides against SARS-CoV-2,
the peptides were first incubated with SARS-CoV-2 and then added to
a Vero E6 cell monolayer. All tested peptides (0.15 and 0.30 mg mL–1) significantly reduced the plaque formation in Vero
E6 cell co-cultures with SARS-CoV-2 (Figure ). The least effective peptide was Mo-CBP3-PepII, which inhibited plaque formation
only by about 35 at both concentrations (Figure A). In turn, the most effective peptide was
PepKAA, which achieved inhibition of around 60% at both concentrations
(Figure D).
Figure 2
Plaque reduction
neutralization effect of synthetic peptides against
SARS-CoV-2. (A) Mo-CBP3-PepII, (B) RcAlb-PepI, (C) PepGAT, and (D) PepKAA inhibiting SARS-CoV-2 plaque
formation on Vero E6 cells. One hundred PFU of SARS-CoV-2 were incubated
with the peptides at different at 37 °C for 1 h. Then, they were
added to pre-seeded Vero E6 cells at 90–100% confluence. After
fixation for 1 h, the overlay was removed, and cells were stained
with 0.5% crystal violet. Vero E6 cells with and without SARS-CoV-2
were used as negative controls for SARS-CoV-2 neutralization. Data
are shown as mean ± standard deviation of three independent experiments.
***P < 0.001.
Plaque reduction
neutralization effect of synthetic peptides against
SARS-CoV-2. (A) Mo-CBP3-PepII, (B) RcAlb-PepI, (C) PepGAT, and (D) PepKAA inhibiting SARS-CoV-2 plaque
formation on Vero E6 cells. One hundred PFU of SARS-CoV-2 were incubated
with the peptides at different at 37 °C for 1 h. Then, they were
added to pre-seeded Vero E6 cells at 90–100% confluence. After
fixation for 1 h, the overlay was removed, and cells were stained
with 0.5% crystal violet. Vero E6 cells with and without SARS-CoV-2
were used as negative controls for SARS-CoV-2 neutralization. Data
are shown as mean ± standard deviation of three independent experiments.
***P < 0.001.
Molecular Docking of Peptides toward the D614G Mutant of S Protein
The results of a recent study[8] revealed
that all peptides tested interacted with the S protein from the Wuhan
isolate of SARS-CoV-2. Of these, Mo-CBP3-PepII and PepKAA presented the strongest interaction energies.[8] In our study, sequencing revealed that the SARS-CoV-2
isolate contained D614G on the S protein (Figure S1). Thus, we performed a new docking analysis to see if the
peptides would still interact with this mutant protein.The
molecular docking analysis predicted that Mo-CBP3-PepII, RcAlb-PepI, PepGAT, and PepKAA would
interact in the S1 region of the mutant protein S (D614G) of SARS-CoV-2
(Figure A,C,E,G).
We observed alterations in the atomic positions of the S protein after
interaction with peptides, revealed by variations in RMSD values of
0.663, 0.661, 0.662, and 0.661 Å, respectively, to Mo-CBP3-PepII, RcAlb-PepI, PepGAT, and
PepKAA.
Figure 3
Peptides Mo-CBP3-PepII, RcAlb-PepI, PepGAT, and PepKAA interacting with the mutated S protein
(SD614G) of SARS-CoV-2. (A,C,E,G) General view of the interaction
between SD614G with Mo-CBP3-PepII, RcAlb-PepI, PepGAT, and PepKAA, respectively.
(B,D,F,H) Zoomed overview of the interaction and the hydrogen bonds
between SD614G and peptides Mo-CBP3-PepII, RcAlb-PepI, PepGAT, and PepKAA, respectively.
Toxicity of peptides to human cells.
Peptides Mo-CBP3-PepII, RcAlb-PepI, PepGAT, and PepKAA interacting with the mutated S protein
(SD614G) of SARS-CoV-2. (A,C,E,G) General view of the interaction
between SD614G with Mo-CBP3-PepII, RcAlb-PepI, PepGAT, and PepKAA, respectively.
(B,D,F,H) Zoomed overview of the interaction and the hydrogen bonds
between SD614G and peptides Mo-CBP3-PepII, RcAlb-PepI, PepGAT, and PepKAA, respectively.
Toxicity of peptides to human cells.Mo-CBP3-PepII presented the lowest
binding energy of interaction (LBEI), −676.1 kJ·mol–1 with SD614G. The interaction was supported
by hydrogen bonds among amino acid residues of SD614G Lys386, Asp389, Asn544, Arg567, Asn969, Ile973 and Arg983, with
respective distances of 2.8, 2.7, 2.6, 2.4, 3.3, 2.8, and 2.4 Å
(Figure B). Mo-CBP3-PepII also exhibited hydrophobic interactions
with the residues Leu518, Ser974, Leu517, Thr430, Phe565, His519, Cys391, Ala522, Gly545, Leu546, Ser982, and Leu390 of SD614G (Figure S2A).RcAlb-PepI
had a LBEI value with SD614G of −646.9 kJ·mol–1, supported by six
hydrogen bonds and 14 hydrophobic interactions with SD614G. The hydrogen bonds were with the residues Asp428, Leu546, Thr547, and Asp979, with distances
of 2.6, 2.6, 2.6, 2.7, 2.8, and 3.0 Å (Figure D). The hydrophobic interactions were formed
by Ile973, Thr430, Ser974, Leu518, Ser982, His519, Arg983, Leu517, Gln564, Phe565, Gly545, Leu390, Cys391, and Asn544 of SD614G (Figure S2B).The LBEI between PepGAT and SD614G was −605.5
kJ·mol–1. Ten hydrogen bonds between the PepGAT
and the amino acid residues Thr430, Asp979,
Asp428, Arg567, Ser975, and Ala520 of SD614G and 10 hydrophobic interactions between
Arg983, Ile973, Cys391, Leu390, Leu518, Leu517, Ser974, Asp40,
Phe565, Val42, and His519 supported
the PepGAT-SD614G complex (Figures F and S2C).PepKAA (LBEI, −779.4 kJ·mol–1) interacted
with SD614G by hydrogen bonds with the amino acid residues
Asp198, Glu516, Arg567, Asp571, Thr547, Thr573, and Asn544, with
distances of 2.7, 2.6, 2.8, 2.9, 3.1, 2.7, and 2.7 Å, respectively
(Figure H). PepKAA
also exhibited hydrophobic interactions with Tyr200, Leu517, Leu518, Ser974, Asn969, His519, Arg983, Ser975, Val976, Phe565, Ala522, Leu390, Gly545, Leu546, and Cys391 of
SD614G with PepKAA (Figure S2D).To find possible clinical applications of peptides, their
toxicity
to human cells was assessed (Figure ). The MTT assay revealed that the peptides were not
toxic to the human cells tested. All cells treated with peptides presented
100% viability (Figure S3). Additionally,
we performed two other experiments to evaluate the peptides’
safety on human cells. In the first experiment, we tested whether
the peptides would induce DNA damage (Figure A–C) to three cell lines: L929 fibroblast
cells from mice and two human lines—human fetal lung fibroblast
(MRC-5 line) and human keratinocytes (HaCaT line). At a concentration
of 1 mg mL–1, the comet assay revealed that all
peptides caused no damage to DNA (Figure A–C). In contrast, in the positive
control for DNA damage, the methyl methanesulfonate (MMS) (4 ×
10–5 M) agent caused severe damages to the cells’
DNA.
Figure 4
Assessment of toxicity of Mo-CBP3-PepII, RcAlb-PepI, PepGAT, and PepKAA to human cell lines. (A)
L929, (B) HaCat, and (C) MRC-5 lines were incubated with synthetic
peptides at a concentration of 1 mg mL to evaluate the damage to DNA
by comet assay. (D–F) Cell lines were incubated with peptides
as described and evaluated for viable cells and cells in apoptosis.
MMS (4 × 10–5 M) was employed as a positive
control for cell toxicity and healthy cells as a negative control
for toxicity. Data are shown as mean ± standard deviation of
three independent experiments. ***P < 0.001.
Assessment of toxicity of Mo-CBP3-PepII, RcAlb-PepI, PepGAT, and PepKAA to human cell lines. (A)
L929, (B) HaCat, and (C) MRC-5 lines were incubated with synthetic
peptides at a concentration of 1 mg mL to evaluate the damage to DNA
by comet assay. (D–F) Cell lines were incubated with peptides
as described and evaluated for viable cells and cells in apoptosis.
MMS (4 × 10–5 M) was employed as a positive
control for cell toxicity and healthy cells as a negative control
for toxicity. Data are shown as mean ± standard deviation of
three independent experiments. ***P < 0.001.The second experiment analyzed whether peptides
can induce apoptosis
in the same cell line at a concentration of 1 mg mL–1 (Figure D–F).
At that concentration of peptides, the treated cells presented no
characteristics of apoptosis (Figure D,E). In contrast, cells treated with 4 × 10–5 M presented all aspects of apoptosis, such as small
cell volume, fragmented nucleus, peripheral condensation of chromatin,
and apoptotic bodies (Figure D,E).
Toxicity of Peptides to Zebrafish Embryos
To assess
in depth the safety of peptides, the toxicology to zebrafish embryos
was evaluated (Figure ). The survival rates of zebrafish larvae and embryos after exposure
to 1 mg mL–1 of the peptides for 96 h were ≥90%
(Figure A). After
96 h of treatment, morphological analysis revealed no alterations
(nonlethal effects) in the embryos exposed to the control (Figure B) and peptides (Figure C–F). The
embryo coagulation rates were ≤20% in the control and peptide-tested
embryos (Figure B–F).
This is an expected and spontaneous natural process that happens in
zebrafish embryos, leading to a mortality rate of 5–25%.
Figure 5
Assessment
of toxicity of Mo-CBP3-PepII, RcAlb-PepI, PepGAT, and PepKAA to zebrafish embryos. (A)
Survival rate (%) of zebrafish embryos and larvae exposed to 1 mg
mL of each synthetic peptide and control (E3 medium) samples after
96 h. (B–F) Zebrafish larvae exposed to, respectively, control
(E3 medium), Mo-CBP3-PepII, RcAlb-PepI, PepGAT, and PepKAA 1 mg mL–1 for 96 h.
All organisms presented normal development even exposed to peptides.
Data are shown as mean ± standard deviation of three independent
experiments.
Assessment
of toxicity of Mo-CBP3-PepII, RcAlb-PepI, PepGAT, and PepKAA to zebrafish embryos. (A)
Survival rate (%) of zebrafish embryos and larvae exposed to 1 mg
mL of each synthetic peptide and control (E3 medium) samples after
96 h. (B–F) Zebrafish larvae exposed to, respectively, control
(E3 medium), Mo-CBP3-PepII, RcAlb-PepI, PepGAT, and PepKAA 1 mg mL–1 for 96 h.
All organisms presented normal development even exposed to peptides.
Data are shown as mean ± standard deviation of three independent
experiments.
Discussion
At
the beginning of the current outbreak, it was thought that vaccines
would be the only way to fight back SARS-CoV-2. Thus, an unprecedented
collaboration worldwide led to the development of vaccines in record
times.[14,15] However, together with the beginning of
vaccination came SARS-CoV-2 variants not affected or weakly affected
by the immune response produced by vaccines.[16] The first vaccines applied brought widespread hope that vaccination
would end the pandemic. However, what nobody expected was the emergence
of many SARS-CoV-2 variants due to mutations, reducing the vaccines’
efficiency. For example, a mutation in the RBD of the S protein at
position E484 reduced the SARS-CoV-2 neutralization by monoclonal
antibodies and convalescent sera.[16] Thus
far, there are more than 3.5 billion people vaccinated worldwide.[17,18] Even though this number has been reached in short a time, the WHO
and American Centers for Disease Control (CDC) have reported that
some fully vaccinated people have still been infected by SARS-CoV-2.
Thus, they are advising even fully vaccinated people to continue using
masks to prevent infection by variants.[19−21]The emergence
of these SARS-CoV-2 variants and the existence of
people at greater risk such as immunosuppressed patients and those
who cannot be vaccinated, such as young children, makes it important
to develop drug/treatments specific for SARS-CoV-2 that can abolish
or alleviate the symptoms of infected patients. To find such drugs
quickly, many groups have examined drug repositioning, so far without
success. Many antiviral drugs such as arbidol, an anti-influenza drug
targeting the S protein, and galidesivir, remdesivir, tenofovir, sofosbuvir,
and ribavirin, which target the RdRp, have been submitted to in silico
assays for use against SARS-CoV-2.[22,23] Of these drugs,
the studies with arbidol have gone the furthest. Wang et al.[24] reported that at a concentration of 1.9 mg mL–1, arbidol achieved EC50 to SARS-CoV-2.
They also suggested a dosage of 200 mg 3 times/day or even higher
to alleviate COVID-19 symptoms. In another study, Yang et al.[25] revealed that in a group of 82 health professionals
treated prophylactically with arbidol, 48 people (58.5%) were infected
by SARS-CoV-2 and hospitalized and 34 (41.5%) developed mild symptoms
of COVID-19. Altogether, these results suggest that arbidol is not
very effective, so new drugs need to be developed. Nevertheless, it
was approved for the treatment of patients with COVID-19 by the National
Medical Products Administration of China.In our previous study,[8] it was shown
that Mo-CBP3-PepII and PepKAA strongly
bind to the S protein, leading to changes in three-dimensional (3D)
conformational structures and misplaced interactions with the ACE2
receptor.[8] Here, we have shown that all
four peptides interact with SD614G, altering its conformational
structure (Figures and S2). Thus, it is feasible to suggest
that these peptides block the entrance of SARS-CoV-2 in cells by interacting
with the S protein. In silico studies predicted that PepKAA, by strongly
interacting with both S8 and SD614G (Figures and S2) and changing their 3D, could be the best
peptide to block SARS-CoV-2 entrance in cells. Those results were
confirmed by in vitro experiments (Figures and 2). Based on
all our results, PepKAA is the best peptide to inhibit SARS-CoV-2
SD614G in cells, by reducing the virus plaque number by
60% at the lowest concentration (0.15 mg mL–1).
This study is very important because the active concentration of PepKAA
is 47.5-fold lower than that of arbidol, the most widely studied drug
against SARS-CoV-2.[24,25]Studying the repositioning
of other antiviral drugs, Sacramento
et al.[26] employed a combination of daclatasvir
and sofosbuvir, two anti-HCV drugs (Hepatitis C virus), against SARS-CoV-2
at concentrations higher than that employed to treat HCV. Wang et
al.[27] reported that other antiviral drugs
such as penciclovir, ribavirin, faviparin, and remdesivir had EC50 values in vitro against SARS-CoV-2, respectively, of 0.1,
0.98, 0.062, and 0.062 mg mL–1. These results show
that PepKAA is a good candidate to be used as a source to develop
a new drug against SARS-CoV-2. Of these drugs, only arbidol and remdesivir
have been investigated beyond in vitro tests. However, as reported
above, arbidol was not very efficient.[25] Although approved by the FDA (U.S. Food and Drug Administration),
the results of remdesivir in clinical trials presented no significant
differences between the placebo and the drug-treated group of patients.[28] Therefore, the search for new drugs to block
SARS-CoV-2 is still necessary.Based on our data, of all peptides
tested here, PepKAA is likely
the best one to neutralize SARS-CoV-2 SD614G, by preventing
cell infection (Figures and 2). However, other peptides are also
promising to inhibit SARS-CoV-2, such as RcAlb-PepI
and PepGAT (Figures and 2). Here, data from in vitro (Figures and 2) and in silico tests (Figures and S2) of PepKAA are in
harmony and may explain its high efficiency against SARS-CoV-2. Souza
et al.[8] reported that PepKAA interacted
with the S protein with an LBEI of −715.6 kJ·mol–1, by inducing a conformational change in the S protein. The authors
showed that PepKAA led to the incorrect interaction of S protein and
ACE2 receptor, suggesting that the formation of the PepKAA-S protein
complex inhibits the entrance of SARS-CoV-2 in cells.[8]In another study, Amaral et al.[9] predicted
that PepKAA has a strong interaction with Mpro of SARS-CoV-2,
leading to conformational changes and reduction of the active site.
These findings suggest that PepKAA reduced the Mpro activity.
Mpro is vital for SARS-CoV-2 replication because SARS-CoV-2
is an RNA virus producing a polyprotein that is cleaved by Mpro, releasing
SARS-CoV-2 proteins.[4,9] Here, the molecular docking study
predicted that PepKKA has the highest LBEI for SD614G compared
to other peptides (Figure ). The LBEI of PepKAA to SD614G was −779.4
kJ·mol–1, which is very similar to that presented
for the wild-type S protein, −715.6 kJ·mol–1.[8] PepKAA did not interact close to the
G614 mutation, but the interaction was still important to induce conformational
changes in the SD614G protein, which could induce the wrong
interaction with the ACE2 receptor.Here, in all assays, PepKAA
was first incubated with SARS-CoV-2
for 30 min to interact with the SD614G protein before infecting
cells. Interestingly, even when the virus, peptides, and cells were
incubated at the same time, PepKAA could not prevent cell infection
by SARS-CoV-2 (data not shown). These results strongly suggest that
the mechanism behind PepKAA’s anti-SARS-CoV-2 activity is by
interacting with SD614G, inducing conformational changes
(as revealed by docking), wrongly interacting with the ACE2 receptor,[8] and thus inhibiting SARS-CoV-2 from invading
the cell. Also, it is feasible to suggest that by a so-far-unclear
mechanism, PepKAA inhibits Mpro activity, which also prevents
SARS-CoV-2 replication and infection. We hypothesize this because
PepKAA is a membrane-penetrating peptide[12] and SARS-CoV-2 has a lipid envelope.[4] Thus, PepKAA could also target the SARS-CoV-2 membrane, inactivating
it. In that case, PepKAA could inhibit SARS-CoV-2 by different mechanisms.Recently, other synthetic peptides have been tested against SARS-CoV-2.[29−31] Curreli et al.[29] have analyzed the interaction
of RBD from SARS-CoV-2 with human ACE2. From this analysis, the authors
designed and synthesized four peptides. Among those, the synthetic
peptide NYBSP-4 presented an IC50 of 1.97 μM against
SARS-CoV-2. This result is better than that showed for PepKAA, which
reached the same inhibition at a concentration of 12.1 μM. However,
compared to results reported by Larue et al.,[30] PepKAA seems to be more effective. Larue et al.[30] reported that synthetic peptides SAP1, SAP2, and SAP6 derived
from human ACE2 receptor displayed an IC50 toward SARS-CoV-2
at concentrations of, 2.39, 3.72, and 1.90 mM, respectively, which
are much higher than the concentration of 12.1 μM presented
by PepKAA. In another study, Han et al.[31] reported that the peptide GK-7, also derived from human ACE2, presented
an IC50 toward SARS-CoV-2 at a concentration of 3.8 μM,
which is lower than the concentration presented by PepKAA (12.1 μM).
It is important to notice that all studies performed the antiviral
assay using a pseudovirus expressing the S protein.[29−31] In our case,
we employed the entire natural virus isolated from a patient with
full fitness to infect cells in the assay, which is closer than what
occurs during infection. This could be an explanation for the elevated
concentration required for PepKAA to reach the same concentration
presented by other peptides. This result still highlights the efficiency
and the potential of PepKAA toward SARS-CoV-2.One important
feature of a candidate drug is safety. New drugs
must cause no or very low side effects. In the case of arbidol and
remdesivir, this is not true. Both have considerable side effects.[32−34] Yet despite this, they were approved for treatment given the emergency
faced by the population from the SARS-CoV-2 outbreak. Arbidol was
hastily approved to treat COVID-19 in China, even without satisfactory
results. The use of arbidol was associated in rats with loss of body
weight, loss of organ weight (mainly liver), and piloerection in females.[34] In the case of remdesivir, the side effects
reported were from patients with COVID-19 using the drug to treat
symptoms. The side effects reported were increased nausea, diarrhea,
vomiting, gastroparesis, atrial fibrillation, cardiac arrest, and
acute kidney injury.[35] It is important
to state that we are not criticizing or judging the use of these drugs.
Because they were approved, they should be used. However, it is also
urgently necessary to seek a new type of drug that has lower or no
toxic effects.In this context, PepKAA is a strong candidate.
PepKAA was meticulously
designed to prevent any kind of toxic effect.[12] Indeed, the bioinformatics analysis revealed no allergic or toxic
potential. The hemolytic analysis revealed only 5% chance to cause
hemolysis on erythrocytes. However, the in vitro tests against erythrocytes
revealed no hemolysis. PepKAA did not present any toxicity to Vero
cells.[12] Here, looking toward clinical
trials with PepKAA and/or other peptides, we performed additional
toxicity tests. Not only PepKAA but also all synthetic peptides presented
no genotoxicity or pro-apoptotic effects on human cells L929, MRC-5,
and HaCaT (Figure ) at a concentration of 1 mg mL–1, which is threefold
higher than the highest concentration tested (Figure ). To shed more light on the toxicity of
peptides, we employed an important tool for drug development: testing
by the zebrafish model.[36] At a concentration
of 1 mg mL, none of the peptides presented toxicity to zebrafish larvae
and embryos (Figure ). That concentration (1 mg mL–1) is 25-fold higher
than the EC50 value of PepKAA against SARS-CoV-2. In this assay, we
used 20 zebrafish embryos, and the embryos incubated with the peptides
had a survival of ≥95%. In the case of PepKAA, the survival
was 100% (Figure A),
and no damage was found in the zebrafish (Figure C–H). Altogether, the results of efficacy
against SARS-CoV-2 and safety indicated PepKAA as a potential substance
for testing the development of new drugs against SARS-CoV-2.PepKAA is a synthetic peptide. One question always arises when
working with synthetic peptides: are they cost-effective for commercial
use? During the 1990s, the employment of synthetic peptides was impossible
given the high cost of chemicals used in the synthesis combined with
a very low yield. However, new technologies allowing the recovery
and recycling of solvents used during synthesis have made it feasible
to produce with a kilogram scale, leading to a dramatic reduction
in the cost of peptide synthesis.[37] The
synthetic peptide Fuzeon is an example. It is a peptide used in the
treatment of HIV requiring kilogram-scale production (≥100
kg), which is possible due to these new technologies.[38] In 2018, the FDA approved a glucagon-like synthetic peptide,
Rybelsus, used in the treatment of type II diabetes.[39] Therefore, if the pharmaceutical industry demonstrates
interest, the application of synthetic peptides is surely practicable.
Based on the potential of our peptides, we filed a patent application
in Brazil with the National Industrial Property Institute, under number
BR 10 2020 023728 4.
Conclusions
Here, all synthetic
peptides were active against SARS-CoV-2 to
some extent. The ability to interact with SD614G provided
a clue about how peptides act to inhibit SARS-CoV-2. PepKAA was the
most prominent peptide to inhibit SARS-CoV-2 while showing no toxicity
to human cells and zebrafish embryos. PepKAA thus has a higher potential
to develop new anti-SARS-CoV-2 drugs that are effective without adverse
effects.
Methodology
Ethical Statement
This experiment
conducted with a
human patient in this study was approved by the Research Ethics Committee
involving human beings on the use of humans in experiments of the
Federal Univerisity of Ceará, with authorization documented
by protocol no. 4.029.490. All experiments with SARS-CoV-2 were done
in accordance with relevant guidelines and regulations. Additionally,
the informed consent was obtained from all participants and/or their
legal guardian(s).The experiments conducted with zebrafish
in this study were approved by the Ethics Committee in the Use of
Animals of the Federal University of Paraíba, with authorization
documented by protocol no. 4460140920. In addition to this, animal
use methods were carried out in compliance with the ARRIVE guidelines
and in accordance with relevant guidelines and regulations.
Sample
Collection
The clinical sample was collected
from a patient with positive real-time (RT)-qPCR result and presenting
symptoms of SARS-CoV-2 infection. A nasopharyngeal swab was used to
collect the sample, and the sample was placed into a 3 mL tube containing
a viral transportation solution, as described by Holshue et al.[40] All experiments were carried out in the biosafety
level 3 facility (NB-3) of the Laboratory of Emerging and Reemerging
Pathogens of the Federal University of Ceará (Fortaleza, Brazil).
Viral Isolation and Titration
The SARS-CoV-2 isolation
was performed following the protocol described by Harcourt et al.[41] with modifications. Vero E6 cells (ATCC number
CCL-81) were cultured in a Leibovitz medium (L-15) supplemented with
2% heat-inactivated fetal bovine serum (FBS) and 1% penicillin/streptomycin
solution (GIBCO). The frozen sample was thawed and passed through
a 22 μm syringe filter. The viral isolation was carried out
in a 96-well plate containing a Vero cell monolayer with 90–100%
confluence. Then, 50 μL of the L-15 medium without FBS was added
to 100 μL of the clinical material. Then, the plate was incubated
for 2 h at 37 °C, with shaking every 15 min to facilitate the
infection of cells by the virus. After the incubation, the medium
was removed and added to all wells containing 100 μL of L-15
with 2% FCS and 1% penicillin/streptomycin, incubated at 37 °C
and observed daily for the presence of cytopathic effects. The material
in the wells in which cytopathic effects were observed was submitted
to confirmatory testing using RT reverse transcription PCR. The supernatant
from the infected Vero cells was collected, placed into cryotubes,
and stored at −80 °C. The virus titration was done following
the method described by Mendoza et al.[42]
RT PCR and Sequencing
SARS-CoV-2 RNA was extracted
using the QIAamp Viral RNA kit (Qiagen) following the manufacturer’s
instructions for SARS-CoV-2 detection using a one-step procedure.
For RT PCR (qPCR), the CDC 2019-nCoV qPCR diagnostic panel was followed,
using specific primers to confirm the presence of SARS-CoV-2 in all
cell cultures. In this kit, the primer–probe mixes target two
regions of the nucleocapsid gene (N1 and N2) as well as the human
endogenous control (RNase P gene), a control for sample integrity.Thereafter, a specific RBD region of the spike gene was amplified
from the cDNA sample using the paired primers (F-AATCTATCAGGCCGGTAGCAC
and R-CACCAATGGGTATGTCACACT) and Platinum Taq DNA Polymerase High-Fidelity
kit (Invitrogen). The PCR product was analyzed by electrophoresis
through 1.5% agarose gel and purified using the PureLink PCR Purification
Kit (Thermo Fisher). The purified product was sequenced using the
BigDye Terminator v1.1 Cycle Sequencing Kit (Thermo Fisher) according
to the manufacturer’s instructions. The basic local alignment
search tool software was used for the computer analysis of sequence
data with the reference sequence from Wuhan, China (NC_045512.2).By amplifying the N gene as a target, we confirmed the presence
of SARS-CoV-2. After the confirmation, we performed partial sequencing
of the RBD region of the S protein, revealing the presence of mutation
A23231G, which corresponds to the D614G mutation in the spike protein
(SD614G).
Anti-SARS-CoV-2 Activity of Peptides by the
MTT Assay
The peptides were diluted to a concentration ranging
from 0 to 2.5
mg mL–1 in the L-15 medium (Cultilab, Brazil) without
fetal serum and filtered through 22 μm filters. The peptides
were mixed with an equal volume of viral solution and incubated for
30 min at 37 °C. After incubation, the mixture from each peptide
dilution and a multiplicity of infection (MOI) of 1.85 were added
in triplicate to a 96-well plate containing 2.5 × 105 cells per well. After 2 h of incubation at 37 °C, the virus-containing
mixtures were removed from the wells and replaced by fresh L-15 medium
containing 2% FBS and 1% antibiotic. The plate was incubated for 4
days at 37 °C. The negative control wells received only the culture
medium, and the positive control wells received the virus. After this
period, the medium was removed, and 50 μL of 5 mg mL–1 MTT (3-4,5-dimethylthiazol-2-yl)-25 diphenyl tetrazolium bromide
(Life Technologies, USA) was added. The plate remained for 4 h at
37 °C, after which the solution was removed, and the formazan
crystals were diluted with 50 μL DMSO. The plate was read at
540 nm, and the percentage of protection (PP) was calculated by the
following formula: PP = [(AB)/(CB) × 100], where A, B, and C
indicate the absorbance of the peptide, virus, and control cells,
respectively.[43]
Plaque Reduction Neutralization
Tests
The neutralization
assay was performed as described by Muruato et al.[44] Vero E6 cells (6 × 105 per well) were seeded
in 12-well plates and left at rest overnight. Then, SARS-CoV-2 samples
at an MOI of 1.85 were incubated with the peptides at concentrations
of 0.15 and 0.30 mg mL–1 at 37 °C for 1 h.
The virus–peptide mixture was added in triplicate to pre-seeded
Vero E6 cells at 90–100% confluence. After 2 h of incubation
at 37 °C, 1 mL of the overlay containing 1.5% carboxymethylcellulose
in L-15 containing 2% FBS and 1% penicillin/streptomycin antibiotics
(GIBCO) was added to the infected cells. After 3 days of incubation,
1 mL of 3.65% formaldehyde in phosphate-buffered saline (PBS) was
added to the overlay-covered cells. After fixation for 1 h, the overlay
was removed, and the contents of the cells were stained with 0.5%
crystal violet. The plates were washed with water to remove excess
dye, photographed, and submitted for counting using the ImageJ program.
Percentage plaque reduction was calculated using the following formula:
[(sample × 100)/positive control] −100.
Molecular Docking
Analysis
The crystallographic data
of the mutant (D614G) SARS-CoV-2 protein S was obtained from the Protein
Data Bank, with accession number PDB ID: 7DX1. The 3D structures of
peptides Mo-CBP3-PepII, RcAlb-PepI, PepGAT, and PepKAA were predicted using the PEP-FOLD 3.5
software.[45] The protein and peptide structures
were determined and the protonation states were adjusted using the
ProteinPrepare software.[46]To carry
out molecular docking between the peptides and the mutated protein
S of SARS-CoV-2, we used the ClusPro 2.0 server,[47] which showed the best results in the CAPRI challenge.[48] The results were analyzed using the number of
members in each cluster and the lowest energies calculated through
the Balanced software method, which considers the energies obtained
from electrostatic and hydrophobic interactions.To analyze
the interactions between the peptides and the mutated
S protein, the software LigPlot+ v. 2.2.4 was used.[49] The preparation of figures and measurement of
the variation in rmsd were performed with the Pymol software.
Assessment
of Cytotoxicity
The cytotoxicity was quantified
by the ability of live cells to reduce the yellow dye 3-(4,5-dimethyl-2-thiozolyl)-2,5-diphenyl-2H-tetrazolium
bromide (MTT) to formazan.[50] Cytotoxicity
was checked against L929 (murine fibroblasts, ATCC number CCL-1),
MRC-5 (human lung fibroblasts, ATCC number CCL-171), and HaCat (human
keratinocytes), provided by the Rio de Janeiro Cell Bank (BCRJ, Brazil).
All cell lines were washed and resuspended in the DMEM medium supplemented
with 10% FBS, 2 mM of glutamine, 100 U/mL of penicillin, and 100 μg/mL
of streptomycin, at 37 °C under 5% CO2. For the experiments,
cells were plated in 96-well plates (0.1 × 106 cells/mL
for HaCat cells and 0.1 × 104 cells/mL for L929 and
MRC-5 cell lines). After 24 h, the tested peptides (1 mg mL–1 in the culture medium) were added to each well, and the cells were
incubated for 72 h. MMS (4 × 10–5 M) was used
as the positive control. Thereafter, the plates were centrifuged,
and the medium was replaced with a fresh medium (150 μL) containing
0.5 mg mL of MTT. Three hours later, the MTT formazan product was
dissolved in 150 μL DMSO, and the absorbance was measured using
a multiplate reader (Spectra Count, Packard, Ontario, Canada). Drug
effect was quantified as the percentage of control absorbance of the
reduced dye at 595 nm.
Comet Assay
For this assay, the
concentration of peptides
was 1 mg mL–1. DMSO–NaCl was the negative
control for damage, and MMS (4 × 10–5 M) was
used as the positive control for DNA damage. The standard alkaline
comet assay (single-cell gel electrophoresis) was performed as previously
described.[51] After treatment (24 h), cells
were washed with ice-cold PBS, trypsinized, and resuspended in the
complete medium. Then, 20 μL of cell suspension (0.7 ×
105 cells/mL) was dissolved in 0.75% low-melting-point
agarose and immediately spread onto a glass microscope slide pre-coated
with a layer of 1% agarose with a normal melting point. The agarose
was allowed to set at 4 °C for 5 min. Slides were incubated in
an ice-cold lysis solution (2.5 M NaCl, 0.01 M Tris, 0.1 M EDTA, 1%
Triton X-100, and 10% DMSO, pH 10.0) at 4 °C for at least 1 h
to remove cell membranes, leaving DNA as “nucleoids”.After that, the slides were placed in a horizontal electrophoresis
unit and incubated with a fresh buffer solution (0.3 M NaOH, 0.001
M EDTA, pH 13.0) at 4 °C for 20 min to allow DNA unwinding and
the expression of alkali-labile sites. Electrophoresis was conducted
for 20 min at 25 V (94 V/cm). All the above steps were performed in
the dark to prevent additional DNA damage. Slides were neutralized
(0.4 M Tris, pH 7.5) and stained using 20 μg mL ethidium bromide
(EB). One hundred and fifty cells (50 cells from each of the three
replicate slides for each treatment) were selected, coded, and blindly
analyzed for DNA migration. These cells were visually scored according
to the tail length into five classes: (1) class 0: undamaged, without
a tail; (2) class 1: with a tail shorter than the diameter of the
head nucleus; (3) class 2: with a tail length 1–2× the
diameter of the head; (4) class 3: with a tail longer than 2×
the diameter of the head; and (5) class 4: comets with no heads. The
damage index (DI) value was assigned to each sample. DI is an arbitrary
score based on the number of cells in the different damage classes,
which are visually scored by measuring the DNA migration length and
the amount of DNA in the tail. DI ranges from 0 (no tail: 100 cells
× 0) to 400 (with maximum migration: 100 cells × 4).[52]
Morphological Characterization of Apoptotic
PBLs
For
this assay, the concentration of peptides was 1 mg mL–1. DMSO–NaCl was the negative control for damage, and MMS (4
× 10–5 M) was used as the positive control
for DNA damage. The peptide and control solutions were incubated as
described above. Then, cells with morphological characteristics of
apoptosis (i.e., small cell volume, peripheral condensation of chromatin,
fragmented nucleus, and apoptotic bodies) were determined after each
treatment (24 h) by the acridine orange (AO)/EB staining assay: 25
μL of the cell suspension was mixed with 1 μL of the staining
solution (100 μg/mL AO + 100 μg/mL EB in PBS) and spread
on a slide, where 300 cells were counted per data point. The percentage
of apoptotic cells was then calculated.[53]
Zebrafish Toxicity
Zebrafish Embryos
Zebrafish embryos
(AB wild-type strain)
with approximately 1 HPF (hour post-fertilization) were provided by
the Production Unit for Alternative Model Organisms (UniPOM), Federal
University of Paraiba, João Pessoa, Brazil. The parents were
maintained in a recirculation system with regular monitoring of water
quality parameters (pH, ammonia, and nitrite levels) in a room with
controlled temperature (26 ± 1 °C) and photoperiod (14:10
light/dark cycle). Fish were fed daily with commercial feed (Color
Bits Tetra, Melle, Germany) and freeze-dried spirulina (Fazenda Tamanduá,
Patos, Brazil) and were monitored for abnormal behavior or disease
development.Before the experiment, adult zebrafish (male-to-female
ratio of 2:1) were transferred to a 7 L spawning tank with a bottom
mesh and a quick-opening valve for embryo collection. Embryos were
collected on the day of the experiment and cultured in an adapted
embryonic medium E3 (5.0 mM NaCl, 0.17 mM KCl, 0.33 mM CaCl, and 0.33
mM MgSO4) containing 0.005% methylene blue. Only spawning samples
with a fertilization rate ≥90% were used. Viable embryos (normal
cleavage pattern and without morphological changes) were selected
under an inverted light microscope (Televal 31, Zeiss, Germany), at
50× magnification.
Acute Toxicity Test
The fish embryo
acute toxicity
test was independently conducted with four peptides according to the
OECD Guideline number 236[54] adapted for
96-well plates by Muniz et al.[55] Zebrafish
embryos with up to 3 hpf were exposed to 1 mg mL–1 of each sample. For each test and control sample, 20 wells were
filled with 0.3 mL of solution and 1 embryo.Additionally, 20
embryos were exposed only to the E3 medium (the solvent control).
Lethal and non-lethal effects were observed daily for 96 h. Embryos
showing lethality endpoints (coagulation, no formation of somites,
no detachment of tail, or absence of heartbeat) were considered dead.
This number was used to determine the survival percentage (number
of live organisms/total number of organisms × 100) per tested
sample. The exposures were under static conditions (without the renovation
of the exposure solution). Observations were performed with a stereomicroscope
(Olympus SZX7, Japan) at 56× magnification and photographed (Moticam
5+, China). After 96 h, surviving larvae were euthanized with eugenol
and properly discarded.
Statistical Analysis
The assays were performed in three
independent experiments. The statistics were expressed as the mean
± standard deviation. The data were submitted to ANOVA followed
by the Tukey test, using GraphPad Prisma 5.01, with a significance
of p < 0.05.
Authors: Meagan L Adamsick; Ronak G Gandhi; Monique R Bidell; Ramy H Elshaboury; Roby P Bhattacharyya; Arthur Y Kim; Sagar Nigwekar; Eugene P Rhee; Meghan E Sise Journal: J Am Soc Nephrol Date: 2020-06-08 Impact factor: 10.121
Authors: Michelle L Holshue; Chas DeBolt; Scott Lindquist; Kathy H Lofy; John Wiesman; Hollianne Bruce; Christopher Spitters; Keith Ericson; Sara Wilkerson; Ahmet Tural; George Diaz; Amanda Cohn; LeAnne Fox; Anita Patel; Susan I Gerber; Lindsay Kim; Suxiang Tong; Xiaoyan Lu; Steve Lindstrom; Mark A Pallansch; William C Weldon; Holly M Biggs; Timothy M Uyeki; Satish K Pillai Journal: N Engl J Med Date: 2020-01-31 Impact factor: 91.245
Authors: Paul T Heath; Eva P Galiza; David N Baxter; Marta Boffito; Duncan Browne; Fiona Burns; David R Chadwick; Rebecca Clark; Catherine Cosgrove; James Galloway; Anna L Goodman; Amardeep Heer; Andrew Higham; Shalini Iyengar; Arham Jamal; Christopher Jeanes; Philip A Kalra; Christina Kyriakidou; Daniel F McAuley; Agnieszka Meyrick; Angela M Minassian; Jane Minton; Patrick Moore; Imrozia Munsoor; Helen Nicholls; Orod Osanlou; Jonathan Packham; Carol H Pretswell; Alberto San Francisco Ramos; Dinesh Saralaya; Ray P Sheridan; Richard Smith; Roy L Soiza; Pauline A Swift; Emma C Thomson; Jeremy Turner; Marianne E Viljoen; Gary Albert; Iksung Cho; Filip Dubovsky; Greg Glenn; Joy Rivers; Andreana Robertson; Kathy Smith; Seth Toback Journal: N Engl J Med Date: 2021-06-30 Impact factor: 91.245
Authors: Alexis Lamiable; Pierre Thévenet; Julien Rey; Marek Vavrusa; Philippe Derreumaux; Pierre Tufféry Journal: Nucleic Acids Res Date: 2016-04-29 Impact factor: 16.971
Authors: Jeffrey V Lazarus; Scott C Ratzan; Adam Palayew; Lawrence O Gostin; Heidi J Larson; Kenneth Rabin; Spencer Kimball; Ayman El-Mohandes Journal: Nat Med Date: 2020-10-20 Impact factor: 53.440
Authors: Tawanny K B Aguiar; Nilton A S Neto; Cleverson D T Freitas; Ayrles F B Silva; Leandro P Bezerra; Ellen A Malveira; Levi A C Branco; Felipe P Mesquita; Gustavo H Goldman; Luciana M R Alencar; Jose T A Oliveira; Ralph Santos-Oliveira; Pedro F N Souza Journal: Pharmaceutics Date: 2022-08-12 Impact factor: 6.525