Progress toward chemical accuracy in the computer simulation of condensed phase reactions.

We describe a procedure for the generation of chemically accurate computer-simulation models to study chemical reactions in the condensed phase. The process involves (i) the use of a coupled semiempirical quantum and classical molecular mechanics method to represent solutes and solvent, respectively; (ii) the optimization of semiempirical quantum mechanics (QM) parameters to produce a computationally efficient and chemically accurate QM model; (iii) the calibration of a quantum/classical microsolvation model using ab initio quantum theory; and (iv) the use of statistical mechanical principles and methods to simulate, on massively parallel computers, the thermodynamic properties of chemical reactions in aqueous solution. The utility of this process is demonstrated by the calculation of the enthalpy of reaction in vacuum and free energy change in aqueous solution for a proton transfer involving methanol, methoxide, imidazole, and imidazolium, which are functional groups involved with proton transfers in many biochemical systems. An optimized semiempirical QM model is produced, which results in the calculation of heats of formation of the above chemical species to within 1.0 kcal/mol (1 kcal = 4.18 kJ) of experimental values. The use of the calibrated QM and microsolvation QM/MM (molecular mechanics) models for the simulation of a proton transfer in aqueous solution gives a calculated free energy that is within 1.0 kcal/mol (12.2 calculated vs. 12.8 experimental) of a value estimated from experimental pKa values of the reacting species.

Fluorescence and photobleaching dynamics of single light-harvesting complexes

Single light-harvesting complexes LH-2 from Rhodopseudomonas acidophila were immobilized on various charged surfaces under physiological conditions. Polarized light experiments showed that the complexes were situated on the surface as nearly upright cylinders. Their fluorescence lifetimes and photobleaching properties were obtained by using a confocal fluorescence microscope with picosecond time resolution. Initially all molecules fluoresced with a lifetime of 1 ± 0.2 ns, similar to the bulk value. The photobleaching of one bacteriochlorophyll molecule from the 18-member assembly caused the fluorescence to switch off completely, because of trapping of the mobile excitations by energy transfer. This process was linear in light intensity. On continued irradiation the fluorescence often reappeared, but all molecules did not show the same behavior. Some LH-2 complexes displayed a variation of their quantum yields that was attributed to photoinduced confinement of the excited states and thereby a diminution of the superradiance. Others showed much shorter lifetimes caused by excitation energy traps that are only ≈3% efficient. On repeated excitation some molecules entered a noisy state where the fluorescence switched on and off with a correlation time of ≈0.1 s. About 490 molecules were examined.

Stability of nonrelativistic quantum mechanical matter coupled to the (ultraviolet cutoff) radiation field

We announce a proof of H-stability for the quantized radiation field, with ultraviolet cutoff, coupled to arbitrarily many non-relativistic quantized electrons and static nuclei. Our result holds for arbitrary atomic numbers and fine structure constant. We also announce bounds for the energy of many electrons and nuclei in a classical vector potential and for the eigenvalue sum of a one-electron Pauli Hamiltonian with magnetic field.

Fluorescence correlation spectroscopy reveals fast optical excitation-driven intramolecular dynamics of yellow fluorescent proteins

Fast excitation-driven fluctuations in the fluorescence emission of yellow-shifted green fluorescent protein mutants T203Y and T203F, with S65G/S72A, are discovered in the 10−6–10−3-s time range, by using fluorescence correlation spectroscopy at 10−8 M. This intensity-dependent flickering is conspicuous at high pH, with rate constants independent of pH and viscosity with a minor temperature effect. The mean flicker rate increases linearly with excitation intensity for at least three decades, but the mean dark fraction of the molecules undergoing these dynamics is independent of illumination intensity over ≈6 × 102 to 5 × 106 W/cm2. These results suggest that optical excitation establishes an equilibration between two molecular states of different spectroscopic properties that are coupled only via the excited state as a gateway. This reversible excitation-driven transition has a quantum efficiency of ≈10−3. Dynamics of external protonation, reversibly quenching the fluorescence, are also observed at low pH in the 10- to 100-μs time range. The independence of these two bright–dark flicker processes implies the existence of at least two separate dark states of these green fluorescent protein mutants. Time-resolved fluorescence measurements reveal a single exponential decay of the excited state population with 3.8-ns lifetime, after 500-nm excitation, that is pH independent. Our fluorescence correlation spectroscopy results are discussed in terms of recent theoretical studies that invoke isomerization of the chromophore as a nonradiative channel of the excited state relaxation.

The importance of reactant positioning in enzyme catalysis: A hybrid quantum mechanics/molecular mechanics study of a haloalkane dehalogenase

Hybrid quantum mechanics/molecular mechanics calculations using Austin Model 1 system-specific parameters were performed to study the SN2 displacement reaction of chloride from 1,2-dichloroethane (DCE) by nucleophilic attack of the carboxylate of acetate in the gas phase and by Asp-124 in the active site of haloalkane dehalogenase from Xanthobacter autotrophicus GJ10. The activation barrier for nucleophilic attack of acetate on DCE depends greatly on the reactants having a geometry resembling that in the enzyme or an optimized gas-phase structure. It was found in the gas-phase calculations that the activation barrier is 9 kcal/mol lower when dihedral constraints are used to restrict the carboxylate nucleophile geometry to that in the enzyme relative to the geometries for the reactants without dihedral constraints. The calculated quantum mechanics/molecular mechanics activation barriers for the enzymatic reaction are 16.2 and 19.4 kcal/mol when the geometry of the reactants is in a near attack conformer from molecular dynamics and in a conformer similar to the crystal structure (DCE is gauche), respectively. This haloalkane dehalogenase lowers the activation barrier for dehalogenation of DCE by 2–4 kcal/mol relative to the single point energies of the enzyme's quantum mechanics atoms in the gas phase. SN2 displacements of this sort in water are infinitely slower than in the gas phase. The modest lowering of the activation barrier by the enzyme relative to the reaction in the gas phase is consistent with mutation experiments.

Molecular dynamics of a grid-mounted molecular dipolar rotor in a rotating electric field

Classical molecular dynamics is applied to the rotation of a dipolar molecular rotor mounted on a square grid and driven by rotating electric field E(ν) at T ≃ 150 K. The rotor is a complex of Re with two substituted o-phenanthrolines, one positively and one negatively charged, attached to an axial position of Rh\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{2}^{4+}}}\end{equation*}\end{document} in a [2]staffanedicarboxylate grid through 2-(3-cyanobicyclo[1.1.1]pent-1-yl)malonic dialdehyde. Four regimes are characterized by a, the average lag per turn: (i) synchronous (a < 1/e) at E(ν) = |E(ν)| > Ec(ν) [Ec(ν) is the critical field strength], (ii) asynchronous (1/e < a < 1) at Ec(ν) > E(ν) > Ebo(ν) > kT/μ, [Ebo(ν) is the break-off field strength], (iii) random driven (a ≃ 1) at Ebo(ν) > E(ν) > kT/μ, and (iv) random thermal (a ≃ 1) at kT/μ > E(ν). A fifth regime, (v) strongly hindered, W > kT, Eμ, (W is the rotational barrier), has not been examined. We find Ebo(ν)/kVcm−1 ≃ (kT/μ)/kVcm−1 + 0.13(ν/GHz)1.9 and Ec(ν)/kVcm−1 ≃ (2.3kT/μ)/kVcm−1 + 0.87(ν/GHz)1.6. For ν > 40 GHz, the rotor behaves as a macroscopic body with a friction constant proportional to frequency, η/eVps ≃ 1.14 ν/THz, and for ν < 20 GHz, it exhibits a uniquely molecular behavior.