Background/Objectives: Gabapentin has variable pharmacokinetics (PK), which contributes to difficulty in dosing and increased risk of adverse events. The objective of this study was to leverage gabapentin concentrations from therapeutic drug monitoring (TDM) to develop a population PK (popPK) model and characterize significant covariates that impact gabapentin PK. Methods: Data were retrospectively collected from 82 hospitalized adult patients with TDM gabapentin concentrations. Renal function indicators (i.e., estimated glomerular filtration rate, creatinine clearance, acute kidney injury), body weight parameters (i.e., actual body weight, ideal body weight, adjusted body weight, lean body weight, body mass index, obesity status), fasting plasma glucose levels, and diagnosis of type 2 diabetes were tested as potential covariates. A popPK model was developed in MONOLIX (2020R1, Lixoft, France). Results: A one-compartment model best described gabapentin PK with first-order absorption, dose-dependent bioavailability, first-order elimination, and no lag time. Population parameter estimates for the volume of distribution (V_d), and clearance (Cl) were 44.61 L, and 5.73 L/h, respectively. Serum creatinine was a significant covariate on Cl. Conclusions: The popPK model highlights the importance of renal function in the interindividual variability of gabapentin PK and suggests that diabetes and body weight parameters have no impact on gabapentin PK. Moreover, our study supports the utility of leveraging data obtained from clinical TDM for popPK model development.