Cleavage of Menger’s aliphatic amide: A model for peptidase enzyme solely explained by proximity orientation in intramolecular proton transfer

Enzyme Models-From Catalysis to Prodrugs

Enzymes are highly specific biological catalysts that accelerate the rate of chemical reactions within the cell. Our knowledge of how enzymes work remains incomplete. Computational methodologies such as molecular mechanics (MM) and quantum mechanical (QM) methods play an important role in elucidating the detailed mechanisms of enzymatic reactions where experimental research measurements are not possible. Theories invoked by a variety of scientists indicate that enzymes work as structural scaffolds that serve to bring together and orient the reactants so that the reaction can proceed with minimum energy. Enzyme models can be utilized for mimicking enzyme catalysis and the development of novel prodrugs. Prodrugs are used to enhance the pharmacokinetics of drugs; classical prodrug approaches focus on alternating the physicochemical properties, while chemical modern approaches are based on the knowledge gained from the chemistry of enzyme models and correlations between experimental and calculated rate values of intramolecular processes (enzyme models). A large number of prodrugs have been designed and developed to improve the effectiveness and pharmacokinetics of commonly used drugs, such as anti-Parkinson (dopamine), antiviral (acyclovir), antimalarial (atovaquone), anticancer (azanucleosides), antifibrinolytic (tranexamic acid), antihyperlipidemia (statins), vasoconstrictors (phenylephrine), antihypertension (atenolol), antibacterial agents (amoxicillin, cephalexin, and cefuroxime axetil), paracetamol, and guaifenesin. This article describes the works done on enzyme models and the computational methods used to understand enzyme catalysis and to help in the development of efficient prodrugs.

Design, synthesis and in vitro kinetic study of tranexamic acid prodrugs for the treatment of bleeding conditions

Based on density functional theory (DFT) calculations for the acid-catalyzed hydrolysis of several maleamic acid amide derivatives four tranexamic acid prodrugs were designed. The DFT results on the acid catalyzed hydrolysis revealed that the reaction rate-limiting step is determined on the nature of the amine leaving group. When the amine leaving group was a primary amine or tranexamic acid moiety, the tetrahedral intermediate collapse was the rate-limiting step, whereas in the cases by which the amine leaving group was aciclovir or cefuroxime the rate-limiting step was the tetrahedral intermediate formation. The linear correlation between the calculated DFT and experimental rates for N-methylmaleamic acids 1-7 provided a credible basis for designing tranexamic acid prodrugs that have the potential to release the parent drug in a sustained release fashion. For example, based on the calculated B3LYP/6-31G(d,p) rates the predicted t1/2 (a time needed for 50 % of the prodrug to be converted into drug) values for tranexamic acid prodrugs ProD 1-ProD 4 at pH 2 were 556 h [50.5 h as calculated by B3LYP/311+G(d,p)] and 6.2 h as calculated by GGA: MPW1K), 253 h, 70 s and 1.7 h, respectively. Kinetic study on the interconversion of the newly synthesized tranexamic acid prodrug ProD 1 revealed that the t1/2 for its conversion to the parent drug was largely affected by the pH of the medium. The experimental t1/2 values in 1 N HCl, buffer pH 2 and buffer pH 5 were 54 min, 23.9 and 270 h, respectively.

Computationally designed prodrugs of statins based on Kirby's enzyme model

DFT calculations at B3LYP/6-31G(d,p) for intramolecular proton transfer in Kirby's enzyme models 1-7 demonstrated that the reaction rate is dependent on the distance between the two reacting centers, rGM, and the hydrogen bonding angle, α, and the rate of the reaction is linearly correlated with rGM and α. Based on these calculation results three simvastatin prodrugs were designed with the potential to provide simvastatin with higher bioavailability. For example, based on the calculated log EM for the three proposed prodrugs, the interconversion of simvastatin prodrug ProD 3 to simvastatin is predicted to be about 10 times faster than that of either simvastatin prodrug ProD 1 or simvastatin ProD 2. Hence, the rate by which the prodrug releases the statin drug can be determined according to the structural features of the promoiety (Kirby's enzyme model).

Targeted prodrugs in oral drug delivery: the modern molecular biopharmaceutical approach

INTRODUCTION

The molecular revolution greatly impacted the field of drug design and delivery in general, and the utilization of the prodrug approach in particular. The increasing understanding of membrane transporters has promoted a novel 'targeted-prodrug' approach utilizing carrier-mediated transport to increase intestinal permeability, as well as specific enzymes to promote activation to the parent drug.

AREAS COVERED

This article provides the reader with a concise overview of this modern approach to prodrug design. Targeting the oligopeptide transporter PEPT1 for absorption and the serine hydrolase valacyclovirase for activation will be presented as examples for the successful utilization of this approach. Additionally, the use of computational approaches, such as DFT and ab initio molecular orbital methods, in modern prodrugs design will be discussed.

EXPERT OPINION

Overall, in the coming years, more and more information will undoubtedly become available regarding intestinal transporters and potential enzymes that may be exploited for the targeted modern prodrug approach. Hence, the concept of prodrug design can no longer be viewed as merely a chemical modification to solve problems associated with parent compounds. Rather, it opens promising opportunities for precise and efficient drug delivery, as well as enhancement of treatment options and therapeutic efficacy.

Prodrugs of acyclovir--a computational approach

Density functional theory calculation results demonstrated that the efficiency of the acid-catalyzed hydrolysis of Kirby's acid amides 1-15 is strongly dependent on the substitution on the C-C double bond and the nature of the amide N-alkyl group. Further, the results established that while in the gas phase the hydrolysis rate-limiting step is the tetrahedral intermediate formation in polar solvents such as water, the rate-limiting step could be either the formation or the collapse of the tetrahedral intermediate depending on the substitution on the C-C double bond and on the amide nitrogen substituent. Based on a linear correlation between the calculated and experimental effective molarities, the study on the systems reported herein could provide a good basis for designing prodrug systems that are less hydrophilic than their parental drugs and can be used, in different dosage forms, to release the parent drug in a controlled manner. For example, based on the calculated log effective molarities values, the predicted t(1/2) (a time needed for 50% of the reactant to be hydrolyzed to products) for acyclovir prodrugs, ProD 1-4, was 29.2 h, 6097 days, 4.6 min, and 8.34 h, respectively. Hence, the rate by which acyclovir prodrug releases acyclovir can be determined according to the structural features of the linker (Kirby's acid amide moiety).

Computational-aided design for dopamine prodrugs based on novel chemical approach

The goal of this project was to design novel dopamine prodrugs for the treatment of Parkinson's disease that can improve the overall biopharmaceutical profile of the current medications to enhance effectiveness and to ease the use of the medications. Based on Menger's study and other's kinetic studies on the cleavage reactions of some Kemp's acid amides to the corresponding amines and anhydrides, DFT calculations were made to find a candidate to be used as an efficient dopamine prodrug. The proposed dopamine prodrugs have a carboxylic group as a hydrophilic moiety and a hydrocarbon skeleton as a lipophilic moiety, where the combination of both groups ensures a moderate hydrophilic lipophilc balance value. The potential prodrugs are expected to give better bioavailability than the parental drug owing to improved absorption. Furthermore, the proposed prodrugs are believed to be more effective than l-dopa because the latter undergoes decarboxylation in the periphery before reaching the blood-brain barrier. Additionally, the predicted prodrugs can be used in different dosage forms (I.V., S.C., tablets, and others) because of their potential solubility in organic and aqueous media. It is expected that the proposed prodrugs will undergo cleavage reactions to release the parental dopamine drug with half-life time (t(½) ) of 12-20 h.

Computer-assisted design for atenolol prodrugs for the use in aqueous formulations

Based on stability studies on the drugs atenolol and propranolol and some of their derivatives it is believed that increasing the lipophilicity of the drug will lead to an increase in the stability of its aqueous solutions and will provide a prodrug system with the potential for releasing atenolol in a controlled manner. Using DFT theoretical calculations we have calculated an intramolecular acid catalyzed hydrolysis in nine maleamic (4-amino-4-oxo-2butenoic) acids (Kirby's N-alkylmaleamic acids), 1-9. The DFT calculations confirmed that the acid-catalyzed hydrolysis mechanism in these systems involves: (1) a proton transfer from the hydroxyl of the carboxyl group to the adjacent amide carbonyl carbon, (2) an approach of the carboxylate anion toward the protonated amide carbonyl carbon to form a tetrahedral intermediate; and (3) a collapse of the tetrahedral intermediate into products. Furthermore, DFT calculations in different media revealed that the reaction rate-limiting step depends on the reaction medium. In aqueous medium the rate-limiting step is the collapse of the tetrahedral intermediate whereas in the gas phase the formation of the tetrahedral intermediate is the rate-limiting step. Furthermore, the calculations establish that the acid-catalyzed hydrolysis efficiency is largely sensitive to the pattern of substitution on the carbon-carbon double bond. Based on the experimental t(1/2) (the time needed for the conversion of 50% of the reactants to products) and EM (effective molarity) values for processes 1-9 we have calculated the t(1/2) values for the conversion of the two prodrugs to the parental drug, atenolol. The calculated t(1/2) values for ProD 1-2 are predicted to be 65.3 hours and 11.8 minutes, respectively. Thus, the rate by which atenolol prodrug undergoes cleavage to release atenolol can be determined according to the nature of the linker of the prodrug (Kirby's N-alkylmaleamic acids 1-9).

Computer-assisted design for paracetamol masking bitter taste prodrugs

It is believed that the bitter taste of paracetamol, a pain killer drug, is due to its hydroxyl group. Hence, it is expected that blocking the hydroxy group with a suitable linker could inhibit the interaction of paracetamol with its bitter taste receptor/s and hence masking its bitterness. Using DFT theoretical calculations we calculated proton transfers in ten different Kirby's enzyme models, 1-10. The calculation results revealed that the reaction rate is linearly correlated with the distance between the two reactive centers (r(GM)) and the angle of the hydrogen bonding (α) formed along the reaction pathway. Based on these results three novel tasteless paracetamol prodrugs were designed and the thermodynamic and kinetic parameters for their proton transfers were calculated. Based on the experimental t(1/2) (the time needed for the conversion of 50% of the reactants to products) and EM (effective molarity) values for processes 1-10 we have calculated the t(1/2) values for the conversion of the three prodrugs to the parental drug, paracetamol. The calculated t(1/2) values for ProD 1-3 were found to be 21.3 hours, 4.7 hours and 8 minutes, respectively. Thus, the rate by which the paracetamol prodrug undergoes cleavage to release paracetamol can be determined according to the nature of the linker of the prodrug (Kirby's enzyme model 1-10). Further, blocking the phenolic hydroxyl group by a linker moiety is believed to hinder the paracetamol bitterness.

Prodrugs of aza nucleosides based on proton transfer reaction

DFT calculation results for intramolecular proton transfer reactions in Kirby's enzyme models 1-7 reveal that the reaction rate is quite responsive to geometric disposition, especially to distance between the two reactive centers, r (GM), and the angle of attack, α (the hydrogen bonding angle). Hence, the study on the systems reported herein could provide a good basis for designing aza nucleoside prodrug systems that are less hydrophilic than their parental drugs and can be used, in different dosage forms, to release the parent drug in a controlled manner. For example, based on the calculated log EM, the cleavage process for prodrug 1ProD is predicted to be about 10¹⁰ times faster than that for prodrug 7ProD and about 10⁴ times faster than prodrug 3ProD: rate( 1ProD ) > rate( 3ProD ) > rate( 7ProD ). Hence, the rate by which the prodrug releases the aza nucleoside drug can be determined according to the structural features of the linker (Kirby's enzyme model).

A computational analysis of intramolecularity in proton transfer reactions

Hydrogen bonds along which most proton transfers take place are analyzed regarding their contribution to the efficiency of intramolecular processes through the assessment of effective molarities by DFT calculations on rigid intra-molecular model systems in comparison with their intermolecular counterparts. A few cyclic systems in which there is no possibility of hydrogen bonding between reactants neither in the ground state nor in the product are identified by theoretical calculations as leading to important rate increases in intramolecular general-base catalysis constituting a new exception to the rule of low effective molarities (EMs) for proton transfers. The recovery of high efficiency in these systems is attributed to the absence of special features facilitating the corresponding intermolecular reactions. This work confirms that any explanation of EMs found in intramolecular systems also requires a careful analysis of the corresponding bimolecular systems.

Computer-assisted design of pro-drugs for antimalarial atovaquone

Density Functional Theory (DFT) and ab initio calculation results for the proton transfer reaction in Kirby's enzyme models 1-6 reveal that the reaction rate is largely dependent on the existence of a hydrogen bonding net in the reactants and the corresponding transition states. Further, the distance between the two reacting centers and the angle of the hydrogen bonding formed along the reaction path has profound effects on the rate. Hence, the study on the systems reported herein could provide a good basis for designing antimalarial (atovaquone) pro-drug systems that can be used to release the parent drug in a controlled manner. For example, based on the calculated log EM, the cleavage process for pro-drug 1Pro may be predicted to be about 10¹¹ times faster than that for a pro-drug 4Pro and about 10⁴ times faster than pro-drug 2Pro: rate (1Pro) > rate (2Pro > rate (4Pro). Thus, the rate by which the pro-drug releases the antimalarial drug can be determined according to the nature of the linker (Kirby's enzyme model 1-6).

The effective molarity (EM)--a computational approach

The effective molarity (EM) for 12 intramolecular S(N)2 processes involving the formation of substituted aziridines and substituted epoxides were computed using ab initio and DFT calculation methods. Strong correlation was found between the calculated effective molarity and the experimentally determined values. This result could open a door for obtaining EM values for intramolecular processes that are difficult to be experimentally provided. Furthermore, the calculation results reveal that the driving forces for ring-closing reactions in the two different systems are proximity orientation of the nucleophile to the electrophile and the ground strain energies of the products and the reactants.