Semiempirical Study of the Bergman Reaction: Towards a Computationally Efficient and Accurate Method for Modeling Enediyne Anticancer Antibiotics

Over the past 7 years, the enediyne anticancer antibiotics have been widely studied due to their DNA cleaving ability. The focus of these antibiotics, represented by kedarcidin chromophore, neocarzinostatin chromophore, calicheamicin, esperamicin A, and dynemicin A, is on the enediyne moiety contained within each of these antibiotics. In its inactive form, the moiety is benign to its environment. Upon suitable activation, the system undergoes a Bergman cycloaromatization proceeding through a 1,4-dehydrobenzene diradical intermediate. It is this diradical intermediate that is thought to cleave double-stranded dna through hydrogen atom abstraction. Semiempirical, semiempiricalci, Hartree–Fock ab initio, and mp2 electron correlation methods have been used to investigate the inactive hex-3-ene-1,5-diyne reactant, the 1,4-dehydrobenzene diradical, and a transition state structure of the Bergman reaction. Geometries calculated with different basis sets and by semiempirical methods have been used for single-point calculations using electron correlation methods. These results are compared with the best experimental and theoretical results reported in the literature. Implications of these results for computational studies of the enediyne anticancer antibiotics are discussed.

Ability of the PM3 Quantum Mechanical Method to Model Intermolecular Hydrogen Bonding between Neutral Molecules

The PM3 semiempirical quantum-mechanical method was found to systematically describe intermolecular hydrogen bonding in small polar molecules. PM3 shows charge transfer from the donor to acceptor molecules on the order of 0.02-0.06 units of charge when strong hydrogen bonds are formed. The PM3 method is predictive; calculated hydrogen bond energies with an absolute magnitude greater than 2 kcal mol-' suggest that the global minimum is a hydrogen bonded complex; absolute energies less than 2 kcal mol-' imply that other van der Waals complexes are more stable. The geometries of the PM3 hydrogen bonded complexes agree with high-resolution spectroscopic observations, gas electron diffraction data, and high-level ab initio calculations. The main limitations in the PM3 method are the underestimation of hydrogen bond lengths by 0.1-0.2 for some systems and the underestimation of reliable experimental hydrogen bond energies by approximately 1-2 kcal mol-l. The PM3 method predicts that ammonia is a good hydrogen bond acceptor and a poor hydrogen donor when interacting with neutral molecules. Electronegativity differences between F, N, and 0 predict that donor strength follows the order F > 0 > N and acceptor strength follows the order N > 0 > F. In the calculations presented in this article, the PM3 method mirrors these electronegativity differences, predicting the F-H- - -N bond to be the strongest and the N-H- - -F bond the weakest. It appears that the PM3 Hamiltonian is able to model hydrogen bonding because of the reduction of two-center repulsive forces brought about by the parameterization of the Gaussian core-core interactions. The ability of the PM3 method to model intermolecular hydrogen bonding means reasonably accurate quantum-mechanical calculations can be applied to small biologic systems.

A Microfluidic Method to Measure Small Molecule Diffusion in Hydrogels

Drug release from a fluid-contacting biomaterial is simulated using a microfluidic device with a channel defined by solute-loaded hydrogel; as water is pumped through the channel, solute transfers from the hydrogel into the water. Optical analysis of in-situ hydrogels, characterization of the microfluidic device effluent, and NMR methods were used to find diffusion coefficients of several dyes (model drugs) in poly( ethylene glycol) diacrylate (PEG-DA) hydrogels. Diffusion coefficients for methylene blue and sulforhodamine 101 in PEG-DA calculated using the three methods are in good agreement; both dyes are mobile in the hydrogel and elute from the hydrogel at the aqueous channel interface. However, the dye acid blue 22 deviates from typical diffusion behavior and does not release as expected from the hydrogel. Importantly, only the microfluidic method is capable of detecting this behavior. Characterizing solute diffusion with a combination of NMR, optical and effluent methods offer greater insight into molecular diffusion in hydrogels than employing each technique individually. The NMR method made precise measurements for solute diffusion in all cases. The microfluidic optical method was effective for visualizing diffusion of the optically active solutes. The optical and effluent methods show potential to be used to screen solutes to determine if they elute from a hydrogel in contact with flowing fluid. Our data suggest that when designing a drug delivery device, analyzing the diffusion from the molecular level to the device level is important to establish a complete picture of drug elution, and microfluidic methods to study such diffusion can play a key role. (C) 2013 Elsevier B.V. All rights reserved.