We present a molecular simulation protocol to compute free energies of binding, which combines a QM/MM correction term with rigorous classical free energy techniques, thereby accounting for electronic polarization effects. Relative free energies of binding are first computed using classical force fields, Monte Carlo sampling, and replica exchange thermodynamic integration. Snapshots of the configurations at the end points of the perturbation are then subjected to DFT-QM/MM single-point calculations using the B3LYP functional and a range of basis sets. The resulting quantum mechanical energies are then processed using the Zwanzig equation to give free energies incorporating electronic polarization. Our approach is conceptually simple and does not require tightly coupled QM and MM software. The method has been validated by calculating the relative free energies of hydration of methane and water and the relative free energy of binding of two inhibitors of cyclooxygenase-2. Closed thermodynamic cycles are obtained across different pathways, demonstrating the correctness of the technique, although significantly more sampling is required for the protein-ligand system. Our method offers a simple and effective way to incorporate quantum mechanical effects into computed free energies of binding.

We present a computational model study designed to simulate the results of time-resolved fluorescence spectra of tryptophan in proteins. In such measurements, the occurrence of more than one fluorescence lifetime is generally attributed to the existence of several tryptophan rotamers and/or structural conformations of the protein structure. The protein system we chose for this initial study is the tetracycline repressor (TetR), an interesting model system for the investigation of the mechanisms of transcriptional regulation. Fluorescence resonance energy transfer (FRET) from tryptophan to tetracycline is frequently observed in complexes of the TetR with the antibiotic tetracycline. We use a combined classical/quantum mechanical approach to model the structure and the spectroscopic properties of the TetR-tetracycline complex. A classical molecular dynamics simulation provides input geometries for semiempirical quantum mechanical/molecular mechanical (QM/MM) single-point configuration interaction (CI) calculations, which are used to calculate tryptophan vertical absorption and fluorescence energies and intensities as well as relative FRET rate constants. These rate constants together with the Einstein coefficients for spontaneous emission and an assumed rate for nonradiative deactivation allow us to simulate fluorescence decay curves with and without FRET and for the entire ensemble as well as for individual rotamers. Our results indicate that the classical "rotamer model", used to explain the multiexponential fluorescence-decay curves of time-resolved tryptophan emission spectra, can be extended to systems with FRET acceptors present in the protein matrix but that the interpretation of the fitted lifetimes is different to that usually used.

Density functional theory (DFT) has been used to investigate the conformations and tautomeric forms of neutral tetracycline in aqueous solution. The results suggest that the extended conformation is 3−3.5 kcal mol-1 more stable than the twisted one for equivalent tautomers and that as many as six different tautomeric forms lie within 10 kcal mol-1 of the most stable one. The energetic preference for the extended conformation is a solvent effect. Calculated infrared, NMR, and UV/vis spectra are used to suggest ways of differentiating between tautomers and conformations experimentally. It is shown that the 13C chemical shifts of C4a, C5, and C6 can be used to distinguish between the twisted and extended conformations.

No abstract

N-Alkylaminoferrocene (NAAF)-based prodrugs are activated in the presence of elevated amounts of reactive oxygen species (ROS), which corresponds to cancer specific conditions, with formation of NAAF and p-quinone methide. Both products act synergistically by increasing oxidative stress in cancer cells that causes their death. Though it has already been demonstrated that the best prodrugs of this type retain their antitumor activity in vivo, the effects were found to be substantially weaker than those observed in cell cultures. Moreover, the mechanistic studies of these compounds in vivo are missing. For clarification of these important questions, labeling of the prodrugs with radioactive moieties would be necessary. In this paper, we first observed that the representative NAAF-based prodrugs are hydrolyzed in dilute aqueous solutions to the corresponding arylboronic acids. We confirmed that these products are responsible for ROS amplification and anticancer properties of the parent prodrugs. Next, we developed the efficient synthetic protocol for radiolabeling the hydrolyzed NAAF-based prodrugs by [ 18 F]fluoroglucosylation under the conditions of the copper(I)-catalyzed azide−alkyne 1,3-dipolar cycloaddition and used this protocol to prepare one representative hydrolyzed NAAF-based prodrug radiolabeled with 18 F. Finally, we studied the stability of the 18 F-labeled compound in human serum in vitro and in rat blood in vivo and obtained preliminary data on its biodistribution in vivo in mice carrying pancreatic (AR42J) and prostate (PC3) tumors by applying PET imaging studies. The compound described in this paper will help to understand in vivo effects (e.g., pharmacokinetics, accumulation in organs, the nature of side effects) of these prodrugs that will strongly contribute to their advancement to clinical trials.