key: cord-0874121-20t2zpo3
authors: Jeong, Yoo-Seong; Jusko, William J.
title: Meta-Assessment of Metformin Absorption and Disposition Pharmacokinetics in Nine Species
date: 2021-06-07
journal: Pharmaceuticals (Basel)
DOI: 10.3390/ph14060545
sha: 3e00258baddd19d79447e247d2464bbd0d0b3c1d
doc_id: 874121
cord_uid: 20t2zpo3

The objective of this study was to systematically assess literature datasets and quantitatively analyze metformin PK in plasma and some tissues of nine species. The pharmacokinetic (PK) parameters and profiles of metformin in nine species were collected from the literature. Based on a simple allometric scaling, the systemic clearances ([Formula: see text]) of metformin in these species highly correlate with body weight ([Formula: see text]) (R(2) = 0.85) and are comparable to renal plasma flow in most species except for rabbit and cat. Reported volumes of distribution ([Formula: see text]) varied appreciably (0.32 to 10.1 L/kg) among species. Using the physiological and anatomical variables for each species, a minimal physiologically based pharmacokinetic (mPBPK) model consisting of blood and two tissue compartments (Tissues 1 and 2) was used for modeling metformin PK in the nine species. Permeability-limited distribution (low [Formula: see text] and [Formula: see text]) and a single tissue-to-plasma partition coefficient ([Formula: see text]) value for Tissues 1 and 2 were applied in the joint mPBPK fitting. Nonlinear regression analysis for common tissue distribution parameters along with species-specific [Formula: see text] values reasonably captured the plasma PK profiles of metformin across most species, except for rat and horse with later time deviations. In separate fittings of the mPBPK model to each species, Tissue 2 was considered as slowly-equilibrating compartment consisting of muscle and skin based on in silico calculations of the mean transit times through tissues. The well-fitted mPBPK model parameters for absorption and disposition PK of metformin for each species were compared with in vitro/in vivo results found in the literature with regard to the physiological details and physicochemical properties of metformin. Bioavailability and absorption rates decreased with the increased [Formula: see text] among the species. Tissues such as muscle dominate metformin distribution with low permeability and partitioning while actual tissue concentrations found in rats and mice show likely transporter-mediated uptake in liver, kidney, and gastrointestinal tissues. Metformin has diverse pharmacologic actions, and this assessment revealed allometric relationships in its absorption and renal clearance but considerable variability in actual and modeled tissue distribution probably caused by transporter differences.

Metformin is a biguanide compound that was derived from Galega officinalis (French lilac) having an abundant content of the hypoglycemic substance guanidine [1, 2] . Although the specific compound metformin was first synthesized in 1920s [3, 4] and its use was later replaced by the breakthrough discovery of insulin in the same era [5] , experience with guanidine and diguanides for treating diabetes led to development of phenformin [6] , buformin [7] , and metformin (Glucophage was named by Jean Sterne) in the 1950s [8, 9] . However, after phenformin and buformin were withdrawn from the market in most countries due to their high risks of lactic acidosis [10] , metformin has been recommended as the first-line oral medication for treatment of type 2 diabetes by the American Diabetes Association (ADA) and the European Association for the Study of Diabetes (EASD) since 2009 [11] . In addition to its antihyperglycemic effect, the United Kingdom Prospective Diabetes Study (UKPDS) reported long-term metabolic effects and reduced cardiovascular risk with metformin [12] . Metformin is increasingly recognized as a potential anticancer agent due to a reduced cancer incidence in diabetic patients treated with the drug [13] [14] [15] , and recently, patients taking metformin were associated with a reduced risk of COVID-19related mortality [16, 17] . The history of metformin from its discovery in the 18th century to clinical uses for over 60 years is well-summarized in the literature [18] , and the drug is now widely used with more than 500 million prescriptions in United States with growth trends since 2000 [19] .

Although its mechanism of action remains unclear even with the long-term experiences of clinical use, the blood glucose-lowering effect of metformin is attributed to several tissues. The major mode of its pharmacological action is decreased hepatic gluconeogenesis by modulation of mitochondrial enzymes and by the antagonistic regulation of the glucagon signaling pathway [20, 21] . In vitro studies in rodents showed that metformin improved insulin sensitivity of muscle by increasing the expression and activity of insulin receptors, leading to increased insulin-dependent glucose uptake in cells [22] [23] [24] . In addition, increased glucose uptake into the intestinal mucosa [25] and reduced intestinal glucose absorption [26] were observed in metformin-treated rats. A compartmental model accounting for amounts of metformin in plasma, liver, muscle, and gastrointestinal (GI) tract provided an insightful quantitative assessment of glucose-lowering effects that differ with route of administration in diabetic rats [27] . An alternative model for glucose production inhibited by metformin distributed to a biophase compartment described its blood-glucose lowering effect in diabetic rats at different doses [28] . Pharmacodynamic studies of metformin in rat or man describing its effects over a range of doses are few.

Oral absorption of metformin from the immediate-release dosage forms is incomplete in man, with an estimated population mean of 55% for bioavailability [29] . Metformin, a strong base with pKa = 11.8 [30] , is present mostly as a cationic form at biologically relevant pH ranges and has a high water solubility larger than 100 mg/mL [31] . Its hydrophilicity (log P = −1.43) [32] is associated with the low intestinal [33] and cell membrane permeability [34] , which is recognized as a primary limiting step for metformin oral absorption and accounts for its insignificant metabolism in the liver. Metformin is mainly eliminated by renal excretion mediated by transporters, including organic cation transporters (OCT1 and OCT2) and multidrug and toxin extrusion proteins (MATE1 and MATE2) [35] . Dosage adjustments are thus necessary for patients with severe renal dysfunction [36] [37] [38] . In healthy subjects, a longer half-life of metformin in erythrocytes (23.4 h) was observed compared to plasma (2.7 h) [39] , but this may be attributed to the time frames assessed. Similarly, a longer half-life of metformin was also observed in following urinary excretion rates (9.3 to 19.0 h) [37, 40] . Metformin concentrations in tissues have been measured using radiolabeled compounds (i.e., 14 C-label for mouse and 11 C-label for man), as well as use of PET (Positron Emission Tomography) scanning [41, 42] . However, a limited number of time points/tissues, potential involvement of its metabolite(s), and indistinguishability of intra-and extracellular and the intravascular fractions using PET are limiting factors. The drug exhibits either poly-exponential or power function disposition with uncertainties regarding its distribution kinetics, the terminal half-life, and recoveries in urine.

When only pharmacokinetic (PK) profiles in plasma are available, minimal physiologically based pharmacokinetic (mPBPK) models allow a rational interpretation of plasma and expected tissue PK by the incorporation of physiological and anatomical elements into PK analyses [43] . In line with this, allometric scaling is often utilized to relate drug disposition processes to body size [44] . The renal clearance of metformin was reported to have a simple allometric relationship with body weight (BW) with a common exponent of 0.75 [45] . However, traditional allometric scaling assessing PK parameters in relation to BW (i.e., Y = a·BW b ) among species is considered empirical [46] . Thus, integrating allometric features into mPBPK modeling is useful for interspecies scaling of PK characteristics. Previously, mPBPK models were successfully applied for across-species fitting and scaling of full PK profiles of moxifloxacin [43] , dexamethasone [47] , and several monoclonal antibodies [48] . Despite abundant PK data for metformin in various species from its wide use in clinical/veterinary and experimental settings for over 60 years, there has been only limited quantitative analysis assessing interspecies relationships regarding PK properties of metformin [49, 50] .

In reviewing the extensive literature on metformin PK, we first assessed its clearance (CL) in relation to BW and renal plasma flow (RPF) across species using traditional allometry. We then employed mPBPK models to quantitatively analyze metformin PK in all nine species. Joint fittings using mPBPK models that adopted the physiological and anatomical properties of each species were applied to obtain common PK parameters with conserved features across the species seeking a minimal number of parameters that could describe all data. As a third step, a mPBPK model structure based on in silico-predicted tissue constituents was used to further determine and compare the species-dependent PK parameters related to absorption and disposition of metformin for each of the nine species. This systematic review of literature datasets and quantitative analyses of metformin PK in plasma and some tissues of all available species should be useful for further evaluation of the pleiotropic pharmacological effects of metformin.

The PK parameters of metformin after the intravenous (IV) and oral doses found in nine species are summarized in Table 1 . The raw data for metformin PK in dog were provided by Johnston et al. [51] , and all others were digitized for the current analysis. When PK parameters in any species were found to be variable among multiple studies (e.g., CL values ranging from 12.7 to 81.7 mL/min/kg in mouse and from 6.0 to 42.8 mL/min/kg in rat; Supplementary Table S1), PK data were selected on the basis of the richness of sampling time points, duration of sampling periods, details of study design, and alignment of CL values with basic allometric scaling. It is noted that, due to variable PK profiles depending on the oral doses in rats in previous literature that included the same IV doses (50, 100, and 200 mg/kg) [52] , another oral PK profile at a similar experimental setting [53] was considered for the individualized analysis. As shown in Figure 1 , the reported CL values for metformin showed a simple allometric relationship with BW of the nine species (R 2 = 0.85). Despite this reasonable correlation, the systemic CL values of metformin in rabbit and cat deviate from the fitted relationship (with a factor of 6. 33 and 4.49) . Except for these two species, CL values were comparable with RPF, indicating a high extraction ratio of metformin in the kidney, and were consistent with evidence that metformin is mainly eliminated via renal excretion [40, 73] . ) values (red) and body weights of 9 species (Table 3, see below) . Except for rabbit and cat, the systemic clearance of metformin was comparable to . Lines were fitted to the indicated allometric relationships ( = • , where is the body weight).

The metformin PK data were first jointly analyzed using the common parameters across the species adopted in the mPBPK model ( , , , , , and ), along with the individual optimized for each species. Figure 2 shows the mPBPK model structure used for metformin PK and the differential equations for the current analysis are described in Methods section. As shown in Figure 3 , the plasma concentrations across species ranged from about 1000 to 0.001 µg/mL, with rapid and slower phases often extending to 72 h. A wide range of IV doses were given (Table 1) , and the plasma concentrations were not normalized for the dose. The mPBPK model anticipates triexponential disposition profiles that appear in all species. Consistent with its hydrophilicity (log P = −1.43) [32] and low passive permeability [74] , the assumption of perfusion-limited distribution (i.e., + = 1) did not adequately capture the PK data in our preliminary mPBPK modeling. Thus, it was reasoned that the sum of the two fractional distribution parameters would be less than 1 and conserved across the species in the joint fitting approach (Equation (7)). In addition, a consideration of different values for each tissue compartment ( and ) did not work well, leading to unreliable estimates (extremely high CV%) for parameters , , and . As a result, the PK profiles were reasonably captured by the current mPBPK model, except for rat and horse, with later time deviations. However, the late PK data in horses were difficult to digitize and the mouse data lacked a terminal phase. The three IV dose profiles for man came from three different studies (Table 1 ) and all were very well-captured. Despite these divergent fittings, CV% values for the distribution parameters of the joint mPBPK model were all less than 10% ( Table 2 ). The fitted for each species from the joint fitting are listed in Table 3 (all CV% less than 15%). Based on the single assumption, the physical volume of Tissue 1 (0.16• ) was smaller than Tissue 2 (0.77• ). The , is 0.0457, which implies that 4.57% of the drug is subjected to extraction from the blood during each passage through tissues, as the total distributional clearance is • , (Equation (2)). The distributional clearance to Tissue 1 or 2 can be calculated by multiplying either or by the cardiac output ( ) of each species, indicating that Tissue 1, which has a shorter , is more rapidly equilibrated than Tissue 2 ( = /( ) and = /( )). The overall value is 1.29, which is midrange to many values found by noncompartmental analysis (NCA) (0.321 to 2.26 L/kg, except for dogs), as listed in Table 1 . The value for dogs calculated by NCA using the averaged concentration-time data was 10.1 L/kg, which (Table 3 , see below). Except for rabbit and cat, the systemic clearance of metformin was comparable to RPF. Lines were fitted to the indicated allometric relationships (y = a·BW b , where BW is the body weight).

The metformin PK data were first jointly analyzed using the common parameters across the species adopted in the mPBPK model (K p , f d,total , f d1 , f d2 , and f t ), along with the individual optimized CL for each species. Figure 2 shows the mPBPK model structure used for metformin PK and the differential equations for the current analysis are described in Methods section. As shown in Figure 3 , the plasma concentrations across species ranged from about 1000 to 0.001 µg/mL, with rapid and slower phases often extending to 72 h. A wide range of IV doses were given (Table 1) , and the plasma concentrations were not normalized for the dose. The mPBPK model anticipates triexponential disposition profiles that appear in all species. Consistent with its hydrophilicity (log P = −1.43) [32] and low passive permeability [74] , the assumption of perfusion-limited distribution (i.e., f d1 + f d2 = 1) did not adequately capture the PK data in our preliminary mPBPK modeling. Thus, it was reasoned that the sum of the two fractional distribution parameters would be less than 1 and conserved across the species in the joint fitting approach (Equation (7)). In addition, a consideration of different K p values for each tissue compartment (K p1 and K p2 ) did not work well, leading to unreliable estimates (extremely high CV%) for parameters K p1 , K p2 , and f t . As a result, the PK profiles were reasonably captured by the current mPBPK model, except for rat and horse, with later time deviations. However, the late PK data in horses were difficult to digitize and the mouse data lacked a terminal phase. The three IV dose profiles for man came from three different studies (Table 1 ) and all were very well-captured. Despite these divergent fittings, CV% values for the distribution parameters of the joint mPBPK model were all less than 10% ( Table 2 ). The fitted CL for each species from the joint fitting are listed in Table 3 (all CV% less than 15%). Based on the single K p assumption, the physical volume of Tissue 1 (0.16·BW) was smaller than Tissue 2 (0.77·BW). The f d,total is 0.0457, which implies that 4.57% of the drug is subjected to extraction from the blood during each passage through tissues, as the total distributional clearance CL D is Q CO · f d,total (Equation (2)). The distributional clearance to Tissue 1 or 2 can be calculated by multiplying either f d1 or f d2 by the cardiac output (Q CO ) of each species, indicating that Tissue 1, which has a shorter MTT, is more rapidly equilibrated than Tissue 2 (MTT 1 = V 1 K p /(Q CO f d1 R b ) and MTT 2 = V 2 K p /(Q CO f d2 R b )). The overall K p value is 1.29, which is midrange to many V SS values found by noncompartmental analysis (NCA) (0.321 to 2.26 L/kg, except for dogs), as listed in Table 1 . The V SS value for dogs calculated by NCA using the averaged concentration-time data was 10.1 L/kg, which was accompanied by the underestimation of the systemic CL for dogs in the joint fitting (9.01 mL/min/kg, compared to the reported CL of 24.1 mL/min/kg) within a restricted range of distribution parameters. Due to the model prediction showing the terminal phase with a longer half-life than the limited observations in mouse (Figure 3 ), divergent NCA and jointly fitted CL values were observed (Table 3) . was accompanied by the underestimation of the systemic for dogs in the joint fitting (9.01 mL/min/kg, compared to the reported of 24.1 mL/min/kg) within a restricted range of distribution parameters. Due to the model prediction showing the terminal phase with a longer half-life than the limited observations in mouse ( Figure 3 ), divergent NCA and jointly fitted values were observed (Table 3 ). A single K p value (K p1 = K p2 ) was used in the joint fitting, while the cardiac output (Q CO ) was divided into tissue blood flow to each compartment (Q 1 + Q 2 = Q CO ) in the separate fittings. Symbols are defined in the text and in Table 2 and Table 5 (see below) . 

In order to investigate potential differences in PK variables for metformin a species, we carried out further model fittings of both IV and oral PK profiles se for each of the nine species. Having multiple PK profiles in each species adds gre tainty in resolving model parameters. As shown in Figure 4 , all PK profiles we captured. Tissue lumping theories [75] [76] [77] [78] state that the major organ(s) compo rapidly and slowly equilibrating compartments (Tissues 1 and 2) that, in turn, go exponential components/shapes of the systemic PK may not be identifiable witho 

In order to investigate potential differences in PK variables for metformin across the species, we carried out further model fittings of both IV and oral PK profiles separately for each of the nine species. Having multiple PK profiles in each species adds greater certainty in resolving model parameters. As shown in Figure 4 , all PK profiles were well-captured. Tissue lumping theories [75] [76] [77] [78] state that the major organ(s) composing the rapidly and slowly equilibrating compartments (Tissues 1 and 2) that, in turn, govern the exponential components/shapes of the systemic PK may not be identifiable without information on tissue PK. However, muscle was included in Tissue 2, since the f t estimate of 0.172 from the joint fitting (Table 2 ) was small. Using S e f f values in a WB-PBPK structure for rats [79] , the tissue permeability coefficient P was estimated to be 0.0102 × 10 −6 cm/s with a CV% of 29.4% (R 2 = 0.98). As a result, in silico calculations ( Table 5 ) consistently indicate that skin may be the kinetically largest tissue with the longest MTT for a variable range of permeability coefficients. A rule of thumb for lumping kinetically large tissues suggests that tissues within a two-fold difference in MTT may be lumped into a single compartment [75] , which was satisfied by MTT values calculated from the in vitro PAMPA P for muscle and skin ( Table 5 ). The brain, the second-largest tissue based on the case of fitted P that satisfies the acceptance condition with a factor of 1.62, was not included in Tissue 2, considering the variable K p values for rat brains (cf. in vivo K p ranged from 0.2 to 1.48 ( Table 6 , see below)), compared to an in silico value of 0.86 (Table 5) , where the involvement of the efflux transport of metformin by P-glycoprotein [80] may reduce K p and, thus, MTT. Based on these considerations, the mPBPK model in our separate fitting approach utilizes the blood flow rate and anatomical volume of skin and muscle to define Tissue 2, with the fractional tissue volume ( f t = V 2 /(BW − V B )) ranging from 0.472 to 0.678 and blood flow (Q 2 /Q CO ) from 0.217 to 0.563 (Supplementary Table S2 ). range of permeability coefficients. A rule of thumb for lumping kinetically large tissues suggests that tissues within a two-fold difference in may be lumped into a single compartment [75] , which was satisfied by values calculated from the in vitro PAMPA for muscle and skin ( Table 4 ). The brain, the second-largest tissue based on the case of fitted that satisfies the acceptance condition with a factor of 1.62, was not included in Tissue 2, considering the variable values for rat brains (cf. in vivo ranged from 0.2 to 1.48 ( Table 6 , see below)), compared to an in silico value of 0.86 (Table  4) , where the involvement of the efflux transport of metformin by P-glycoprotein [80] may reduce and, thus, . Based on these considerations, the mPBPK model in our separate fitting approach utilizes the blood flow rate and anatomical volume of skin and muscle to define Tissue 2, with the fractional tissue volume ( = /( − )) ranging from 0.472 to 0.678 and blood flow ( / ) from 0.217 to 0.563 (Supplementary Table S2 ). Table 5 . Figure 4 . Fitting of intravenous and oral pharmacokinetics of metformin individually for 9 species to the mPBPK model. The slowly-equilibrating compartment (Tissue 2) was assumed to consist of skin and muscle (i.e.,

. Fitted parameters are listed in Table 4 .

The IV and oral PK profiles were reasonably captured ( Figure 4 ) in all species with the fitted parameters listed in Table 4 . It is noted that the 48-h time point in the oral PK profile in minipig was excluded in the current analysis as an obvious outlier as inclusion produced parameter estimates with an unreliable CV% and an extremely shallow terminal phase. It should be noted that the terminal phase slope in flip-flop kinetics is dependent on the total escaping rate from the absorption site (k a ), not F·k a . The half-lives of k a in mouse, rat, and man are 1.85, 2.82, and 1.74 h (Table 4) , which are longer than model-predicted terminal phase half-lives, leading to flip-flop kinetics in these species ( Figure 4 ). As shown in Table 4 , the f d1 values were smaller than f d2 for each species, while K p1 was smaller than K p2 except for rat, indicating that MTT 2 is longer than MTT 1 . Except for f d1 for mouse and K p2 for horse, the largest CV% of the primary estimates was 93.5%, but the others were quite small (Table 4 ), indicating reasonable model performances. The estimated PK parameters related to oral absorption of metformin from the separate fittings are shown in Figure 5 . The model-fitted bioavailability F for the nine species showed a decreasing trend with increasing BW. While the k a values were statistically independent of BW (i.e., allometric exponent with statistically insignificant difference from 0), the effective absorption rate constant (F·k a ) was highly correlated with BW with an allometric exponent of 0.257 (R 2 = 0.54). This is consistent with a common allometric exponent of 0.25 found for the rate constants [44] . According to the Biopharmaceutics Classification System (BCS) [83] , metformin is considered a Class III drug based on its high water solubility as a hydrochloride salt (i.e., >100 mg/mL) [31] and low intestinal effective permeability (e.g., P e f f ,rat = 0.27 × 10 −4 cm/s) [33] . Based on the assumption of instantaneous solubilization of metformin into the fluid volumes of the stomach and small intestine (V Lumen values summarized in Supplementary Table S2) , Dose/V Lumen was estimated to be far below 100 mg/mL for the maximum dose of 175 mg/kg in mouse (i.e., 10.9 mg/mL). Therefore, expecting that intestinal permeability is the major determinant governing metformin absorption, the effective permeability (P e f f ) was estimated by considering the GI tract as a cylinder (i.e., F·k a = 2P e f f R , where R is the gut radius). These were particularly higher in minipig, man, and dog than other species.

Part of this literature review included the assessment of actual tissue-to-plasma partition coefficients of metformin. All available measurements are listed in Supplementary  Table S3 , with means and ranges summarized in Table 6 . In vivo data could be mostly found for rats and mice, except for protein binding and erythrocyte distribution. The drug appeared to sequester most highly in the GI tract, kidney, and liver, with tissue-to-plasma ratios of 3 or greater. Large values are expected in kidney owing to the robust transport mechanism and probable inclusion of tubular fluid with excreted drug in the measurements. The high ratios in the liver and kidney were not anticipated when using the tissue K p prediction methods (Table 5) for metformin in rats, as these methods do not account for the role of transporters. Muscle is the dominant tissue in the body by weight and has low tissue-to-plasma ratios of 0.6 to 1.0 that are in accordance with the predicted value (0.79). This likely explains the overall low K p and V SS values found when fitting mouse and rat data using the mPBPK model (Table 4 ). Using the in vivo tissue K p values listed in Table 6 , along with the muscle K p assumed to be the same for the unavailable tissues, the calculated V SS (V SS = V B + ∑ V T,i K p,i ) [84] was 2140 mL/kg for mouse and 777 mL/kg for rat (Supplementary Table S4 ), comparable to the observed V SS of 1840 mL/kg for mouse and 586-693 mL/kg for rat (Table 1) . maceuticals 2021, 14, x FOR PEER REVIEW Figure 5 . Fitted pharmacokinetic parameters related to oral absorption of metformin in re body weights, including bioavailability ( ), absorption rate constants ( and · ), and i effective permeability ( ). The fitted parameters for each species are listed in Table 5 .

Part of this literature review included the assessment of actual tissue-to-plasm tition coefficients of metformin. All available measurements are listed in Supplem Table S3 , with means and ranges summarized in Table 6 . In vivo data could be Table 4 . Supplementary Table S3 . Data from normal (nondiabetic) species were included after an outlier test based on the interquartile range rule. ** For data presented as a range from a single study, the average value of the maximum and minimum was adopted.

Metformin is widely used as a first-line agent for the treatment of type 2 diabetes. Its antihyperglycemic effect is mainly achieved by decreasing hepatic gluconeogenesis and glycogenolysis and by increasing the insulin sensitivity of muscle tissue [85, 86] . Its advantages include insignificant risks of clinical hypoglycemia, which is a serious side effect caused by antidiabetic agents such as insulin and sulfonylureas [87] . In addition, metformin reduces plasma triglycerides and low-density lipoproteins (LDL) [88] [89] [90] [91] , which may decrease the risk of cardiovascular disease. The drug has received considerable attention due to a reduced cancer incidence in diabetic patients treated with metformin [13] [14] [15] . Inhibition of mitochondrial complex I activity has been proposed to be the mechanism of the antitumor effect [92, 93] that could be mostly observed at supra-pharmacological concentrations of metformin (i.e., IC50 of 19-66 mM) [94, 95] , while the blood glucose lowering effect was mediated by increased mitochondrial respiration and fission at clinically relevant concentrations (e.g., 75 µM) [96] . Other beneficial effects of metformin reported are weight reduction [97] , delay of aging [98, 99] , and reduced risk of COVID-19-related mortality [16, 17] . In contrast, metformin use may cause adverse effects such as lactic acidosis, hepatotoxicity, acute pancreatitis, and vitamin B12 deficiency [100] . Due to these pleiotropic pharmacological and toxicological effects of metformin, this modeling analysis was pursued to gain improved understanding of its PK and tissue distribution properties among various species. Metformin PK were found in nine species owing to its wide use in the treatment of diabetes in clinical and veterinary settings for over 60 years. The PubMed search (last accessed 1 May 2021) provided 1733 citations since 1969 using the key words 'metformin' and 'pharmacokinetics' and 717 citations since 2010 with the topic of 'the pharmacokinetics of metformin'.

Metformin is well-known to be mostly or entirely eliminated via urinary excretion [40, 73] . Though there has been some evidence of incomplete urinary recovery of metformin (e.g., approximately 80%) after IV doses in man [36, 37] and rats [52, 53, 101] , the systemic CL was comparable with RPF in most species, except for rabbits and cats (Figure 1 ). Since the glomerular filtration rate (essentially f up ·GFR) is much smaller than RPF (Supplementary  Table S2 ), the concordance between CL and RPF implies significant active secretion of metformin by the kidney. While hyperthermia (e.g., over 41 • C) is associated with decreased GFR and RPF [102] , the RPF of relatively hyperthermic rabbits (39.1 • C) [103] was aligned with the fitted trendline between BW and RPF (Figure 1 ). This is consistent with the assessment that body temperature may be not a crucial term for allometric scaling across the mammalian species [104, 105] . It is evident in the literature that renal transporters including organic cation transporters (OCT1 and OCT2) and multidrug and toxin extrusion proteins (MATE1 and MATE2) are involved in the active secretion of metformin [35] . Based on a quantitative proteomics approach, related transporters (e.g., OCT2 and MATE1) are highly expressed in the kidney, which were generally comparable among several species (e.g., mouse, rat, monkey, and man), arguing for the case that transporters are conserved [106] . While rabbits exhibited higher CL R particularly for anionic drugs [45] , which may be related to a higher extent of organic anion transporters (e.g., rbOAT1 and rbOAT3) compared to OCT transporters [107] , discordant CL and RPF in rabbits ( Figure 1 ) may be related to a lower affinity of metformin to rbOCT2 (e.g., IC50 value of 808 µM) than hOCT2 (339 µM) [108] . However, a previous PBPK model incorporating an electrogenic component for OCT1 and OCT2 transport suggested the plasma concentration profiles of metformin may be sensitive to transport activity mainly by MATE1/2 in the kidney [109, 110] . Therefore, further studies to explain the systemic CL smaller than RPF in rabbit and cat with regard to the differential expression/activities in the apical membrane transport of metformin in their kidneys would be interesting.

As shown in Figure 3 , the joint fitting of the mPBPK model reasonably captured the observed PK profiles for metformin in the nine species. One advantage of this approach is that fitting all profiles together allows the data to be shared in resolving the model parameters. The data from mice were rather limited, while most other species had more extensive measurements over longer times. Based on the in vitro PAMPA P of metformin (0.5 × 10 −6 cm/s) [74] , the apparent PAMPA P coefficient (i.e., f up P/R b = 0.5 × 10 −6 cm/s) supported the finding of permeability-limited distribution in mPBPK modeling [79] . Another in vitro study [111] showed that PAMPA permeability increased from 0.15 to 0.75 × 10 −6 cm/s with the increasing pH (i.e., pH ranging from 4.6 to 9.32), which indicates its low permeability over a range of physiologically relevant pH values. In addition to in vitro physicochemical properties, we alternatively applied a typical triexponential function for describing metformin PK (Supplementary Table S5 ). The total distributional clearance (CL D ) calculated from these in vivo observations [112, 113] (i.e., with initial condition restricted to be Dose/V B ) was significantly lower than Q CO , which also indicates that such modeling reveals that tissue permeability governs the rate of metformin distribution. It is noted that model-dependent CL values from the triexponential fitting were generally comparable with CL values from the separate fitting of mPBPK models in most species, except for dogs. This may be due to the overestimation of AUC in the triexponential function by restricting the initial condition to Dose/V B while the drug was an infused IV over 5 min (as we applied in the mPBPK models). Since f d1 (0.0390) is higher than f d2 (0.00677) ( Table 2) , Tissue 1 could be considered as rapidly perfused, while Tissue 2 as slowly perfused, when making a comparison to a traditional three-compartment model.

In this study, the apparent K p values observed for tissues in mouse and rat were collected from the literature as listed in Table 6 (for detailed information, refer to Supplementary Table S3 ). The K p1 values for mouse (0.921) and rat (0.575) obtained from the separate fitting approach appeared to be somewhat underestimated considering the in vivo tissue K p values corresponding to Tissue 1 (e.g., for tissues other than muscle and skin, 0.213-11.3 in mouse and 0.956-4.63 in rat). Nevertheless, the K p2 values in the mPBPK modeling for mouse (1.44) and rat (0.466) were consistent with their muscle K p (i.e., 1.03 for mouse and 0.597 for rat, averaged from multiple studies), supporting our kinetic assumption that Tissue 2 included muscle. For the larger species, however, the K p2 values were estimated to range from 0.675 to 7.53, which may imply species-dependent metformin partitioning to muscle and skin and/or different tissue components that are kinetically necessary for Tissue 2 in these species. Since metformin is a strong base (pKa = 11.8) and mainly exists in a cationic form in biological matrices, its binding to acidic phospholipids [114] and possible sequestration in the lysosomal/mitochondrial sub-compartments [115, 116] may complicate its tissue distribution kinetics with regard to both the extent (K p ) and rate ( f d ). Interestingly, the metformin distribution within cells showed that concentrations in the mitochondrial fraction (4.6 and 64.5 µM, after 22 h of metformin exposure in Hepa1-6 cells at 75 and 1000 µM) were significantly lower than in the cytosolic fraction (38.2 and 996 µM) [96] .

OCT1 is expressed in the basolateral membrane of hepatocytes, which is a major site of pharmacological action for metformin inhibiting glucose production [117] . In the literature [118] , the uptake rate of metformin into hepatocytes was highest in monkeys among four species (rats, monkeys, dogs, and man), which is consistent with the highest protein expression of hepatic OCT1 in monkeys [119] . These properties could contribute to the higher estimate of K p1 in monkeys (0.615) than rats (0.575) and man (0.484). Other distribution mechanisms may be involved, as V SS is markedly higher in dogs (10.1 L/kg in Table 1 and 6.0 L/kg in Table 4 ). Similarly, the 4.4-fold higher metformin uptake by mOCT1 than by hOCT1, along with the comparable expression level of the transporter in the liver [120] , was consistent with our K p1 values for mouse (0.921) and man (0.484). The higher uptake clearance normalized by the mRNA expression levels of hOCT1 compared to rOCT1 (with a factor of 13.1) [121] may explain the difference in our mPBPK model parameter f d1 between man (0.993) and rat (0.0694).

Based on the IV PK of metformin reasonably capturing the observations in nine species, the kinetic parameters for absorption could be estimated as shown in Figure 5 . There was a decreasing bioavailability of metformin with increasing BW, which was associated with the decrease in its effective absorption rate constant F·k a . From a single-pass intestinal perfusion study in rats [33] , P e f f ,rat was lower at a higher concentration of metformin (e.g., 47.1 × 10 −6 at 10 µg/mL, and 27.0 × 10 −6 cm/s at 200 µg/mL).

The P e f f ,rat calculated from our mPBPK modeling (4.61 × 10 −6 cm/s) was much lower than from this in situ study, probably because the 100 mg/kg oral dose given to rats may result in a luminal concentration of 4 mg/mL even under an assumption of instantaneous dissolution of the drug in the total fluid content of the stomach (2.29 g) and small intestine (3.89 g) [122] . This is much higher than 200 µg/mL, showing a probable saturable permeability. In addition, that a saturable oral absorption may be involved is also based on K m values (1.15 to 4 mM) for uptake transporters in the apical membrane of human intestine (e.g., thiamine transporter 2 (THTR-2), plasma membrane monoamine transporter (PMAT), and serotonin reuptake transporter (SERT)) [35] . The differential apical uptake of metformin in the intestinal barrier between hTHTR-2 and mTHTR-2 transporters [123] may be related to the P e f f ,man (42.6 × 10 −6 cm/s) being higher than P e f f ,mouse (7.0 × 10 −6 cm/s). In addition, the decrease or knockout of OCT1 activity in man and mouse led to the increase of oral absorption of metformin [117] . Since mOCT1/hOCT1 transporters were suggested to play a crucial role in basolateral uptake of metformin into intestinal tissue [124] , higher P e f f ,man may be involved with the species-dependent transport of metformin by intestinal OCT1. From our mPBPK modeling, the P e f f ,in vivo was estimated to be relatively higher in minipig, dog, and human than other species ( Figure 5 ). Considering the probable absence of solubility limitations of metformin in gut lumen, species differences in P e f f ,in vivo may result, at least in part, from the water content per gut length being higher in pig, dog, and human [122] , leading to enhanced active uptake across the gut wall with lower effective concentrations in the lumen. The presence of THTR-2 and OCT1 transporters in the GI tract undoubtedly account for the high tissue-to-plasma ratios observed in mice and rats (Table 6) .

Proctor et al. [34] examined metformin transport in Caco-2 cells to demonstrate carriermediated apical uptake that was seven-fold greater than basolateral efflux. The latter is limiting for transcellular transport and appears to account for the high metformin K p values in GI tissues ( Table 6 ).

The P e f f ,man predicted from this in vitro Caco-2 permeability (e.g., log P e f f = 0.4926 × log P app − 0.1454, where P app = 500 nm/s for metformin) [125] was 15.3 µm/s, which is comparable to P e f f ,man from our modeling analysis (42.6 µm/s). Based on the mPBPK model, metformin exhibits partial bioavailability in man (F = 0.485), compared to the fraction absorbed from the enteric compartment (F a = 0.185) calculated from a physiological absorption model (e.g., F a = 1 − 1 + 0.54P e f f −7 , using the P e f f of 15.3 µm/s in the Advanced Compartmental Absorption and Transit (ACAT) model) [126] . Using T SI and R (gut radius) values collected from the literature (Supplementary Table S2 ), the coefficient corresponding to P e f f term (i.e., 2T SI /7R) generally decreased with increasing body size. Thus, if there is no species difference in P e f f , which can be assumed to be the same with P e f f ,rat (e.g., 27.0 × 10 −6 cm/s) or P e f f ,man predicted from Caco-2 cell permeability (e.g., 15.3 × 10 −6 cm/s), the ACAT model consisting of seven enteric compartments argues for a decreasing F with increasing BW (Supplementary Figure S2) . This assessment suggests that the joint role of metformin permeability and gut physiology accounts for the decreasing bioavailability with increasing BW in these species. Species differences in bioavailability were judged to be considerable for a variety of compounds in a literature review [127] . However, the general applicability of this conclusion for F a needs further evaluation by factoring in the first-pass hepatic and gut wall extraction from the overall oral bioavailability.

Despite the abundant evidence for the involvement of carrier-mediated transport, the assumption of linear PK appeared to be adequate for describing the plasma PK of metformin in the current meta-analysis. In single-dose clinical studies [37, 40] , the terminal half-life of metformin based on urinary excretion rates was 4.22 to 5.11 times longer than that in plasma. However, this appears due to the longer sampling time in urine, which may better capture a late terminal phase. Metformin has a reported longer PK half-life when followed in erythrocytes (23.4 h) compared to plasma (2.7 h) after single doses [39] . This nonequilibrium between plasma and erythrocytes [37] is ascribed to a slow uptake (t 1/2 > 80 h) and slow efflux (t 1/2 > 32 h) from the cells [128] . While the PK profiles in these publications may not have been monitored over the same or for a sufficient period of time, Kajbaf et al. [129] followed metformin washout kinetics in patients with regard to both plasma and erythrocytes over 5 to 14 days. Although these patients had clinical complications, such as acute renal failure, the corresponding half-lives were 51.9 and 43.4 h. It is likely that many single-dose studies failed to assess metformin disposition over sufficiently long times to account for late phases and small amounts of the retained drug.

The marked difference between V SS (6-10 L/kg for dogs and 0.432-0.856 L/kg for man) and V d,area (44.8 L/kg for dogs and 0.859-3.83 L/kg for man) [36, 37, 40, 71] is likely related to differences in transporters and related tissue uptake. From the collected literature information (Supplementary Table S3 ), tissue K p values at early times (e.g., <24 h) ranged 0.359 to 0.771 for muscle and 0.0354 for brain in mouse [55, 130, 131] , while the prolonged treatment of metformin (i.e., from 7 to 30 days) produced K p values of 1.00-2.06 for muscle and 0.174-0.257 for the brain [132] . Similarly in rats, the brain K p increased from 0.2 (at 1 h) to 1.48 (at 24 h) [133] , while the K p values for muscle were 0.455 at 2 h of day 1 and 0.738 at 2 h of day 7 [134] . This time-dependent increase in tissue uptake of metformin may be related to differing terminal phase half-lives between the plasma and erythrocytes/tissues or may reflect drug-induced changes in the distribution. Collectively, therefore, carriermediated transport and/or continuous tissue redistribution may vary within or among species and alter the parameters in the PBPK model. More extensive studies that follow plasma, tissue, and urinary excretion of metformin during multiple dosing and washout are needed.

An insightful assessment of metformin PK in dogs was carried out by Wesolowski et al. [71] . Their fitting algorithm accounts for the early distribution kinetics during and after the IV infusion (much like seen in PET scans by Gormsen et al. [135] ) and handles the late disposition phase as a power function that captures the straight line phase on a log-log plot. This method applies a gamma-Pareto convolution that predicts slow continuing distribution consistent with the slow kinetics associated with red blood cells, the very low P-values (Table 5) , and the small f d values (Tables 2 and 4 ) for metformin. This approach yields a lower CL than the NCA analysis, and their mean CL was 15.1 compared to our mPBPK value of 19.3 mL/min/kg. It is also likely that V SS values are underestimated when applying an exponential to a power function curve. While this elegant method may better capture the terminal phase of metformin PK in dogs, a species with unusually high V SS (Table 1 ), it appears to be only applicable to single PK profiles and cannot, as yet, be utilized for tandem (oral and IV) PK profiles, as assessed here.

The present meta-assessment utilized the mPBPK modeling approach to compare the PK properties of metformin across nine species and supported the distribution model parameters with established calculation methods [81, 82] in Table 5 and extensive tissueto-plasma ratios for rodents summarized from the literature (Table 6 ). Typically, full PBPK models are built using both types of information. Burt et al. [109] utilized the Simcyp Simulator (Simcyp Ltd., Certara Co., Sheffield, UK) to predict metformin plasma concentrations in man, assuming F = 0.7, calculated K p [114] values for various tissues, the k a re-estimated using oral dose PK data, and renal and hepatic uptake related to transporter function with optimization using observed plasma and urinary excretion data. With such manipulations, single oral dose profiles of metformin were reasonably 'predicted'. Hanke et al. [42] carried out a similar PBPK modeling effort using PK-Sim ® (Bayer Technology Services, Leverkusen, Germany) to capture the metformin profiles from 22 studies with doses ranging 0.001 to 2550 mg for up to 696 h during multiple dosing. The analysis included K p values based on PET data obtained in plasma, blood, kidney, liver, intestines, and muscle over 1.5 h by Gormsen et al. [135] . Most of the assessment was based on matching plasma concentration profiles, which did not lend insight into longer time tissue distribution/accumulation kinetics of metformin.

The present meta-assessment has some limitations. While we sought all available PK studies of metformin that could be found in PubMed and by reference tracking, the drug has been subject to immense scrutiny, and undoubtedly, some relevant material may be missed. Diverse analytical methods ranging from early HPLC-UV to more recent LC-MS/MS may have a range of specificities and sensitivities ( Table 1) . Most of the data were digitized and dependent on published graphs with limitations based on single doses and in the duration of sampling (Figures 3 and 4) . While metformin is purported to have linear kinetics (as we assumed), this depends on substrate exposures being lower than transporter K m values. The washout PK in some species (like dogs) may better follow a power function, indicating complexities in distribution mechanisms that need tissue PK data and longer-term studies for better resolution. Nevertheless, the mPBPK models represent a state-of-the-art approach to modeling, summarizing, and comparing available single-dose profiles of metformin in the various species. Extending this modeling approach to include and compare oral and IV profiles across species provided unique insights into the absorption properties of metformin.

Upon completion of our analysis, Morse et al. [136] published a report demonstrating the interspecies comparison of the expression and activities of OCT transporters in the liver, which was investigated in relation to metformin PK. Using a single IV dose of 5-mg/kg metformin, wild-type rats showed comparable PK properties with the datasets (50, 100, and 200 mg/kg) used in the current analysis [52, 53] , in terms of CL (29.4 versus 23.6 to 26.4 mL/min/kg) and CL R (i.e., 21.8 versus 17.8 to 19.5 mL/min/kg). Consistent with the previous PBPK model suggesting the importance of apical efflux process of metformin by MATE1/2 transporters [109, 110] , single knockouts of the rOCT1 or rOCT2 transporters did not affect the plasma PK and tissue partitioning of metformin in rats. Different from man and mouse [117] , however, knockout of the rOCT1 transporter exhibited no significant effect on the oral absorption of metformin in rats. It is noted that their measurements of metformin concentrations in plasma up to 24 h after IV administrations led to observing a longer half-life (5.12 h) compared to 2.1-2.5 h determined using the data for 10 h [52] . Thus, mean residence time in the body increased by using the longer half-life values resulted in higher V SS (2.25 L/kg). The tissue-to-plasma ratios observed for various tissues at 4 h after the IV dose may govern the determination of V d,area rather than V SS [137] . Tissue partitioning properties mostly in a range far above the V SS value in the literature (2.25 L/kg) or in our current analysis (0.536 L/kg) may imply the involvement of low tissue permeability (Tables 2 and 4 ) and/or slow/continuous redistribution of metformin between plasma and tissues in rats, as discussed above.

In conclusion, we collected PK profiles from nine species for metformin, an effective antidiabetic agent widely used for various clinical and veterinary purposes. Utilizing a mPBPK model, the common features of metformin distribution that appear conserved among species were assessed in the joint fitting approach, while potential differences in PK parameters were evaluated in the individual fittings, with the results compared with the related literature data. Despite the lack of integrated information for tissue PK of metformin, a priori in silico calculations of physicochemical properties adopted in mPBPK models enabled the physiological interpretation of the absorption and disposition PK of metformin in various species. Unique perspectives were provided regarding the bioavailability that was found to be reduced with increased BW. This study provided a systematic review of the literature datasets and a quantitative analysis of metformin PK in all the available species and demonstrated the advantages of the simultaneous fitting of data from multiple species using system-specific physiological and anatomical variables. Further experimental studies are needed to resolve the mysteries of metformin absorption and distribution with the implementation of full and extended PBPK models utilizing longer term data to characterize its kinetics more completely.

The PK parameters for metformin including the systemic clearance (CL), oral clearance (CL/F), renal clearance (CL R ), and steady-state volume of distribution (V SS ) were collected from the literature for nine species (Table 1) based on a PubMed search and reference tracking. If not directly obtainable in the literature, these PK parameters were either calculated from available descriptors or obtained by noncompartmental analysis (NCA) of the plasma concentration-time profiles that were digitized from the published graphs using GetData Graph Digitizer version 2.26 (http://getdata-graph-digitizer.com/, accessed on 1 May 2021).

The physiological and anatomical variables necessary for mPBPK modeling were obtained for the species in which the metformin PK was found in the literature (mouse, rat, rabbit, cat, monkey, minipig, dog, man, and horse). The blood volumes were from one literature source [138] , except for man [139] . The cardiac output for each species was determined by an allometric relationship [66] :

In addition, the renal plasma flow (RPF) [65] [66] [67] [68] [69] [70] [71] [72] , glomerular filtration rate (GFR) [70, 105, 140, 141] , gut radius (R) [70, [142] [143] [144] [145] [146] , luminal volume of the GI tract (V Lumen for the stomach and small intestine) [122] , and the small intestinal transit time (T SI ) [122, [147] [148] [149] [150] [151] were obtained for each species. The detailed values including the anatomical volumes and blood perfusion rates to muscle (V MU , Q MU ) and skin (V SK , Q SK ) for the species-specific input variables are summarized in Supplementary Table S2 . Allo-metric relationships of these variables among the 9 species are provided in Supplementary Figure S1.

In the mPBPK model structure (Figure 2 ), blood and two peripheral tissue compartments (Tissues 1 and 2) were assumed. The differential equations used for the joint fitting are:

where C p is the plasma concentration of metformin, C 1 and C 2 are the metformin concentrations in Tissues 1 and 2, f d1 and f d2 are the fractional distribution parameters for Tissues 1 and 2, V 1 and V 2 are the anatomical volumes of Tissues 1 and 2, Q CO is cardiac output, K p is the tissue-to-plasma partition coefficient, R b is the blood-to-plasma ratio, V B is the blood volume, and CL is the systemic clearance. The IV Input was defined in accordance with the methodology of the literature source (e.g., bolus or short infusion, where Input = Dose/Duration of infusion). The initial conditions for C 1 and C 2 were set as zero. The blood partitioning properties of metformin (e.g., R b and f up (free fraction in the plasma)) were found to be in unity [39, 109, 128] . The BW (in kg) for each species was obtained from the publications with the metformin PK data ( Table 1 ). The tissue volume fractions for Tissue 1 ( f t ) and Tissue 2 (1 − f t ) were used to calculate V 1 and V 2 , assuming the density of all tissues to be in unity:

Due to the low permeability characteristics of metformin (e.g., in vitro PAMPA permeability of 0.5 × 10 −6 cm/s) [74] , the permeability-limited distribution [79] was assumed by using f d,total conserved across all species:

The primary estimated parameters, along with CL, were K p , f t , f d,total , and f d1 , while f d2 was a secondary parameter. It can be noted that, except for CL, the other parameters are dimensionless.

Based on the collected PK data ( Table 1) , species differences were observed in metformin PK (e.g., bioavailability from 7.1% to 100% and V SS from 0.321 to 10.1 L/kg). Since the involvement of various transporters is evident for metformin PK [35] , there is a possibility that species-dependent transport activity may lead to different PK variables for each species. Therefore, we further conducted the mPBPK model fitting to oral and IV profiles for each species separately with the addition of more PK data.

A priori in silico calculation of the tissue distribution parameters (K p ) was carried out for 11 typical tissues of rats as applied in whole-body PBPK (WB-PBPK) models. Physicochemical (e.g., log P = −1.43, pKa = 11.8) [30, 32] and blood partitioning properties of metformin ( f up = 1, R b = 1) [39, 109, 128] were used. The K p values from two calculation methods (i.e., Poulin and Theil [81] , and Berezhkovskiy [82] ) were applied for the prediction of V SS amongst the four methods found in the GastroPlus PBPK simulator (version 9.8.0002; Simulations Plus, Inc., Lancaster, CA, USA). For the calculation of f d , two approaches were considered for the tissue permeability coefficient P: (i) in vitro PAMPA P (0.5 × 10 −6 cm/s) [74] and (ii) the P-value fitted to the rat PK data (0.0102 × 10 −6 cm/s). For the latter case, a WB-PBPK model structure employing S e f f (i.e., the effective surface area across the interface between the systemic circulation and tissues) [79] was used. For the calculation of f d,tissue :

fup S e f f P Q T R b (8) where Q T is the blood flow to each tissue. Mean transit time (MTT) through the tissues was calculated as:

where V T is the anatomical tissue volume. Compared to traditional compartment models, mPBPK models enable a more physiological interpretation of the systemic PK. Since f t determined by the joint fitting was 0.172 with a single K p assumption (see Results), it was reasoned that muscle, the anatomically largest organ, is included in Tissue 2 (e.g., V MU /(BW − V B ) ranges from 0.411 to 0.574 (>0.172) in the 9 species; Supplementary Table S2 ). Since our in silico calculations suggested skin as the kinetically largest tissue for variable P (see Results), then Tissue 2, the slowly equilibrating compartment, also included skin (i.e., calculated f t value ranging from 0.472 to 0.678 for Tissue 2 in the 9 species). The differential equations used for the separate fittings are:

where Q 2 is the blood flow rate to Tissue 2 (Q MU + Q SK ), Q 1 is the blood flow rate to Tissue 1 (Q CO − Q 2 ), V 2 is the anatomical volume of Tissue 2 (V MU + V SK ), V 1 is the anatomical volume of Tissue 1 (BW − V B − V 2 ), and K p1 and K p2 are the tissue-to-plasma partition coefficients for Tissues 1 and 2. The IV doses had bolus or infusion inputs, while, for the oral administration, the drug Input from the absorption compartment is:

where X a is the drug amount at the absorption site with the initial condition of X a (0) = Dose, k a is the first-order absorption rate constant, and F is the bioavailability.

For the nonlinear regression analysis, the maximum likelihood method in ADAPT 5 was used [152] . The variance model was:

where V i is the variance of the ith data point, Y i is the ith model prediction, σ 1 and σ 2 are the variance model parameters that were estimated together with the system parameters during the model fittings. The goodness-of-fit of the model was monitored by visual inspection, Akaike Information Criterion (AIC), and Coefficient of Variation (CV%) of the estimates. The ADAPT model code for the joint mPBPK is provided in the Supplementary Materials. The NCA were conducted using WinNonlin Professional 5.0.1. software (Pharsight Corporation, Mountain View, CA, USA).

The following are available online at https://www.mdpi.com/article/ 10.3390/ph14060545/s1: Figure S1 : Physiological and anatomical information (i.e., V SK , V MU , Q SK , Q MU , RPF, and GFR) collected from different sources (Supplementary Table S2 ) showed an allometric relationship among 9 species (i.e., mouse (0.025 kg), rat (0.25 kg), rabbit (2 kg), cat (5 kg), monkey (7 kg), minipig (14 kg), dog (28 kg), man (70 kg), and horse (530 kg))., Figure S2 : ACAT modelbased estimation of bioavailability using F a = 1 − 1 + 2P e f f T SI /7R −7 (7-enteric compartments), assuming that the effective intestinal permeability (P e f f ) of the 9 species is the same with that of rat (P e f f ,rat ), and that of man (P e f f ,man ) predicted from Caco-2 cell permeability, which were compared with bioavailability determined from our mPBPK modeling (1-enteric compartment as a cylinder). Table S1 : Literature reports and pharmacokinetic parameters of metformin given intravenously in different species. Table S2 : Summary of the physiological input variables for the meta-analysis of metformin pharmacokinetics in various species. Table S3 : Literature information collected for tissue distribution and blood partitioning of metformin. Table S4 : Calculation of the steady-state volume of distribution (V SS ) using in vivo tissue K p values based on the equation V SS = V B + ∑ V T,i K p,i [84] . Tissue K p values were assumed to be muscle K p (1.03 for mouse and 0.597 for rat) if unavailable in Table 6 . Table S5 : Model fitting results using a triexponential function (C p (t) = C 1 e −λ 1 t + C 2 e −λ 2 t + C 3 e −λ 3 t ) for metformin PK in various species (CV% obtained by fitting). It is noted that five parameters (C 1 , λ 1 , C 2 , λ 2 , and λ 3 ) were optimized, while C 3 values were estimated as a secondary parameter (C 3 = C 0 − C 1 − C 2 ; C 0 = Dose/V B ). In addition, the ADAPT model code for the joint mPBPK is provided in the Supplementary Materials.

Author Contributions: Conceptualization, methodology, validation, formal analysis, investigation, resources, data curation, writing-original draft preparation and review and editing, and visualization, Y.-S.J. and W.J.J. and supervision, project administration, and funding acquisition, W.J.J. All authors have read and agreed to the published version of the manuscript.

Funding: This research was supported by the NIH, grant R35 GM131800.

Data Availability Statement: All the data generated during this study are included in this article. Further datasets digitized and/or analyzed during the study are available from the corresponding author upon reasonable request.

The Vegetable System, or, the Internal Structure and the Life of Plants: Their Classes, Orders, Genera, and Species Ascertained and Described in a Method Altogether New: Comprehending an Artificial Index, and a Natural System: With Figures of all the Plants Designed and Engraved by the Author: The Whole from Nature Only

Traditional Plant Medicines as Treatments for Diabetes

Metformin; a review of its history and future: From lilac to longevity

CCXIV-The preparation of methylguanidine, and of ββ-dimethylguanidine by the interaction of dicyanodiamide, and methylammonium and dimethylammonium chlorides respectively

The Discovery of Insulin: An Important Milestone in the History of

Pharmacological Studies of a New Oral Hypoglycemic Drug

Zur Behandlung der Zuckerkrankheit mit Biguaniden

Du nouveau dans les antidiabetiques. La NN dimethylamine guanyl guanide (NNDG)

Diabetes and the Biguanides: The Mystery of Each; Pasik, C., Ed.; Glucophage: Serving diabetology for 40 years

Medical Management of Hyperglycemia in Type 2 Diabetes: A Consensus Algorithm for the Initiation and Adjustment of Therapy: A consensus statement of the American Diabetes Association and the European Association for the Study of Diabetes

Effect of intensive blood-glucose control with metformin on complications in overweight patients with type 2 diabetes (UKPDS 34)

Metformin and reduced risk of cancer in diabetic patients

Antidiabetic Therapies Affect Risk of Pancreatic Cancer

Metformin Associated with Lower Cancer Mortality in Type 2 Diabetes: ZODIAC-16

Metformin Use Is Associated with Reduced Mortality in a Diverse Population With COVID-19 and Diabetes

Metformin and risk of mortality in patients hospitalised with COVID-19: A retrospective cohort analysis

Metformin: Historical overview

Emerging Trends in Metformin Prescribing in the United States from

Biguanides suppress hepatic glucagon signalling by decreasing production of cyclic AMP

Metformin suppresses gluconeogenesis by inhibiting mitochondrial glycerophosphate dehydrogenase

Effect of metformin on glucose metabolism in mouse soleus muscle

Effect of metformin treatment on insulin action in diabetic rats: In vivo and in vitro correlations

Metformin: An Update

Importance of the intestine as a site of metformin-stimulated glucose utilization

Inhibitory effect of metformin on intestinal glucose absorption in the perfused rat intestine

Pharmacokinetic-Pharmacodynamic Analysis of the Glucose-Lowering Effect of Metformin in Diabetic Rats Reveals First-Pass Pharmacodynamic Effect

Pharmacokinetic/Pharmacodynamic Analysis of Metformin using Different Models in Diabetic Rats

Complex Compounds of Biguanides and Guanylureas with Metallic Elements

Biowaiver extension potential to BCS Class III high solubility-low permeability drugs: Bridging evidence for metformin immediate-release tablet

Interaction of Biguanides with Mitochondrial and Synthetic Membranes

Intestinal permeability of metformin using single-pass intestinal perfusion in rats

Mechanisms Underlying Saturable Intestinal Absorption of Metformin

Transporters Involved in Metformin Pharmacokinetics and Treatment Response

Disposition of metformin (N, N-dimethylbiguanide) in man

Metformin kinetics in healthy subjects and in patients with diabetes mellitus

Kidney Function and Age Are Both Predictors of Pharmacokinetics of Metformin

Kinetics of plasma and erythrocyte metformin after acute administration in healthy subjects

Pharmacokinetics of metformin after intravenous and oral administration to man

Accumulation of metformin by tissues of the normal and diabetic mouse

A Comprehensive Whole-Body Physiologically Based Pharmacokinetic Drug-Drug-Gene Interaction Model of Metformin and Cimetidine in Healthy Adults and Renally Impaired Individuals

Applications of minimal physiologically-based pharmacokinetic models

Interspecies scaling, allometry, physiological time, and the ground plan of pharmacokinetics

Humans are animals, but are animals human enough? A systematic review and meta-analysis on interspecies differences in renal drug clearance

Interspecies Scaling in Pharmacokinetics: A Novel Whole-Body Physiologically Based Modeling Framework to Discover Drug Biodistribution Mechanisms in vivo

Across-species meta-analysis of dexamethasone pharmacokinetics utilizing allometric and scaling modeling approaches

Across-Species Scaling of Monoclonal Antibody Pharmacokinetics Using a Minimal PBPK Model

Interspecies allometric meta-analysis of the comparative pharmacokinetics of 85 drugs across veterinary and laboratory animal species

Evaluation of the Allometric Exponents in Prediction of Human Drug Clearance

Pharmacokinetics and oral bioavailability of metformin hydrochloride in healthy mixed-breed dogs

Dose-Independent Pharmacokinetics of Metformin in Rats: Hepatic and Gastrointestinal First-Pass Effects

Pharmacokinetic interaction between itraconazole and metformin in rats: Competitive inhibition of metabolism of each drug by each other via hepatic and intestinal CYP3A1/2

Ablation of Both Organic Cation Transporter (Oct)1 and Oct2 Alters Metformin Pharmacokinetics but Has No Effect on Tissue Drug Exposure and Pharmacodynamics

Involvement of Carnitine/Organic Cation Transporter OCTN1/SLC22A4 in Gastrointestinal Absorption of Metformin

Fernández-Varón, E. Potential of sustained release microparticles of metformin in veterinary medicine: An in vivo pharmacokinetic study of metformin microparticles as oral sustained release formulation in rabbits

Drug interaction of metformin and aspirin in rabbits

Pharmacokinetics of the antihyperglycemic agent metformin in cats

Evaluation of the oral antihyperglycemic drug metformin in normal and diabetic cats

Cynomolgus Monkey as a Clinically Relevant Model to Study Transport Involving Renal Organic Cation Transporters: In Vitro and In Vivo Evaluation

Fast liquid chromatographic-tandem mass spectrometric (LC-MS-MS) determination of metformin in plasma samples

Pharmacokinetics of intravenous and oral metformin and r,s-verapamil in Sinclair, Hanford, Yucatan and Göttingen minipigs

Pharmacokinetic-pharmacodynamic modeling for the relationship between glucose-lowering effect and plasma concentration of metformin in volunteers

Pharmacokinetics and bioavailability of metformin in horses

Differential effect of T-type voltage-gated Ca 2+ channel disruption on renal plasma flow and glomerular filtration rate in vivo

Physiological Parameter Values for Physiologically Based Pharmacokinetic Models

Development of a physiologically based pharmacokinetic (PBPK) model for methyl iodide in rats, rabbits, and humans

Use of allometry in predicting anatomical and physiological parameters of mammals

The Simcyp ® Population-based ADME Simulator

A Physiologically Based Pharmacokinetic Model of the Minipig: Data Compilation and Model Implementation

Comparison of the gamma-Pareto convolution with conventional methods of characterising metformin pharmacokinetics in dogs

Glomerular filtration rate, effective renal plasma flow, blood pressure and pulse rate in the equine neonate during the first 10 days post partum

Absorption, distribution in the organism and elimination of metformin

Current industrial practices of assessing permeability and P-glycoprotein interaction

Lumping of whole-body physiologically based pharmacokinetic models

Reduction and Lumping of Physiologically Based Pharmacokinetic Models: Prediction of the Disposition of Fentanyl and Pethidine in Humans by Successively Simplified Models

Reducing Whole Body Physiologically Based Pharmacokinetic Models Using Global Sensitivity Analysis: Diazepam Case Study

Lumping of physiologically-based pharmacokinetic models and a mechanistic derivation of classical compartmental models

Estimation of the minimum permeability coefficient in rats for perfusion-limited tissue distribution in whole-body physiologically-based pharmacokinetics

Role of human placental apical membrane transporters in the efflux of glyburide, rosiglitazone, and metformin

Prediction of Pharmacokinetics Prior to In Vivo Studies. 1. Mechanism-Based Prediction of Volume of Distribution

Volume of Distribution at Steady State for a Linear Pharmacokinetic System with Peripheral Elimination

A theoretical basis for a biopharmaceutic drug classification: The correlation of in vitro drug product dissolution and in vivo bioavailability

Effect of Altered Plasma Protein Binding on Apparent Volume of Distribution

Pharmacologic Therapy for Type 2 Diabetes Mellitus

Drugs inducing hypoglycemia

Effect of Metformin on Postprandial Lipemia in Patients with Fairly to Poorly Controlled NIDDM

Efficacy of Metformin in Patients with Non-Insulin-Dependent Diabetes Mellitus

Mechanism of Metformin Action in Obese and Lean Noninsulin-Dependent Diabetic Subjects

Combined metforminsulfonylurea treatment of patients with noninsulin-dependent diabetes in fair to poor glycemic control

Effective breast cancer combination therapy targeting BACH1 and mitochondrial metabolism

Metformin as an anti-cancer agent: Actions and mechanisms targeting cancer stem cells

Effects of metformin and other biguanides on oxidative phosphorylation in mitochondria

Biguanide-induced mitochondrial dysfunction yields increased lactate production and cytotoxicity of aerobically-poised HepG2 cells and human hepatocytes in vitro

Effects of metformin on weight loss: Potential mechanisms

Metformin as a Tool to Target Aging

Benefits of Metformin in Attenuating the Hallmarks of Aging

Metformin: A review of its therapeutic efficacy and adverse effects

Effects of enzyme inducers and inhibitors on the pharmacokinetics of metformin in rats: Involvement of CYP2C11, 2D1 and 3A1/2 for the metabolism of metformin

Glomerular filtration and renal plasma flow during hyperthermia

Scaling of antipyrine intrinsic clearance of unbound drug in 15 mammalian species

Similarity Principles and Intrinsic Geometries: Contrasting Approaches to Interspecies Scaling

Species similarities and differences in pharmacokinetics

Kidney Cortical Transporter Expression across Species Using Quantitative Proteomics

Sex Differences in the mRNA, Protein, and Functional Expression of Organic Anion Transporter (Oat) 1, Oat3, and Organic Cation Transporter (Oct) 2 in Rabbit Renal Proximal Tubules

Molecular Determinants of Substrate/Inhibitor Binding to the Human and Rabbit Renal Organic Cation Transporters hOCT2 and rbOCT2

Metformin and cimetidine: Physiologically based pharmacokinetic modelling to investigate transporter mediated drug-drug interactions

Physiologically-Based Pharmacokinetic Modeling Analysis for Quantitative Prediction of Renal Transporter-Mediated Interactions between Metformin and Cimetidine

A New PAMPA Model Proposed on the Basis of a Synthetic Phospholipid Membrane

Theorems and implications of a model independent elimination/distribution function decomposition of linear and some nonlinear drug dispositions. I. Derivations and theoretical analysis

Theorems and implications of a model-independent elimination/distribution function decomposition of linear and some nonlinear drug dispositions. II. Clearance concepts applied to the evaluation of distribution kinetics

Physiologically Based Pharmacokinetic Modeling 1: Predicting the Tissue Distribution of Moderate-to-Strong Bases

Incorporation of lysosomal sequestration in the mechanistic model for prediction of tissue distribution of basic drugs

Quantitative modeling of selective lysosomal targeting for drug design

Effect of genetic variation in the organic cation transporter 1 (OCT1) on metformin action

Comparison of uptake transporter functions in hepatocytes in different species to determine the optimal model for evaluating drug transporter activities in humans

Interspecies Variability in Expression of Hepatobiliary Transporters across Human, Dog, Monkey, and Rat as Determined by Quantitative Proteomics

Differences in Metformin and Thiamine Uptake between Human and Mouse Organic Cation Transporter 1: Structural Determinants and Potential Consequences for Intrahepatic Concentrations

Metformin is a Superior Substrate for Renal Organic Cation Transporter OCT2 rather than Hepatic OCT1

Animal Farm: Considerations in Animal Gastrointestinal Physiology and Relevance to Drug Delivery in Humans

Metformin Is a Substrate and Inhibitor of the Human Thiamine Transporter, THTR-2 (SLC19A3)

Involvement of Organic Cation Transporter 1 in Hepatic and Intestinal Distribution of Metformin

Prediction of Human Intestinal First-Pass Metabolism of 25 CYP3A Substrates from In Vitro Clearance and Permeability Data

A compartmental absorption and transit model for estimating oral drug absorption

Animal versus human oral drug bioavailability: Do they correlate?

Metformin's intrinsic blood-to-plasma partition ratio (B/P): Reconciling the perceived high in vivo B/P> 10 with the in vitro equilibrium value of unity

Unexpectedly long half-life of metformin elimination in cases of metformin accumulation

Taste of a Pill

Loss of multidrug and toxin extrusion 1 (MATE1) is associated with metformin-induced lactic acidosis

Determination of metformin bio-distribution by LC-MS/MS in mice treated with a clinically relevant paradigm

Quantification of metformin by the HPLC method in brain regions, cerebrospinal fluid and plasma of rats treated with lipopolysaccharide

Metoprolol decreases the plasma exposure of metformin via the induction of liver, kidney and muscle uptake in rats

In Vivo Imaging of Human 11C-Metformin in Peripheral Organs: Dosimetry, Biodistribution, and Kinetic Analyses

Comparison of Hepatic Transporter Tissue Expression in Rodents and Interspecies Hepatic OCT1 Activity

The Connection Between the Steady State (Vss) and Terminal (Vβ) Volumes of Distribution in Linear Pharmacokinetics and The General Proof That Vβ ≥ Vss

Handbook of Laboratory Animal Management and Welfare

Physiological Parameters in Laboratory Animals and Humans

An overview of glomerular filtration rate testing in dogs and cats

Evaluation of a single injection of 99mTc-labeled diethylenetriaminepentaacetic acid for measuring glomerular filtration rate in horses

Comparison of the gastrointestinal anatomy, physiology, and biochemistry of humans and commonly used laboratory animals

Origin of regional and species differences in intestinal glucose uptake

Assessment of gastrointestinal pH, fluid and lymphoid tissue in the guinea pig, rabbit and pig, and implications for their use in drug development

Hypertrophy of intestinal smooth muscle in cats

The maximum attainable body size of herbivorous mammals: Morphophysiological constraints on foregut, and adaptations of hindgut fermenters

Establishment of a Protocol for Determining Gastrointestinal Transit Time in Mice Using Barium and Radiopaque Markers

Employment of a noninvasive magnetic method for evaluation of gastrointestinal transit in rats

Rabbit gastrointestinal physiology

Radiopaque Markers to Evaluate Gastric Emptying and Small Intestinal Transit Time in Healthy Cats

A wireless endoscopy capsule suitable for imaging of the equine stomach and small intestine

ADAPT 5 User's Guide: Pharmacokinetic/Pharmacodynamic Systems Analysis Software

Biomedical Simulations Resource

The authors greatly appreciate Charlotte A. Johnston, who provided the raw data of metformin PK in dogs for the current modeling analysis. We appreciated the review and comments of Geoffrey T. Tucker of the Medical and Biological Sciences, University of Sheffield, and Carl A. Wesolowski of the College of Medicine, University of Saskatchewan.

The authors have declared that there are no conflicts of interest.