Comparison of Experimental and Theoretical Structures of a Transition State Analogue Used for the Induction of Anti-Cocaine Catalytic Antibodies

Comparison of the crystal structure of a transition state analogue that was used to raise catalytic antibodies for the benzoyl ester hydrolysis of cocaine with structures calculated by ab initio, semiempirical, and solvation semiempirical methods reveals that modeling of solvation is crucial for replicating the crystal structure geometry. Both SM3 and SM2 calculations, starting from the crystal structure TSA I, converged on structures similar to the crystal structure. The 3-21G(*)/HF, 6-31G*/HF, PM3, and AM1 calculations converged on structures similar to each other, but these gas-phase structures were significantly extended relative to the condensed phase structures. Two transition states for the hydrolysis of the benzoyl ester of cocaine were located with the SM3 method. The gas phase calculations failed to locate reasonable transition state structures for this reaction. These results imply that accurate modeling of the potential energy surfaces for the hydrolysis of cocaine requires solvation methods.

Quantum Mechnical Investion of Cyclic 3',5'- Adenosine Monophosphate, the Second Hormonal Messenger

Full geometry optimizations using the PM3, AM1, 3-21G∗/HF and 6-31G∗/HF levels of theory were conducted on the syn and anti conformations of cyclic3′,5′-adenosine monophosphate (cAMP). Comparison of the anti crystal structures with the semiempirical and ab initio results revealed that the ab initio results agree well with the experimental results. The results of semiempirical calculations are in qualitative agreement with experimental and ab initio values, with the exception of the glycosyl torsion angle for the anti conformer. Sugar puckering, which is not handled properly by semiempirical methods for unconstrained sugars, nucleosides, nucleotides and nucleotide base pairs, is modeled reasonably well by the semiempirical methods for cAMP. This improvement results from the constraints introduced by the cyclization of AMP to form the phosphodiester.

Hydrogen bonding of nucleotide base pairs: Application of the PM3 method

The ability of the pm3 semiempirical quantum mechanical method to reproduce hydrogen bonding in nucleotide base pairs was assessed. Results of pm3 calculations on the nucleotides 2′-deoxyadenosine 5′-monophosphate (pdA), 2′-deoxyguanosine 5′-monophosphate (pdG), 2′-deoxycytidine 5′-monophosphate (pdC), and 2′-deoxythymidine 5′-monophosphate (pdT) and the base pairs pdA–pdT, pdG–pdC, and pdG(syn)–pdC are presented and discussed. The pm3 method is the first of the parameterized nddo quantum mechanical models with any ability to reproduce hydrogen bonding between nucleotide base pairs. Intermolecular hydrogen bond lengths between nucleotides displaying Watson–Crick base pairing are 0.1–0.2 Å less than experimental results. Nucleotide bond distances, bond angles, and torsion angles about the glycosyl bond (χ), the C4′C5′ bond (γ), and the C5′O5′ bond (β) agree with experimental results. There are many possible conformations of nucleotides. pm3 calculations reveal that many of the most stable conformations are stabilized by intramolecular CHO hydrogen bonds. These interactions disrupt the usual sugar puckering. The stacking interactions of a dT–pdA duplex are examined at different levels of gradient optimization. The intramolecular hydrogen bonds found in the nucleotide base pairs disappear in the duplex, as a result of the additional constraints on the phosphate group when part of a DNA backbone. Sugar puckering is reproduced by the pm3 method for the four bases in the dT–pdA duplex. pm3 underestimates the attractive stacking interactions of base pairs in a B-DNA helical conformation. The performance of the pm3 method implemented in SPARTAN is contrasted with that implemented in MOPAC. At present, accurate ab initio calculations are too timeconsuming to be of practical use, and molecular mechanics methods cannot be used to determine quantum mechanical properties such as reaction-path calculations, transition-state structures, and activation energies. The pm3 method should be used with extreme caution for examination of small DNA systems. Future parameterizations of semiempirical methods should incorporate base stacking interactions into the parameterization data set to enhance the ability of these methods.

A Computationally Efficient Procedure for Modeling the First Step in the Alkaline Hydrolysis of Esters

A computationally efficient procedure for modeling the alkaline hydrolysis of esters is proposed based on calculations performed on methyl acetate and methyl benzoate systems. Extensive geometry and energy comparisons were performed on the simple ester methyl acetate. The effectiveness of performing high level single point ab initio energy calculations on the geometries obtained from semiempirical and ab initio methods was determined. The AM1 and PM3 semiempirical methods are evaluated for their ability to model the transition states and intermediates for ester hydrolysis. The Cramer/Truhlar SM3 solvation method was used to determine activation energies. The most computationally efficient way to model the transition states of large esters is to use the PM3 method. The PM3 transition structure can then be used as a template for the design of haptens capable of inducing catalytic antibodies.

The reaction of 8-amino-p-menthene derivatives with electrophiles

Attempts to ring-close the nitrogen atom of 8-amino-p-menth-1-ene and of N-substituted 8-amino-p-menth-1-enes onto the C1 - C2 double-bond carbons has led to a range of bicyclo[2.2.2] and bicyclo[3.2.1] products, together with the novel bicyclo[4.3.1]-1,3-oxazepine 9.

Relative basicities of free base porphyrins; understanding the role of macrocyclic distortion

Porphyrins have been the center of numerous investigations in different areas of chemistry, geochemistry, and the life sciences. In nature the conformation of the porphyrin macrocycle varies, depending on the function of its apoenzyme. It is believed that the conformation of the porphyrin ring is necessary for the enzyme to achieve its function and modify its reactivity. It is important to understand how the conformation of the porphyrin ring will influence its properties. ^ In synthetic porphyrins particular conformations and ring deformations can be achieved by peripheral substitution, metallation, core substitution, and core protonation among other alterations of the macrocycle. The macrocyclic distortions will affect the ring current, the ability of pyrroles to intramolecularly hydrogen bond and the relative basicity of each of the porphyrins. To understand these effects different theoretical models are used. The ground state structure of each of 19 free base porphyrins is determined using molecular mechanics (MM+) and semiempirical methods (PM3). The energetics of deformation of the macrocyclic core is calculated by carrying out single point energy calculations for the conformation achieved by each synthetic compound. Enthalpies of solution and enthalpies of protonation of 10 porphyrins with varying degrees of macrocyclic deformation and varying electron withdrawing groups in the periphery are determined using solution calorimetry. Using Hess's Law, the relative basicity of each of the different free base porphyrins is calculated. NMR results are described, including the determination of free energies of activation of ring tautomerization and hydrogen bonding for several compounds. It was found that in the absence of electronic effects, the greater macrocyclic deformation, the greater the basicity of the porphyrins. This basicity is attenuated by the presence of electron withdrawing groups and ability to of the macrocycle to intramolecularly hydrogen bond. ^

Relative basicities of free base porphyrins : understanding the role of macrocyclic distortion

Porphyrins have been the center of numerous investigations in different areas of chemistry, geochemistry, and the life sciences. In nature the conformation of the porphyrin macrocycle varies, depending on the function of its apoenzyme. It is believed that the conformation of the porphyrin ring is necessary for the enzyme to achieve its function and modify its reactivity. It is important to understand how the conformation of the porphyrin ring will influence its properties. In synthetic porphyrins particular conformations and ring deformations can be achieved by peripheral substitution, metallation, core substitution, and core protonation among other alterations of the macrocycle. The macrocyclic distortions will affect the ring current, the ability of pyrroles to intramolecularly hydrogen bond and the relative basicity of each of the porphyrins. To understand these effects different theoretical models are used. The ground state structure of each of 19 free base porphyrins is determined using molecular mechanics (MM+) and semiempirical methods (PM3). The energetics of deformation of the macrocyclic core is calculated by carrying out single point energy calculations for the conformation achieved by each synthetic compound. Enthalpies of solution and enthalpies of protonation of 10 porphyrins with varying degrees of macrocyclic deformation and varying electron withdrawing groups in the periphery are determined using solution calorimetry. Using Hess's Law, the relative basicity of each of the different free base porphyrins is calculated. NMR results are described, including the determination of free energies of activation of ring tautomerization and hydrogen bonding for several compounds. It was found that in the absence of electronic effects, the greater macrocyclic deformation, the greater the basicity of the porphyrins. This basicity is attenuated by the presence of electron withdrawing groups and ability to of the macrocycle to intramolecularly hydrogen bond.

Hydrolysis of genotoxic methyl-substituted oxiranes : Experimental kinetic and semiempirical studies

The kinetics of acid-catalyzed hydrolysis of seven methylated aliphatic epoxides - R1R2C(O)CR3R4 (A: R1=R2=R3=R4=H; B: R1=R2=R3=H, R4=Me; C: R1=R2=H, R3=R4=Me; D: R1=R3=H, R2=R4=Me(trans); E: R1=R3=H, R2=R4=Me(cis); F: R1=R3=R4=Me, R2=H; G: R1=R2=R3=R4=Me) - has been studied at 36 ± 1.5°C. Compounds with two methyl groups at the same carbon atom of the oxirane ring exhibit highest rate constants (k(eff) in reciprocal molar concentration per second: 11.0 ± 1.3 for C, 10.7 ± 2.1 for F, and 8.7 ± 0.7 for G as opposed to 0.124 ± 0.003 for B, 0.305 ± 0.003 for D, and 0.635 ± 0.036 for E). Ethylene oxide (A) displays the lowest rate of hydrolysis (0.027 M-1 s-1). The results are consistent with literature data available for compounds A, B, and C. To model the reactivities we have employed quantum chemical calculations (MNDO, AM1, PM3, and MINDO/3) of the main reaction species. There is a correlation of the logarithm k(eff) with the total energy of epoxide ring opening. The best correlation coefficients (r) were obtained using the AM1 and MNDO methods (0.966 and 0.957, respectively). However, unlike MNDO, AM1 predicts approximately zero energy barriers for the oxirane ring opening of compounds B, C, E and G, which is not consistent with published kinetic data. Thus, the MNDO method provides a preferential means of modeling the acidic hydrolysis of the series of methylated oxiranes. The general ranking of mutagenicity in vitro, A > B > C, is in line with the concept that this sequence also gradually leaves the expoxide reactivity optimal for genotoxicity toward reactivities leading to higher biological detoxifications.

A comparison of the molecular orbital (ab initio, CNDO and EHT) methods applied to the conformational study of thiocarbonohydrazide in its basic and diprotonated forms

Conformational preferences of thiocarbonohydrazide (H2NNHCSNHNH2) in its basic and N,N′-diprotonated forms are examined by calculating the barrier to internal rotation around the C---N bonds, using the theoretical LCAO—MO (ab initio and semiempirical CNDO and EHT) methods. The calculated and experimental results are compared with each other and also with values for N,N′-dimethylthiourea which is isoelectronic with thiocarbonohydrazide. The suitability of these methods for studying rotational isomerism seems suspect when lone pair interactions are present.