RESUMEN
Methotrexate is the gold standard treatment in rheumatoid arthritis. Once absorbed, it is internalized in cells, where glutamate residues are added to produce polyglutamated forms, which are responsible for the effect of methotrexate. The aim of the current study is to determine the relationship between methotrexate triglutamate concentrations and the clinical evolution in rheumatoid arthritis patients, as well as to characterize the variability in both features to propose strategies for low-dose methotrexate optimization. The quantification of methotrexate triglutamate concentration in red blood cells was performed through ultra-performance liquid chromatography coupled with mass spectrometry. Polymorphisms of genes involved in the formation of polyglutamates were determined by real-time polymerase chain reaction. A multivariate regression was performed to determine the covariates involved in the variability of methotrexate triglutamate concentrations and a population pharmacokinetics model was developed through nonlinear mixed-effects modeling. Disease activity score changed according to methotrexate triglutamate concentrations; patients with good response to treatment had higher concentrations than moderate or nonresponding patients. The methotrexate triglutamate concentrations were related to time under treatment, dose, red blood cells, and body mass index. A 1-compartment open model was selected to estimate the pharmacokinetic parameters; the typical total clearance (L/day) was determined as 1.45 * (body mass index/28 kg/m2 ) * (red blood cells/4.6 × 106 cells/µL) and the volume of distribution was 52.4 L, with an absorption rate of 0.0346/day and a fraction metabolized of 1.03%. Through the application of the model, the initial dose of methotrexate is proposed on the basis of stochastic simulations and considering methotrexate triglutamate concentrations found in responders patients.