Studies on the Structures and Hydrogen-bonding Interactions in Aqueous Glycine Solutions Using All-atom Molecular Dynamics Simulations and NMR Spectroscopy

All-atom molecular dynamics （MD） simulations and chemical shifts were used to study interactions and structures in the glycine-water system. Radial distribution functions and the hydrogen-bond network were applied in MD simulations. Aggregates in the aqueous glycine solution could be classified into different regions by analysis of the hydrogen-bonding network. Temperature-dependent NMR spectra and the viscosity of glycine in aqueous solutions were measured to compare with the results of MD simulations. The variation tendencies of the hydrogen atom chemical shifts and viscosity with concentration of glycine agree with the statistical results of hydrogen bonds in the MD simulations.     

Structures and Hydrogen Bonding Interactions in Urea-water System Studied by All-atom MD Simulation and Chemical Shifts in NMR Spectrum

The interactions and structures of the urea-water system are studied by an all-atom molecular dynamics （MD） simulation. The hydrogen-bonding network and the radial distribution functions are adopted in MD simulations. The structures of urea-water mixtures can be classified into different regions from the analysis of the hydrogen-bonding network. The urea molecule shows the certain tendency to the self-aggregate with the mole fraction of urea increasing. Moreover, the results of the MD simulations are also compare with the chemical shifts and viscosities of the urea aqueous solutions, and the statistical results of the average number hydrogen bonds in the MD simulations are in agreement with the experiment data such as chemical shifts of the hydrogen atom and viscosity.     

Elucidating the origin of diastereoselectivity in a self-replicating system: selfishness versus altruism

We have investigated a diastereoselective self‐replicating system based on a cycloaddition of a fulvene derivative and a maleimide using a two‐pronged approach of combining NMR spectroscopy with computational modelling. Two diastereomers are formed with identical rates in the absence of replication. When replication is enabled, one diastereomer takes over the resources as a “selfish” autocatalyst, while exploiting the competitor as a weak “altruist”, resulting in a diastereoselectivity of 16:1. We applied 1D and 2D NMR spectroscopic techniques supported by ab initio chemical shifts as well as ab initio molecular dynamics simulations to study the structure and dynamics of the underlying network. This powerful combination allowed us to decipher the energetic and structural rationale behind the observed behaviour, while static computational methods currently used in the field did not.     

Prediction of the 1H and 13C NMR spectra of alpha-D-glucose in water by DFT methods and MD simulations

We have applied computational protocols based on DFT and molecular dynamics simulations to the prediction of the alkyl 1H and 13C chemical shifts of alpha-d-glucose in water. Computed data have been compared with accurate experimental chemical shifts obtained in our laboratory. 13C chemical shifts do not show a marked solvent effect. In contrast, the results for 1H chemical shifts provided by structures optimized in the gas phase are only fair and point out that it is necessary to take into account both the flexibility of the glucose structure and the strong effect exerted by solvent water thereupon. Thus, molecular dynamics simulations were carried out to model both the internal geometry as well as the influence of solvent molecules on the conformational distribution of the solute. Snapshots from the simulation were used as input to DFT NMR calculations with varying degrees of sophistication. The most important factor that affects the accuracy of computed 1H chemical shifts is the solute geometry; the effect of the solvent on the shielding constants can be reasonably accounted for by self-consistent reaction field models without the need of explicitly including solvent molecules in the NMR property calculation.     

特殊缔合体系TFE水溶液分子动力学模拟

三氟乙醇(TFE)水溶液是一类特殊的缔合体系. 采用分子动力学模拟方法结合核磁共振化学位移研究了TFE水溶液体系全浓度范围的氢键网络, 并对动力学模拟结果和核磁共振化学位移进行了比较. 从径向分布函数(RDF)发现, TFE水溶液中存在着强氢键, 而体系中的C—H…O弱相互作用较为明显, 也不能忽略. 氢键网络分析发现TFE 水溶液体系的氢键大致分为以下三个区域: 在水富集区域, 水分子倾向于自身缔合形成稳定的簇结构, 随着TFE 浓度的增加, 水的有序结构受到破坏, 水分子和TFE分子发生交叉缔合作用形成氢键; 在TFE富集区域, 水分子较少, TFE分子自身通过氢键形成多缔体结构. 此外, 分子动力学统计的平均氢键数的变化和文献报导的核磁共振化学位移变化趋势相同, 实验和理论的结果吻合较好.     

All‐atom Molecular Dynamic Simulations Combined with the Chemical Shifts Study on the Weak Interactions of Ethanol‐water System

All‐atom molecular dynamics (MD) simulation combined with chemical shifts was performed to investigate the interactions over the entire concentration range of the ethanol (EtOH)‐water system. The results of the simulation were adopted to explain the NMR experiments by hydrogen bonding analysis. The strong hydrogen bonds and weak C–H···O contacts coexist in the mixtures through the analysis of the radial distribution functions. And the liquid structures in the whole concentration of EtOH‐water mixtures can be classified into three regions by the statistic analysis of the hydrogen‐bonding network in the MD simulations. Moreover, the chemical shifts of the hydrogen atom are in agreement with the statistical results of the average number hydrogen bonds in the MD simulations. Interestingly, the excess relative extent of η_rel^E calculated by the MD simulations and chemical shifts in the EtOH aqueous solutions shows the largest deviation at x_EtOH≈0.18. The excess properties present good agreement with the excess enthalpy in the concentration dependence.     

Long‐Lasting Oscillations in the Electro‐Oxidation of Formic Acid on PtSn Intermetallic Surfaces

Even when in contact with virtually infinite reservoirs, natural and manmade oscillators typically drift in phase space on a time‐scale considerably slower than that of the intrinsic oscillator. A ubiquitous example is the inexorable aging process experienced by all living systems. Typical electrocatalytic reactions under oscillatory conditions oscillate for only a few dozen stable cycles due to slow surface poisoning that ultimately results in destruction of the limit cycle. We report the observation of unprecedented long‐lasting temporal oscillations in the electro‐oxidation of formic acid on an ordered intermetallic PtSn phase. The introduction of Sn substantially increases the catalytic activity and retards the irreversible surface oxidation, which results in the stabilization of more than 2200 oscillatory cycles in about 40 h; a 30–40‐fold stabilization with respect to the behavior of pure Pt surfaces. The dynamics were modeled and numerical simulations point to the surface processes underlying the high stability.     

Simulations of the dynamics of thermal undulations in lipid bilayers in the tensionless state and under stress

The relaxation processes of height undulations and density fluctuations in a membrane have been studied by molecular dynamics simulations of a coarse grained amphiphilic bilayer model. We observe a double exponential decay in their time correlations, with relaxation rates in good quantitative agreement with the theory by Seifert and Langer [Europhys. Lett. 23, 71 (1993)]. Intermonolayer friction due to slippage between the two monolayers is shown to be the dominant dissipative mechanism at the high wave numbers, q>10 mum(-1), typically encountered in computer simulations. We briefly discuss the ramifications of the slow undulatory relaxation process for the calculation of bending rigidities from the static undulation structure factors. The relaxation rates are sensitive to the surface tension, and at high elongations an oscillatory contribution is observed in the time correlation of the undulations.     

Non-Newtonian behavior in simple fluids

Using nonequilibrium molecular dynamics simulations, we study the non-Newtonian rheology of a microscopic sample of simple fluid. The calculations were performed using a configurational thermostat which unlike previous nonequilibrium molecular dynamics or nonequilibrium Brownian dynamics methods does not exert any additional constraint on the flow profile. Our findings are in agreement with experimental results on concentrated "hard sphere"-like colloidal suspensions. We observe: (i) a shear thickening regime under steady shear; (ii) a strain thickening regime under oscillatory shear at low frequencies; and (iii) shear-induced ordering under oscillatory shear at higher frequencies. These results significantly differ from previous simulation results which showed systematically a strong ordering for all frequencies. They also indicate that shear thickening can occur even in the absence of a solvent.     

Protonation of the chromophore in the photoactive yellow protein

The photoactive yellow protein (PYP) acts as a light sensor to its bacterial host: it responds to light by changing shape. After excitation by blue light, PYP undergoes several transformations, to partially unfold into its signaling state. One of the crucial steps in this photocycle is the protonation of p-coumaric acid after excitation and isomerization of this chromophore. Experimentalists still debate on the nature of the proton donor and on whether it donates the hydrogen directly or indirectly. To obtain better knowledge of the mechanism, we studied this proton transfer using Car-Parrinello molecular dynamics, classical molecular dynamics, and computer simulations combining these two methods (quantum mechanics/molecular mechanics, QMMM). The simulations reproduce the chromophore structure and hydrogen-bond network of the protein measured by X-ray crystallography and NMR. When the chromophore is protonated, it leaves the assumed proton donor, glutamic acid 46, with a negative charge in a hydrophobic environment. We show that the stabilization of this charge is a very important factor in the mechanism of protonation. Protonation frequently occurs in simplified ab initio simulations of the chromophore binding pocket in vacuum, where amino acids can easily hydrogen bond to Glu46. When the complete protein environment is incorporated in a QMMM simulation on the complete protein, no proton transfer is observed within 14 ps. The hydrogen-bond rearrangements in this time span are not sufficient to stabilize the new protonation state. Force field molecular dynamics simulations on a much longer time scale have shown which internal rearrangements of the protein are needed. Combining these simulations with more QMMM calculations enabled us to check the stability of protonation states and clarify the initial requirements for the proton transfer in PYP.     

Stokesian Dynamics Simulations of Ferromagnetic Colloidal Dispersions Subjected to a Sinusoidal Shear Flow

We have conducted Stokesian dynamics simulations to investigate the dynamic properties of ferromagnetic colloidal dispersions subjected to a sinusoidal shear flow. Thick chain-like cluster formation is significantly influenced by an oscillatory shear flow even if the amplitude is relatively small, since the internal structures of thick chain-like clusters are highly sensitive to the change in the direction of the shear flow. The motion of thick chain-like clusters is out of phase to a sinusoidal shear rate, and the phase difference is strongly correlated with that of the viscosity and normal stress coefficients. The viscoelastic properties become more apparent with decreasing frequency of the oscillatory shear flow, since such properties have a strong relationship with the thick chain-like cluster formation. In other words, since thick chain-like clusters are more stable for the case of a smaller frequency shear flow, such stable clusters induce significant viscoelastic properties of ferromagnetic colloidal dispersions in a strong, applied magnetic field. Copyright 2000 Academic Press.     

Liquid-vapor phase diagram and cluster formation of two-dimensional ionic fluids

Direct molecular dynamics simulations on interfaces at constant temperature are performed to obtain the liquid-vapor phase diagram of the two-dimensional soft primitive model, an equimolar mixture of equal size spheres carrying opposite charges. Constant temperature and pressure simulations are also carried out to check consistency with interface simulations results. In addition, an analysis of the cluster formation of mixtures of particles with charge asymmetry in the range 1:1 to 1:36 at low and high densities is performed. The number of free ions, when plotted as a function of the positive ion charge, Z(+), has an oscillatory behavior and is independent of the density. The formation of aggregates is analyzed in terms of the attraction and repulsion between ions.     

All-atom Molecular Dynamics Simulationsand NMR Spectroscopy Study on Interactions and Structures in N-Glycylglycine Aqueous Solution

All-atom molecular dynamics (MD) simulation and the NMR spectra are used to investi-gate the interactions in N-glycylglycine aqueous solution. Different types of atoms exhibit different capability in forming hydrogen bonds by the radial distribution function analysis. Some typical dominant aggregates are found in different types of hydrogen bonds by the statistical hydrogen-bonding network. Moreover, temperature-dependent NMR are used to compare with the results of the MD simulations. The chemical shifts of the three hydrogen atoms all decrease with the temperature increasing which reveals that the hydrogen bonds are dominant in the glycylglycine aqueous solution. And the NMR results show agreement with the MD simulations. All-atom MD simulations and NMR spectra are successful in revealing the structures and interactions in the N-glycylglycine-water mixtures.     

Oscillatory dynamics protect enzymes and possibly cells against toxic substances

We have used the oscillating peroxidase-oxidase (PO) reaction as a model system to study how oscillatory dynamics may affect the influence of toxic reaction intermediates on enzyme stability. In the peroxidase-oxidase reaction reactive intermediates, such as hydrogen peroxide, superoxide, and hydroxyl radical are formed. Such intermediates inactivate many cellular macromolecules such as proteins and nucleic acids. These reaction intermediates also react with peroxidase itself to form an inactive enzyme. The fact that the PO reaction shows bistability between an oscillatory and a steady state gives us a unique possibility to compare such inactivation when the system is in one of these two states. We show that inactivation of peroxidase is slower when the system is in an oscillatory state, and using numerical simulations we provide evidence that oscillatory dynamics lower the average concentration of the reactive intermediates.     

Interaction of water with the model ionic liquid [bmim][BF4]: molecular dynamics simulations and comparison with NMR data

We report on molecular dynamics simulations of the ionic liquid [bmim][BF 4] and its mixtures with water, from zero up to 0.5 mol fraction of water. All of the simulations are carried out with two published force fields. The results are compared with each other and with published as well as new NMR data on the same mixtures, whenever possible. We perform extensive analyses of structural quantities, such as pair correlation functions, nearest-neighbor analysis and size distribution of the water clusters formed at higher concentrations. We show that the water clusters are formed almost exclusively by linear chains of hydrogen-bonded molecules. There is a nanoscale structuring of the mixtures but no macroscopic phase separation among the components, in agreement with experiment. Roughly, we identify two solvation regimes. At low water content, the ions are selectively coordinated by individual water molecules, but their ionic network is largely unperturbed. At high water content, the ionic network is somewhat disrupted or swollen in a nonspecific way by the water clusters.     

Interpreting protein structural dynamics from NMR chemical shifts

In this investigation, semiempirical NMR chemical shift prediction methods are used to evaluate the dynamically averaged values of backbone chemical shifts obtained from unbiased molecular dynamics (MD) simulations of proteins. MD-averaged chemical shift predictions generally improve agreement with experimental values when compared to predictions made from static X-ray structures. Improved chemical shift predictions result from population-weighted sampling of multiple conformational states and from sampling smaller fluctuations within conformational basins. Improved chemical shift predictions also result from discrete changes to conformations observed in X-ray structures, which may result from crystal contacts, and are not always reflective of conformational dynamics in solution. Chemical shifts are sensitive reporters of fluctuations in backbone and side chain torsional angles, and averaged (1)H chemical shifts are particularly sensitive reporters of fluctuations in aromatic ring positions and geometries of hydrogen bonds. In addition, poor predictions of MD-averaged chemical shifts can identify spurious conformations and motions observed in MD simulations that may result from force field deficiencies or insufficient sampling and can also suggest subsets of conformational space that are more consistent with experimental data. These results suggest that the analysis of dynamically averaged NMR chemical shifts from MD simulations can serve as a powerful approach for characterizing protein motions in atomistic detail.     

Finite-size effects in dissipative particle dynamics simulations

We have performed dissipative particle dynamics (DPD) simulations to evaluate the effect that finite size of transversal area has on stress anisotropy and interfacial tension. The simulations were carried out in one phase and two phases in parallelepiped cells. In one-phase simulations there is no finite-size effect on stress anisotropy when the simulation is performed using repulsive forces. However, an oscillatory function of stress anisotropy is found for attractive-repulsive interactions. In the case of liquid-liquid interfaces with repulsive interaction between molecules, there is only a small effect of surface area on interfacial tension when the simulations are performed using the Monte Carlo method at constant temperature and normal pressure. An important but artificial finite-size effect of interfacial area on surface tension is found in simulations in the canonical ensemble. Reliable results of interfacial tension from DPD simulations can be obtained using small systems, less than 2000 particles, when they interact exclusively with repulsive forces.     

Simulating rare events in equilibrium or nonequilibrium stochastic systems

We present three algorithms for calculating rate constants and sampling transition paths for rare events in simulations with stochastic dynamics. The methods do not require a priori knowledge of the phase-space density and are suitable for equilibrium or nonequilibrium systems in stationary state. All the methods use a series of interfaces in phase space, between the initial and final states, to generate transition paths as chains of connected partial paths, in a ratchetlike manner. No assumptions are made about the distribution of paths at the interfaces. The three methods differ in the way that the transition path ensemble is generated. We apply the algorithms to kinetic Monte Carlo simulations of a genetic switch and to Langevin dynamics simulations of intermittently driven polymer translocation through a pore. We find that the three methods are all of comparable efficiency, and that all the methods are much more efficient than brute-force simulation.     

Excitability in chemical and biochemical pH-autocatalytic systems

Using two different kinds of pH systems--the papain catalyzed hydrolysis of N-benzoyl-L-arginine ethyl ester in a membrane reactor and the bromate-sulfite-ferrocyanide (BSF) reaction in the CSTR--we study the relation among excitability, oscillations and bistability, and the ability of the system to respond to external periodic perturbations. Excitable properties of dynamical systems are examined in terms of a threshold set which is used to characterise dynamics in the reactor subject to external periodic stimuli. A precise definition and a method of calculating the threshold set are formulated. Two kinds of excitability distinguished by either direct or indirect initiation of the activatory process are found in both pH systems. Periodic pulsed perturbations of the BSF system display a nontrivial dependence of an excitation number on the forcing period. We examined this system also in oscillatory mode by looking at the phase shifts caused by single-pulse perturbations and constructing the phase transition curves (PTCs).     

Study on Non-Newtonian Behaviors of Lennard-Jones Fluids via Molecular Dynamics Simulations

Using nonequilibrium molecular dynamics simulations, we study the non-Newtonian rheological behaviors of a monoatomic fluid governed by the Lennard-Jones potential. Both steady Couette and oscillatory shear flows are investigated. Shear thinning and normal stress effects are observed in the steady Couette flow simulations. The radial distribution function is calculated at different shear rates to exhibit the change of the microscopic structure of molecules due to shear. We observe that for a larger shear rate the repulsion between molecules is more powerful while the attraction is weaker, and the above phenomena can also be confirmed by the analyses of the potential energy. By applying an oscillatory shear to the system, several findings are worth mentioning here:First, the phase difference between the shear stress and shear rate increases with the frequency. Second, the real part of complex viscosity first increases and then decreases while the imaginary part tends to increase monotonically, which results in the increase of the proportion of the imaginary part to the real part with the increasing frequency. Third, the ratio of the elastic modulus to the viscous modulus also increases with the frequency. These phenomena all indicate the appearance of viscoelasticity and the domination of elasticity over viscosity at high oscillation frequency for Lennard-Jones fluids.