Research Article
Molecular Modelling Analysis of the Metabolism of Ceftiofur
School of Biomedical Sciences, Faculty of Health Sciences, C42, The University of Sydney, P.O. Box 1 70, Lidcombe, NSW 1 825, Australia
Ceftiofur sodium (CF; (6R,7R)-((2-amino-4-thiazolyl)-Z-(methoxyimino)acetyl)amino-3-((2-+++furanylcarbonyl)thio)methyl-8-oxo-5-thia-1-azabicyclo(4.2.0)oct-2-ene-2-carboxylic acid) is a third generation broad-spectrum cephalosphorin that is active against both Gram-positive and Gram-negative pathogenic bacteria of veterinary importance (Becker et al. 2003). It has been approved in the EU for subcutaneous treatment of certain respiratory diseases in cattle, horses, pigs, poultry and dogs. In horses, CF has been found to be effective in the treatment and control of respiratory infections caused by Streptococcus zooepidemiincus (Folz et al., 1992) and also in the treatment against septic arthritis (Mills et al., 2000). CF can be administered as an i.m. or s.c. injection. It is minimally bioavailable after oral administration.
The metabolism of CF has been described in rats (Jaglan et al., 1989) and swine (Gilbertson et al., 1995). CF is rapidly metabolized to the active metabolite desfuroylceftiofur (DFC) and Furoic Acid (FA) after parenteral administration. The half-life of the parent drug is less than 10 min, because of rapid cleavage of the thioseter bond. DFC which has a free thiol group is further metabolized to disulfides such as desfuroylceftiofur dimer (DFC-D), desfuroylceftiofur cysteine disulfide (DFC-CYS) and desfuroylceftiofur glutathione (DFC-GS) and desfuroylcetiofur protein conjugate (Olson et al., 1998; Beconi-Barker et al., 1995 a,b; Glilbertson et al., 1995; Beconi-Barker et al., 1995) which may be playing a role in the activity and efficacy of CF. CF has been shown to produce chromosomal aberrations following extended treatment of Chinese hamster ovary cells in culture at high doses (Aaron et al., 1995). The clastogenic effect due to CF is found to be due to prolongation of cell cycle. In this study molecular modelling analyses have been carried out to obtain information on the relative toxicity of CF and its metabolites. The study was done in The School of Biomedical Sciences, The University of Sydney during February to July 2006.
COMPUTATIONAL METHODS
The geometries of CF and its metabolites DFC, FA, DFC-D and DFC-CYS have been optimized based on molecular mechanics (Fig. 1) semi-empirical and DFT calculations, using the molecular modelling program Spartan 02 (Spartan 2002).
Fig. 1: | Metabolic pathway for ceftiofur (Based on De Baere et al., 2004) |
Fig. 2: | Structure of CF giving in: (a) the electrostatic potential (greyish envelope denotes negative electrostatic potential), (b) the HOMOs, (where red indicates HOMOs with high electron density) (c) the LUMOs (where blue indicates LUMOs) and in (d) surface electric charges (where red indicates negative, blue indicates positive and green indicates neutral) |
Fig. 3: | Structure of DFC giving in: (a) the electrostatic potential (greyish envelope denotes negative electrostatic potential), (b) the HOMOs, (where red indicates HOMOs with high electron density) (c) the LUMOs (where blue indicates LUMOs) and in (d) surface electric charges (where red indicates negative, blue indicates positive and green indicates |
Fig. 4: | Structure of FA giving in: (a) the electrostatic potential (greyish envelope denotes negative electrostatic potential), (b) the HOMOs, (where red indicates HOMOs with high electron density) (c) the LUMOs (where blue indicates LUMOs) and in (d) surface electric charges (where red indicates negative, blue indicates positive and green indicates neutral) |
No calculations were done for DFC-GS. Molecular mechanics calculations were carried out using MMFF force field. Semi-empirical calculations were carried out using the routine PM3. DFT calculations were carried at B3LYP/6-31G* level. In optimization calculations, a RMS gradient of 0.001 was set as the terminating condition. For the optimized structures, single point calculations were carried out to give heat of formation, enthalpy, entropy, free energy, dipole moment, solvation energy, energies for HOMO and LUMO. The order of calculations: molecular mechanics followed by semi-empirical followed by DFT ensured that the structure was not embedded in a local minimum. To further check whether the global minimum was reached, some calculations were carried out with improvable structures. It was found that when the stated order was followed, structure corresponding to the global minimum or close to that could ultimately be reached in all cases.
Fig. 5: | Structure of DFC-D giving in: (a) the electrostatic potential (greyish envelope denotes negative electrostatic potential), (b) the HOMOs, (where red indicates HOMOs with high electron density) (c) the LUMOs (where blue indicates LUMOs) and in (d) surface electric charges (where red indicates negative, blue indicates positive and green indicates neutral) |
Although RMS gradient of 0.001 may not be sufficiently low for vibrational analysis, it is believed to be sufficient for calculations associated with electronic energy levels.
Fig. 6: | Structure of DFC-CYS giving in: (a) the electrostatic potential (greyish envelope denotes negative electrostatic potential), (b) the HOMOs, (where red indicates HOMOs with high electron density) (c) the LUMOs (where blue indicates LUMOs) and in (d) surface electric charges (where red indicates negative, blue indicates positive and green indicates neutral) |
Table 1 gives the total energy, heat of formation as per PM3 calculation, enthalpy, entropy, free energy, surface area, volume, dipole moment, energies of HOMO and LUMO as per both PM3 and DFT calculations for CF and its metabolites DFC, FA, DFC-D and DFC-CYS.
Table 1: | Calculated thermodynamic and other parameters of CF and its metabolites |
* in atomic units from DFT calculations |
Figure 2-6 give the regions of negative electrostatic potential (greyish-white envelopes) in (a), HOMOs (where red indicates HOMOs with high electron density) in (b), LUMOs in (c) and surface charges (where red indicates negative, blue indicates positive and green indicates neutral) in (d) as applied to the optimized structures of CF and its metabolites DFC, FA, DFC-D and DFC-CYS.
The calculated solvation energies of CF and its metabolites DFC, FA, DFC-CYS and DFC-D from PM3 calculations in kcal mol-1 are -30.34, -27.71, -9.79, -44.14 and -36.69 and their dipole moments from DFT calculations are 6.8, 7.7, 5.1, 5.3 and 11.6 , respectively. The values suggest that CF and its metabolites would vary in their solubility in water with the terminal metabolites DFC-D and DFC-CYS having much greater solubility.
CF and its metabolites have LUMO-HOMO energy differences ranging from 4.0 to 5.6 eV from DFT calculations, indicating that the compounds would differ in their kinetic lability. DFC-D, CF and DFC are expected to be more labile as they have smaller LUMO-HOMO energy differences, the most labile one being DFC-D.
In the case of CF, DFC, DFC-D and DFC-CYS, electrostatic potential is found to be more negative around the various oxygen centers so that they may be more likely subject to electrophilic attack. In the case of DFC, electrostatic potential is also found to be more negative around one of the sulfur centers so that it may also be subject to electrophilic attack. In the case of FA, electrostatic potential is found to be more negative around the two oxygen centers of the carboxyl group so that they may be more likely subject to electrophilic attack. In the case of DFC-CYS, electrostatic potential is found to be more negative around the two oxygen centers of the carboxyl group so that they may be more likely subject to electrophilic attack.
In the case of CF and DFC, HOMOs with high electron density are found close to the most of the non-hydrogen atoms of the thiazolyl ring whereas LUMOs are found close to the non-hydrogen atoms of the linking chain between thiazolyl ring and the four-membered and six-membered fused rings. In the case of DFC-D, the positions of HOMOs with high electron density and LUMOs are found to be quite different but not in the case of DFC-CYS. In the case FA, both HOMOs with high electron density and LUMOs are found close to the most of the non-hydrogen atoms.
The overlap of HOMO with high electron density and region of negative electrostatic potential close to sulfur, gives further support to the idea that the position may be subject to electrophilic attack.
When the surface areas and volumes are considered, it is found that the active metabolite DFC has distinctly different values from those of CF and other metabolites, suggesting that CF and its other metabolites may not be able to act substrates for the enzyme(s) to which the active metabolite DFC associates.
It can be seen that the molecular surfaces of CF and all its metabolites have some electron-deficient regions (blue) so that they may be subject to nucleophilic attacks such as those by glutathione and nucleobases in DNA. This means that DFC-D, CF and DFC can react more readily with cellular glutathione and nucleobases in DNA resulting into glutathione depletion and oxidation of nucleobases , respectively. Oxidative stress induces cellular toxicity whereas oxidation of nucleobases causes DNA damage.
Molecular modelling analyses based on semi-empirical and DFT calculations show that CF and its metabolites would differ in their LUMO-HOMO energy differences such that DFC-D, CF and DFC would be more kinetically labile. The high kinetic lability and the presence of electron-deficient regions on the molecular surface mean that DFC-D, CF and DFC may react more readily with glutathione and nucleobases in DNA resulting into glutathione depletion (thus causing oxidative stress) and DNA damage, respectively. The effects of such adverse reactions may however be more significant in the case of CF and DFC than the more reactive metabolite DFC-D as DFC-D as it is likely to be more readily excreted because of much greater solubility in water.
Abbreviations | |
CF: | Eftiofur; (6R,7R)-((2-amino-4-thiazolyl)-Z-(methoxyimino)acetyl) amino-3-((2-+++furanylcarbonyl) thio) methyl-8-oxo-5-thia-1-azabicyclo (4.2.0)oct-2-ene-2-carboxylic acid |
DFC-CYS: | Desfuroylceftiofur cysteine disulfide |
LUMO: | Lowest unoccupied molecular orbital |
HOMO: | Highest occupied molecular orbital |
Fazlul Huq is grateful to the School of Biomedical Sciences, The University of Sydney for the time release from teaching.