- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
The Scientific World Journal
Volume 2014 (2014), Article ID 248651, 13 pages
Design, Synthesis, and In Vitro Kinetics Study of Atenolol Prodrugs for the Use in Aqueous Formulations
1Faculty of Pharmacy, Al-Quds University, P.O. Box 20002, Jerusalem, Palestine
2Department of Sciences, University of Basilicata, Via dell’Ateneo Lucano 10, 85100 Potenza, Italy
3Faculty of Public Health Sciences, Al-Quds University, P.O. Box 20002, Jerusalem, Palestine
Received 10 October 2013; Accepted 1 December 2013; Published 12 January 2014
Academic Editors: R. I. Cukier, E. Gomez-Bengoa, R. Luo, D. Quintanar-Guerrero, and T. Takayanagi
Copyright © 2014 Rafik Karaman et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
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 (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 was largely affected by the pH of the medium. The determined values in 1N HCl, buffer pH 2, and buffer pH 5 were 2.53, 3.82, and 133 hours, respectively.
Atenolol, 4-[2-hydroxy-3-[(1-methylethyl)amino]propoxy] benzene acetamide, is a relatively polar hydrophilic compound with a water solubility of 26.5 mg/mL at 37°C and a log partition coefficient (octanol/water) of 0.23. Lipid insoluble hydrophilic compounds such as atenolol, sotalol, and nadolol are excreted only by the kidneys and have low brain penetration. Atenolol is a selective -adrenoceptor antagonist, applied in the treatment of numerous cardiovascular disorders including hypertension, angina, acute myocardial infarction, supraventricular tachycardia, ventricular tachycardia, and the symptoms of alcohol withdrawal via restricting certain nerve impulses, thereby controlling the rate and force of contraction and consequently reducing blood pressure in addition to its treatment of angina pectoris. On the outset, atenolol actively reduces the heart rate by decreasing systolic and diastolic blood pressure. The net effect of atenolol on controlling both the heart rate and blood pressure is the reduction in myocardial work and oxygen requirement which reduces cardiovascular stress, thereby preventing arrhythmia and angina attacks.
Atenolol is marketed as tablets of 25 mg, 50 mg, and 100 mg and an injectable formulation of 5 mg/mL supplied in 10 mL ampoules [1, 2]. Atenolol has a of 9.6; it undergoes ionization in the stomach and intestine, thus its oral bioavailability is low due to inefficient absorption through membranes. The bioavailability of atenolol is between 45% and 55% of the given dose and is not increased by administration of the drug in a solution form [3–5]. About 50% of administered atenolol is absorbed; however, most of the absorbed quantity reaches the systemic circulation. Atenolol peak blood levels are reached within two to four hours after ingestion. Differently from propranolol or metoprolol, atenolol is resistant to metabolism by the liver and the absorbed dose is eliminated by renal excretion. More than 85% of IV dose is excreted in urine within 24 hours compared with 50% for an oral dose. Only 6–16% is protein bound resulting in relatively consistent plasma drug levels with about a fourfold interpatient variation. The elimination half-life of atenolol is between 6 and 7 hours and there is no alteration of kinetic profile of a drug by chronic administration. Atenolol is one of the most important medicines used for prevention of several types of arrhythmias in childhood, but unfortunately it is still unlicensed . On the other hand, atenolol is indicated as a first-step therapy for hypertension in elderly patients, who have difficulty in swallowing and, thus, tablets and capsules are frequently avoided. The ease of administration makes a liquid formulation an ideal dosage form for such patients . Therefore, extemporaneous compounding (off-label), involving preparation of an oral liquid from a pure drug powder, is required. However, formulations compounded from tablets and pure active drug suffer instability and are only stable for less than one week [8–10]. Furthermore, atenolol bitterness is considered as a great challenge to health sector when used among children and geriatrics . The main problem in oral administration of bitter drugs such as atenolol is incompliance by the patients  and this can be overcome by masking the bitterness of a drug either by decreasing its oral solubility on ingestion or eliminating the interaction of drug particles to taste buds .
There are various techniques available which are commonly used for masking drug’s bitterness: (1) taste masking with flavors, sweeteners, and amino acids ; (2) taste masking with lipophilic vehicles such as lipids, lecithin, and lecithin-like substances ; (3) coating which is classified based on the type of coating material, coating solvent system, and the number of coating layers ; (4) microencapsulation based on the principle of solvent extraction or evaporation ; (5) sweeteners are generally used in combination with other taste masking technologies ; (6) taste suppressants and potentiators, such as Linguagen’s bitter blockers (e.g., adenosine monophosphate), are used for masking the bitter taste of various compounds by competing with binding to the G-protein coupled receptor sites (GPCR) ; (7) resins are utilized to mask pharmaceuticals bitterness by forming insoluble resonates [20, 21]; (8) inclusion complex by which the drug molecule fits into the cavity of a complexing agent and forms a stable complex that masks the bitter taste of a drug by decreasing its oral solubility ; (9) pH modifiers are capable of generating a specific pH microenvironment in aqueous media that has the ability to facilitate in situ precipitation of the bitter drug compound in saliva thus reducing the overall taste sensation for liquid dosage forms ; (10) adsorbates, by which the compound may be adsorbed or entrapped in the matrix of the adsorbate pore, which may result in a delayed release of the bitter tastant during passage through the oral cavity and mask the taste ; and (11) prodrug approach by which the functional group/s bind to the bitter taste receptor is blocked by a promoiety .
Previous studies on stability of atenolol ester prodrugs for the use in transdermal preparations have shown that these ester derivatives are much more stable than the corresponding alcohol, atenolol, when they are formulated in aqueous solutions [26, 27]. On the other hand, the only atenolol prodrug intended to be used for oral dosage form was atenolol aspirinate prodrug (aspirin); it is described for antihypertensive therapy to reduce cardiovascular death, stroke, and myocardial infarction (MI); however, recent studies reported that the coupling of atenolol with acetyl salicylic acid by means of an ester linkage did not produce any efficient pharmacological profile, neither in vitro nor in vivo .
Many therapeutic drugs have undesirable properties that may become pharmacological, pharmaceutical, or pharmacokinetic barriers in clinical drug application; one of the most promising strategies for precise and efficient drug delivery and enhancement of a drug therapeutic effect is the prodrug design approach which is becoming more elaborate in the development of efficient and selective drug delivery systems . Approximately, tenth of all worldwide marketed medicines can be categorized as prodrugs, and, in 2008 alone, third of all approved small-molecular-weight drugs were prodrug . This data signifies the importance of the prodrug approach and encourages the scientific community to consider it in the early stages of preclinical development and not as a last resort .
Prodrugs are substances administered in an inactive form that is then metabolized in the body into the corresponding active drug. The rationale behind administering prodrugs is to optimize the drug absorption, distribution, metabolism, and excretion. Since first described in the 1950s, prodrugs continue to be a fertile area of research. There are a number of small pharmaceutical/biotech companies dedicated to use prodrugs for the delivery of older but problematic drugs as well as to develop broad-based prodrug technologies for application to new and future drugs .
The use of the term prodrug usually implies a covalent link between a drug and a chemical promoiety. Generally, prodrugs can be enzymatically or nonenzymatically converted in vivo to provide the parent active drug which exerts the desired therapeutic effect. Ideally, the prodrug should be converted to the parent drug as soon as its goal is achieved, followed by the subsequent rapid elimination of the released promoiety group [33, 34].
The novel prodrug design approach that was utilized for the synthesis and kinetic studies described in this paper is based on intramolecular process using molecular orbital methods and correlations between experimental and calculated values using density functional theory (DFT)  and molecular mechanics (MM)  methods which are increasingly being exploited and widely recognized as tools for providing structure-energy calculations for predictions of potential drugs and prodrugs . In this approach, no enzyme is needed for the catalysis of the intraconversion of a prodrug to its parent drug. The rate of the drug release is dependent solely on the rate-limiting step of the intraconversion process. Therefore, the reaction mechanism must be investigated in order to assign the factors playing role in the reaction rate. This information is then utilized to design an efficient chemical device to be exploited as a prodrug linker capable of liberating the parent drug in a programmable manner (slow or fast release).
As a result, the intended prodrugs can be designed to undergo cleavage reactions in physiological environments such as stomach (pH 1-2), intestine (pH 5.5-7), and/or blood circulation (pH 7.4), with rates that are solely dependent on the structural features of the pharmacologically inactive promoiety.
Computational chemistry methods based on quantum molecular (QM) and molecular mechanics (MM) theories could be useful for the design of innovative prodrugs for hydroxyl, phenol, or amine containing drugs. For instance, mechanisms of intramolecular processes for a respected number of enzyme models that were previously investigated by others to understand enzyme catalysis have been recently explored by us and exploited for a design of some new novel prodrugs. Using DFT, molecular mechanics, and ab initio at different levels, numerous enzyme model processes were calculated for determining the factors governed their rates. These processes include (i) intramolecular acid-catalyzed hydrolysis of N-alkylmaleamic acid derivatives [38, 39]; (ii) proton transfers between two oxygen atoms and between oxygen and nitrogen in Kirby’s acetals [40–48]; (iii) proton transfer between two oxygen atoms in nonflexible (rigid) systems as studied by Menger and Ladika [49–53]; (iv) acid-catalyzed ring-closing reactions of hydroxy-acids as investigated by Milstien and Cohen [54–57] and Menger and Ladika [49–53]; and (v) SN2-based cyclization reactions as researched by Brown and van Gulick , Bruice and Pandit [59, 60], and Galli and Mandolini .
These studies have demonstrated that there is a must to further explore the mechanisms for the above-mentioned processes for assigning the factors affecting the nature and the mode of the reaction. Unraveling the reaction mechanism would allow for an accurate design of an efficient chemical device to be utilized as a prodrug promoiety that can be covalently linked to a parent drug to provide chemically and not enzymatically the parent drug in a programmable manner upon exposure to physiological environments. For example, exploring the mechanism for proton transfer in Kirby’s acetals has led to design and synthesis of novel prodrugs of aza-nucleosides to treat myelodysplastic syndromes , statins to treat high cholesterol blood levels , bitterless paracetamol prodrugs to be administered to children and elderly as antipyretic and pain killer , and prodrugs of phenylephrine as decongestants . In the above-mentioned examples, the prodrug moiety was attached to the hydroxyl group of the active drug such that the drug promoiety (prodrug) has a potential to degrade upon exposure to physiological environment such as stomach, intestine, and/or blood circulation, with rates that are solely dependent on the structural features of the pharmacologically inactive promoiety (Kirby’s enzyme model). Other different linkers such as Kirby’s N-alkylmaleamic acids (enzyme model) were also investigated for the design of some prodrugs such as those of tranexamic acid to treat bleeding conditions  and acyclovir (antiviral drug) to treat Herpes Simplex . Further, prodrugs for masking the bitterness of antibacterial drugs such as cefuroxime were designed and made as well . The role of the promoiety in the antibacterial (cefuroxime) and paracetamol prodrugs was to block the free amine (cefuroxime) or phenol (paracetamol), which is believed to be responsible for the drug bitterness, and to enable the release of the drug in a programmable manner. Menger’s Kemp acid enzyme model was also exploited for the design of dopamine prodrugs for the treatment of Parkinson’s disease .
The goals of this work were (1) design of atenolol prodrugs that can be (a) formulated in aqueous solutions and maintain stability over a long period of time and (b) bitterless and capable of undergoing intraconversion in physiological environment to provide the parent active drug, atenolol, in a programmable manner and (2) synthesis, characterization, and in vitro kinetic study of the interconversion of the designed prodrugs in different media: 1 N HCl and at buffers of pH 2, pH 5, and pH 7.4.
2. Materials and Methods
The melting points were determined in open capillaries on electrothermal stuart SMP3 advanced melting point apparatus; IR spectra were obtained from a KBr matrix (4000–400 cm−1) using a Perkin-Elmer, Spectrum 100, FT-IR spectrometer.
The LC/ESI-MS/MS system used was Agilent 1200 series liquid chromatography coupled with a 6520 accurate mass quadruple-time of flight mass spectrometer (Q-TOF LC/MS). The analysis was performed in the positive electrospray ionization mode. The m/z range for MS scans was 100–1600 Da. The high pressure liquid chromatography (HPLC) system consisted of an Alliance 2695 module coupled with 2996 Photodiode array detector from Waters (Germany). Data acquisition and control were carried out using Empower 2 software (Waters, Germany). Analytes were separated on a 4.6 mm × 150 mm XBridge C18 column (5 μm particle size) used in conjunction with a 4.6 × 20 mm, XBridge C18 guard column. Microfilters of 0.45 μm porosity were normally used (Acrodisc GHP, Waters).
1H-NMR experiments were performed with a Bruker AvanceII 500 spectrometer equipped with a 5 mm BBO probe. pH values were recorded on pH meter model HM-30G: TOA electronics; thin-layer chromatography (TLC) was carried out on TLC plastic sheets silica gel, 20 × 20 cm, layer thickness 0.2 mm, and the spots on the plates were localized by exposure to UV light.
2.2. Chemicals and Reagents
Pure atenolol standard was commercially available; maleic anhydride, anhydrous sodium dihydrogen phosphate, sodium hydroxide, concentrated hydrochloric acid, and anhydrous magnesium sulfate were commercially obtained from sigma Aldrich (Israel). HPLC grade solvents of methanol, acetonitrile, and water were purchased from Sigma Aldrich (Israel). High purity chloroform, DMF, and diethyl ether (>99%) were purchased from Biolab (Israel).
2.3. Preparation of Atenolol ProD 1
Synthesis of the atenolol ProD 1 was accomplished using Kirby’s procedure  (Figure 1); in a 250 mL round bottom flask equipped with a stirrer and a condenser, standard atenolol (10 mmol) was dissolved in DMF (10 mL), then maleic anhydride (10 mmol) was added, the reaction mixture was refluxed for four hours and cooled and the DMF was removed by evaporation under vacuum, the resulting residue was dissolved in 20 mL 1 N NaOH, extracted with 3 × 30 mL of ether, the organic layer was dried over anhydrous MgSO4, and the dry residue was crystallized (yield 81%). The melting point of the resulting white precipitate was 105°C, and UV () was 276 nm. 1H NMR δ (ppm) (500 MHz, MeOD) 1.16–1.23 (d, 6H, CH3CHCH3), 3.45–3.47 (d, 2H, Hz, NH2-CO-CH2–), 3.49-3.50 (m, 2H, N–CH2–) 3.52–3.55 (m, 1H, –NCH(CH3)2), 4.00-4.01 (m, 1H, CH2CH(OH)CH2), 4.24–4.27 (m, 1H, –O–CH2–), 4.41–4.43 (m, 1H, –OCH2–), 5.96 (s, 2H, –NH2), 6.04–6.07 (d, 1H, Hz, –CH=CH–), 6.28–6.32 (d, 1H, Hz, –CH=CH–), 6.88–6.93 (m, 2H, aromatic), 7.18–7.22 (m, 2H, aromatic). IR (KBr/ cm−1) 1712 (C=O), 1643 (C=O), 1587 (C=C), 1420, 1280, 1200, 1132, 1058. LC-ESi-MS (positive mode) m/z (M+1)+ found 365.1705 Da, calculated 365.1712 Da.
2.4. Kinetic Method
2.4.1. Buffer Preparation
Buffer pH 2 : 6.8 g potassium dihydrogen phosphate were dissolved in 900 mL water for HPLC and the pH was adjusted by diluted o-phosphoric acid and water was added to a final volume of 1000 mL (0.05 M). The same procedure was done for the preparation of buffers pH 5 and 7.4; however, in these two cases, the pH was adjusted using 1 N NaOH.
2.4.2. Calibration Curve for Atenolol and Atenolol ProD 1
To construct a calibration curve for atenolol prodrug and the parent drug, atenolol, several concentrations of 10, 50, 100, 200, 300, and 500 ppm were prepared. All samples were injected into HPLC-PDA. The optimal HPLC conditions used for the analysis of atenolol were 4.6 mm × 150 mm, 5 μm, XBridge C18 column, a mobile phase consisting of binary solvent mixture of water: acetonitrile (water pH adjusted to 3.4 using dilute o-phosphoric acid) (55 : 45 v/v), a flow rate of 0.5 mL/minute, and a UV detection at a wavelength of 230 nm. Peak area versus concentration of the pharmaceutical (ppm) was then plotted, and the value of the plot was recorded.
2.4.3. Preparation of Standard and Sample Solutions
500 ppm of standard atenolol was prepared by dissolving 50 mg of the drug in 100 mL 1 N HCL, buffer pH 3, buffer pH 5, or buffer pH 7.4. Ten μL of the sample was injected into the HPLC to detect the retention time of atenolol. The same procedure was followed for the preparation of 500 ppm of atenolol prodrug.
2.4.4. Hydrolysis of Atenolol ProD 1
Atenolol ProD 1 rate hydrolysis was studied at 37°C in buffer solutions at different pHs (1 N HCl, pH 2, pH 5, and pH 7.4); samples of the reaction mixtures were analyzed directly by HPLC; quantitative analysis of atenolol prodrug and its hydrolysis product, atenolol, was obtained using a calibration curve.
2.5. Calculations Methods
The Becke three-parameter, hybrid functional combined with the Lee, Yang, and Parr correlation functional, denoted by B3LYP, were employed in the calculations using density functional theory (DFT). The DFT calculations at B3LYP/6-31G (d,p) and B3LYP/311+G (d,p) levels, MP2 calculations, and the density functional from Truhlar group (hybrid GGA: MPW1k) [69–71] were carried out using the quantum chemical package Gaussian-2009 . Calculations were carried out based on the restricted Hartree-Fock method . The starting geometries of all calculated molecules were obtained using the Argus Lab program  and were initially optimized at the HF/6-31G level of theory, followed by optimization at the B3LYP/6-31G (d,p), B3LYP/6-311+G (d,p), MPW1k (mpwpw91/6-31+G (d,p)), or MP2/6-31G (d,p) levels. Total geometry optimizations included all internal rotations. Second derivatives were estimated for all 3N-6 geometrical parameters during optimization. In the DFT calculations for 1–9 and atenolol ProD 1-2, two types of conformations in particular were considered: one in which the amide carbonyl is syn to the carboxyl group and another in which it is anti. An energy minimum (a stable compound or a reactive intermediate) has no negative vibrational force constant. A transition state is a saddle point which has only one negative vibrational force constant . Transition states were located first by the normal reaction coordinate method  where the enthalpy changes were monitored by stepwise changing the interatomic distance between two specific atoms. The geometry at the highest point on the energy profile was reoptimized by using the energy gradient method at the B3LYP/6-31G (d,p) level of theory . The “reaction coordinate method” was used to calculate the activation energy in maleamic (4-amino-4-oxo-2butenoic) acids (Kirby’s N-alkylmaleamic acids), 1–9, and atenolol ProD 1-2. In this method, one bond length is constrained for the appropriate degree of freedom while all other variables are freely optimized. The activation energy values for the approach and dissociation steps were calculated from the difference in energies of the global minimum structures (GM) and the derived transition states. Verification of the desired reactants and products was accomplished using the “intrinsic coordinate method” . The transition state structures were verified by their only one negative frequency. Full optimization of the transition states was accomplished after removing any constrains imposed while executing the energy profile. The activation energies obtained from DFT at B3LYP/6-31G (d,p) level of theory for 1–9 and atenolol ProD 1-2 were calculated with and without the inclusion of solvent (water and ether). The calculations with the incorporation of a solvent were performed using the integral equation formalism model of the Polarizable Continuum Model (PCM) [76–79].
3. Results and Discussion
Our proposed atenolol prodrugs that were designed based on the acid-catalyzed hydrolysis reactions of N-alkyl maleamic acids 1–9 (Figure 2) are depicted in Figure 3. As illustrated in Figure 3, the only difference between the proposed atenolol prodrugs and the parent drug, atenolol, is that the amine group in atenolol was replaced with an amide moiety. This chemical change is expected to increase the stability of the alcohol derivative which is formed (prodrug) compared to the corresponding amine alcohol, atenolol, due to general chemical stability for tertiary alcohols over amine alcohols. In addition, stability studies on atenolol ester derivatives have showed the ester derivatives to be much more stable than their corresponding alcohol, atenolol, upon formulating in aqueous solutions. On the other hand, kinetic study of atenolol and propranolol revealed that increasing the lipophilicity of the drug leads to an increase in the stability of its aqueous solutions. Based on that, it is expected that atenolol prodrugs (atenolol amide derivatives) shown in Figure 3 will be more resistant to heat or/oxidation when standing in aqueous solutions [7, 11, 80–83].
The drug’s bitter taste can be masked by utilizing the prodrug chemical approach. For example, paracetamol, a widely used pain killer and fever reducer found in the urine of patients who had taken phenacetin, has a very unpleasant bitter taste . Phenacetin, on the other hand, lacks or has very slight bitterness . The difference in the structural features of both drugs is only in the nature of the group in the para position of the benzene ring. While in the case of paracetamol the group is hydroxyl, in phenacetin, it is ethoxy. Acetanilide has a chemical structure similar to that of paracetamol and phenacetin but lacks the group in the para position of the benzene ring, making it lack the bitter taste characteristic of paracetamol . These combined facts suggest that the presence of the hydroxy group on the para position is the major contributor for the bitter taste of paracetamol. It is believed that paracetamol interacts with the bitter taste receptors via hydrogen bonding which involves its phenolic group. In fact, our recent study on binding of paracetamol and its prodrugs to bitter taste receptors has confirmed that prodrugs lacking the hydroxyl group have no bitterness. Similarly, it is expected that blocking the amine group in atenolol with a suitable linker might inhibit the interaction between the amine group in atenolol and its bitter taste receptors and hence masks its bitterness. The nature of the bitter taste receptors with either paracetamol (via the phenolic group) or atenolol (via the amine group) is likely to be as a result of hydrogen bonding between the substrate and the receptor .
As shown in Figure 3, the proposed atenolol prodrugs, atenolol ProD 1 and atenolol ProD 2, have a hydroxyl and carboxylic acid groups (hydrophilic moiety) and the rest of the prodrug molecule is a lipophilic moiety, where the combination of both groups ensures a molecule with a moderate hydrophilic lipophilic balance (HLB).
It is worthy to note that the HLB value of atenolol prodrug moiety will be largely determined on the medium (physiologic environment) by which the prodrug is dissolved. For instance, in the stomach (pH 1-2), atenolol prodrugs will exist in the free carboxylic acid form, whereas in the blood circulation (pH 7.4) the carboxylate anion form will be dominant. It is planned that atenolol ProD 1-ProD 2 (Figure 3) will be formulated as sodium carboxylate salt since this form is expected to be stable in neutral aqueous medium. However, upon exposure to stomach (pH less than 3), the prodrugs will exist mainly as free carboxylic acid forms thus enabling the acid-catalyzed hydrolysis to proceed.
3.1. Atenolol Prodrugs Design
3.1.1. DFT Mechanistic Study of the Acid-Catalyzed Hydrolysis of 1–9, Atenolol ProD 1, and Atenolol ProD 2
The acid-catalyzed hydrolysis of 1–9 (Figure 2) was extensively studied by Kirby and Lancaster. The study findings demonstrate that the amide bond cleavage is due to intramolecular nucleophilic catalysis by the adjacent carboxylic acid group and the rate-limiting step is the dissociation of the tetrahedral intermediate (Figure 4) . Two decades later, the reaction was theoretically investigated by Katagi using AM1 semiempirical calculations. In contrast to what was proposed by Kirby, the AM1 study revealed that the rate-limiting step is the tetrahedral intermediate formation and not its collapse . Later on, Kluger and Chin have studied the intramolecular hydrolysis mechanism of a number of N-alkylmaleamic acids derived from aliphatic amines having a wide range of basicity . Their findings revealed that the identity of the rate-limiting step is a function of both the basicity of the leaving group and the acidity of the solution.
Based on the above-mentioned experimental findings we believe that replacing the N-methylamide group in 1–9 (Figure 2) with atenolol molecule, as shown for atenolol ProD 1 and ProD 2 in Figure 3, will not have any significant effect on the relative rates of these processes, if they proceed in the same mechanism. Therefore, DFT calculations of the kinetic and thermodynamic properties for processes 1–9 might fairly predict the rates for the chemical intraconversion of atenolol ProD 1 and ProD 2.
3.1.2. The Factors Determining the Reaction Rate-Limiting Step
Our previous calculations of the mechanism for the acid-catalyzed hydrolysis of 1–7 demonstrate that the process involves at first a proton transfer from the carboxylic group into the adjacent amide carbonyl carbon followed by a nucleophilic attack of the carboxylate anion thus formed onto the protonated carbonyl carbon to form a tetrahedral intermediate which in further step dissociates to give products (Figure 4).
In addition, the DFT calculation results were found to be in accordance with the reports by Kirby and Lancaster.  and Kluger and Chin  that indicate that the N-alkylmaleamic acids hydrolysis occurs by one mechanism where the rate-limiting step is largely dependent on the nature of the amine leaving group and the medium solvent. Two different rate-limiting steps were suggested: (1) formation of atetrahedral intermediate  and (2) a tetrahedral intermediate dissociation (Figure 4) [38, 88]. Our gas phase DFT calculated barriers revealed that the rate-limiting step for the hydrolysis of all maleamic acid amides studied regardless of the nature of the amine leaving group is a tetrahedral formation, whereas when the calculations were run in a dielectric constant of 78.39 (water), the rate-limiting step for the hydrolysis of acid amides with primary amine was the dissociation of the tetrahedral intermediate; however, when an acyclovir (antiviral drug) or a cefuroxime (antibacterial drug) moieties (for prodrug structures see Figure 5) were a leaving group, the tetrahedral intermediate formation was found to be involved in the hydrolysis rate-limiting step.
In order to assign the factors creating the unusual accelerations in rate and to explore whether the discrepancy in processes 1–9 rates stems from strain energy or other effects using Allinger’s MM2 method , the strain energies for the reactants (GM), intermediates (INT), and products (P) in 1–9 and atenolol ProD 1-2 were calculated and their values were examined for correlation with the corresponding experimental relative rates, . Strong correlations were obtained with a correlation coefficient . Similarly, fair correlations were obtained between the experimental and calculated free activation energy values and strain energy values of the intermediates and reactants (Es INT-GM) (). The combined correlation results imply that the rate of the reaction for systems having less strained intermediates or products such as 2 and 5 are higher than that having more strained intermediates or products such as 1 and 4. This might be attributed to the fact that the transition state structures resemble that of the corresponding intermediate.
3.1.3. Calculation of the Values for the Cleavage Reactions of Prodrugs ProD 1-2
Intramolecular efficiency is widely accepted to be measured by its effective molarity (EM) value. EM is defined as a ratio of the intramolecular and its corresponding intermolecular rate, where both processes are driven by identical mechanisms. The major factors affecting the EM parameter are the size of a ring, reaction medium, and reaction type. Cyclization processes which proceed via intramolecular nucleophilic addition are much more efficient than intramolecular proton transfers. Values in the order of 109–1013 M have been reported for the EM parameter in intramolecular processes occurring through nucleophilic addition. On the other hand, EM values less than 10 M were measured for proton transfer processes until recently where values of 1010 were reported by Kirby on the hydrolysis of some enzyme models [40–46, 89–92].
For obtaining credibility to our calculation results, we introduce our computation rational for calculating the EM values for processes 1–9 and ProD 1-2 based on the DFT calculated activation energies () of 1–9 and ProD 1-2, and the corresponding intermolecular process Inter (Figure 5 and Table 1).
Using (1)–(4), (5) which describes the EM parameter as a function of the difference in the activation energies of the intramolecular and the corresponding intermolecular processes was derived. The calculated EM values for processes 1–9 and ProD 1-2 were obtained as where is the temperature in Kelvin and is the gas constant.
The calculated log EM values for 1–5 were examined for correlation with the log EM experimental values  (Figure 6(a)). The correlation results revealed that 2 and 5 were the most efficient processes among 1–5, whereas process 4 was the least. The discrepancy in rates between 2 and 5 on one hand and 4 on the other hand is attributed to strain effects.
Figure 6(a) indicates that although the calculated and experimental EM values are comparable, their absolute values slightly differ. This might be due to the fact that the experimental EM values for 1–5 were measured in the presence of aqueous acid, whereas the DFT calculations were run in plain water. The dielectric constant value for a mixture of acid/water is expected to be different from pure water (78.39) and hence the discrepancy in the calculated and experimental EM values.
In addition, for further support to the credibility of our DFT calculations, the calculated free activation energies in water () were correlated with the experimental free activation energies (). Strong correlation was obtained with value of 0.96 (Figure 6(b)).
Using (6) obtained from the correlation of log versus and the value for process 2 ( second)38, the values for ProD 1 and ProD 2 at pH 2 were calculated and their values were 65.3 hours and 11.8 minutes, respectively. Consider
3.1.4. Hydrolysis Studies
The kinetics of the acid-catalyzed hydrolysis study for atenolol ProD 1 was carried out in an aqueous buffer in a similar manner to that done by Kirby on N-alkylmaleamic acids 1–7. This is for investigating whether atenolol prodrug hydrolyzes in aqueous medium and to what extent, suggesting its fate in the system. Acid-catalyzed hydrolysis kinetics of the synthesized atenolol ProD 1 was studied in four different aqueous media: 1 N HCl and buffers pH 2, pH 5, and pH 7.4. Under the experimental conditions, the target prodrug (atenolol ProD 1) was hydrolyzed to release the parent drug, atenolol, (Figure 7) as was evident by HPLC determination. At constant pH and temperature, the reaction displayed strict first-order kinetics as the was fairly constant and a straight line was obtained on plotting log concentration of residual prodrug versus time. The rate constant () and the corresponding half-lives () for atenolol prodrug ProD 1 in the different media were calculated from the linear regression equation correlating the log concentration of the residual prodrug versus time. The kinetic data, and values, are listed in Table 2. One N HCl, pH 2, and pH 5 were selected to examine the interconversion of atenolol ProD 1 in pH as of stomach, because the mean fasting stomach pH of adult is approximately 1-2 and increases up to 5 following ingestion of food. In addition, buffer pH 5 mimics the beginning of the small intestine pathway. The medium at pH 7.4 was selected to examine the interconversion of the tested atenolol prodrug in the blood circulation system. Acid-catalyzed hydrolysis of the atenolol ProD 1 was found to be higher in 1 N HCl than at pH 2 and 5 (Figure 7). At 1 N HCl, the atenolol ProD 1 was hydrolyzed to release the parent drug in 2.53 hours. On the other hand, at pH 7.4, the prodrug was entirely stable and no release of the parent drug was observed. Since the of atenolol ProD 1 carboxylic acid is in the range of 3-4, it is expected at pH 5 that the anionic form of the prodrug will be dominant and the percentage of the free acidic form that undergoes the acid-catalyzed hydrolysis will be relatively low. At 1 N HCl and pH 2, most of the prodrug will exist as the free acid form, whereas at pH 7.4 most of the prodrug will be in the anionic form. Thus, the difference in rates at the different pH buffers.
The QM calculations at different levels demonstrated that the efficiency of Kirby’s N-alkylmaleamic acids and atenolol ProD 1-ProD 2 is largely sensitive to the pattern of substitution on the carbon-carbon double bond and nature of the alkyl group on amide nitrogen. The linear correlation found between the acid-catalyzed hydrolysis strain energy difference between the intermediate and the reactant (Es INT-GM) suggests that the reaction is governed by strain effect. Moreover, a strong correlation obtained from plotting the calculated and experimental EM values (effective molarity) reinforces the credibility of using such methods for energy and rate predictions of the kind of processes reported herein.
Comparison of the calculated DFT properties for 1–7 atenolol prodrugs ProD 1-ProD 2 with previously calculated properties for the acid-catalyzed hydrolysis of acyclovir prodrugs and cefuroxime ProD 1–ProD 4 (Figure 5) indicates that while, for systems 1–7 and atenolol prodrugs ProD 1-ProD 2, the rate-limiting step was the breakdown of the tetrahedral intermediate in the reactions of cefuroxime prodrugs ProD 1–ProD 4 and acyclovir prodrugs ProD 1–ProD 4, the rate-limiting step was the formation of the tetrahedral intermediate (Figure 8). This might be due to the nature of the amine leaving group involved in the tetrahedral intermediate dissociation step.
Using the correlation equation obtained from the plot of the calculated and experimental EM values, the values of two different atenolol prodrugs (ProD 1-ProD 2) were estimated.
Comparison between the calculated values 63.2 hours) for atenolol ProD 1 to the experimental value (3.82 hours) indicates that while the value obtained by B3LYP/6-31G (d,p) is overestimated (about 17 times larger than the experimental), the values obtained by mpwpw91/6-31+G (d,p) were much more closer, 6.3 hours. This discrepancy between the calculated and experimental values might be attributed to the following. (1) B3LYP/6-31 G (d,p) is a DFT method without dispersion corrections and (2) PCM solvation model (calculations in presence of solvent) is not capable of handling calculations in acidic aqueous solvent (medium), since the dielectric constant for pH 2 aqueous solutions is not known. In the study, calculations of the value of 78.39 (dielectric constant for pure water) were used instead (it should be noted that DFT calculations in the presence of a mixture of acid and water are not feasible).
The experimental value at pH 5 was 133 hours and at pH 7.4 no interconversion was observed. The lack of the reaction at the latter pH might be due to the fact that at this pH atenolol ProD 1 exists solely in the ionized form ( about 3-4). As mentioned before, the free acid form is a mandatory requirement for the reaction to proceed.
Future strategy to achieve more efficient atenolol prodrugs capable of increasing the liquid formulation stability, eliminating atenolol bitterness, and releasing the parent drug in a programmable manner is (a) synthesis of atenolol prodrugs having around 6 (intestine pH) such as atenolol ProD 3; (ii) in vitro kinetic studies of atenolol ProD 3 performed at pH 6.5 (intestine) and pH 7.4 (blood circulation system); and (iii) in vivo pharmacokinetic studies done in order to determine the bioavailability and the duration of action of the tested prodrug. Furthermore, based on the in vivo pharmacokinetics characteristics of atenolol ProD 3, new prodrugs may be design and synthesized.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Rafik Karaman would like to acknowledge funding by the German Research Foundation (DFG, ME 1024/8-1). Special thanks are also given to Nardene Karaman, Angi Karaman, Donia Karaman, and Rowan Karaman for the technical assistance.
- G. K. McEvoy, AHFS Drug Information 93, American Society of Health System Pharmacists, 1993.
- A. Melander, P. Stenberg, and H. Liedholm, “Food-induced reduction in bioavailability of atenolol,” European Journal of Clinical Pharmacology, vol. 16, no. 5, pp. 327–330, 1979.
- J. McAinsh, W. T. Simpson, and B. F. Holmes, “Bioavailability of atenolol formulations,” Biopharmaceutics and Drug Disposition, vol. 1, no. 6, pp. 323–332, 1980.
- J. F. Standing and C. Tuleu, “Paediatric formulations—getting to the heart of the problem,” International Journal of Pharmaceutics, vol. 300, no. 1-2, pp. 56–66, 2005.
- S. S. Garner, D. B. Wiest, and E. R. Reynolds Jr., “Stability of atenolol in an extemporaneously compounded oral liquid,” American Journal of Hospital Pharmacy, vol. 51, no. 4, pp. 508–511, 1994.
- T. Foppa, F. S. Murakami, and M. A. S. Silva, “Development, validation and stability study of pediatric atenolol syrup,” Pharmazie, vol. 62, no. 7, pp. 519–521, 2007.
- D. Patel, D. H. Doshi, and A. Desia, “Short term stability of atenolol in oral liquid formulations,” International Journal of Pharmaceutical Compounding, vol. 1, no. 6, pp. 437–439, 1997.
- K. Chan and J. Swenden, “Pilot study of the short-term physico-chemical stability of atenolol tablets stored in a multi-compartment compliance aid,” European Journal of Hospital Pharmacy, vol. 13, pp. 60–66, 2007.
- R. Karaman, K. Dajani, and H. Hallak, “Computer-assisted design for atenolol prodrugs for the use in aqueous formulations,” Journal of Molecular Modeling, vol. 18, pp. 1523–1540, 2012.
- B. Anroop, B. Ghosh, V. Parcha, and J. Khanam, “Transdermal delivery of atenolol: effect of prodrugs and iontophoresis,” Current Drug Delivery, vol. 6, no. 3, pp. 280–290, 2009.
- H. Sohi, Y. Sultana, and R. K. Khar, “Taste masking technologies in oral pharmaceuticals: recent developments and approaches,” Drug Development and Industrial Pharmacy, vol. 30, no. 5, pp. 429–448, 2004.
- W. J. Reilly, “Pharmaceutical necessities,” in Remington: The Science and Practice of Pharmacy, pp. 1018–1020, Mack Publishing Company, 2002.
- A. A. Fawzy, “Pleasant-tasting aqueous liquid composition of a bitter-tasting drug,” PCT International Applications, WO 1998005312 A1, 1998.
- W. G. Gowan, “Aliphatic esters as a solventless coating for pharmaceuticals,” Patent CA 2082137 C, 1993.
- K. Gowthamarajan, G. T. Kulkarni, and M. N. Kumar, “Pop the pills without bitterness: taste-masking technologies for bitter drugs,” Resonance, vol. 9, no. 12, pp. 25–32, 2004.
- J. A. Bakan, “Microencapsulation,” in Theory and Practice of Industrial Pharmacy, pp. 412–429, 3rd edition, 1986.
- V. S. Iyer and S. C. Srinivas, “Effervescent granular formulations of antiretroviral drugs,” Patent WO2007060682, 2007.
- L. Bush, “Bitter taste bypass need for sugar spoon,” Pharm Technol, 2004, http://www.pharmtech.com/pharmtech/data/articlestandard//pharmtech/072004/84521/article.pdf.
- “Ion exchange resin complexes: an approach to mask the taste of bitter drugs,” http://www.pharmainfo.net/reviews/ion-exchange-resin-complexes-approach-mask-taste-bitter-drugs.
- W. S. Bress, N. Kulkarni, S. Ambike, and M. P. Ramsay, “Fast dissolving orally consumable filsm containing an ion exchange resin as a taste masking agent,” Patent EP1674078, 2006.
- W. R. Mendes, Theory and Practice of Industrial Pharmacy, 3rd edition, 1976.
- A. M. J. Redondo and L. B. Abanades, “Aqueous base liquid pharmaceutical compositions in suspension form for the oral administration of ibuprofen,” Patent WO 047550, 2003.
- N. Kashid, P. Chouhan, and G. Mukherji, “Taste masked pharmaceutical composition for oral solid dosage form and process for preparing the same using magnesium aluminium silicate,” Patent WO2007108010, 2007.
- H. Hejaz, R. Karaman, and M. Khamis, “Computer-assisted design for paracetamol masking bitter taste prodrugs,” Journal of Molecular Modeling, vol. 18, no. 1, pp. 103–114, 2012.
- S. C. Shin and J. S. Choi, “Enhanced bioavailability of atenolol by transdermal administration of the ethylene-vinyl acetate matrix in rabbits,” European Journal of Pharmaceutics and Biopharmaceutics, vol. 56, no. 3, pp. 439–443, 2003.
- A. C. Montes-Gil, M. Zanfolin, C. E. Okuyama et al., “Pharmacokinetic profile of atenolol aspirinate,” Archiv der Pharmazie, vol. 340, no. 9, pp. 445–455, 2007.
- V. J. Stella, “Prodrugs as therapeutics,” Expert Opinion on Therapeutic Patents, vol. 14, no. 3, pp. 277–280, 2004.
- V. J. Stella, “Prodrugs: some thoughts and current issues,” Journal of Pharmaceutical Sciences, vol. 99, no. 12, pp. 4755–4765, 2010.
- K. M. Huttunen, H. Raunio, and J. Rautio, “Prodrugs-from serendipity to rational design,” Pharmacological Reviews, vol. 63, no. 3, pp. 750–771, 2011.
- V. J. Stella, R. T. Borchardt, M. Hageman, R. Oliyai, H. Maag, and J. Tilley, Prodrugs: challenges and Rewards Part 1 and Part 2, Springer, New York, NY, USA, 2007.
- V. J. Stella, W. N. A. Charman, and V. H. Naringrekar, “Prodrugs. Do they have advantages in clinical practice?” Drugs, vol. 29, no. 5, pp. 455–473, 1985.
- P. K. Banerjee and G. L. Amidon, “Design of prodrugs based on enzymes substrate specificity,” in Design of Prodrugs, pp. 93–133, 1985.
- C. E. Müller, “Prodrug approaches for enhancing the bioavailability of drugs with low solubility,” Chemistry and Biodiversity, vol. 6, no. 11, pp. 2071–2083, 2009.
- R. G. Parr and W. Yang, Density Functional Theory of Atoms and Molecules, Oxford University Press, Oxford, UK, 1989.
- U. Burker and N. L. Allinger, Molecular Mechanics, American Chemical Society, Washington, DC, USA, 1982.
- A. Warshel and M. Levitt, “Theoretical studies of enzymatic reactions: dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme,” Journal of Molecular Biology, vol. 103, no. 2, pp. 227–249, 1976.
- A. J. Kirby and P. W. Lancaster, “Structure and efficiency in intramolecular and enzymic catalysis. Catalysis of amide hydrolysis by the carboxy-group of substituted maleamic acids,” Journal of the Chemical Society, Perkin Transactions 2, no. 9, pp. 1206–1214, 1972.
- R. Karaman, “Analyzing the efficiency in intramolecular amide hydrolysis of Kirby's N-alkylmaleamic acids—a computational approach,” Computational and Theoretical Chemistry, vol. 974, no. 1–3, pp. 133–142, 2011.
- S. E. Barber, K. E. S. Dean, and A. J. Kirby, “A mechanism for efficient proton-transfer catalysis. Intramolecular general acid catalysis of the hydrolysis of 1-arylethyl ethers of salicylic acid,” Canadian Journal of Chemistry, vol. 77, no. 5-6, pp. 792–801, 1999.
- A. J. Kirby, M. F. Lima, D. da Silva, C. D. Roussev, and F. Nome, “Efficient intramolecular general acid catalysis of nucleophilic attack on a phosphodiester,” Journal of the American Chemical Society, vol. 128, no. 51, pp. 16944–16952, 2006.
- A. J. Kirby and N. H. Williams, “Efficient intramolecular general acid catalysis of enol ether hydrolysis. Hydrogen-bonding stabilization of the transition state for proton transfer to carbon,” Journal of the Chemical Society, Perkin Transactions 2, vol. 2, pp. 643–648, 1994.
- A. J. Kirby and N. H. Williams, “Efficient intramolecular general acid catalysis of vinyl ether hydrolysis by the neighbouring carboxylic acid group,” Journal of the Chemical Society D, no. 22, pp. 1643–1644, 1991.
- A. J. Kirby, “Enzyme Mechanisms, Models, and Mimics,” Angewandte Chemie, vol. 35, no. 7, pp. 707–724, 1996.
- T. H. Fife and T. J. Przystas, “Intramolecular general acid catalysis in the hydrolysis of acetals with aliphatic alcohol leaving groups,” Journal of the American Chemical Society, vol. 101, no. 5, pp. 1202–1210, 1979.
- A. J. Kirby, “Efficiency of proton transfer catalysis in models and enzymes,” Accounts of Chemical Research, vol. 30, no. 7, pp. 290–296, 1997.
- R. Karaman, “The efficiency of proton transfer in Kirby's enzyme model, a computational approach,” Tetrahedron Letters, vol. 51, no. 16, pp. 2130–2135, 2010.
- R. Karaman, “The effective molarity (EM) puzzle in proton transfer reactions,” Bioorganic Chemistry, vol. 37, no. 4, pp. 106–110, 2009.
- F. M. Menger and M. Ladika, “Fast hydrolysis of an aliphatic amide at neutral pH and ambient temperature. A peptidase model,” Journal of the American Chemical Society, vol. 110, no. 20, pp. 6794–6796, 1988.
- F. M. Menger, “On the source of intramolecular and enzymatic reactivity,” Accounts of Chemical Research, vol. 18, no. 5, pp. 128–134, 1985.
- F. M. Menger, J. F. Chow, H. Kaiserman, and P. C. Vasquez, “Directionality of proton transfer in solution. Three systems of known angularity,” Journal of the American Chemical Society, vol. 105, no. 15, pp. 4996–5002, 1983.
- F. M. Menger, A. L. Galloway, and D. G. Musaev, “Relationship between rate and distance,” Chemical Communications, vol. 9, no. 18, pp. 2370–2371, 2003.
- R. Karaman, “Analysis of Menger's 'spatiotemporal hypothesis',” Tetrahedron Letters, vol. 49, no. 41, pp. 5998–6002, 2008.
- S. Milstien and L. A. Cohen, “Concurrent general-acid and general-base catalysis of esterification,” Journal of the American Chemical Society, vol. 92, no. 14, pp. 4377–4382, 1970.
- S. Milstien and L. A. Cohen, “Rate acceleration by stereopopulation control: models for enzyme action,” Proceedings of the National Academy of Sciences of the United States of America, vol. 67, no. 3, pp. 1143–1147, 1970.
- S. Milstien and L. A. Cohen, “Stereopopulation control. I. Rate enhancement in the lactonizations of o-hydroxyhydrocinnamic acids,” Journal of the American Chemical Society, vol. 94, no. 26, pp. 9158–9165, 1972.
- R. Karaman, “Proximity vs. strain in intramolecular ring-closing reactions,” Molecular Physics, vol. 108, no. 13, pp. 1723–1730, 2010.
- R. F. Brown and N. M. van Gulick, “The geminal alkyl effect on the rates of ring closure of bromobutylamines,” Journal of Organic Chemistry, vol. 21, no. 9, pp. 1046–1049, 1956.
- T. C. Bruice and U. K. Pandit, “The effect of geminal substitution ring size and rotamer distribution on the intra molecular nucleophilic catalysis of the hydrolysis of monophenyl esters of dibasic acids and the solvolysis of the intermediate anhydrides,” Journal of the American Chemical Society, vol. 82, no. 22, pp. 5858–5865, 1960.
- T. C. Bruice and U. K. Pandit, “Intramolecular models depicting the kinetic importance of “Fit” in enzymatic catalysis,” Proceedings of the National Academy of Sciences of the United States of America, vol. 46, pp. 402–404, 1960.
- C. Galli and L. Mandolini, “The role of ring strain on the ease of ring closure of bifunctional chain molecules,” European Journal of Organic Chemistry, no. 18, pp. 3117–3125, 2000.
- R. Karaman, “Prodrugs of aza nucleosides based on proton transfer reaction,” Journal of Computer-Aided Molecular Design, vol. 24, no. 12, pp. 961–970, 2010.
- R. Karaman, W. Amly, S. A. Bufo, and L. Scrano, “Computationally designed prodrugs of statins based on Kirby's enzyme model,” Journal of Molecular Modeling, vol. 19, no. 9, pp. 3969–3982, 2013.
- R. Karaman, D. Karaman, and I. Zeiadeh, “Computationally-designed phenylephrine prodrugs—a model for enhancing bioavailability,” Molecular Physics, vol. 111, no. 21, pp. 3249–3264, 2013.
- R. Karaman, H. Ghareeb, K. K. Dajani et al., “Design, synthesis and in-vitro kinetic study of tranexamic acid prodrugs for the treatment of bleeding conditions,” Journal of Computer-Aided Molecular Design, vol. 27, no. 7, pp. 615–635, 2013.
- R. Karaman, K. K. Dajani, A. Qtait, and M. Khamis, “Prodrugs of acyclovir—a computational approach,” Chemical Biology and Drug Design, vol. 79, no. 5, pp. 819–834, 2012.
- R. Karaman, “Prodrugs for masking bitter taste of antibacterial drugs—a computational approach,” Journal of Molecular Modeling, vol. 19, no. 6, pp. 2399–2412, 2013.
- R. Karaman, “Computational-aided design for dopamine prodrugs based on novel chemical approach,” Chemical Biology and Drug Design, vol. 78, no. 5, pp. 853–863, 2011.
- Y. Zhao and D. G. Truhlar, “Density functionals with broad applicability in chemistry,” Accounts of Chemical Research, vol. 41, no. 2, pp. 157–167, 2008.
- Y. Zhao and D. G. Truhlar, “Exploring the limit of accuracy of the global hybrid meta density functional for main-group thermochemistry, kinetics, and noncovalent interactions,” Journal of Chemical Theory and Computation, vol. 4, no. 11, pp. 1849–1868, 2008.
- J. Zheng, Y. Zhao, and D. G. Truhlar, “The DBH24/08 database and its use to assess electronic structure model chemistries for chemical reaction barrier heights,” Journal of Chemical Theory and Computation, vol. 5, no. 4, pp. 808–821, 2009.
- M. J. Frisch, G. W. Trucks, H. B. Schlegel, et al., Gaussian 09, Revision A.1, Gaussian Inc., Wallingford, Conn, USA, 2009.
- C. J. Casewit, K. S. Colwell, and A. K. Rappé, “Application of a universal force field to main group compounds,” Journal of the American Chemical Society, vol. 114, no. 25, pp. 10046–10053, 1992.
- J. N. Murrell and K. J. Laidler, “Symmetries of activated complexes,” Transactions of the Faraday Society, vol. 64, pp. 371–377, 1968.
- K. Muller, “Reaction paths on multidimensional energy hypersurfaces,” Angewandte Chemie, vol. 19, no. 1, pp. 1–13, 1980.
- E. Cancès, B. Mennucci, and J. Tomasi, “A new integral equation formalism for the polarizable continuum model: theoretical background and applications to Isotropic and anisotropic dielectrics,” Journal of Chemical Physics, vol. 107, no. 8, pp. 3032–3041, 1997.
- B. Mennucci and J. Tomasi, “Continuum solvation models: a new approach to the problem of solute's charge distribution and cavity boundaries,” Journal of Chemical Physics, vol. 106, no. 12, pp. 5151–5158, 1997.
- B. Mennucci, E. Cancès, and J. Tomasi, “Evaluation of solvent effects in isotropic and anisotropic dielectrics and in ionic solutions with a unified integral equation method: theoretical bases, computational implementation, and numerical applications,” Journal of Physical Chemistry B, vol. 101, no. 49, pp. 10506–10517, 1997.
- J. Tomasi, B. Mennucci, and E. Cancès, “The IEF version of the PCM solvation method: an overview of a new method addressed to study molecular solutes at the QM ab initio level,” Journal of Molecular Structure, vol. 464, no. 1–3, pp. 211–226, 1999.
- J. McAinsh, W. T. Simpson, and B. F. Holmes, “Bioavailability of atenolol formulations,” Biopharmaceutics and Drug Disposition, vol. 1, no. 6, pp. 323–332, 1980.
- H. Vergin and V. Nitsche, “Oral bioavailability of atenolol,” Journal of International Medical Research, vol. 17, no. 5, pp. 417–425, 1989.
- R. Hipple and M. Nahahta, “Atenolol oral suspension,” in Pediatric Drug Formulations, 4th edition, 2000.
- T. Katagi, “AM1 study of acid-catalyzed hydrolysis of maleamic (4-amino-4-oxo-2-butenoic) acids,” Journal of Computational Chemistry, vol. 11, no. 9, pp. 1094–1100, 1990.
- R. Kluger and J. Chin, “Carboxylic acid participation in amide hydrolysis. Evidence that separation of a nonbonded complex can be rate determining,” Journal of the American Chemical Society, vol. 104, no. 10, pp. 2891–2897, 1982.
- R. Karaman, “Effects of substitution on the effective molarity (EM) for five membered ring-closure reactions—a computational approach,” Journal of Molecular Structure, vol. 939, no. 1–3, pp. 69–74, 2010.
- R. Karaman, “The effective molarity (EM) puzzle in intramolecular ring-closing reactions,” Journal of Molecular Structure, vol. 940, no. 1–3, pp. 70–75, 2010.
- R. Karaman, “The effective molarity (EM)-A computational approach,” Bioorganic Chemistry, vol. 38, no. 4, pp. 165–172, 2010.
- A. J. Kirby, “Effective molarities for intramolecular reactions,” Journal of Physical Organic Chemistry, vol. 18, pp. 101–278, 2005.