Electrostatics of phosphatidic acid monolayers: Insights from computer simulations

In this paper, we argue that many of the fascinating electrostatic effects that take place in amphiphilic systems are strongly related to the particular organization of the oxygen atoms within each individual molecule. In particular, we focus on two effects: charge inversion and dielectric overscreening. For that purpose, we present molecular dynamics simulations of phosphatidic acid (DMPA²⁻) in the presence of divalent counterions. Our results show that the many oxygens present in DMPA²⁻ cooperatively create strong binding sites for counterions, which in some cases lead to charge inversion. We also present an analysis of the role of interfacial water and relate our analysis to the phenomenon of dielectric overscreening. Several experimental implications are discussed in the conclusions.     

On the relationship betweenT 1 andT 2 for stochastic relaxation models

We consider the relaxation dynamics of two quantum levels coupled to a stochastic bath. We emphasize that even if the matrix elements of the fluctuating Hamiltonian are Gaussian, a second-order cumulant truncation is not exact. For various stochastic models, including the case of a spin-1/2 particle in a fluctuating magnetic field, we calculate 1/T ₁, the population relaxation rate, and 1/T ₂, the phase relaxation rate, up to fourth order in perturbation theory. We show that unlike the commonly accepted second-order result that 1/T ₂1/2T ₁, when fourth-order terms are included, in some instances 1/T ₂<1/2T ₁.     

Spin interaction with an ideal Fermi gas

We consider the equilibrium dynamics of a system consisting of a spin interacting with an ideal Fermi gas on the lattice , 3. We present two examples: when this system is unitarily equivalent to an ideal Fermi gas or to a spin in an ideal Fermi gas without interaction between them.     

Heat conductivity of three periodic Hard Disks via nonequilibrium molecular dynamics

We use a driving field, of the type first suggested by Evans, to generate a steady heat current in the simplest possible system, a two-dimensional periodic fluid of three hard disks. Hard-disk motion equations can be conveniently derived from repulsive constant-force or linear-force potentials by considering the infinitely repulsive limit of these potentials. We show that the isoenergetic and isokinetic forms of the nonequilibrium equations of motion generate steady-state heat conductivities differing by terms of order 1/N, whereN is the number of particles. The resulting conductivities appear to vary as the logarithm of the driving field strength. Even at low fields, the three-body periodic-system results lie well below Enskog's infinite-system prediction.     

室温熔盐二(三氟甲基磺酸酰)亚胺锂-尿素体系的分子动力学模拟

应用分子动力学模拟了25 ℃和50 ℃时新型室温熔盐二(三氟甲基磺酸酰)亚胺锂[LiN(SO2CF3)2, LiTFSI]与尿素(摩尔比为1:3.6)体系的结构与动力学性质. 在两个温度下体系的微观结构基本相同, Li+的配位数约为5, 且都是与溶剂和阴离子中的氧原子发生配位. 对TFSI-的研究表明, 每个TFSI-只提供四个氧中的一个与Li+配位; 而且在Li+的配位层中, TFSI-具有顺、反和gauche 等不同的构象, 并且不同构象出现的几率会随着温度的改变而改变.     

Binding Ensembles of p53-MDM2 Peptide Inhibitors by Combining Bayesian Inference and Atomistic Simulations

Designing peptide inhibitors of the p53-MDM2 interaction against cancer is of wide interest. Computational modeling and virtual screening are a well established step in the rational design of small molecules. But they face challenges for binding flexible peptide molecules that fold upon binding. We look at the ability of five different peptides, three of which are intrinsically disordered, to bind to MDM2 with a new Bayesian inference approach (MELD × MD). The method is able to capture the folding upon binding mechanism and differentiate binding preferences between the five peptides. Processing the ensembles with statistical mechanics tools depicts the most likely bound conformations and hints at differences in the binding mechanism. Finally, the study shows the importance of capturing two driving forces to binding in this system: the ability of peptides to adopt bound conformations (ΔGconformation) and the interaction between interface residues (ΔGinteraction).     

A computational model of the nicotinic acetylcholine binding site

We have derived a model of the nicotinic acetylcholine binding site. This was accomplished by using three known agonists (acetylcholine, nicotine and epibatidine) as templates around which polypeptide side chains, found to be part of the receptor cavity from published molecular biology studies, are allowed to flow freely in molecular dynamics simulations and mold themselves around these templates. The resulting supramolecular complex should thus be a complement, both in terms of steric effects as well as electronic effects, to the agonists and it should be a good estimation of the true receptor cavity structure. The shapes of those minireceptor cavities equilibrated rapidly on the simulation time scale and their structural congruence is very high, implying that a satisfactory model of the nicotinic acetylcholine binding site has been achieved. The computational methodology was internally tested against two rigid and specific antagonists (dihydro--erytroidine and erysoidine), that are expected to give rise to a somewhat differently shaped binding site compared to that derived from the agonists. Using these antagonists as templates there were structural reorganizations of the initial receptor cavities leading to distinctly different cavities compared to agonists. This indicates that adequate times and temperatures were used in our computational protocols to achieve equilibrium structures for the agonists. Overall, both minireceptor geometries for agonists and antagonists are similar with the exception of one amino acid (ARG209).     

Lithium bromide in acetonitrile: Thermodynamics,theory, and simulation

A variety of methods has been used for the study of lithium bromide solutions in acetonitrile yielding by their combination reliable information on different levels of approximation. Osmotic coefficients based on precise vapor pressure measurements are reproduced by CM (chemical model) and HNC (hypernetted chain) calculations and by BD (brownian dynamics) simulations. The results of neutron scattering experiments are treated with the help of HNC and BD methods. Hartree-Fock calculations on isolated LiBr pairs and solvated lithium ions yield reliable particle distances and reveal the geometry of the lithium solvation sphere.     

A force field study of diacetylene molecule by kinetic constant method

The potential constants of diacetylene molecule has been evaluated using kinetic constants. The other molecular constants such as the generalised vibrational mean amplitudes, shrinkage constants, Coriolis coupling constants and centrifugal distortion constants are also calculated using the vibrational frequencies and the results discussed.