Proximity vs. strain in intramolecular ring-closing reactions

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.

Semisynthesis and Inhibitory Effects of Solidagenone Derivatives on TLR-Mediated Inflammatory Responses

A series of nine derivatives (2⁻10) were prepared from the diterpene solidagenone (1) and their structures were elucidated by means of spectroscopic studies. Their ability to inhibit inflammatory responses elicited in peritoneal macrophages by TLR ligands was investigated. Compounds 5 and 6 showed significant anti-inflammatory effects, as they inhibited the protein expression of nitric oxide synthase (NOS-2), cyclooxygenase-2 (COX-2), and cytokine production (TNF-α, IL-6, and IL-12) induced by the ligand of TLR4, lipopolysaccharide (LPS), acting at the transcriptional level. Some structure⁻activity relationships were outlined. Compound 5 was selected as a representative compound and molecular mechanisms involved in its biological activity were investigated. Inhibition of NF-κB and p38 signaling seems to be involved in the mechanism of action of compound 5. In addition, this compound also inhibited inflammatory responses mediated by ligands of TLR2 and TLR3 receptors. To rationalize the obtained results, molecular docking and molecular dynamic studies were carried out on TLR4. All these data indicate that solidagenone derivative 5 might be used for the design of new anti-inflammatory agents.

Accessing pyrrolones and pyridinones: controlling 5-exo and 6-endo ring closures in heterocyclic alkynylamides

Competitive 5-exo and 6-endo anionic intramolecular cyclization reactions in heterocyclic alkynylamides were explored via experimental and computational analysis. The 5-exo-dig cyclization pathway is usually disfavoured in heterocyclic systems, and 6-endo products are often both the kinetic and thermodynamic products. However, we’ve found that it is possible to shift selectivity toward the 5-exo-dig pyrrolone products away from the less strained pyridinone products that are produced via the 6-endo-dig cyclization. Parameters such as identity of heteroatom, heteroatom positioning within the heterocycle, and functionality on the alkyne were investigated and, in many cases, were found to strongly influence product ratios. A series of computational studies was performed to provide further insight into the 5-exo-dig and 6-endo-dig pathways in these heterocyclic systems. Theoretical predictions were found to reproduce experimental results, highlighting the predictive capabilities of the computations in determining preferred products.

Design, synthesis, and in vitro kinetics study of atenolol prodrugs for the use in aqueous formulations

Based on DFT, MP2, and the density functional from Truhlar group (hybrid GGA: MPW1k) calculations for an acid-catalyzed hydrolysis of nine Kirby's N-alkylmaleamic acids and two atenolol prodrugs were designed. The calculations demonstrated that the amide bond cleavage is due to intramolecular nucleophilic catalysis by the adjacent carboxylic acid group and the rate-limiting step is determined based on the nature of the amine leaving group. In addition, a linear correlation of the calculated and experimental rate values has drawn credible basis for designing atenolol prodrugs that are bitterless, are stable in neutral aqueous solutions, and have the potential to release the parent drug in a sustained release manner. For example, based on the calculated B3LYP/6-31 G (d,p) rates, the predicted t_1/2 (a time needed for 50% of the prodrug to be converted into drug) values for atenolol prodrugs ProD 1-ProD 2 at pH 2 were 65.3 hours (6.3 hours as calculated by GGA: MPW1K) and 11.8 minutes, respectively. In vitro kinetic study of atenolol prodrug ProD 1 demonstrated that the t_1/2 was largely affected by the pH of the medium. The determined t_1/2 values in 1N HCl, buffer pH 2, and buffer pH 5 were 2.53, 3.82, and 133 hours, respectively.

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).

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).