Photoisomerization and photodissociation of aniline and 4-methylpyridine

Photoisomerization and photodissociation of aniline and 4-methylpyridine at 193 nm were studied separately using multimass ion imaging techniques. Photofragment translational energy distributions and dissociation rates were measured. Our results demonstrate that more than 23% of the ground electronic state aniline and 10% of 4-methylpyridine produced from the excitation by 193 nm photons after internal conversion isomerize to seven-membered ring isomers, followed by the H atom migration in the seven-membered ring, and then rearomatize to both methylpyridine and aniline prior to dissociation. The significance of this isomerization is that the carbon, nitrogen, and hydrogen atoms belonging to the alkyl or amino groups are involved in the exchange with those atoms in the aromatic ring during the isomerization.     

Approach to nonadiabatic transitions by density matrix evolution and molecular dynamics simulations

Many biological processes are characterized by an essentially quantum dynamical event, such as electron or proton transfer, in a complex classical environment. To treat such processes properly by computer simulation, allowing nonadiabatic transitions involving excited states, we recently developed a density matrix evolution (DME ) method [H. J. C. Berendsen and J. Mavri, J. Phys. Chem, 97 , 13464 (1993)] which simulates the dynamics of quantum systems embedded in a classical environment. The formalism of the method is presented and an overview of the applications ranging from collisions of a quantum harmonic oscillator with noble gas atoms to proton tunneling in a double-well hydrogen bond is given. The methodology for treatment of proton-transfer processes with inclusion of excited states is presented. Future application of the method on biologically interesting processes, such as proton transfer in enzymatic reactions, is discussed. © 1996 John Wiley & Sons, Inc.     

Matching-pursuit/split-operator Fourier-transform simulations of nonadiabatic quantum dynamics

A rigorous and practical approach for simulations of nonadiabatic quantum dynamics is introduced. The algorithm involves a natural extension of the matching-pursuitsplit-operator Fourier-transform (MPSOFT) method [Y. Wu and V. S. Batista, J. Chem. Phys. 121, 1676 (2004)] recently developed for simulations of adiabatic quantum dynamics in multidimensional systems. The MPSOFT propagation scheme, extended to nonadiabatic dynamics, recursively applies the time-evolution operator as defined by the standard perturbation expansion to first-, or second-order, accuracy. The expansion is implemented in dynamically adaptive coherent-state representations, generated by an approach that combines the matching-pursuit algorithm with a gradient-based optimization method. The accuracy and efficiency of the resulting propagation method are demonstrated as applied to the canonical model systems introduced by Tully for testing simulations of dual curve-crossing nonadiabatic dynamics.     

Methyaluminoxane (MAO) polymerization mechanism and kinetic model from ab initio molecular dynamics and electronic structure calculations

MAO is the co-catalyst in metallocene catalytic systems, which are widely used in single-site olefin polymerization due to their high stereoselectivity. To date, the structures of the catalytically active compound or compounds in MAO have eluded researchers. Although many structural models have been proposed, none are generally accepted. In this study, aspects of the formation mechanism of MAO are addressed. Molecular dynamics simulations at the MP2 level of theory were carried out for presumed elementary steps in MAO formation via hydrolysis of trimethylaluminum (TMA). Methane production was observed, in agreement with experiment, as well as intermediate species that are consistent with the known structural features of MAO and similar to isolated and structurally characterized aluminoxanes. A (CH3)3Al-OH2 species, which we denote as TMA-OH2, containing a stable Al-O single bond emerged as the building block molecule. From this species, a hexameric cage was formed and activation barriers for the various reactions were calculated. Three distinct channels were identified for growth beyond the hexameric cage. It was concluded that MAO formation is a step polymerization through a bifunctional monomer, with [(CH3)Al-O] as the structural unit and a kinetic model was proposed. The structures that emerged were in agreement with the crystallographic evidence for aluminoxanes and support the experimental data regarding the MAO chemical composition.     

Interpreting ultrafast molecular fragmentation dynamics with ab initio electronic structure calculations

High-level ab initio electronic structure calculations are used to interpret the fragmentation dynamics of CHBr(2)COCF(3), following excitation with an intense ultrafast laser pulse. The potential energy surfaces of the ground and excited cationic states along the dissociative C-CF(3) bond have been calculated using multireference second order perturbation theory methods. The calculations confirm the existence of a charge transfer resonance during the evolution of a dissociative wave packet on the ground state potential energy surface of the molecular cation and yield a detailed picture of the dissociation dynamics observed in earlier work. Comparisons of the ionic spectrum for two similar molecules support a general picture in which molecules are influenced by dynamic resonances in the cation during dissociation.     

Trajectory-based solution of the nonadiabatic quantum dynamics equations: an on-the-fly approach for molecular dynamics simulations

The non-relativistic quantum dynamics of nuclei and electrons is solved within the framework of quantum hydrodynamics using the adiabatic representation of the electronic states. An on-the-fly trajectory-based nonadiabatic molecular dynamics algorithm is derived, which is also able to capture nuclear quantum effects that are missing in the traditional trajectory surface hopping approach based on the independent trajectory approximation. The use of correlated trajectories produces quantum dynamics, which is in principle exact and computationally very efficient. The method is first tested on a series of model potentials and then applied to study the molecular collision of H with H(2) using on-the-fly TDDFT potential energy surfaces and nonadiabatic coupling vectors.     

The decay mechanism of photoexcited guanine - a nonadiabatic dynamics study

Ab initio surface hopping dynamics calculations were performed for the biologically relevant tautomer of guanine in gas phase excited into the first ππ? state. The results show that the complete population of UV-excited molecules returns to the ground state following an exponential decay within ～220 fs. This value is in good agreement with the experimentally obtained decay times of 148 and 360 fs. No fraction of the population remains trapped in the excited states. The internal conversion occurs in the ππ? state at two related types of conical intersections strongly puckered at the C2 atom. Only a small population of about 5% following an alternative pathway via a nπ? state was found in the dynamics.     

Insights into the phosphoryl-transfer mechanism of cAMP-dependent protein kinase from quantum chemical calculations and molecular dynamics simulations

To investigate the molecular details of the phosphoryl-transfer mechanism catalyzed by cAMP-dependent protein kinase, we performed quantum mechanical (QM) calculations on a cluster model of the active site and molecular dynamics (MD) simulations of a ternary complex of the protein with Mg(2)ATP and a 20-residue peptide substrate. Overall, our theoretical results confirm the participation of the conserved aspartic acid, Asp(166), as an acid/base catalyst in the reaction mechanism catalyzed by protein kinases. The MD simulation shows that the contact between the nucleophilic serine side chain and the carboxylate group of Asp(166) is short and dynamically stable, whereas the QM study indicates that an Asp(166)-assisted pathway is structurally and energetically feasible and is in agreement with previous experimental results.     

Aggregation mechanism of polyglutamine diseases revealed using quantum chemical calculations, fragment molecular orbital calculations, molecular dynamics simulations, and binding free energy calculations

Polyglutamine (polyQ) diseases, including Huntington’s disease (HD), are caused by expansion of polyQ-encoding repeats within otherwise unrelated gene products. The aggregation mechanism of polyQ diseases, the inhibition mechanism of Congo red, and the alleviation mechanism of trehalose were proposed here based on quantum chemical calculations and molecular dynamics simulations. The calculations and simulations revealed the following. The effective molecular bonding is between glutamine (Gln) and Gln (Gln + Gln), between Gln and Congo red (Gln + Congo red), and between Gln and trehalose (Gln + trehalose). The bonding strength is −13.1 kcal/mol for Gln + Gln, −24.4 kcal/mol for Gln + Congo red, and −12.0 kcal/mol for Gln + trehalose. In the polyQ region, both the number of intermolecular Gln + Gln formations and the total calories generated by the Gln + Gln formation are proportional to the number of repetitions of Gln. We propose an aggregation mechanism whose heat generated by the intermolecular Gln + Gln formation causes the pathogeny of polyQ disease. In our aggregation mechanism, this generated heat collapses the host protein and promotes fibrillogenesis. Without contradiction, our mechanism can explain all the experimental results reported to date. Our mechanism can also explain the inhibition mechanism by Congo red as an inhibitor of polyglutamine-induced protein aggregation and the alleviation mechanism by trehalose as an alleviator of that aggregation. The inhibition mechanism by Congo red is explained by the strong interaction with Gln and by the characteristic structure of Congo red.     

Symmetry and equivalence restrictions in electronic structure calculations

A simple method for obtaining MCSCF orbitals and CI natural orbitals adapted to degenerate point groups, with full symmetry and equivalence restrictions, is described. Among several advantages accruing from this method are the ability to perform atomic SCF calculations on states for which the SCF energy expression cannot be written in terms of Coulomb and exchange integrals over real orbitals, and the generation of symmetry-adapted atomic natural orbitals for use in a recently proposed method for basis set contraction.     

The photoisomerization mechanism of azobenzene: a semiclassical simulation of nonadiabatic dynamics

We have simulated the photoisomerization dynamics of azobenzene, taking into account internal conversion and geometrical relaxation processes, by means of a semiclassical surface hopping approach. Both n-->pi* and pi-->pi* excitations and both cis-->trans and trans-->cis conversions have been considered. We show that in all cases the torsion around the N==N double bond is the preferred mechanism. The quantum yields measured are correctly reproduced and the observed differences are explained as a result of the competition between the inertia of the torsional motion and the premature deactivation of the excited state. Recent time-resolved spectroscopic experiments are interpreted in the light of the simulated dynamics.     

Molecular dynamics simulations and free energy calculations on the enzyme 4‐hydroxyphenylpyruvate dioxygenase

4‐Hydroxyphenylpyruvate dioxygenase is a relevant target in both pharmaceutical and agricultural research. We report on molecular dynamics simulations and free energy calculations on this enzyme, in complex with 12 inhibitors for which experimental affinities were determined. We applied the thermodynamic integration approach and the more efficient one‐step perturbation. Even though simulations seem well converged and both methods show excellent agreement between them, the correlation with the experimental values remains poor. We investigate the effect of slight modifications on the charge distribution of these highly conjugated systems and find that accurate models can be obtained when using improved force field parameters. This study gives insight into the applicability of free energy methods and current limitations in force field parameterization. © 2011 Wiley Periodicals, Inc. J Comput Chem 2011     

Conformational manifold of alpha-aminoisobutyric acid (Aib) containing alanine-based tripeptides in aqueous solution explored by vibrational spectroscopy, electronic circular dichroism spectroscopy, and molecular dynamics simulations

Replacement of the alpha-proton of an alanine residue to generate alpha-aminoisobutyric acid (Aib) in alanine-based oligopeptides favors the formation of a 3(10) helix when the length of the oligopeptide is about four to six residues. This research was aimed at experimentally identifying the structural impact of an individual Aib residue in an alanine context of short peptides in water and Aib's influence on the conformation of nearest-neighbor residues. The amide I band profile of the IR, isotropic and anisotropic Raman, and vibrational circular dichroism (VCD) spectra of Ac-Ala-Ala-Aib-OMe, Ac-Ala-Aib-Ala-OMe, and Ac-Aib-Ala-Ala-OMe were measured and analyzed in terms of different structural models by utilizing an algorithm that exploits the excitonic coupling between amide I' modes. The conformational search was guided by the respective 1H NMR and electronic circular dichroism spectra of the respective peptides, which were also recorded. From these analyses, all peptides adopted multiple conformations. Aib predominantly sampled the right-handed and left-handed 3(10)-helix region and to a minor extent the bridge region between the polyproline (PPII) and the helical regions of the Ramachandran plot. Generally, alanine showed the anticipated PPII propensity, but its conformational equilibrium was shifted towards helical conformations in Ac-Aib-Ala-Ala-OMe, indicating that Aib can induce helical conformations of neighboring residues positioned towards the C-terminal direction of the peptide. An energy landscape exploration by molecular dynamics simulations corroborated the results of the spectroscopic studies. They also revealed the dynamics and pathways of potential conformational transitions of the corresponding Aib residues.     

Rotational dynamics of benzene and water in an ionic liquid explored via molecular dynamics simulations and NMR T1 measurements

The rotational dynamics of benzene and water in the ionic liquid (IL) 1-butyl-3-methylimidazolium chloride are studied using molecular dynamics (MD) simulation and NMR T(1) measurements. MD trajectories based on an effective potential are used to calculate the (2)H NMR relaxation time, T(1) via Fourier transform of the relevant rotational time correlation function, C(2R)(t). To compensate for the lack of polarization in the standard fixed-charge modeling of the IL, an effective ionic charge, which is smaller than the elementary charge is employed. The simulation results are in closest agreement with NMR experiments with respect to the temperature and Larmor frequency dependencies of T(1) when an effective charge of ±0.5e is used for the anion and the cation, respectively. The computed C(2R)(t) of both solutes shows a bi-modal nature, comprised of an initial non-diffusive ps relaxation plus a long-time ns tail extending to the diffusive regime. Due to the latter component, the solute dynamics is not under the motional narrowing condition with respect to the prevalent Larmor frequency. It is shown that the diffusive tail of the C(2R)(t) is most important to understand frequency and temperature dependencies of T(1) in ILs. On the other hand, the effect of the initial ps relaxation is an increase of T(1) by a constant factor. This is equivalent to an "effective" reduction of the quadrupolar coupling constant (QCC). Thus, in the NMR T(1) analysis, the rotational time correlation function can be modeled analytically in the form of aexp (-t/τ) (Lipari-Szabo model), where the constant a, the Lipari-Szabo factor, contains the integrated contribution of the short-time relaxation and τ represents the relaxation time of the exponential (diffusive) tail. The Debye model is a special case of the Lipari-Szabo model with a = 1, and turns out to be inappropriate to represent benzene and water dynamics in ILs since a is as small as 0.1. The use of the Debye model would result in an underestimation of the QCC by a factor of 2-3 as a compensation for the neglect of the Lipari-Szabo factor.     

Insights into mechanistic photodissociation of acetyl chloride by ab initio calculations and molecular dynamics simulations

The potential energy surfaces of dissociation and elimination reactions for CH(3)COCl in the ground (S0) and first excited singlet (S1) states have been mapped with the different ab inito calculations. Mechanistic photodissociation of CH(3)COCl has been characterized through the intrinsic reaction coordinate and ab initio molecular dynamics calculations. The alpha-C-C bond cleavage along the S1 pathway leads to the fragments of COCl((2)A' ') and CH(3) ((2)A') in an excited electronic state and a high barrier exists on the pathway. This channel is inaccessible in energy upon photoexcitation of the CH(3)COCl molecules at 236 nm. The S1 alpha-C-Cl bond cleavage yields the Cl((2)P) and CH(3)CO(X(2)A') fragments in the ground state and there is very small or no barrier on the pathway. The S1 alpha-C-Cl bond cleavage proceeds in a time scale of picosecond in the gas phase, followed by CH(3)CO decomposition to CH(3) and CO. The barrier to the C-Cl bond cleavage on the S1 surface is significantly increased by effects of the argon matrix. The S1 alpha-C-Cl bond cleavage in the argon matrix becomes inaccessible in energy upon photoexcitation of CH(3)COCl at 266 nm. In this case, the excited CH(3)COCl(S1) molecules cannot undergo the C-Cl bond cleavage in a short period. The internal conversion from S1 to S0 becomes the dominant process for the CH(3)COCl(S1) molecules in the condensed phase. As a result, the direct HCl elimination in the ground state becomes the exclusive channel upon 266 nm photodissociation of CH(3)COCl in the argon matrix at 11 K.     

Molecular dynamics simulations of structure and dynamics of organic molecular crystals

A set of model compounds covering a range of polarity and flexibility have been simulated using GAFF, CHARMM22, OPLS and MM3 force fields to examine how well classical molecular dynamics simulations can reproduce structural and dynamic aspects of organic molecular crystals. Molecular structure, crystal structure and thermal motion, including molecular reorientations and internal rotations, found from the simulations have been compared between force fields and with experimental data. The MM3 force field does not perform well in condensed phase simulations, while GAFF, CHARMM and OPLS perform very similarly. Generally molecular and crystal structure are reproduced well, with a few exceptions. The atomic displacement parameters (ADPs) are mostly underestimated in the simulations with a relative error of up to 70%. Examples of molecular reorientation and internal rotation, observed in the simulations, include in-plane reorientations of benzene, methyl rotations in alanine, decane, isopropylcyclohexane, pyramidal inversion of nitrogen in amino group and rotation of the whole group around the C-N bond. Frequencies of such dynamic processes were calculated, as well as thermodynamic properties for reorientations in benzene and alanine. We conclude that MD simulations can be used for qualitative analysis, while quantitative results should be taken with caution. It is important to compare the outcomes from simulations with as many experimental quantities as available before using them to study or quantify crystal properties not available from experiment.     

The NMR structure of [Xd(C2)]4 investigated by molecular dynamics simulations

The i‐motif tetrameric structure is built up from two parallel duplexes intercalated in a head‐to‐tail orientation, and held together by hemiprotonated cytosine pairs. Two topologies exist for the i‐motif structure, one with outermost 3′ extremities and the other with outermost 5′ extremities, called the 3′E and 5′E topology, respectively. Since the comparison of sugar and phosphate group interactions between the two topologies is independent of the length of the intercalation motif, the relative stability of the 3′E and 5′E topologies therefore should not depend on this length. Nevertheless, it has been shown that the 3′E topology of the [d(C2)]4 is much more stable than the 5′E topology, and that the former is the only species observed in solution. In order to understand the reason for this atypical behavior, the NMR structure of the [Xd(C2)]4 was determined and analyzed by molecular dynamics simulations. In the NMR structure, the width of the narrow groove is slightly smaller than in previously determined i‐motif structures, which supports the importance of phosphodiester backbone interactions in the structure stability. The simulations show that the stacking of cytosines, essential for the i‐motif stability, is produced by a similar and non‐negative twisting of the phosphodiester backbones. The twisting is induced by an interaction between the backbones; the [Xd(C2)]4 in 5′E topology, exhibiting very limited interaction between the phosphodiester backbones, is thus unstable. Copyright © 2002 John Wiley & Sons, Ltd.