Vibrational Spectroscopy and Proton Transfer Dynamics in Protonated Oxalate.

The dynamics and infrared spectroscopic signatures of proton transfer in protonated oxalate (p-Oxa) are studied using classical and quantum dynamics. The intermolecular interactions are described by a force field suitable to follow proton transfer. This allows to carry out multiple extended classical molecular dynamics (MD) and ring polymer MD simulations from which the infrared spectrum is determined. Simulations at 600 K sample the quantum mechanical ground state probability distribution and best reproduce the experimentally observed maximum absorption wavelength and part of the line shape. Comparison with the experimentally measured spectrum provides an estimate for the barrier height for proton transfer which can not be determined directly from experiment. A barrier of 4.2 kcal/mol is found to best reproduce the position and width of the infrared absorption of the transferring proton in p-Oxa and also leads to an infrared (IR) spectrum in good agreement with experiment for the deuterated species d-Oxa. A novel means to capture the two resonance forms of oxalate depending on the localization of the excess proton on either CO moiety is found to yield improved results for the spectroscopy in the framework region between 1000 and 2000 cm-1.

Spectroscopy and dynamics of double proton transfer in formic acid dimer.

We present the isolated gas phase infrared spectra of formic acid dimer, (HCOOH)2, and its deuterated counterpart formic-d acid, (DCOOH)2, at room temperature. The formic acid dimer spectrum was obtained by spectral subtraction of a spectrum of formic acid vapor recorded at low pressure from that recorded at a higher pressure. The spectra of formic acid vapor contain features from both formic acid monomer and formic acid dimer, but at low and high pressures of formic acid, the equilibrium is pushed towards the monomer and dimer, respectively. A similar approach was used for the formic-d acid dimer. Building on the previous development of the Molecular Mechanics with Proton Transfer (MMPT) force field for simulating proton transfer reactions, molecular dynamics (MD) simulations were carried out to interpret the experimental spectra in the OH-stretching region. Within the framework of MMPT, a combination of symmetric single and double minimum potential energy surfaces (PESs) provides a good description of the double proton transfer PES. In a next step, potential morphing together with electronic structure calculations at the B3LYP and MP2 level of theory was used to align the computed and experimentally observed spectral features in the OH-stretching region. From this analysis, a barrier for double proton transfer between 5 and 7 kcal mol(-1) was derived, which compares with a CCSD(T)/aug-cc-pVTZ calculated barrier of 7.9 kcal mol(-1). Such a combination of experimental and computational techniques for estimating barriers for proton transfer in gas phase systems is generic and holds promise for further improved PESs and energetics of these important systems. Additional MD simulations at the semi-empirical DFTB level of theory agree quite well for the center band position but underestimate the width of the OH-stretching band.

Multicomponent kinetic analysis and theoretical studies on the phenolic intermediates in the oxidation of eugenol and isoeugenol catalyzed by laccase.

Laccase catalyzes the oxidation of natural phenols and thereby is believed to initialize reactions in lignification and delignification. Numerous phenolic mediators have also been applied in laccase-mediator systems. However, reaction details after the primary O-H rupture of phenols remain obscure. In this work two types of isomeric phenols, EUG (eugenol) and ISO (trans-/cis-isoeugenol), were used as chemical probes to explore the enzymatic reaction pathways, with the combined methods of time-resolved UV-Vis absorption spectra, MCR-ALS, HPLC-MS, and quantum mechanical (QM) calculations. It has been found that the EUG-consuming rate is linear to its concentration, while the ISO not. Besides, an o-methoxy quinone methide intermediate, (E/Z)-4-allylidene-2-methoxycyclohexa-2,5-dienone, was evidenced in the case of EUG with the UV-Vis measurement, mass spectra and TD-DFT calculations; in contrast, an ISO-generating phenoxyl radical, a (E/Z)-2-methoxy-4-(prop-1-en-1-yl) phenoxyl radical, was identified in the case of ISO. Furthermore, QM calculations indicated that the EUG-generating phenoxyl radical (an O-centered radical) can easily transform into an allylic radical (a C-centered radical) by hydrogen atom transfer (HAT) with a calculated activation enthalpy of 5.3 kcal mol(-1) and then be fast oxidized to the observed eugenol quinone methide, rather than an O-radical alkene addition with barriers above 12.8 kcal mol(-1). In contrast, the ISO-generating phenoxyl radical directly undergoes a radical coupling (RC) process, with a barrier of 4.8 kcal mol(-1), while the HAT isomerization between O- and C-centered radicals has a higher reaction barrier of 8.0 kcal mol(-1). The electronic conjugation of the benzyl-type radical and the aromatic allylic radical leads to differentiation of the two pathways. These results imply that competitive reaction pathways exist for the nascent reactive intermediates generated in the laccase-catalyzed oxidation of natural phenols, which is important for understanding the lignin polymerization and may shed some light on the development of efficient laccase-mediator systems.

Mechanistic insights into the effects of diuron exposure on Alexandrium pacificum.

Diuron (N-(3,4-dichlorophenyl)-N,N­dimethylurea, DCMU), a ureic herbicide, is extensively used in agriculture to boost crop productivity; however, its extensive application culminates in notable environmental pollution, especially in aquatic habitats. Therefore, the present study investigated the effect of diuron on the dinoflagellate Alexandrium pacificum, which is known to induce harmful algal blooms (HAB), and its potential to biodegrade DCMU. Following a four-day DCMU exposure, our results revealed that A. pacificum proficiently assimilated DCMU at concentrations of 0.05 mg/L and 0.1 mg/L in seawater, attaining a complete reduction (100 % efficiency) after 96 h for both concentrations. Moreover, evaluations of paralytic shellfish toxins content indicated that cells subjected to higher DCMU concentrations (0.1 mg/L) exhibited reductions of 73.4 %, 86.7 %, and 75 % in GTX1, GTX4, and NEO, respectively. Exposure to DCMU led to a notable decrease in A. pacificum's photosynthetic efficacy, accompanied by increased levels of reactive oxygen species (ROS) and suppressed cell growth, with a growth inhibition rate of 41.1 % at 72 h. Proteomic investigations pinpointed the diminished expression levels of specific proteins like SxtV and SxtW, linked to paralytic shellfish toxins (PSTs) synthesis, as well as key proteins associated with Photosystem II, namely PsbA, PsbD, PsbO, and PsbU. Conversely, proteins central to the cysteine biosynthesis pathways exhibited enhanced expression. In summary, our results preliminarily resolved the molecular mechanisms underlying the response of A. pacificum to DCMU and revealed that DCMU affected the synthesis of PSTs. Meanwhile, our data suggested that A. pacificum has great potential in scavenging DCMU.

Identification and characterization of an anti-fibrotic benzopyran compound isolated from mangrove-derived Streptomyces xiamenensis.

An anti-fibrotic compound produced by Streptomycesn xiamenensis, found in mangrove sediments, was investigated for possible therapeutic effects against fibrosis. The compound, N-[[3,4-dihydro-3S-hydroxy-2S-methyl-2-(4'R-methyl-3'S-pentenyl)-2H-1-benzopyran-6-yl]carbonyl]-threonine (1), was isolated from crude extracts and its structure, including the absolute configuration was determined by extensive spectroscopic data analyses, Mosher's method, Marfey's reagent and quantum mechanical calculations. In terms of biological effects, this compound inhibits the proliferation of human lung fibroblasts (WI26), blocks adhesion of human acute monocytic leukemia cells (THP-1) to a monolayer of WI26 cells, and reduces the contractile capacity of WI26 cells in three-dimensional free-floating collagen gels. Altogether, these data indicate that we have identified a bioactive alkaloid (1) with multiple inhibitory biological effects on lung excessive fibrotic characteristics, that are likely involved in fibrosis, suggesting that this molecule might indeed have therapeutic potential against fibrosis.

PTEN Lipid Phosphatase Activity Enhances Dengue Virus Production through Akt/FoxO1/Maf1 Signaling.

Dengue virus (DENV) is an arthropod-borne viral pathogen and a global health burden. Knowledge of the DENV-host interactions that mediate virus pathogenicity remains limited. Host lipid metabolism is hijacked by DENV for virus replication in which lipid droplets (LDs) play a key role during the virus lifecycle. In this study, we reveal a novel role for phosphatase and tensin homolog deleted on chromosome 10 (PTEN) in LDs-mediated DENV infection. We demonstrate that PTEN expression is downregulated upon DENV infection through post-transcriptional regulation and, in turn, PTEN overexpression enhances DENV replication. PTEN lipid phosphatase activity was found to decrease cellular LDs area and number through Akt/FoxO1/Maf1 signaling, which, together with autophagy, enhanced DENV replication and virus production. We therefore provide mechanistic insight into the interaction between lipid metabolism and the DENV replication cycle.

Multistate Reactive Molecular Dynamics Simulations of Proton Diffusion in Water Clusters and in the Bulk.

The molecular mechanics with proton transfer (MMPT) force field is combined with multistate adiabatic reactive molecular dynamics (MS-ARMD) to describe proton transport in the condensed phase. Parametrization for small protonated water clusters based on electronic structure calculations at the MP2/6-311+G(2d,2p) level of theory and refinement by comparing with infrared spectra for a protonated water tetramer yields a force field which faithfully describes the minimum energy structures of small protonated water clusters. In protonated water clusters up to (H2O)100H+, the proton hopping rate is around 100 hops/ns. This rate converges for 21 ≤ n ≤ 31, and no further speedup in bulk water is found. This indicates that bulklike behavior requires the solvation of a Zundel motif by ∼25 water molecules, which corresponds to the second solvation sphere. For smaller cluster sizes, the number of available states (i.e., the number of proton acceptors) is too small and slows down proton-transfer rates. The cluster simulations confirm that the excess proton is typically located on the surface. The free-energy surface as a function of the weights of the two lowest states and a configurational parameter suggests that the "special pair" plays a central role in rapid proton transport. The barriers between this minimum-energy structure and the Zundel and Eigen minima are sufficiently low (∼1 kcal/mol, consistent with recent experiments and commensurate with a hopping rate of ∼100/ns or 1 every 10 ps), leading to a highly dynamic environment. These findings are also consistent with recent experiments which find that Zundel-type hydration geometries are prevalent in bulk water.

Kinetic isotope effects and how to describe them.

We review several methods for computing kinetic isotope effects in chemical reactions including semiclassical and quantum instanton theory. These methods describe both the quantization of vibrational modes as well as tunneling and are applied to the ⋅H + H2 and ⋅H + CH4 reactions. The absolute rate constants computed with the semiclassical instanton method both using on-the-fly electronic structure calculations and fitted potential-energy surfaces are also compared directly with exact quantum dynamics results. The error inherent in the instanton approximation is found to be relatively small and similar in magnitude to that introduced by using fitted surfaces. The kinetic isotope effect computed by the quantum instanton is even more accurate, and although it is computationally more expensive, the efficiency can be improved by path-integral acceleration techniques. We also test a simple approach for designing potential-energy surfaces for the example of proton transfer in malonaldehyde. The tunneling splittings are computed, and although they are found to deviate from experimental results, the ratio of the splitting to that of an isotopically substituted form is in much better agreement. We discuss the strengths and limitations of the potential-energy surface and based on our findings suggest ways in which it can be improved.

Ultrafast dynamics induced by the interaction of molecules with electromagnetic fields: Several quantum, semiclassical, and classical approaches.

Several strategies for simulating the ultrafast dynamics of molecules induced by interactions with electromagnetic fields are presented. After a brief overview of the theory of molecule-field interaction, we present several representative examples of quantum, semiclassical, and classical approaches to describe the ultrafast molecular dynamics, including the multiconfiguration time-dependent Hartree method, Bohmian dynamics, local control theory, semiclassical thawed Gaussian approximation, phase averaging, dephasing representation, molecular mechanics with proton transfer, and multipolar force fields. In addition to the general overview, some focus is given to the description of nuclear quantum effects and to the direct dynamics, in which the ab initio energies and forces acting on the nuclei are evaluated on the fly. Several practical applications, performed within the framework of the Swiss National Center of Competence in Research "Molecular Ultrafast Science and Technology," are presented: These include Bohmian dynamics description of the collision of H with H2, local control theory applied to the photoinduced ultrafast intramolecular proton transfer, semiclassical evaluation of vibrationally resolved electronic absorption, emission, photoelectron, and time-resolved stimulated emission spectra, infrared spectroscopy of H-bonding systems, and multipolar force fields applications in the condensed phase.

Vibrational relaxation and energy migration of N-methylacetamide in water: the role of nonbonded interactions.

Nonequilibrium molecular dynamics (MD) simulations together with physics-based force fields are used to follow energy flow between vibrationally excited N-methylacetamide (NMA) and water. The simulations are carried out with a previously validated force field for NMA, based on a multipolar representation of the electrostatics, and with a new fluctuating point charge model. For the water solvent, a flexible and a rigid model was employed to distinguish between the role of inter- and intramolecular degrees of freedom. On a 10 ps time scale about 90% of the available energy goes into the solvent. The remaining energy resides within internal NMA-degrees of freedom from where energy flow takes place on longer time scales. The total amount of energy transferred to the solvent on the 10 ps time scale does not depend on whether the water molecules are rigid or flexible during the simulations. Vibrational energy relaxation time scales include two regimes: one on the several 100 fs time scale and a longer one, ranging from 6 to 10 ps. This longer time scale agrees with previous simulations but overestimates the experimentally determined relaxation time by a factor of 2, which can be explained by the classical treatment of the vibrations. Including a previously determined quantum correction factor brings the long time scale into quite favorable agreement with experiment. Coupling to the bending vibration of the water molecules in H-bonding contact with the excited C═O chromophore is substantial. The equilibrium and nonequilibrium distribution of the bending angles of the water molecules in contact with the local oscillator are non-Gaussian, and one approaches the other on the subpicosecond time scale. Analysis of the water velocity distribution suggests that the C═O vibrational energy relaxes into the solvent water shells in an impulsive fashion on a picosecond time scale.