Next Article in Journal
The Impact of YRNAs on HNSCC and HPV Infection
Next Article in Special Issue
Identification of Natural Lead Compounds against Hemagglutinin-Esterase Surface Glycoprotein in Human Coronaviruses Investigated via MD Simulation, Principal Component Analysis, Cross-Correlation, H-Bond Plot and MMGBSA
Previous Article in Journal
Overcoming EGFR Resistance in Metastatic Colorectal Cancer Using Vitamin C: A Review
Previous Article in Special Issue
Reversal of Azole Resistance in Candida albicans by Human Neutrophil Peptide
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biomedicines
Volume 11
Issue 3
10.3390/biomedicines11030682
biomedicines-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

In Vitro Pharmacokinetic Behavior of Antiviral 3-Amidinophenylalanine Derivatives in Rat, Dog and Monkey Hepatocytes

by
Katalin Lányi
1,
Katalin Monostory
2,
Torsten Steinmetzer
3,
Ákos Jerzsele
4,5 and
Erzsébet Pászti-Gere
4,*
1
Department of Food Hygiene, University of Veterinary Medicine, István utca 2, H-1078 Budapest, Hungary
2
Research Centre for Natural Sciences, Institute of Enzymology, Magyar Tudósok 2, H-1117 Budapest, Hungary
3
Department of Pharmacy, Institute of Pharmaceutical Chemistry, Philipps University Marburg, Marbacher Weg 6, G-35032 Marburg, Germany
4
Department of Pharmacology and Toxicology, University of Veterinary Medicine, István utca 2, H-1078 Budapest, Hungary
5
National Laboratory of Infectious Animal Diseases, Antimicrobial Resistance, Veterinary Public Health and Food Chain Safety, University of Veterinary Medicine, István utca 2, H-1078 Budapest, Hungary
*
Author to whom correspondence should be addressed.
Biomedicines 2023, 11(3), 682; https://doi.org/10.3390/biomedicines11030682
Submission received: 31 January 2023 / Revised: 17 February 2023 / Accepted: 21 February 2023 / Published: 23 February 2023
(This article belongs to the Special Issue Drug Discovery for Infectious Diseases)
Browse Figures
Versions Notes

Abstract

:
Type II transmembrane serine proteases represent pharmacological targets for blocking entry and spread of influenza or coronaviruses. In this study, the depletion rates of the 3-amidinophenylalanine (3-APhA)-derived matriptase/TMPRSS2 inhibitors MI-463, MI-472, MI-485 or MI-1900 were determined by LC-MS/MS measurements over a period of 300 min using suspensions of rat, dog and cynomolgus monkey primary hepatocytes. From these in vitro pharmacokinetic (PK) experiments, intrinsic clearance values (Clint) were evaluated, and in vivo pharmacokinetic parameters (hepatic clearance, hepatic extraction ratio and bioavailability) were predicted. It was found that rat hepatocytes were the most active in the metabolism of 3-APhA derivatives (Clint 31.9–37.8 mL/min/kg), whereas dog and monkey cells displayed somewhat lower clearance of these compounds (Clint 6.6–26.7 mL/min/kg). These data support elucidation of important PK properties of anti-TMPRSS2/anti-matriptase 3-APhAs using mammalian hepatocyte models and thus contribute to the optimization of lead compounds.

1. Introduction

Both matriptase and transmembrane serine protease 2 (TMPRSS2) belong to the type II transmembrane serine proteases (TTSPs), and these membrane-anchored enzymes are involved in several physiological and pathological processes [1]. Excessive TTSP activity can be correlated with the formation of several disorders, thus specific inhibition of matriptase- and TMPRSS2-caused proteolysis with 3-amidinophenylalanine (3-APhA) derivatives can be an effective pharmacological option to alleviate tumor formation and metastatic potential [2], to counteract cartilage damage [3] and to restore hepcidin-mediated tipped iron homeostasis [4].
Moreover, these enzymes appear to be responsible for proteolytic activation of viral glycoproteins, thus facilitating certain viruses’ replication and spread. TMPRSS2 was reported to cleave hemagglutinin (HA) of avian and mammalian influenza A viruses (IAV) with a monobasic cleavage site in mice [5]. Other viruses such as Middle East respiratory syndrome- (MERS-) and severe acute respiratory syndrome (SARS) coronaviruses (CoV) also require activation of spike (S) glycoprotein for early entry into the cells via host TMPRSS2-mediated proteolytic cleavage [6,7]. Matriptase, which is highly expressed in a membrane-bound form, can also cleave HA precursor [8].
A new type of coronavirus, SARS-CoV-2 emerged in 2019 in China and often caused serious pneumonia (COVID-19) in humans, becoming rapidly pandemic. Vaccination is one of the main therapeutic measures to limit the spread of viral infection; however, the appearance of mutated variants of concern (VOCs) necessitates the clinical application of antivirals such as the direct-acting nucleotide prodrug remdesivir [9,10] and the nucleoside molnupiravir as well as the orally available inhibitor of the viral M pro protease inhibitor Paxlovid (combination of nirmatrelvir and ritonavir) [11,12]. The non-selective serine protease inhibitors, camostate mesylate originally used in the therapy of postoperative reflux esophagitis and acute exacerbations of chronic pancreatitis and antifibrinolytic aprotinin, were repurposed as drug candidates for COVID-19; however, camostate failed during clinical testing [13]. So far, host-directed small molecule peptidomimetics have been developed, such as the sulfonylated 3-APhA derivatives [14,15] used in this study as well as several ketobenzothiazole (ktb)-based inhibitors [16,17] to impede the proteolytic activation of the SARS-CoV-2 S protein. For improving the stability and potency, the structure of the latter compound type was optimized. The most promising inhibitor, N-0385 (Ms-Gln-Phe-Arg-kbt), possesses an IC50 value of 1.9 ± 1.4 nM. It was also found that N-0385 was effective in inhibition of SARS-CoV-2 infection in human lung epithelium-derived human airway Calu-2 cells and in a K18-human ACE2 transgenic mouse model [18].
Several animal species have been identified to be susceptible to SARS-CoV-2 infections [19]. Although the disease seems to be of zoonotic origin, transmission of the pathogen from humans to animals was also reported. Dogs show some susceptibility to infection, possibly due to the presence of angiotensin-converting enzyme 2 (ACE2), as in humans with similar clinical signs of the disease, and it was reported that dogs could also produce positive virological tests and that virus transmission was possible from humans to dogs under natural conditions [20]. Non-human primates (NHP) have high susceptibility to SARS-CoV-2 infection, and their symptoms resemble human pathological conditions to a great extent. Type II pneumocyte hyperplasia and alveolar fibrosis were mainly found in lung tissue from cynomolgus macaques, proving that these animals can be successfully infected by airborne SARS-CoV-2 [21]. On the other hand, animal models enable in-depth characterization of SARS-CoV-2 pathophysiology and can contribute to the development and preclinical investigations of antiviral therapies. Severity of respiratory disorders in cynomolgus macaques was correlated well with the ages and underlying diseases (diabetes and hyperlipidemia) of the animals infected experimentally with coronavirus [22]. Other mammalian models including hamsters, ferrets, and mice with hACE2 Tg were found to be also susceptible to SARS-CoV-2, showing respiratory transmission; however, sex- and age-related severity differences can be mainly seen in cynomolgus monkey models [23].
The tested 3-APhA-derived inhibitors of trypsin-like serine proteases can serve as start structures for the development of potential antiviral drugs. The suppression of host cell proteases can be favorable versus classical antiviral therapies due to the low incidence of drug resistance [24,25]. The application of these drug candidates can hinder the entry of SARS-CoV-2 via preventing TMPRSS2-mediated activation of S protein necessary for the attachment of virus particles to the host cells through surface ACE2 receptors [26]. Combinative administration of furin and TMPRSS2 inhibitors due to their different cleavage sites of S (S1/S2 and S2′) was proven to have excellent antiviral properties in airway epithelial cells infected with SARS-CoV-2; therefore, their simultaneous application can lead to a more effective synergistic effect. This was the first study in which peptidomimetic sulfonylated 3-APhA-based MI compounds including MI-432 and MI-1900 were used to effectively decrease SARS-CoV-2 virus titer [27].
Until now, no predicted pharmacokinetic (PK) data have been available for 3-APhA-based inhibitors after their interaction with hepatocytes of different animal species. In this study, depletion % of the applied inhibitors (MI-463, MI-472, MI-485 and MI-1900) was determined in primary rat, dog and cynomolgus monkey hepatocytes using LC-MS/MS measurements. Interspecies comparisons of depletion of inhibitors over the chosen time window (0–300 min) were also performed to assess the impact of certain chemical structures and the effect of species on biotransformation rates of these 3-APhA derivatives. In addition, in vitro intrinsic clearance (Clint) values were evaluated for each MI compound, and in vivo PK parameters, such as predicted Cl (ClH), hepatic extraction ratio (EH) and bioavailability (F), were predicted.

2. Materials and Methods

2.1. Preparation of Primary Hepatocyte Suspensions

Cryopreserved rat, dog and cynomolgus monkey primary hepatocytes were obtained from Primacyt (Primacyt, Schwerin, Germany). Hepatocytes were thawed using cryopreserved hepatocyte thawing medium (Primacyt, Schwerin, Germany), and cells were centrifuged (10 min, 100× g, 20 °C) to remove supernatants. The viability of hepatocytes was examined by the Trypan blue exclusion method, and it was confirmed that more than 90% of the cells were alive. Cell suspensions were then diluted with serum-free maintenance medium (Primacyt, Schwerin, Germany) after cell counting in Bürker chamber to the appropriate concentration (original concentration of the viable hepatocytes was 107 cells/mL).

2.2. Incubation of Primary Hepatocytes with Inhibitors

Stock solutions of MI-463, MI-472, MI-485 and MI-1900 (10 mM) in dimethyl-sulfoxide were kept at −20 °C. The structures of the inhibitors are shown in Figure 1. Time courses of the unchanged amount of the inhibitors (in %) after incubation with hepatocytes were measured. Each compound was incubated in cell suspension (at the concentration of 5 × 105 cells/mL in 6-well plates) at constant and low-speed shaking at 37 °C in a humid atmosphere containing 5% CO2. The solutions of the inhibitors at 10 µM were prepared with maintenance medium freshly from the stock prior to each experiment (DMSO concentration did not exceed 0.1% in the final solutions). Cells were incubated with serum-free maintenance medium either without or with each selected MI for 5 h. At various time points (0, 30, 60, 120, 180, and 300 min), the incubation mixtures were sampled (aliquots: 0.15 mL) and terminated by the addition of 0.1 mL ice-cold acetonitrile. The cell debris was separated by centrifugation, and the supernatants were analyzed by LC-MS/MS for quantitation of the remaining amount of the inhibitors.

2.3. Determination of Unchanged Parent Compound Concentrations

For the LC-MS/MS investigations, 120 μL from the incubation mixture samples obtained as described above was removed and mixed with 315 μL chilled acetonitrile and 15 μL triphenylphosphate solution (TPP; 12000 ng/mL in acetonitrile) used as internal standard for the quantitation. After vortex mixing, the samples were centrifuged at 4 °C and 13300 RPM for 10 min. The supernatant was filtered through a PVDF syringe filter of 0.22 μm pore diameter, then transferred to chromatographic vials.
The LC-MS/MS instrument used for analyzing the MI compounds was a Shimadzu LCMS 8030+ system operated with an electrospray ion source (ESI) in positive multiple reaction monitoring (MRM) mode. A Phenomenex Kinetex C18 EVO 100 × 4.6 mm ID (2.6 µm particle size) column equipped with a 40 × 2 mm C18 guard column was used for the chromatographic separation. Gradient elution was used with eluents “A”: 50 mM ammonium acetate in water (pH 5 set by acetic acid) and “B”: 0.1 v/v% formic acid in acetonitrile. The flow rate was 0.4 mL/min, a chromatographic run took 8 min. The column oven was set to 30 °C, and the autoinjector’s temperature was at 6 °C. The injected volume was 10 µL. Other MS parameters were as follows: interface 4.5 kV; interface temperature 350 °C; desolvation line 300 °C; heat block 400 °C; detector 1.78 kV, nebulizing gas (N2) 3 liters/min, drying gas (N2) 15 liters/min; collision gas (Ar) 230 kPa.

2.4. Assessment of Pharmacokinetics Parameters

The intrinsic clearance for hepatocytes (Clint) [mL (min 5 × 105 cells) −1] was calculated from the decrease in the concentration of MI compounds as follows: Clint = β = ln 2 t 1 / 2 [28] For scaling up, the Clint value to obtain Clint per whole liver (g)/bw(kg), the cell concentration in liver, the average liver weight and average body weight parameters were used (Table 1) [29]. The value for hepatic clearance (ClH) was calculated as follows [28]:
C l H = C l int l i v e r / b w ∗ f u ∗ Q p l a s m a ( C l int l i v e r / b w ∗ f u ) + Q p l a s m a
where the flow rate of Qplasma = QH* plasma/blood ratio. To calculate ClH, the hepatic flow rate, plasma/blood ratio (Table 1) and fu = 1 values were used [30]. (QH is the hepatic blood flow, while fu is the unbound fraction of the compound.) The hepatic extraction ratio was defined as E = ClH/QH, and the bioavailability (%) was F = (1 − E) * 100 [28,31,32,33].

2.5. Statistical Analysis

Data of normal distribution and homogeneity of variances were evaluated with one-way ANOVA and Tukey’s post hoc test using the R version 4.0.4 software package (Vienna, Austria, 2021). The results were expressed as means ± SD, and differences between obtained results were determined to be significant if p was <0.05 (* p  <  0.05, ** p  <  0.01, *** p  <  0.001). Every measurement was repeated at least in triplicate with 5% interday and intraday precision.

3. Results

In rat hepatocytes, all inhibitors were significantly depleted even at 60 min sampling time (pMI-463, 60 min = 0.02755; pMI-472, 60 min, pMI-485, 60 min and pMI-1900, 60 min < 0.001). It was also ascertained that conversions of the inhibitors were elevated to a greater extent independently of the structural differences in the 3-APhA derivatives (p < 0.001 in each case) at 120, 180 and 300 min. Comparing the depletion rates of 3-APhA derivatives with various structural moieties showed that they did not differ from each other significantly at 300 min sampling time (p values were in the range of 0.3367 and 1); however, they were depleted significantly after 5 h incubation with hepatocytes (p < 0.001) (Figure 2). It is also interesting to note that although not much difference in hepatic extraction ratios (E) for MI compounds was predicted, the 0.3 cut-off value between low and intermediate extraction drugs assigned MI-472 and MI-485 as intermediate, while MI-463 and MI-1900 were assigned as low-extraction drugs.
In dog hepatocytes, MI-472, MI-485 and MI-1900 were significantly depleted even at 30 min and 60 min sampling times (pMI-472, 30 min = 0.0412; pMI-485, 30 min and pMI-1900, 30 min < 0.001; pMI-472, 60 min = 0.0014; both pMI-485, 60 min and pMI-1900, 60 min < 0.001). MI-463 content slowly decreased from 180 min (pMI-463, 180 min = 0.0054 and pMI-463, 300 min = 0.0042), in contrast to accelerated transformation rates of MI-485 and MI-1900, which could suggest that minor structural variations present in the 3-APhA derivatives could cause different depletion rates observed in dog hepatocytes (Figure 3). Due to the slow decomposition of the MI-463, significant differences in depletion rates could be seen at 300 min sampling point between MI-463 and the other inhibitors (MI-472, MI-485 and MI-1900, <0.001 in each case).
Cynomolgus monkey hepatocytes could convert all inhibitors significantly after 120 min treatments (pMI-463, 120 min < 0.001; pMI-472, 120 min = 0.0109; pMI-485, 120 min = 0.0015; pMI-1900, 120 min < 0.001). It was also observed that MI-1900 was depleted at the highest speed, and even from 60 min the decomposition was much higher compared to initial values (Figure 4). No chemical structure-related differences were observed in cynomolgus macaque hepatocytes at 300 min sampling time in the case of MI-463, MI-472, MI-485 and MI-1900 (p values ranged between 0.06 and 1).
Interspecies comparisons of depletion rates of 3-APhA derivatives at 300 min proved that MI-463 decomposed differently in rat hepatocytes compared to that in dog and cynomolgus macaque (p < 0.001); however, there were no differences between dog and cynomolgus macaque samples (p = 0.964). This interspecies variation was also seen in the case of MI-472 (p rat-dog = 0.0106; p rat-monkey < 0.001; p dog-cynomolgus monkey = 0.9713) and MI-1900 (p rat-dog < 0.001; p rat-monkey = 0.4974; p dog-cynomolgus monkey = 0.0128); however, MI-485 did not decompose differently in hepatocytes of the used mammalian species (p rat-dog = 0.1341; p dog-monkey = 1.000; p rat-monkey = 0.2269).
Some interspecies differences in the intrinsic clearance of inhibitors were observed between rat, dog and monkey hepatocytes (Table 2). Rat cells were the most active in biotransformation of all of these compounds (Clint: 31.9–37.8 mL/min/kg), whereas the metabolizing capacity of monkey hepatocytes was the lowest (Clint: 6.6–12.0 mL/min/kg), except for MI-1900, which was eliminated somewhat faster than in dog cells (Clint: 16.1 vs. 10.5 mL/min/kg, respectively). On the basis of hepatic extraction ratios, all four 3-APhA derivatives were designated as low-extraction drugs in monkey cells, while in rat and dog hepatocytes, MI-463 and MI-1900 were low-extraction drugs, and MI-472 and MI-485 were considered to be intermediate-extraction ones. The bioavailability values ranged between 65.7% and 87.9%, predicting low hepatic first-pass effect.

4. Discussion

Infections with SARS-CoV-2 caused the COVID-19 pandemic, leading to a significant global health threat. Emerging findings support that entry of SARS-CoV-2 into host cells also highly depends on serine protease activities and that the suppression of TTSP-related cleavage of viral surface glycoprotein hinders fusion of the virus with host cell membranes.
In vitro studies also confirmed that the 3-APhA-derived inhibitors such as MI-432 and MI-1900 could successfully lower viral titers in Calu-3 cells infected with SARS-CoV-2 [27] in a dose-dependent manner, most likely through TMPRSS2-mediated cleavage occurring at the S2′ site. Inhibitor MI-485, with improved selectivity against thrombin and factor Xa, showed inhibition of matriptase-dependent H9N2 influenza virus replication in Madin-Darby canine kidney (MDCK) cells [15].
Four structurally related 3-APhA-derived inhibitors with an N-terminal dimethoxy-substituted biphenyl-3-sulfonyl group were investigated in this work. In our previous studies, we proved that these 3-APhA-derived inhibitors, such as MI-1900 and its structurally related analogues, did not exert any cytotoxic properties, and it was also found that the investigated inhibitors MI-432, MI-463, MI-482 and MI-1900 did not bind to human serum albumin to a great extent and did not influence the activity of cytochrome P 450 isoenzymes (CYP)1A2, 2C9 and 2C19. However, they could modulate human CYP3A4 function and cynomolgus monkey CYP1A2 [34,35].
The present study examined the pharmacokinetic behavior of 3-APhA-derived protease inhibitors MI-463, MI-472, MI-485 and MI-1900 with potential antiviral effects in primary hepatocytes derived from rat, dog and cynomolgus monkey. Based on previous enzyme kinetic measurements, MI-463 and MI-1900 acted as competitive inhibitors of TMPRSS2 and matriptase with very low, nanomolar Ki values, and MI-485 could also suppress matriptase activity with Ki lower than 5 nM. To characterize the pharmaco-toxicological properties of these compounds, cell-based assays were used without excessive sacrifice of experimental animals, in good agreement with animal welfare policies, and at the same time, extrapolatable in vitro PK data were obtained. Primary hepatocytes are widely accepted as a suitable model for studying in vitro pharmacokinetic properties of newly synthesized or repurposed compounds, including CYP inducibility and biotransformation rates in the liver, compared to tumorigenic cellular systems such as human hepatoma cell line HepG2 [36,37,38]. Hence, they could serve as suitable tools for research related to the hepatic metabolism of MI compounds in rat, dog and cynomolgus monkey.
It was ascertained that in rat hepatocytes, the depletion rates of the investigated MI compounds were higher than in dog or monkey cells and were similar to each other independently of their chemical structural differences. Nevertheless, significant interspecies differences in depletion rates of MI compounds were demonstrated. The piperidine-derived MI-463 with C-terminal cyclohexyl urea moiety showed pronounced stability against biotransformation in dog samples. In contrast, elevated depletion rates of this inhibitor were observed in rat and macaque hepatocytes.
The fact that the intrinsic clearance values for these compounds were lower than the rates of hepatic blood flow means that drug metabolizing capacity of the liver rather than the liver blood flow is expected to primarily determine the elimination of MI compounds. The interspecies differences in hepatic extraction properties of MI derivatives anticipate different systemic exposures of rat, dog and monkey at the same dose level. Furthermore, 3-APhA derivatives were proven to be low- or intermediate-extraction drugs with E < 0.34, which means that slight change in drug metabolizing ability is expected to highly influence the hepatic clearance, whereas variations in hepatic blood flow will poorly affect the clearance of these compounds. The bioavailability values were predicted to be higher than 65% for all 3-APhA derivatives, anticipating that hepatic metabolism will poorly decrease the systemic exposure of these drugs.
Several cellular and animal models have been developed for investigation of virus biology under experimental conditions. Due to their low cost, small size and rapid breeding, mouse models (transgenic mouse expressing hACE2, hACE2-transduced mouse, mouse adapted SARS-CoV-2 strain and hACE2 humanized mouse) and hamster models have been used with some limitations, such as lack of severity of the respiratory symptoms. The cynomolgus macaque model was reported to be suitable for studying the pathology and immune responses of the SARS-CoV-2 infection due to the close genetic similarity of cynomolgus macaques to humans; however, their limited availability, the ethical concerns, and the high cost related to the need for complex husbandry restrict their frequent usage in preclinical COVID-19 investigations [39,40,41,42,43]. It was also reported that liver metabolism in cynomolgus monkeys shows strong resemblance to that in humans, including the CYP expressions and drug uptake of hepatocytes, thereby allowing better prediction of human in vivo intrinsic hepatic clearance [44,45,46]. Cynomolgus monkey hepatocytes could metabolize the piperazine-derived MI-1900 rapidly. That is in good agreement with our previous experimental findings pointing out that previously tested TTSP inhibitors could be metabolized significantly even upon 60 min exposure to mammalian hepatic microsomal CYP isoenzymes [34,35].
The One Health movement includes the screening and development of novel drugs for the prevention and treatment of human and animal diseases (www.onehealthinitiative.com). The findings of these in vitro studies can contribute to the practical application of integrated multidisciplinary approaches in both human medicine and veterinary practice. In vitro studies such as our primary hepatocytes-based depletion assays facilitate the selection of suitable drug candidates and thus provide PK support prior to in vivo studies and indirectly contribute to preventing severe adverse effects in laboratory animals and later in clinical trials. In our further experiments, additional PK parameters for these protease inhibitors are planned to be elucidated, including the protein binding properties and potential modulatory effects on isoenzyme activity, in addition to elucidation of their safety profiles, such as the selectivity of these inhibitors against blood clotting factors.

5. Conclusions

In conclusion, the current study also confirmed that mammalian cellular studies can be applied to predict important PK parameters of 3-APhA-type antiviral drug candidates prior to preclinical animal studies using depletion rates of the inhibitors without preliminary animal sacrifice. Differences in metabolic rates of the selected inhibitors could mainly be seen in dog and cynomolgus monkey hepatocytes, of which the latter possesses very strong similarity in biotransformation pattern to humans, whereas the rat hepatocyte-mediated depletion was accelerated and almost independent of the chemical structures of MI compounds. Among the investigated inhibitors, MI-1900 with piperazine ring was depleted faster in macaque hepatocytes compared to other derivatives having a piperidine unit, suggesting that different types of heterocyclic rings on the C terminal might influence the decomposition rates. These results may suggest that further modification of these inhibitors on the C-terminal ureido moieties could improve their PK parameters, such as Clint, CLH, EH and bioavailability.

Author Contributions

Conceptualization, K.L., K.M. and E.P.-G.; methodology, K.L., K.M. and E.P.-G.; software, K.M.; validation, K.L. and K.M.; formal analysis, K.L., K.M. and E.P.-G.; investigation, K.L. and E.P.-G.; resources, Á.J.; data curation, K.L., K.M. and E.P.-G.; writing—original draft preparation, E.P.-G.; writing—review and editing, K.L., K.M., E.P.-G. and T.S; visualization, K.L.; supervision, E.P.-G.; project administration, K.L., K.M. and E.P.-G.; funding acquisition, Á.J. and E.P.-G., T.S. provided the inhibitors. All authors have read and agreed to the published version of the manuscript.

Funding

Project no. RRF-2.3.1-21-2022-00001 has been implemented with the support provided by the Recovery and Resilience Facility (RRF), financed under the National Recovery Fund budget estimate, RRF-2.3.1-21 funding scheme.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All raw data supporting the results of the present study can be obtained from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Bugge, T.H.; Antalis, T.M.; Wu, Q. Type II Transmembrane Serine Proteases. J. Biol. Chem. 2009, 284, 23177–23181. [Google Scholar] [CrossRef] [Green Version]
  2. Förbs, D.; Thiel, S.; Stella, M.C.; Stürzebecher, A.; Schweinitz, A.; Steinmetzer, T.; Stürzebecher, J.; Uhland, K. In vitro inhibition of matriptase prevents invasive growth of cell lines of prostate and colon carcinoma. Int. J. Oncol. 2005, 27, 1061–1070. [Google Scholar] [CrossRef] [PubMed]
  3. Milner, J.M.; Patel, A.; Davidson, R.K.; Swingler, T.E.; Desilets, A.; Young, D.A.; Kelso, E.B.; Donell, S.T.; Cawston, T.E.; Clark, I.M.; et al. Matriptase is a novel initiator of cartilage matrix degradation in osteoarthritis. Arthritis Rheum. 2010, 62, 1955–1966. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Krijt, J.; Frýdlová, J.; Gurieva, I.; Přikryl, P.; Báječný, M.; Steinbicker, A.U.; Vokurka, M.; Truksa, J. Matriptase-2 and Hemojuvelin in Hepcidin Regulation: In Vivo Immunoblot Studies in Mask Mice. Int. J. Mol. Sci. 2021, 22, 2650. [Google Scholar] [CrossRef] [PubMed]
  5. Hatesuer, B.; Bertram, S.; Mehnert, N.; Bahgat, M.M.; Nelson, P.S.; Pöhlman, S.; Schughart, K. Tmprss2 Is Essential for Influenza H1N1 Virus Pathogenesis in Mice. PLoS Pathog. 2013, 9, e1003774. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Böttcher-Friebertshäuser, E. Membrane-Anchored Serine Proteases: Host Cell Factors in Proteolytic Activation of Viral Glycoproteins. Act. Viruses Host Proteases 2018, 16, 153–203. [Google Scholar] [CrossRef]
  7. Hamilton, B.S.; Gludish, D.W.; Whittaker, G.R. Cleavage Activation of the Human-Adapted Influenza Virus Subtypes by Matriptase Reveals both Subtype and Strain Specificities. J. Virol. 2012, 86, 10579–10586. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Baron, J.; Tarnow, C.; Mayoli-Nüssle, D.; Schilling, E.; Meyer, D.; Hammami, M.; Schwalm, F.; Steinmetzer, T.; Guan, Y.; Garten, W.; et al. Matriptase, HAT, and TMPRSS2 Activate the Hemagglutinin of H9N2 Influenza A Viruses. J. Virol. 2013, 87, 1811–1820. [Google Scholar] [CrossRef] [Green Version]
  9. Kokic, G.; Hillen, H.S.; Tegunov, D.; Dienemann, C.; Seitz, F.; Schmitzova, J.; Farnung, L.; Siewert, A.; Höbartner, C.; Cramer, P. Mechanism of SARS-CoV-2 polymerase stalling by remdesivir. Nat. Commun. 2021, 12, 279. [Google Scholar] [CrossRef]
  10. Beigel, J.H.; Tomashek, K.M.; Dodd, L.E.; Mehta, A.K.; Zingman, B.S.; Kalil, A.C.; Hohmann, E.; Chu, H.Y.; Luetkemeyer, A.; Kline, S.; et al. Remdesivir for the Treatment of COVID-19—Preliminary report. N. Engl. J. Med. 2020, 383, 1813–1826. [Google Scholar] [CrossRef]
  11. Owen, D.R.; Allerton, C.M.N.; Anderson, A.S.; Aschenbrenner, L.; Avery, M.; Berritt, S.; Boras, B.; Cardin, R.D.; Carlo, A.; Coffman, K.J.; et al. An oral SARS-CoV-2 M pro inhibitor clinical candidate for the treatment of COVID-19. Science 2021, 374, 1586–1593. [Google Scholar] [CrossRef]
  12. Wahl, A.; Gralinski, L.E.; Johnson, C.E.; Yao, W.; Kovarova, M.; Dinnon, K.H.; Liu, H.; Madden, V.J.; Krzystek, H.M.; De, C.; et al. SARS-CoV-2 infection is effectively treated and prevented by EIDD-2801. Nature 2021, 591, 451–457. [Google Scholar] [CrossRef] [PubMed]
  13. Tobback, E.; Degroote, S.; Buysse, S.; Delesie, L.; Van Dooren, L.; Vanherrewege, S.; Barbezange, C.; Hutse, V.; Romano, M.; Thomas, I.; et al. Efficacy and safety of camostat mesylate in early COVID-19 disease in an ambulatory setting: A randomized placebo-controlled phase II trial. Int. J. Infect. Dis. 2022, 122, 628–635. [Google Scholar] [CrossRef] [PubMed]
  14. Meyer, D.; Sielaff, F.; Hammami, M.; Böttcher-Friebertshäuser, E.; Garten, W.; Steinmetzer, T. Identification of the first synthetic inhibitors of the type II transmembrane serine protease TMPRSS2 suitable for inhibition of influenza virus activation. Biochem. J. 2013, 452, 331–343. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Pilgram, O.; Keils, A.; Benary, G.E.; Müller, J.; Merkl, S.; Ngaha, S.; Huber, S.; Chevillard, F.; Harbig, A.; Magdolen, V.; et al. Improving the selectivity of 3-amidinophenylalanine-derived matriptase inhibitors. Eur. J. Med. Chem. 2022, 238, 114437. [Google Scholar] [CrossRef] [PubMed]
  16. Colombo, É.; Désilets, A.; Duchêne, D.; Chagnon, F.; Najmanovich, R.; LeDuc, R.; Marsault, E. Design and Synthesis of Potent, Selective Inhibitors of Matriptase. ACS Med. Chem. Lett. 2012, 3, 530–534. [Google Scholar] [CrossRef] [Green Version]
  17. Beaulieu, A.; Gravel, É.; Cloutier, A.; Marois, I.; Colombo, É.; Désilets, A.; Verreault, C.; Leduc, R.; Marsault, É.; Richter, M.V. Matriptase Proteolytically Activates Influenza Virus and Promotes Multicycle Replication in the Human Airway Epithelium. J. Virol. 2013, 87, 4237–4251. [Google Scholar] [CrossRef] [Green Version]
  18. Shapira, T.; Monreal, I.A.; Dion, S.P.; Buchholz, D.W.; Imbiakha, B.; Olmstead, A.D.; Jager, M.; Désilets, A.; Gao, G.; Martins, M.; et al. A TMPRSS2 inhibitor acts as a pan-SARS-CoV-2 prophylactic and therapeutic. Nature 2022, 605, 340–348. [Google Scholar] [CrossRef]
  19. Mahdy, M.A.A.; Younis, W.; Ewaida, Z. An Overview of SARS-CoV-2 and Animal Infection. Front. Vet.—Sci. 2020, 7, 596391. [Google Scholar] [CrossRef]
  20. Sit, T.H.C.; Brackman, C.J.; Ip, S.M.; Tam, K.W.S.; Law, P.Y.T.; To, E.M.W.; Yu, V.Y.T.; Sims, L.D.; Tsang, D.N.C.; Chu, D.K.W.; et al. Infection of dogs with SARS-CoV-2. Nature 2020, 586, 776–778. [Google Scholar] [CrossRef]
  21. Johnston, S.C.; Ricks, K.M.; Jay, A.; Raymond, J.L.; Rossi, F.; Zeng, X.; Scruggs, J.; Dyer, D.; Frick, O.; Koehler, J.W.; et al. Development of a coronavirus disease 2019 nonhuman primate model using airborne exposure. PLoS ONE 2021, 16, e0246366. [Google Scholar] [CrossRef]
  22. Urano, E.; Okamura, T.; Ono, C.; Ueno, S.; Nagata, S.; Kamada, H.; Higuchi, M.; Furukawa, M.; Kamitani, W.; Matsuura, Y.; et al. COVID-19 cynomolgus macaque model reflecting human COVID-19 pathological conditions. Proc. Natl. Acad. Sci. USA 2021, 118, e2104847118. [Google Scholar] [CrossRef] [PubMed]
  23. Zhao, S.; Fan, J.; Liu, E. Animal Models for COVID-19 Therapeutic Development: Where We Are and Where We Need to Go. Front. Microbiol. 2022, 13, 907406. [Google Scholar] [CrossRef] [PubMed]
  24. Hammami, M.; Rühmann, E.; Maurer, E.; Heine, A.; Gütschow, M.; Klebe, G.; Steinmetzer, T. New 3-amidinophenylalanine-derived inhibitors of matriptase. MedChemComm 2012, 3, 807–813. [Google Scholar] [CrossRef]
  25. List, K.H.; Steinmetzer, T.; Hardes, K. The antiviral potential of host protease inhibitors. In Activation of Viruses by Host Proteases, 1st ed.; Böttcher-Friebertshäuser, E., Garten, W., Klenk, H.D., Eds.; Springer: Marburg, Germany, 2018; pp. 279–325. [Google Scholar]
  26. Hoffmann, M.; Kleine-Weber, H.; Schroeder, S.; Krüger, N.; Herrler, T.; Erichsen, S.; Schiergens, T.S.; Herrler, G.; Wu, N.-H.; Nitsche, A.; et al. SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor. Cell 2020, 181, 271–280.e278. [Google Scholar] [CrossRef]
  27. Bestle, D.; Heindl, M.R.; Limburg, H.; Van Lam van, T.; Pilgram, O.; Moulton, H.; Stein, D.A.; Hardes, K.; Eickmann, M.; Dolnik, O.; et al. TMPRSS2 and furin are both essential for proteolytic activation of SARS-CoV-2 in human airway cells. Life Sci. Alliance 2020, 3, e202000786. [Google Scholar] [CrossRef]
  28. Houston, J.B. Utility of in vitro drug metabolism data in predicting in vivo metabolic clearance. Biochem. Pharmacol. 1994, 47, 1469–1479. [Google Scholar] [CrossRef]
  29. Sohlenius-Sternbeck, A.-K. Determination of the hepatocellularity number for human, dog, rabbit, rat and mouse livers from protein concentration measurements. Toxicol. Vitr. 2006, 20, 1582–1586. [Google Scholar] [CrossRef]
  30. Davies, B.; Morris, T. Physiological Parameters in Laboratory Animals and Humans. Pharm. Res. 1993, 10, 1093–1095. [Google Scholar] [CrossRef]
  31. Naritomi, Y.; Terashita, S.; Kagayama, A.; Sugiyama, Y. Utility of Hepatocytes in Predicting Drug Metabolism: Comparison of Hepatic Intrinsic Clearance in Rats and Humans in Vivo and in Vitro. Drug Metab. Dispos. 2003, 31, 580–588. [Google Scholar] [CrossRef] [Green Version]
  32. Chiba, M.; Sugiyama, Y.; Ishii, Y.; Takahashi, H.; Shibata, Y. Prediction of Hepatic Clearance in Humans from Experimental Animals and In Vitro Data. In Drug Metabolizing Enyzmes. Cytochrome P450 and Other Enzymes in Drug Discovery and Development; Lee, J.S., Obach, R.S., Fisher, M.B., Eds.; Fontis Media: Lausanne, Switzerland, 2003; pp. 453–481. [Google Scholar] [CrossRef]
  33. Tóth, K.; Sirok, D.; Kiss, Á.; Mayer, A.; Pátfalusi, M.; Hirka, G.; Monostory, K. Utility of in vitro clearance in primary hepatocyte model for prediction of in vivo hepatic clearance of psychopharmacons. Microchem. J. 2018, 136, 193–199. [Google Scholar] [CrossRef] [Green Version]
  34. Fedor, Z.; Szentkirályi-Tóth, A.; Nagy, G.; Szimrók, Z.; Varga, E.; Pászti, A.; Pászti, Z.; Jerzsele, Á.; Pilgram, O.; Steinmetzer, T.; et al. Interspecies Comparisons of the Effects of Potential Antiviral 3-Amidinophenylalanine Derivatives on Cytochrome P450 1A2 Isoenzyme. Vet.—Sci. 2022, 9, 156. [Google Scholar] [CrossRef] [PubMed]
  35. Pászti-Gere, E.; Szentkirályi, A.; Fedor, Z.; Nagy, G.; Szimrók, Z.; Pászti, Z.; Pászti, A.; Pilgram, O.; Steinmetzer, T.; Bodnárová, S.; et al. In vitro interaction of potential antiviral TMPRSS2 inhibitors with human serum albumin and cytochrome P 450 isoenzymes. Biomed. Pharmacother. 2022, 146, 112513. [Google Scholar] [CrossRef]
  36. Gerets, H.H.J.; Tilmant, K.; Gerin, B.; Chanteux, H.; Depelchin, B.O.; Dhalluin, S.; Atienzar, F.A. Characterization of primary human hepatocytes, HepG2 cells, and HepaRG cells at the mRNA level and CYP activity in response to inducers and their predictivity for the detection of human hepatotoxins. Cell Biol. Toxicol. 2012, 28, 69–87. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Wilkening, S.; Stahl, F.; Bader, A. Comparison of primary human hepatocytes and hepatoma cell line hepg2 with regard to their biotransformation properties. Drug Metab. Dispos. 2003, 31, 1035–1042. [Google Scholar] [CrossRef] [PubMed]
  38. Berger, B.; Donzelli, M.; Maseneni, S.; Boess, F.; Roth, A.; Krähenbühl, S.; Haschke, M. Comparison of Liver Cell Models Using the Basel Phenotyping Cocktail. Front. Pharmacol. 2016, 7, 443. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Fan, C.; Wu, Y.; Rui, X.; Yang, Y.; Ling, C.; Liu, S.; Liu, S.; Wang, Y. Animal models for COVID-19: Advances, gaps and perspectives. Signal Transduct. Target. Ther. 2022, 7, 220. [Google Scholar] [CrossRef] [PubMed]
  40. Caldera-Crespo, L.A.; Paidas, M.J.; Roy, S.; Schulman, C.I.; Kenyon, N.S.; Daunert, S.; Jayakumar, A.R. Experimental Models of COVID-19. Front. Cell Infect. Microbiol. 2022, 11, 792584. [Google Scholar] [CrossRef] [PubMed]
  41. Lee, C.-Y.; Lowen, A.C. Animal models for SARS-CoV-2. Curr. Opin. Virol. 2021, 48, 73–81. [Google Scholar] [CrossRef]
  42. Choudhary, S.; Kanevsky, I.; Tomlinson, L. Animal models for studying COVID-19, prevention, and therapy: Pathology and disease phenotypes. Vet.—Pathol. 2022, 59, 516–527. [Google Scholar] [CrossRef]
  43. Rosa, R.B.; Dantas, W.M.; Do Nascimento, J.C.F.; Da Silva, M.V.; De Oliveira, R.N.; Pena, L.J. In Vitro and In Vivo Models for Studying SARS-CoV-2, the Etiological Agent Responsible for COVID-19 Pandemic. Viruses 2021, 13, 379. [Google Scholar] [CrossRef] [PubMed]
  44. Iwasaki, K.; Uno, Y. Cynomolgus monkey CYPs: A comparison with human CYPs. Xenobiotica 2009, 39, 578–581. [Google Scholar] [CrossRef] [PubMed]
  45. De Bruyn, T.; Ufuk, A.; Cantrill, C.; Kosa, R.E.; Bi, Y.-A.; Niosi, M.; Modi, S.; Rodrigues, A.D.; Tremaine, L.M.; Varma, M.V.S.; et al. Predicting Human Clearance of Organic Anion Transporting Polypeptide Substrates Using Cynomolgus Monkey: In Vitro–In Vivo Scaling of Hepatic Uptake Clearance. Drug Metab. Dispos. 2018, 46, 989–1000. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Shen, H.; Yang, Z.; Rodrigues, A.D. Cynomolgus Monkey as an Emerging Animal Model to Study Drug Transporters: In Vitro, In Vivo, In Vitro-to-In Vivo Translation. Drug Metab. Dispos. 2022, 50, 299–319. [Google Scholar] [CrossRef] [PubMed]
Biomedicines 11 00682 g001 550
Figure 1. Structures of the matriptase/TMPRSS2 inhibitors applied in this study.
Figure 1. Structures of the matriptase/TMPRSS2 inhibitors applied in this study.
Biomedicines 11 00682 g001
Biomedicines 11 00682 g002 550
Figure 2. Time course of depletion of inhibitors MI-463, MI-472, MI-485 and MI-1900 in rat primary hepatocytes (5 × 105/mL) using 0, 30, 60, 120, 180 and 300 min sampling times. Relative percentage values of unchanged parent compounds were calculated by dividing the concentrations of the inhibitor measured at a given time point by mean starting MI concentration in the treatment solution (mean ± SD; n = 3–9/group). Asterisks indicate significant differences compared to starting values of the indicated inhibitors in parentheses at 300 min (*** p < 0.001).
Figure 2. Time course of depletion of inhibitors MI-463, MI-472, MI-485 and MI-1900 in rat primary hepatocytes (5 × 105/mL) using 0, 30, 60, 120, 180 and 300 min sampling times. Relative percentage values of unchanged parent compounds were calculated by dividing the concentrations of the inhibitor measured at a given time point by mean starting MI concentration in the treatment solution (mean ± SD; n = 3–9/group). Asterisks indicate significant differences compared to starting values of the indicated inhibitors in parentheses at 300 min (*** p < 0.001).
Biomedicines 11 00682 g002
Biomedicines 11 00682 g003 550
Figure 3. Time course of depletion of inhibitors MI-463, MI-472, MI-485 and MI-1900 in dog primary hepatocytes (5 × 105/mL) using 0, 30, 60, 120, 180 and 300 min sampling times. Relative percentage values of unchanged parent compounds were calculated by dividing the concentrations of the inhibitor measured at given time point by mean starting MI concentration in the treatment solution (mean ± SD; n = 3–9/group). Asterisks indicate significant differences compared to starting values of the indicated inhibitors in parentheses at 300 min (** p < 0.01, *** p < 0.001).
Figure 3. Time course of depletion of inhibitors MI-463, MI-472, MI-485 and MI-1900 in dog primary hepatocytes (5 × 105/mL) using 0, 30, 60, 120, 180 and 300 min sampling times. Relative percentage values of unchanged parent compounds were calculated by dividing the concentrations of the inhibitor measured at given time point by mean starting MI concentration in the treatment solution (mean ± SD; n = 3–9/group). Asterisks indicate significant differences compared to starting values of the indicated inhibitors in parentheses at 300 min (** p < 0.01, *** p < 0.001).
Biomedicines 11 00682 g003
Biomedicines 11 00682 g004 550
Figure 4. Time course of depletion of inhibitors MI-463, MI-472, MI-485 and MI-1900 in macaque primary hepatocytes (5 × 105/mL) using 0, 30, 60, 120, 180 and 300 min sampling times. Relative percentage values of unchanged parent compounds were calculated by dividing the concentrations of the inhibitor measured at given time point by mean starting MI concentration in the treatment solution (mean ± SD; n = 3–9/group). Asterisks indicate significant differences compared to starting values of the indicated inhibitors in parentheses (*** p < 0.001).
Figure 4. Time course of depletion of inhibitors MI-463, MI-472, MI-485 and MI-1900 in macaque primary hepatocytes (5 × 105/mL) using 0, 30, 60, 120, 180 and 300 min sampling times. Relative percentage values of unchanged parent compounds were calculated by dividing the concentrations of the inhibitor measured at given time point by mean starting MI concentration in the treatment solution (mean ± SD; n = 3–9/group). Asterisks indicate significant differences compared to starting values of the indicated inhibitors in parentheses (*** p < 0.001).
Biomedicines 11 00682 g004
Table
Table 1. Physiological parameters for calculation of clearance values [30].
Table 1. Physiological parameters for calculation of clearance values [30].
ParameterRatDogMonkey
Number of hepatocytes in liver (×106 cells/g liver)117215114
Liver weight (g)10320150
Body weight (kg)0.25105
Liver blood flow (mL/min/kg)55.230.943.6
Plasma/blood ratio0.580.570.61
Table
Table 2. Summary of rat, dog and macaque pharmacokinetic parameters of inhibitors. F means bioavailability. Cl int = intrinsic clearance, ClH = hepatic clearance values, EH = hepatic extraction ratio.
Table 2. Summary of rat, dog and macaque pharmacokinetic parameters of inhibitors. F means bioavailability. Cl int = intrinsic clearance, ClH = hepatic clearance values, EH = hepatic extraction ratio.
Clint (mL/min/kg)Predicted Cl CLHExtraction Ratio EHF (%)Classification
MI-463
Rat31.8715.970.28971.07Low-extraction drug
Dog19.4879.250.29970.06Low-extraction drug
Macaque6.585.270.12187.91Low-extraction drug
MI-472
Rat35.5816.850.30569.47Intermediate-extraction drug
Dog26.6910.610.34365.66Intermediate-extraction drug
Macaque12.008.270.19087.91Low-extraction drug
MI-485
Rat37.7817.330.31468.60Intermediate-extraction drug
Dog22.759.930.32167.87Intermediate-extraction drug
Macaque11.407.980.18381.69Low-extraction drug
MI-1900
Rat33.7616.430.29870.23Low-extraction drug
Dog10.466.560.21278.76Low-extraction drug
Macaque16.1010.030.23076.99Low-extraction drug
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lányi, K.; Monostory, K.; Steinmetzer, T.; Jerzsele, Á.; Pászti-Gere, E. In Vitro Pharmacokinetic Behavior of Antiviral 3-Amidinophenylalanine Derivatives in Rat, Dog and Monkey Hepatocytes. Biomedicines 2023, 11, 682. https://doi.org/10.3390/biomedicines11030682

AMA Style

Lányi K, Monostory K, Steinmetzer T, Jerzsele Á, Pászti-Gere E. In Vitro Pharmacokinetic Behavior of Antiviral 3-Amidinophenylalanine Derivatives in Rat, Dog and Monkey Hepatocytes. Biomedicines. 2023; 11(3):682. https://doi.org/10.3390/biomedicines11030682

Chicago/Turabian Style

Lányi, Katalin, Katalin Monostory, Torsten Steinmetzer, Ákos Jerzsele, and Erzsébet Pászti-Gere. 2023. "In Vitro Pharmacokinetic Behavior of Antiviral 3-Amidinophenylalanine Derivatives in Rat, Dog and Monkey Hepatocytes" Biomedicines 11, no. 3: 682. https://doi.org/10.3390/biomedicines11030682

APA Style

Lányi, K., Monostory, K., Steinmetzer, T., Jerzsele, Á., & Pászti-Gere, E. (2023). In Vitro Pharmacokinetic Behavior of Antiviral 3-Amidinophenylalanine Derivatives in Rat, Dog and Monkey Hepatocytes. Biomedicines, 11(3), 682. https://doi.org/10.3390/biomedicines11030682

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

Article Metrics

Back to TopTop