Theory of Monte Carlo simulations of the magnetic circular dichroism spectra of alkali metal/rare gas systems

The history of magnetic circular dichroism (MCD) spectroscopy in the study of alkali metal/rare gas (M/Rg) cryogenic systems is reviewed in the context of developing a better understanding of alkali metal/hydrogen systems of current interest to the U.S. Air Force as enhanced‐performance cryogenic rocket propellants. A new theory for simulating the MCD spectra of M/Rg systems is presented together with a careful discussion of the theory's implicit and explicit approximations and their implications. This theory uses a classical Monte Carlo (MC) simulation scheme to model the perturbing effects of the Rg environment on the ²S → ²P MCD‐active transition of the M atom. The theory sets up the MC–MCD simulation as a 6 × 6 matrix eigenvalue/eigenvector problem in the ²P manifold in which are included the effects of M–Rg interactions, metal atom spin‐orbit coupling in the ²P manifold, magnetic Zeeman perturbations of the ²S and ²P manifolds, Boltzmann temperature factors, and electric dipole transition moment integrals for left circularly polarized (LCP) and right circularly polarized (RCP) light. The theory may be applied to any type of trapping site of the host M in the guest Rg matrix; a single atom substitutional metal atom trapping site (one host Rg atom is replaced by one guest M atom) is modeled in this study for M = Na and Rg = Ar. Two temperature factors are used in these simulations; a lattice temperature to model the mobility of the Rg lattice and a magnetic temperature to model Boltzmann factors in the ²S ground manifold. The 6 × 6 eigenvalue/eigenvector problem is solved for a number of randomly generated and suitably averaged Rg configurations to yield the simulated MC‐MCD spectrum for the single substitutional Na/Ar system. The MC–MCD simulations of Na/Ar give the characteristic triplet MCD spectrum with the correct Boltzmann temperature dependence. The simulated MC–MCD spectrum correctly inverts when the direction of the applied magnetic field is reversed. Addition of the LCP and RCP absorbances gives rise to a characteristic ²S → ²P triplet absorption feature. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Molecular dynamics simulations of charged dendrimers: low-to-intermediate half-generation PAMAMs

We have performed atomistic molecular dynamics simulations of PAMAM dendrimers of generations 0.5, 1.5, 2.5, 3.5, and 4.5. The simulated systems comprise the charged dendrimer and its counterions embedded in a dielectric continuum (i.e., without explicit solvent). Structural properties of these dendrimers, like the radius of gyration, the principal moments of inertia, and the segment density profiles, were evaluated from the simulations. The average radius of gyration obtained for the intermediate half-generations 2.5, 3.5, and 4.5 follows the same scaling law that was previously inferred from simulations of full-generation PAMAMs, Rg approximately M1/3, and is characteristic of space-filling objects. The low half-generations 0.5 and 1.5 deviate, however, to greater Rg values. The shape of the smaller dendrimers is approximately that of a prolate ellipsoid, which becomes more spherical for higher generations. The segment density profiles show features identical to those obtained in other simulations of flexible-chain dendrimers, like dendron-backfolding. Two slightly different configurations, in terms of size and shape, were identified for generation 2.5. The radial distributions of counterions extracted from the simulations compare well with the solutions of Poisson-Boltzmann cell model, and the dendrimer's effective charge was estimated using the Bjerrum criterion. The influence of electrostatic interactions in the dendrimer's conformation due to repulsion between the charged end-groups and its relation to counterion effects is discussed for the several generations simulated. The form factors calculated from the simulations are compared with the model of a homogeneous ellipsoid of revolution. The overall results are in agreement with the previously established morphological transition of PAMAM dendrimers toward a more spherical and compact conformation above generations 3 or 4.     

Systematic approach to the quantitative voltammetric analysis of the FeIII/FeII component of the [alpha2-Fe(OH2)P2W17O61]7-/8- reduction process in buffered and unbuffered aqueous media

The one-electron reduction of [alpha(2)-Fe(III)(OH(2))P(2)W(17)O(61)](7-) at a glassy carbon electrode was investigated using cyclic and rotating-disk-electrode voltammetry in buffered and unbuffered aqueous solutions over the pH range 3.45-7.50 with an ionic strength of approximately 0.6 M maintained. The behavior is well-described by a square-scheme mechanism P + e(-) <--> Q (E(1)(0/) = -0.275 V, k(1)(0/) = 0.008 cm s(-1), and alpha(1) = 1/2), PH(+) + e(-) <--> QH(+) (E(2)(0/) = -0.036 V, k(2)(0/) = 0.014 cm s(-1), and alpha(2) = 1/2), PH(+) <--> P + H(+) (K(P) = 3.02 x 10(-6) M), and QH(+) <--> Q + H(+) (K(Q) = 2.35 x 10(-10) M), where P, Q, PH(+), and QH(+) correspond to [alpha(2)-Fe(III)(OH)P(2)W(17)O(61)](8-), [alpha(2)-Fe(II)(OH)P(2)W(17)O(61)](9-), [alpha(2)-Fe(III)(OH(2))P(2)W(17)O(61)](7-), and [alpha(2)-Fe(II)(OH(2))P(2)W(17)O(61)](8-), respectively; E(1)(0)' and E(2)(0)' are the formal potentials, k(1)(0)' and k(2)(0)' are the formal (standard) rate constants, and K(P) and K(Q) are the acid dissociation constants for the relevant reactions. The analysis for the buffered media is based on the approach of Laviron who demonstrated that a square scheme with fully reversible protonations, reversible or quasi reversible electron transfers with the assumption that alpha(1) = alpha(2), can be well-described by the behavior of a simple redox couple, ox + e(-) <--> red, whose formal potential, E(app)(0)', and standard rate constant, k(app)(0)', are straightforwardly derived functions of pH, as are the values of E(1)(0)', k(1)(0)', E(2)(0)', k(2)(0)', and K(P) (only three of the four thermodynamic parameters in a square scheme can be specified). It was assumed that alpha(app) = 1/2, and the simulation program DigiSim was used to determine the values of E(app)(0)' and k(app)(0)', which are required to describe the cyclic voltammograms obtained in buffered media in the pH range from 3.45 to 7.52 (buffer-related reactions which effect general acid-base catalysis are included in the simulations). DigiSim simulations of cyclic voltammograms obtained in unbuffered media yielded the values of E(1)(0)' and k(1)(0)'; K(Q) was then directly computed from thermodynamic constraints. These simulations included additional reactions between the redox species and H(2)O. The value of the diffusion coefficient of the [alpha(2)-Fe(III)(OH(2))P(2)W(17)O(61)](7-), 2.92 x 10(-6) cm(2) s(-1), was determined using DigiSim simulations of voltammograms at a rotating disk electrode in buffered and unbuffered media at pH 3.45. The diffusion coefficients of all redox species were assumed to be identical. When the pH is greater than 6, instability of P (i.e., [alpha(2)-Fe(III)(OH)P(2)W(17)O(61)](8-)) led to the loss of the reactant and precluded lengthy experimentation.     

Nanoscale colloids in a freely adsorbing polymer solution: a Monte Carlo simulation study

A key issue in nanoscale materials and chemical processing is the need for thermodynamic and kinetic models covering colloid-polymer systems over the mesoscopic length scale (approximately 1-100 nm). We have applied Monte Carlo simulations to attractive nanoscale colloid-polymer mixtures toward developing a molecular basis for models of these complex systems. The expanded ensemble Monte Carlo simulation method is applied to calculate colloid chemical potentials (micro(c)) and polymer adsorption (gamma) in the presence of freely adsorbing Lennard-Jones (LJ) homopolymers (surface modifiers). gamma and micro(c) are studied as a function of nanoparticle diameter (sigma(c)), modifier chain length (n) and concentration, and colloid-polymer attractive strength over 0.3 < Rg/sigma(c) < 6 (Rg is the polymer radius of gyration). In the attractive regime, nanocolloid chemical potential decreases and adsorbed amount increases as sigma(c), or n is increased. The scaling of gamma with n from the simulations agrees with the theory of Aubouy and Raphael (Macromolecules 1998, 31, 4357) in the extreme limits of Rg/sigma(c). When Rg/sigma(c) is large, the "colloid" approaches a molecular size and interacts only locally with a few polymer segments and gamma approximately n. When Rg/sigma(c) is small, the system approaches the conventional colloid-polymer size regime where multiple chains interact with a single particle, and gamma approximately sigma(c)2, independent of n. In contrast, adsorption in the mesoscopic range of Rg/sigma(c) investigated here is represented well by a power law gamma approximately n(p), with 0 < p < 1 depending on concentration and LJ attractive strength. Likewise, the chemical potential from our results is fitted well with micro(c) approximately n(q)sigma(c)3, where the cubic term results from the sigma(c) dependence of particle surface area (approximately sigma(c)2) and LJ attractive magnitude (approximately sigma(c)). The q-exponent for micro(c) (micro(c) approximately n(q)) varies with composition and LJ attractive strength but is always very close to the power exponent for gamma (gamma approximately n(p)). This result leads to the conclusion that in attractive systems, polymer adsorption (and thus polymer-colloid attraction) dominates the micro(c) dependence on n, providing a molecular interpretation of the effect of adsorbed organic layers on nanoparticle stability and self-assembly.     

Effective simulations of gas diffusion through kinetically accessible tunnels in multisubunit proteins: O2 pathways and escape routes in T-state deoxyhemoglobin

The diffusion of small gases to special binding sites within polypeptide matrices pivotally defines the biochemical specificity and reactivity of proteins. We investigate here explicit O(2) diffusion in adult human hemoglobin (HbA) as a case study employing the recently developed temperature-controlled locally enhanced sampling (TLES) method and vary the parameters to greatly increase the simulation efficiency. The method is carefully validated against standard molecular dynamics (MD) simulations and available experimental structural and kinetic data on ligand diffusion in T-state deoxyHbA. The methodology provides a viable alternative approach to traditional MD simulations and/or potential of mean force calculations for: (i) characterizing kinetically accessible diffusion tunnels and escape routes for light ligands in porous proteins; (ii) very large systems when realistic simulations require the inclusion of multiple subunits of a protein; and (iii) proteins that access short-lived conformations relative to the simulation time. In the case of T-state deoxyHbA, we find distinct ligand diffusion tunnels consistent with the experimentally observed disparate Xe cavities in the α- and β-subunits. We identify two distal barriers including the distal histidine (E7) that control access to the heme. The multiple escape routes uncovered by our simulations call for a review of the current popular hypothesis on ligand escape from hemoglobin. Larger deviations from the crystal structure during simulated diffusion in isolated α- and β-subunits highlight the dampening effects of subunit interactions and the importance of including all subunits of multisubunit proteins to map realistic kinetically accessible diffusion tunnels and escape routes.     

Incremental update of electrostatic interactions in adaptively restrained particle simulations

The computation of long‐range potentials is one of the demanding tasks in Molecular Dynamics. During the last decades, an inventive panoply of methods was developed to reduce the CPU time of this task. In this work, we propose a fast method dedicated to the computation of the electrostatic potential in adaptively restrained systems. We exploit the fact that, in such systems, only some particles are allowed to move at each timestep. We developed an incremental algorithm derived from a multigrid‐based alternative to traditional Fourier‐based methods. Our algorithm was implemented inside LAMMPS, a popular molecular dynamics simulation package. We evaluated the method on different systems. We showed that the new algorithm's computational complexity scales with the number of active particles in the simulated system, and is able to outperform the well‐established Particle Particle Particle Mesh (P3M) for adaptively restrained simulations. © 2018 Wiley Periodicals, Inc.     

Matching-pursuit split-operator Fourier-transform simulations of excited-state intramolecular proton transfer in 2-(2'-hydroxyphenyl)-oxazole

The excited-state intramolecular proton-transfer dynamics associated with the keto-enolic tautomerization reaction in 2-(2(')-hydroxyphenyl)-oxazole is simulated according to a numerically exact quantum-dynamics propagation method and a full-dimensional excited-state potential energy surface, based on an ab initio reaction surface Hamiltonian. The reported simulations involve the propagation of 35-dimensional wave packets according to the recently developed matching-pursuit/split-operator-Fourier-transform (MP/SOFT) method by Wu and Batista. The underlying propagation scheme recursively applies the time-evolution operator as defined by the Trotter expansion to second order accuracy in dynamically adaptive coherent-state expansions. Computations of time-dependent survival amplitudes, photoabsorption cross sections, and time-dependent reactant(product) populations are compared to the corresponding calculations based on semiclassical approaches, including the Herman-Kluk semiclassical initial value representation method. The reported results demonstrate the capabilities of the MP/SOFT method as a valuble computational tool to study ultrafast reaction dynamics in polyatomic systems as well as to validate semiclassical simulations of complex (nonintegrable) quantum dynamics in multidimensional model systems.     

Infrared spectra and UV-induced photochemistry of methyl aziridine-2-carboxylate isolated in argon and xenon matrices

Methyl aziridine-2-carboxylate (MA2C) has been isolated in low temperature argon and xenon matrices and its structure and photochemistry were studied by FTIR spectroscopy. The reactant as well as the main photoproducts were characterized by comparison of their experimental IR spectra with spectra calculated at the DFT(B3LYP)/6-311++G(d,p) level. The theoretical calculations predicted the existence of two low energy MA2C conformers, differing by the orientation of the OCCN dihedral angle. Both conformers were identified in the studied matrices. Both narrowband tunable and broadband UV irradiations of matrix-isolated MA2C yielded isomerization photoproducts resulting from cleavage of the CC and weakest CN bonds of the aziridine ring. Irradiation with UV laser-light at λ = 235 nm resulted in the formation of the E isomer of methyl 2-(methylimino)-acetate (MMIA) and the Z isomer of methyl 3-iminopropanoate (M3IP). Subsequent irradiation at 290 nm led to observation of new bands resulting from E → Z isomerization of MMIA, while bands due to M3IP remained unchanged. The photoproduced Z isomer of MMIA could be subsequently consumed upon higher-wavelength irradiation (λ = 330 nm). The initially produced MMIA conformer was found to obey the nonequilibrium of excited rotamers (NEER) principle. No photoproducts resulting from the cleavage of the strongest CN bond of the MA2C aziridine ring were observed, nor that of methyl 3-aminoacrylate (M3AA), which could in principle be obtained also by cleavage of the weakest CN bond of the MA2C aziridine ring, but would imply a different H-atom migration simultaneous with the ring opening process. These results indicate that both the differential electronic characteristics of the CN bonds of substituted aziridine rings and the type of required H-atom migration are major factors in determining the specific photochemistries of substituted aziridines. Photofragmentation reactions of MA2C were also observed, through identification of various related products, e.g., acetonitrile, methanol, methane, CO and CO₂.     

高效液相色谱-二极管阵列检测/质谱法分析生脉饮煎剂中的人参皂甙类成分

采用高效液相色谱-二极管阵列检测/质谱（HPLC-DAD/MS）联用技术,以10 mmol/L醋酸铵和乙腈混合溶液梯度洗脱 系统为流动相,应用C18色谱柱对生脉饮煎剂中人参皂甙类成分进行分离鉴定。分析结果表明:生脉饮煎剂中主要含有17个 人参皂甙类成分,即20(R)-人参皂甙Rh1、Rh2、Rg3、Rg2,20(S)-人参皂甙Rh1、Rh2、Rg3、Rg2,人参皂甙Rf、Rg6、Rg5 、F4、Rk1、Rk3、Rh4,20(S)-和20(R)-原人参三醇。人参皂甙成分在煎煮过程中发生了很大变化,主要变成了一些中低 极性产物,这是因为煎煮过程中发生了水解、差向异构、脱水等反应。该方法简便、精确、灵敏度高,可以用来分析生脉 饮煎剂中人参皂甙的变化。     

Simulation of vibrational dephasing of I(2) in solid Kr using the semiclassical Liouville method

In this paper, we present simulations of the decay of quantum coherence between vibrational states of I(2) in its ground (X) electronic state embedded in a cryogenic Kr matrix. We employ a numerical method based on the semiclassical limit of the quantum Liouville equation, which allows the simulation of the evolution and decay of quantum vibrational coherence using classical trajectories and ensemble averaging. The vibrational level-dependent interaction of the I(2)(X) oscillator with the rare-gas environment is modeled using a recently developed method for constructing state-dependent many-body potentials for quantum vibrations in a many-body classical environment [J. M. Riga, E. Fredj, and C. C. Martens, J. Chem. Phys. 122, 174107 (2005)]. The vibrational dephasing rates gamma(0n) for coherences prepared between the ground vibrational state mid R:0 and excited vibrational state mid R:n are calculated as a function of n and lattice temperature T. Excellent agreement with recent experiments performed by Karavitis et al. [Phys. Chem. Chem. Phys. 7, 791 (2005)] is obtained.     

I(+1) transfer from diiodomalonic acid to malonic acid and a complete inhibition of the CO and CO2 evolution in the Briggs-Rauscher reaction by resorcinol

A recent report on an intense CO 2 and CO evolution in the Briggs-Rauscher (BR) reaction revealed that iodination of malonic acid (MA) is not the only important organic reaction in the classical BR oscillator. To disclose the source of the gas evolution, iodomalonic (IMA) and diiodomalonic (I2MA) acids were prepared by iodinating MA with nascent iodine in a semibatch reactor. The nascent iodine was generated by an iodide inflow into the reactor, which contained a mixture of MA and acidic iodate. Some CO2 and a minor CO production was observed during these iodinations. It was found that in an aqueous acidic medium the produced I2MA is not stable but decomposes slowly to diiodoacetic acid and CO2. The first-order rate constant of the I 2MA decarboxylation at 20 degrees C was found to be k1 = 9 x 10(-5) s(-1), which is rather close to the rate constant of the analogous decarboxylation of dibromomalonic acid under similar conditions (7 x 10(-5)s(-1)). From the rate of the CO2 evolution, the I2MA concentration can be calculated in a MA-IMA-I2MA mixture as only I2MA decarboxylates spontaneously but MA and IMA are stable. Following CO2 evolution rates, it was proven that I2MA can react with MA in the reversible reaction I2MA + MA <--> 2 IMA. The equilibrium constant of this reaction was calculated as K = 380 together with the rate constants of the forward k 2 = 6.2 x 10 (-2) M (-1)s(-1) and backward k-2 = 1.6 x 10(-4) M(-1)s(-1) reactions. The probable mechanism of the reaction is I(+1) transfer from I2MA to MA. The presence of I(+1) in a I2MA solution is demonstrated by its reduction with ascorbic acid. To estimate the fraction of CO2 coming from the decarboxylation of I2MA in an oscillatory BR reaction, the oscillations were inhibited by resorcinol. Unexpectedly, all CO2 and CO evolution was interrupted for more than one hour after injecting a small amount of resorcinol (10(-5) M initial concentration in the reactor). Finally, some implications of the newly found I(+1) transfer reactions and the surprisingly effective inhibition by resorcinol regarding the mechanism of the oscillatory BR reaction are discussed. The latter is explained by the ability of resorcinol to scavenge free radicals including iodine atoms without producing iodide ions.     

Peptide hydrolysis catalyzed by matrix metalloproteinase 2: a computational study

The MMP-2 reaction mechanism is investigated by using different computational methodologies. First, quantum mechanical (QM) calculations are carried out on a cluster model of the active site bound to an Ace-Gly approximately Ile-Nme peptide. Along the QM reaction path, a Zn-bound water molecule attacks the Gly carbonyl group to give a tetrahedral intermediate. The breaking of the C-N bond is completed thanks to the Glu 404 residue that shuttles a proton from the water molecule to Ile-N atom. The gas-phase QM energy barrier is quite low ( approximately 14 kcal/mol), thus suggesting that the essential catalytic machinery is included in the cluster model. A similar reaction path occurs in the MMP-2 catalytic domain bound to an octapeptide substrate according to hybrid QM and molecular mechanical (QM/MM) geometry optimizations. However, the rupture of the Gly( P 1) approximately Ile( P 1') amide bond is destabilized in the static QM/MM calculations, owing to the positioning of the Ile( P 1') side chain inside the MMP-2 S 1' pocket and to the inability of simple energy miminization methodologies to properly relax complex systems. Molecular dynamics simulations show that these steric limitations are overcome easily through structural fluctuations. The energetic effect of structural fluctuations is taken into account by combining QM energies with average MM Poisson-Boltzmann free energies, resulting in a total free energy barrier of 14.8 kcal/mol in good agreement with experimental data. The rate-determining event in the MMP-2 mechanism corresponds to a H-bond rearrangement involving the Glu 404 residue and/or the Glu 404-COOH --> N-Ile( P 1') proton transfer. Overall, the present computational results and previous experimental data complement each other well in order to provide a detailed view of the MMPs catalytic mechanism.     

Isobaric vapor–liquid phase equilibria for the binary systems of tert-butanol + 2-ethyl-1-hexanol and n-butanol + 2-ethyl-1-hexanol

Isobaric vapor–liquid equilibria for the binary mixtures of tert-butanol (TBA) + 2-ethyl-1-hexanol and n-butanol (NBA) + 2-ethyl-1-hexanol were experimentally investigated at atmospheric pressure in the temperature range of 353.2–458.2 K. The raw experimental data were correlated using the UNIQUAC and NRTL models and used to estimate the interaction parameters between each pair of components in the systems. The experimental activity coefficients were obtained using the gas chromatographic method and compared with the calculated data obtained from these equilibrium models. The results show that UNIQUAC model gives better correlation than NRTL for these binary systems. The liquid–liquid extraction of TBA from aqueous solution using 2-ethyl-1-hexanol was demonstrated by simulation and the variation of separation factor of TBA at several temperatures was reported.     

高效液相色谱-飞行时间质谱联用法分析人参及人参皂苷与山楂配伍过程中的水解行为

利用高效液相色谱-飞行时间质谱联用的方法,分别对人参配伍山楂前后人参皂苷的变化进行分析,同时对人参皂苷Re、Rg1、Rb1、Rd与山楂配伍的水解规律进行系统研究,并与单独煎煮液、仿山楂配伍pH值煎煮液的水解产物进行比较,结果发现人参与山楂配伍后人参皂苷Rg1、Rb1含量明显减少,而人参皂苷Re、Rd、Rg2、Rg3、F2、Rh1含量明显增加,其中人参皂苷Re与山楂配伍后水解产物为人参皂苷20(R)-Rg2、20(S)-Rg2,仿山楂配伍pH值水解产物为人参皂苷20(R)-Rg2、20(S)-Rg2、Rg4、Rg6;人参皂苷Rg1与山楂配伍后水解产物为20(S)-Rh1、20(R)-Rh1,仿山楂pH值水解产物为20(S)-Rh1、20(R)-Rh1、Rh4、Rk3;人参皂苷Rb1与山楂配伍后水解产物为Rd、20(S)-Rg3,仿山楂pH值水解产物为F2、20(S)-Rg3;人参皂苷Rd与山楂配伍后水解产物为F2、20(S)-Rg3、20(R)-Rg3,仿山楂pH值水解产物为20(S)-Rg3、20(R)-Rg3。研究表明,不同人参皂苷和山楂配伍后与仿山楂pH值的水解产物并不相同,人参与山楂配伍改变了人参皂苷成分的种类及含量。本研究为临床方剂中人参与山楂配伍后成分的变化提供物质基础数据。     

Free energy perturbation (FEP) simulation on the transition states of cocaine hydrolysis catalyzed by human butyrylcholinesterase and its mutants

A novel computational protocol based on free energy perturbation (FEP) simulations on both the free enzyme and transition state structures has been developed and tested to predict the mutation-caused shift of the free energy change from the free enzyme to the rate-determining transition state for human butyrylcholinesterase (BChE)-catalyzed hydrolysis of (-)-cocaine. The calculated shift, denoted by DeltaDeltaG(1 --> 2), of such kind of free energy change determines the catalytic efficiency (kcat/KM) change caused by the simulated mutation transforming enzyme 1 to enzyme 2. By using the FEP-based computational protocol, the DeltaDeltaG(1 --> 2) values for the mutations A328W/Y332A --> A328W/Y332G and A328W/Y332G --> A328W/Y332G/A199S were calculated to be -0.22 and -1.94 kcal/mol, respectively. The calculated DeltaDeltaG(1 --> 2) values predict that the change from the A328W/Y332A mutant to the A328W/Y332G mutant should slightly improve the catalytic efficiency and that the change from the A328W/Y332G mutant to the A328W/Y332G/A199S mutant should significantly improve the catalytic efficiency of the enzyme for the (-)-cocaine hydrolysis. The predicted catalytic efficiency increases are supported by the experimental data showing that kcat/KM = 8.5 x 10(6), 1.4 x 10(7), and 7.2 x 10(7) min(-1) M(-1) for the A328W/Y332A, A328W/Y332G, and A328W/Y332G/A199S mutants, respectively. The qualitative agreement between the computational and experimental data suggests that the FEP simulations may provide a promising protocol for rational design of high-activity mutants of an enzyme. The general computational strategy of the FEP simulation on a transition state can be used to study the effects of a mutation on the activation free energy for any enzymatic reaction.     

Persistence lengths and structure factors of wormlike polymers under confinement

The behavior of semiflexible chains modeling wormlike polymers such as DNA and actin in confined spaces was explored by coarse-grained Monte Carlo simulations. The persistence length P, mean end-to-end distance R2, mean radius of gyration Rg2, and the size ratio R2/Rg2 were computed for chains in slits, cylinders, and spheres. It was found that the intrinsic persistence length of a free chain undergoes on confinement substantial alteration into the apparent persistence length. The qualitative differences were found in trends of the apparent persistence lengths between slits and cylinders on one side and spheres on the other side. The quantities P, R2, Rg2, and R2/Rg2 display similar dependences upon squeezing the chains in nanopores. The above quantities change nonmonotonically with confinement in slits and cylinders, whereas they drop smoothly with decreasing radius of a sphere. For elongation of a chain in a cylinder, two regimes corresponding to strong and moderate confinements were found and compared to experiments and predictions of the blob and Odijk theories. In a spherical cavity, the toroidal chain structure with a hole in the center was detected under strong confinements. The scattering form factor S(q) computed for semiflexible confined chains revealed three regimes of behavior in a slit and a cylinder that matched up well with the scaling theory. The complex form of the function S(q) computed for a sphere was interpreted as a sign of the toroidal structure. A reasonable agreement was found between the simulations and measurements of DNA and actin filaments, confined in nano- and microfluidic channels and spherical droplets, pertaining to the changes of the persistence lengths, chain elongation, and toroidal structure formation.     

Atmospheric-pressure Penning ionization of aliphatic hydrocarbons

A study has been made of the atmospheric-pressure Penning ionization (APPeI) of aliphatic hydrocarbons (pentane, hexane, heptane, and octane) with long-lived rare gas atoms (Rg*). The metastable rare gas atoms (He*, Ne*, Ar* and Kr*) were generated by the negative-mode corona discharge of atmospheric-pressure rare gases. In the Rg*APPeI mass spectra for aliphatic hyrocarbons, the relative abundances of fragment ions were found to increase in the order of He* --> Ne* --> Ar* --> Kr*. The order is in the opposite direction to the internal energies of the Rg*. The less fragmentation observed for He* may be because the nascent molecular ions [M(+.)]* formed by Penning ionization have lifetimes long enough for them to be collisionally deactivated in the atmospheric-pressure ion source. It was found that the relative abundances of fragment ions in Ar*APPeI increased when the sample pressure in the ion source was reduced. This is attributed to the collision of Ar* with molecular ions followed by fragmentation.     

Inelastic scattering matrix elements for the nonadiabatic collision B(2P1/2)+H2(1Sigmag+,j)<-->B(2P3/2)+H2(1Sigmag+,j')

Inelastic scattering matrix elements for the nonadiabatic collision B(2P1/2)+H2(1Sigmag+,j)<-->B(2P3/2)+H2(1Sigmag+,j') are calculated using the time dependent channel packet method (CPM). The calculation employs 1 2A', 2 2A', and 1 2A" adiabatic electronic potential energy surfaces determined by numerical computation at the multireference configuration-interaction level [M. H. Alexander, J. Chem. Phys. 99, 6041 (1993)]. The 1 2A' and 2 2A', adiabatic electronic potential energy surfaces are transformed to yield diabatic electronic potential energy surfaces that, when combined with the total B+H2 rotational kinetic energy, yield a set of effective potential energy surfaces [M. H. Alexander et al., J. Chem. Phys. 103, 7956 (1995)]. Within the framework of the CPM, the number of effective potential energy surfaces used for the scattering matrix calculation is then determined by the size of the angular momentum basis used as a representation. Twenty basis vectors are employed for these calculations, and the corresponding effective potential energy surfaces are identified in the asymptotic limit by the H2 rotor quantum numbers j=0, 2, 4, 6 and B electronic states 2Pja, ja=1/2, 3/2. Scattering matrix elements are obtained from the Fourier transform of the correlation function between channel packets evolving in time on these effective potential energy surfaces. For these calculations the H2 bond length is constrained to a constant value of req=1.402 a.u. and state to state scattering matrix elements corresponding to a total angular momentum of J=1/2 are discussed for j=0<-->j'=0,2,4 and 2P1/2<-->2P1/2, 2P3/2 over a range of total energy between 0.0 and 0.01 a.u.     

Classical and quantum mechanical/molecular mechanical molecular dynamics simulations of alanine dipeptide in water: comparisons with IR and vibrational circular dichroism spectra

We have implemented the combined quantum mechanical (QM)/molecular mechanical (MM) molecular dynamics (MD) simulations of alanine dipeptide in water along with the polarizable and nonpolarizable classical MD simulations with different models of water. For the QM/MM MD simulation, the alanine dipeptide is treated with the AM1 or PM3 approximations and the fluctuating solute dipole moment is calculated by the Mulliken population analysis. For the classical MD simulations, the solute is treated with the polarizable or nonpolarizable AMBER and polarizable CHARMM force fields and water is treated with the TIP3P, TIP4P, or TIP5P model. It is found that the relative populations of right-handed alpha-helix and extended beta and P(II) conformations in the simulation trajectory strongly depend on the simulation method. For the QM/MM MD simulations, the PM3/MM shows that the P(II) conformation is dominant, whereas the AM1/MM predicts that the dominant conformation is alpha(R). Polarizable CHARMM force field gives almost exclusively P(II) conformation and other force fields predict that both alpha-helical and extended (beta and P(II)) conformations are populated with varying extents. Solvation environment around the dipeptide is investigated by examining the radial distribution functions and numbers and lifetimes of hydrogen bonds. Comparing the simulated IR and vibrational circular dichroism spectra with experimental results, we concluded that the dipeptide adopts the P(II) conformation and PM3/MM, AMBER03 with TIP4P water, and AMBER polarizable force fields are acceptable for structure determination of the dipeptide considered in this paper.