Molecular dynamics simulation of hydrated DPPC monolayers using charge equilibration force fields

We present results of molecular dynamics simulations of a model DPPC-water monolayer using charge equilibration (CHEQ) force fields, which explicitly account for electronic polarization in a classical treatment of intermolecular interactions. The surface pressure, determined as the difference between the monolayer and pure water surface tensions at 323 K, is predicted to be 22.92 ±1.29 dyne/cm, just slightly below the broad range of experimental values reported for this system. The surface tension for the DPPC-water monolayer is predicted to be 42.35 ±1.16 dyne/cm, in close agreement with the experimentally determined value of 40.9 dyne/cm. This surface tension is also consistent with the value obtained from DPPC monolayer simulations using state-of-the-art nonpolarizable force fields. The current results of simulations predict a monolayer-water potential difference relative to the pure water-air interface of 0.64 ±0.02 Volts, an improved prediction compared to the fixed-charge CHARMM27 force field, yet still overestimating the experimental range of 0.3 to 0.45 Volts. As the charge equilibration model is a purely charge-based model for polarization, the current results suggest that explicitly modeled polarization effects can offer improvements in describing interfacial electrostatics in such systems.     

Identification of Intermediates in Zeolite‐Catalyzed Reactions by In Situ UV/Vis Microspectroscopy and a Complementary Set of Molecular Simulations

The optical absorption properties of (poly)aromatic hydrocarbons occluded in a nanoporous environment were investigated by theoretical and experimental methods. The carbonaceous species are an essential part of a working catalyst for the methanol‐to‐olefins (MTO) process. In situ UV/Vis microscopy measurements on methanol conversion over the acidic solid catalysts H‐SAPO‐34 and H‐SSZ‐13 revealed the growth of various broad absorption bands around 400, 480, and 580 nm. The cationic nature of the involved species was determined by interaction of ammonia with the methanol‐treated samples. To determine which organic species contribute to the various bands, a systematic series of aromatics was analyzed by means of time‐dependent density functional theory (TDDFT) calculations. Static gas‐phase simulations revealed the influence of structurally different hydrocarbons on the absorption spectra, whereas the influence of the zeolitic framework was examined by using supramolecular models within a quantum mechanics/molecular mechanics framework. To fully understand the origin of the main absorption peaks, a molecular dynamics (MD) study on the organic species trapped in the inorganic host was essential. During such simulation the flexibility is fully taken into account and the effect on the UV/Vis spectra is determined by performing TDDFT calculations on various snapshots of the MD run. This procedure allows an energy absorption scale to be provided and the various absorption bands determined from in situ UV/Vis spectra to be assigned to structurally different species.     

A fast parallel clustering algorithm for molecular simulation trajectories

We implemented a GPU‐powered parallel k‐centers algorithm to perform clustering on the conformations of molecular dynamics (MD) simulations. The algorithm is up to two orders of magnitude faster than the CPU implementation. We tested our algorithm on four protein MD simulation datasets ranging from the small Alanine Dipeptide to a 370‐residue Maltose Binding Protein (MBP). It is capable of grouping 250,000 conformations of the MBP into 4000 clusters within 40 seconds. To achieve this, we effectively parallelized the code on the GPU and utilize the triangle inequality of metric spaces. Furthermore, the algorithm's running time is linear with respect to the number of cluster centers. In addition, we found the triangle inequality to be less effective in higher dimensions and provide a mathematical rationale. Finally, using Alanine Dipeptide as an example, we show a strong correlation between cluster populations resulting from the k‐centers algorithm and the underlying density. © 2012 Wiley Periodicals, Inc.     

Rotamer decomposition and protein dynamics: Efficiently analyzing dihedral populations from molecular dynamics

Molecular mechanics methods have matured into powerful methods to understand the dynamics and flexibility of macromolecules and especially proteins. As multinanosecond to microsecond length molecular dynamics (MD) simulations become commonplace, advanced analysis tools are required to generate scientifically useful information from large amounts of data. Some of the key degrees of freedom to understand protein flexibility and dynamics are the amino acid residue side chain dihedral angles. In this work, we present an easily automated way to summarize and understand the relevant dihedral populations. A tremendous reduction in complexity is achieved by describing dihedral timeseries in terms of histograms decomposed into Gaussians. Using the familiar and widely studied protein lysozyme, it is demonstrated that our approach captures essential properties of protein structure and dynamics. A simple classification scheme is proposed that indicates the rotational state population for each dihedral angle of interest and allows a decision if a given side chain or peptide backbone fragment remains rigid during the course of an MD simulation, adopts a converged distribution between conformational substates or has not reached convergence yet. © 2012 Wiley Periodicals, Inc.     

Two‐dimensional replica‐exchange method for predicting protein–ligand binding structures

We have developed a two‐dimensional replica‐exchange method for the prediction of protein–ligand binding structures. The first dimension is the umbrella sampling along the reaction coordinate, which is the distance between a protein binding pocket and a ligand. The second dimension is the solute tempering, in which the interaction between a ligand and a protein and water is weakened. The second dimension is introduced to make a ligand follow the umbrella potential more easily and enhance the binding events, which should improve the sampling efficiency. As test cases, we applied our method to two protein‐ligand complex systems (MDM2 and HSP 90‐alpha). Starting from the configuration in which the protein and the ligand are far away from each other in each system, our method predicted the ligand binding structures in excellent agreement with the experimental data from Protein Data Bank much faster with the improved sampling efficiency than the replica‐exchange umbrella sampling method that we have previously developed. © 2013 Wiley Periodicals, Inc.     

Structure and cooperativity of the hydrogen bonds in sodium dihydrogen triacetate

Geometry and energetics of low energy conformers of sodium dihydrogen triacetate (SDHTA) and its anion are studied using density functional theory (DFT) at the Becke, Lee‐Yang‐Parr hybrid functional (BLYP) and Becke, three‐parameter, Lee‐Yang‐Parr hybrid functional (B3LYP) levels. For both cases, two structures of comparable energy are found, which have different symmetry with respect to the two hydrogen bonds (HBs). DFT‐based Born–Oppenheimer molecular dynamics simulations are performed for SDHTA, which show that both structures are visited at room temperature conditions. The trajectory analysis further reveals that the two HBs behave anticooperative, that is, on average elongation of one HB is accompanied by a compression of the other one. This is in accord with nuclear magnetic resonance (NMR) experimental studies for a similar counter ion–dihydrogen triacetate complex. © 2012 Wiley Periodicals, Inc.     

Dispersion of Graphene Sheets in Aqueous Solution by Oligodeoxynucleotides

Applications of graphene sheets in the fields of biosensors and biomedical devices are limited by their insolubility in water. Consequently, understanding the dispersion mechanism of graphene in water and exploring an effective way to prepare stable dispersions of graphene sheets in water is of vital importance for their application in biomaterials, biosensors, biomedical devices, and drug delivery. Herein, a method for stable dispersion of graphene sheets in water by single‐stranded oligodeoxynucleotides (ssODNs) is studied. Owing to van der Waals interactions between graphene sheets, they undergo layer‐to‐layer (LtL) aggregation in water. Molecular dynamics simulations show that, by disrupting van der Waals interaction of graphene sheets with ssODNs, LtL aggregation of graphene sheets is prevented, and water molecules can be distributed stably between graphene sheets. Thus, graphene sheets are dispersed stably in water in the presence of ssODNs. The effects of size and molarity of ssODNs and noncovalent modification of graphene sheets are also discussed.     

Computational Study of the One‐ and Two‐Dimensional Infrared Spectra of a Proton‐Transfer Mode in a Hydrogen‐Bonded Complex Dissolved in a Polar Nanocluster

The signatures of nanosolvation on the one‐ and two‐dimensional (1D and 2D) IR spectra of a proton‐transfer mode in a hydrogen‐bonded complex dissolved in polar solvent molecule nanoclusters of varying size are elucidated by using mixed quantum–classical molecular dynamics simulations. For this particular system, increasing the number of solvent molecules successively from N=7 to N=9 initiates the transition of the system from a cluster state to a bulk‐like state. Both the 1D and 2D IR spectra reflect this transition through pronounced changes in their peak intensities and numbers, but the time‐resolved 2D IR spectra also manifest spectral features that uniquely identify the onset of the cluster‐to‐bulk transition. In particular, it is observed that in the 1D IR spectra, the relative intensities of the peaks change such that the number of peaks decreases from three to two as the size of the cluster increases from N=7 to N=9. In the 2D IR spectra, off‐diagonal peaks are observed in the N=7 and N=8 cases at zero waiting time, but not in the N=9 case. It is known that there are no off‐diagonal peaks in the 2D IR spectrum of the bulk version of this system at zero waiting time, so the disappearance of these peaks is a unique signature of the onset of bulk‐like behavior. Through an examination of the trajectories of various properties of the complex and solvent, it is possible to relate the emergence of these off‐diagonal peaks to an interplay between the vibrations of the complex and the solvent polarization dynamics.