Coupled molecular dynamics/semiempirical simulation of organic solutes in polar liquids. I. Naphthalene in acetonitrile

Spectral shift mechanisms of chlorophylls in liquids and proteins

Origins of non-excitonic spectral shifts of chlorophylls that can reach -1,000 cm(-1) in pigment-protein complexes are actively debated in literature. We investigate possible shift mechanisms, basing on absorption and fluorescence measurements in large number of liquids. Transition wavelength in solvent-free state was estimated (±2 nm) for chlorophyll a (Chl a, 647 nm), Chl b (624 nm), bacteriochlorophyll a (BChl a, 752 nm), and pheophytines. The dispersive-repulsive shift is a predominating mechanism. It depends on polarizability difference between the ground and the excited state Δα and the Lorenz-Lorentz function of refractive index of solvent (n). The approximate (± 2Å(3)) increase of polarizability Δα is close to 15Å(3) for S(1) bands of Chl a, BChl a, and BPheo a, slightly larger for Chl b (18Å(3)), and less for Pheo a (11Å(3)). The effect of solvent polarity, expressed in terms of static dielectric permittivity (ε) is relatively minor, but characteristic for different pigments and transitions. Remarkably, maximum influence of ε on S(1) band of BChl a is less (-20 ± 10 cm(-1)) than that for Chl a (-50 ± 10 cm(-1)), and not correlated with dipole moment changes on excitation Δμ (∼2D and 0.1 ± 0.1D, respectively). Hydrogen bonding in protic solvents produces red shifts in Chl a (-60 cm(-1)) and BChl a (-100 cm(-1)), but not in Chl b. Second axial ligand of BChl a has no influence on the S(1) band, whereas the S(2) transition suffers a -400 to -600 cm(-1) down shift. Aromatic character of solvent is responsible for a ∼-100 cm(-1) red shift of both Q transitions in BChl a. The S(1) bands in chlorophylls are relatively insensitive with respect to dielectric properties and specific solvation. Therefore, nontrivial mechanisms, yielding large site-energy shifts are expected in photosynthetic chlorophyll-proteins.

CAN SEMIEMPIRICAL QUANTUM MODELS CALCULATE THE BINDING ENERGY OF HYDROGEN BONDING FOR BIOLOGICAL SYSTEMS?

A modified semiempirical model named RM1BH, which is based on RM1 parameterizations, is proposed to simulate varied biological hydrogen-bonded systems. The RM1BH is formulated by adding Gaussian functions to the core–core repulsion items in original RM1 formula to reproduce the binding energies of hydrogen bonding of experimental and high-level computational results. In the parameterizations of our new model, 35 base-pair dimers, 18 amino acid residue dimers, 14 dimers between a base and an amino acid residue, and 20 other multimers were included. The results performed with RM1BH were compared with experimental values and the benchmark density-functional (B3LYP/6-31G**/BSSE) and Möller–Plesset perturbation (MP2/6-31G**/BSSE) calculations on various biological hydrogen-bonded systems. It was demonstrated that RM1BH model outperforms the PM3 and RM1 models in the calculations of the binding energies of biological hydrogen-bonded systems by very close agreement with the values of both high-level calculations and experiments. These results provide insight into the ideas, methods, and views of semiempirical modifications to investigate the weak interactions of biological systems.

RM1: A reparameterization of AM1 for H, C, N, O, P, S, F, Cl, Br, and I

Twenty years ago, the landmark AM1 was introduced, and has since had an increasingly wide following among chemists due to its consistently good results and time-tested reliability--being presently available in countless computational quantum chemistry programs. However, semiempirical molecular orbital models still are of limited accuracy and need to be improved if the full potential of new linear scaling techniques, such as MOZYME and LocalSCF, is to be realized. Accordingly, in this article we present RM1 (Recife Model 1): a reparameterization of AM1. As before, the properties used in the parameterization procedure were: heats of formation, dipole moments, ionization potentials and geometric variables (bond lengths and angles). Considering that the vast majority of molecules of importance to life can be assembled by using only six elements: C, H, N, O, P, and S, and that by adding the halogens we can now build most molecules of importance to pharmaceutical research, our training set consisted of 1736 molecules, representative of organic and biochemistry, containing C, H, N, O, P, S, F, Cl, Br, and I atoms. Unlike AM1, and similar to PM3, all RM1 parameters have been optimized. For enthalpies of formation, dipole moments, ionization potentials, and interatomic distances, the average errors in RM1, for the 1736 molecules, are less than those for AM1, PM3, and PM5. Indeed, the average errors in kcal x mol(-1) of the enthalpies of formation for AM1, PM3, and PM5 are 11.15, 7.98, and 6.03, whereas for RM1 this value is 5.77. The errors, in Debye, of the dipole moments for AM1, PM3, PM5, and RM1 are, respectively, 0.37, 0.38, 0.50, and 0.34. Likewise, the respective errors for the ionization potentials, in eV, are 0.60, 0.55, 0.48, and 0.45, and the respective errors, in angstroms, for the interatomic distances are 0.036, 0.029, 0.037, and 0.027. The RM1 average error in bond angles of 6.82 degrees is only slightly higher than the AM1 figure of 5.88 degrees, and both are much smaller than the PM3 and PM5 figures of 6.98 degrees and 9.83 degrees, respectively. Moreover, a known error in PM3 nitrogen charges is corrected in RM1. Therefore, RM1 represents an improvement over AM1 and its similar successor PM3, and is probably very competitive with PM5, which is a somewhat different model, and not fully disclosed. RM1 possesses the same analytical construct and the same number of parameters for each atom as AM1, and, therefore, can be easily implemented in any software that already has AM1, not requiring any change in any line of code, with the sole exception of the values of the parameters themselves.

ULTRAFAST CHEMISTRY: Using Time-Resolved Vibrational Spectroscopy for Interrogation of Structural Dynamics

Time-resolved infrared (IR) and Raman spectroscopy elucidates molecular structure evolution during ultrafast chemical reactions. Following vibrational marker modes in real time provides direct insight into the structural dynamics, as is evidenced in studies on intramolecular hydrogen transfer, bimolecular proton transfer, electron transfer, hydrogen bonding during solvation dynamics, bond fission in organometallic compounds and heme proteins, cis-trans isomerization in retinal proteins, and transformations in photochromic switch pairs. Femtosecond IR spectroscopy monitors the site-specific interactions in hydrogen bonds. Conversion between excited electronic states can be followed for intramolecular electron transfer by inspection of the fingerprint IR- or Raman-active vibrations in conjunction with quantum chemical calculations. Excess internal vibrational energy, generated either by optical excitation or by internal conversion from the electronic excited state to the ground state, is observable through transient frequency shifts of IR-active vibrations and through nonequilibrium populations as deduced by Raman resonances.

Solvation Dynamics in Acetonitrile: A Study Incorporating Solute Electronic Response and Nuclear Relaxation

The solvent reorganization process after electronic excitation of a polar solute in a polar solvent such as acetonitrile is related mainly to the time evolution of the solute-solvent electrostatic interaction. Modern laser-based techniques have sufficient time resolution to follow this decay in real time, providing information to be confirmed and interpreted by theories and models. We present here a study aimed at the investigation of the different steps involved in the process taking place after a vertical S(0) --> S(1) excitation of a large size chromophore, coumarin 153 (C153), in acetonitrile, from both the solute and the solvent points of view. To do this, we use accurate quantum mechanical calculations for the solute properties within the polarizable continuum model (PCM) and classical molecular dynamics (MD) simulations, both equilibrium and nonequilibrium, for C153 in the presence of the solvent. The geometry of the solute is allowed to change in order to study the role of internal motions in the time-dependent solvation process. The solvent response function has been obtained from the simulation data and compared to experiment, while the comparison between equilibrium and nonequilibrium MD results for the solvation response confirms the validity of the linear response approximation in the C153-acetonitrile system. The MD trajectories have also been used to monitor the structure of the solvation shell and to determine its change in response to the change in the solute partial charges.

The theory of electron transfer reactions: what may be missing?

Molecular dynamics simulations are presented for condensed-phase electron transfer (ET) systems where the electronic polarizability of both the solvent and the solute is incorporated. The solute polarizability is allowed to change with electronic transition. The results display notable deviation from the standard free energy parabolas of traditional ET theories. A new three-parameter ET model is applied, and the theory is shown to accurately model the free energy surfaces. This paper presents conclusive evidence that the traditional theory for the free energy barrier of ET reactions requires modification.