CHARMM‐GUI Membrane Builder toward realistic biological membrane simulations

CHARMM‐GUI Membrane Builder, http://www.charmm‐gui.org/input/membrane , is a web‐based user interface designed to interactively build all‐atom protein/membrane or membrane‐only systems for molecular dynamics simulations through an automated optimized process. In this work, we describe the new features and major improvements in Membrane Builder that allow users to robustly build realistic biological membrane systems, including (1) addition of new lipid types, such as phosphoinositides, cardiolipin (CL), sphingolipids, bacterial lipids, and ergosterol, yielding more than 180 lipid types, (2) enhanced building procedure for lipid packing around protein, (3) reliable algorithm to detect lipid tail penetration to ring structures and protein surface, (4) distance‐based algorithm for faster initial ion displacement, (5) CHARMM inputs for P₂₁ image transformation, and (6) NAMD equilibration and production inputs. The robustness of these new features is illustrated by building and simulating a membrane model of the polar and septal regions of E. coli membrane, which contains five lipid types: CL lipids with two types of acyl chains and phosphatidylethanolamine lipids with three types of acyl chains. It is our hope that CHARMM‐GUI Membrane Builder becomes a useful tool for simulation studies to better understand the structure and dynamics of proteins and lipids in realistic biological membrane environments. © 2014 Wiley Periodicals, Inc.     

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.     

FMO‐MD Simulations on the Hydration of Formaldehyde in Water Solution with Constraint Dynamics

Full‐quantum mechanical fragment molecular orbital‐based molecular dynamics (FMO‐MD) simulations were applied to the hydration reaction of formaldehyde in water solution under neutral conditions. Two mechanisms, a concerted and a stepwise one, were considered with respect to the nucleophilic addition and the proton transfer. Preliminary molecular orbital calculations by means of polarized continuum model reaction field predicted that the hydration prefers a concerted mechanism. Because the calculated activation barriers were too high for free FMO‐MD simulations to give reactive trajectories spontaneously, a More O’Ferrall–Jencks‐type diagram was constructed from the statistical analysis of the FMO‐MD simulations with constraint dynamics. The diagram showed that the hydration proceeds through a zwitterionic‐like (ZW‐like) structure. The free energy changes along the reaction coordinate calculated by means of the blue moon ensemble for the hydration and the amination of formaldehyde indicated that the hydration proceeds through a concerted process through the ZW‐like structure, whereas the amination goes through a stepwise mechanism with a ZW intermediate. In inspection of the FMO‐MD trajectories, water‐mediated cyclic proton transfers were observed in both reactions on the way from the ZW‐like structure to the product. These proton transfers also have an asynchronous character, in which deprotonation from the nucleophilic oxygen atom (or nitrogen atom for amination) precedes the protonation of the carbonyl oxygen atom. The results showed the strong advantage of the FMO‐MD simulations to obtain detailed information at a molecular level for solution reactions.     

CHARMM‐GUI ligand reader and modeler for CHARMM force field generation of small molecules

Reading ligand structures into any simulation program is often nontrivial and time consuming, especially when the force field parameters and/or structure files of the corresponding molecules are not available. To address this problem, we have developed Ligand Reader & Modeler in CHARMM‐GUI. Users can upload ligand structure information in various forms (using PDB ID, ligand ID, SMILES, MOL/MOL2/SDF file, or PDB/mmCIF file), and the uploaded structure is displayed on a sketchpad for verification and further modification. Based on the displayed structure, Ligand Reader & Modeler generates the ligand force field parameters and necessary structure files by searching for the ligand in the CHARMM force field library or using the CHARMM general force field (CGenFF). In addition, users can define chemical substitution sites and draw substituents in each site on the sketchpad to generate a set of combinatorial structure files and corresponding force field parameters for throughput or alchemical free energy simulations. Finally, the output from Ligand Reader & Modeler can be used in other CHARMM‐GUI modules to build a protein‐ligand simulation system for all supported simulation programs, such as CHARMM, NAMD, GROMACS, AMBER, GENESIS, LAMMPS, Desmond, OpenMM, and CHARMM/OpenMM. Ligand Reader & Modeler is available as a functional module of CHARMM‐GUI at http://www.charmm-gui.org/input/ligandrm . © 2017 Wiley Periodicals, Inc.     

GridMAT‐MD: A grid‐based membrane analysis tool for use with molecular dynamics

GridMAT‐MD is a new program developed to aid in the analysis of lipid bilayers from molecular dynamics simulations. It reads a GROMACS coordinate file and generates two types of data: a two‐dimensional contour plot depicting membrane thickness, and a polygon‐based tessellation of the individual lipid headgroups. GridMAT‐MD can also account for proteins or small molecules within the headgroups of the lipids, closely approximating their occupied lateral area. The program requires no installation, is fast, and is freely available. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2009     

Acetato(N‐phenylpyridine‐2‐carboxamidato‐κ2N,N)(N‐phenylpyridine‐2‐carboxamide‐κ2N1,O)copper(II)

The title complex, [Cu(C₁₂H₉N₂O)(C₂H₃O₂)(C₁₂H₁₀N₂O)], is a neutral Cu^II complex with a primary N₃O₂ coordination sphere. The Cu centre coordinates to both a deprotonated and a neutral molecule of N‐phenylpyridine‐2‐carboxamide and also to an acetate anion. The coordination around the metal centre is asymmetric, the deprotonated ligand providing two N donor atoms [Cu—N = 1.995 (2) and 2.013 (2) Å] and the neutral ligand providing one N and one O donor atom to the coordination environment [Cu—N = 2.042 (2) Å and Cu—O = 2.2557 (19) Å], the fifth donor being an O atom of the acetate ion [Cu—O = 1.9534 (19) Å]. The remaining O atom from the acetate ion can be considered as a weak donor atom [Cu—O = 2.789 (2) Å], conferring to the Cu complex an asymmetric octahedral geometry. The crystal structure is stabilized by intermolecular N—H...O, C—H...O and C—H...π interactions.     

CHARMM‐GUI Martini Maker for modeling and simulation of complex bacterial membranes with lipopolysaccharides

A complex cell envelope, composed of a mixture of lipid types including lipopolysaccharides, protects bacteria from the external environment. Clearly, the proteins embedded within the various components of the cell envelope have an intricate relationship with their local environment. Therefore, to obtain meaningful results, molecular simulations need to mimic as far as possible this chemically heterogeneous system. However, setting up such systems for computational studies is far from trivial, and consequently the vast majority of simulations of outer membrane proteins still rely on oversimplified phospholipid membrane models. This work presents an update of CHARMM‐GUI Martini Maker for coarse‐grained modeling and simulation of complex bacterial membranes with lipopolysaccharides. The qualities of the outer membrane systems generated by Martini Maker are validated by simulating them in bilayer, vesicle, nanodisc, and micelle environments (with and without outer membrane proteins) using the Martini force field. We expect this new feature in Martini Maker to be a useful tool for modeling large, complicated bacterial outer membrane systems in a user‐friendly manner. © 2017 Wiley Periodicals, Inc.     

Microsecond simulations of DNA and ion transport in nanopores with novel ion–ion and ion–nucleotides effective potentials

We developed a novel scheme based on the grand‐canonical Monte Carlo/Brownian dynamics simulations and have extended it to studies of ion currents across three nanopores with the potential for single‐stranded DNA (ssDNA) sequencing: solid‐state nanopore Si₃N₄, α‐hemolysin, and E111N/M113Y/K147N mutant. To describe nucleotide‐specific ion dynamics compatible with ssDNA coarse‐grained model, we used the inverse Monte Carlo protocol, which maps the relevant ion–nucleotide distribution functions from all‐atom molecular dynamics (MD) simulations. Combined with the previously developed simulation platform for Brownian dynamics simulations of ion transport, it allows for microsecond‐ and millisecond‐long simulations of ssDNA dynamics in the nanopore with a conductance computation accuracy that equals or exceeds that of all‐atom MD simulations. In spite of the simplifications, the protocol produces results that agree with the results of previous studies on ion conductance across open channels and provide direct correlations with experimentally measured blockade currents and ion conductances that have been estimated from all‐atom MD simulations. © 2014 Wiley Periodicals, Inc.     

Computer-aided drug design platform using PyMOL

The understanding and optimization of protein-ligand interactions are instrumental to medicinal chemists investigating potential drug candidates. Over the past couple of decades, many powerful standalone tools for computer-aided drug discovery have been developed in academia providing insight into protein-ligand interactions. As programs are developed by various research groups, a consistent user-friendly graphical working environment combining computational techniques such as docking, scoring, molecular dynamics simulations, and free energy calculations is needed. Utilizing PyMOL we have developed such a graphical user interface incorporating individual academic packages designed for protein preparation (AMBER package and Reduce), molecular mechanics applications (AMBER package), and docking and scoring (AutoDock Vina and SLIDE). In addition to amassing several computational tools under one interface, the computational platform also provides a user-friendly combination of different programs. For example, utilizing a molecular dynamics (MD) simulation performed with AMBER as input for ensemble docking with AutoDock Vina. The overarching goal of this work was to provide a computational platform that facilitates medicinal chemists, many who are not experts in computational methodologies, to utilize several common computational techniques germane to drug discovery. Furthermore, our software is open source and is aimed to initiate collaborative efforts among computational researchers to combine other open source computational methods under a single, easily understandable graphical user interface.     

(E)‐1‐(2′‐Deoxy‐β‐d‐ribo­furan­osyl)‐2,4‐di­fluoro‐5‐(2‐iodo­vinyl)­benzene

This analysis of the title compound, C₁₃H₁₃F₂IO₃, establishes the orientation of (E)‐5‐(CH=CH—I) as antiperiplanar (ap) to the C—C bond (5–6 position) of the 2,4‐di­fluoro­phenyl ring system, with the (E)‐5‐(CH=CH—I) H atom located in close proximity (2.17 Å) to the F4 atom of the 2,4‐di­fluoro­phenyl moiety.     

Binding structures of tri‐N‐acetyl‐β‐glucosamine in hen egg white lysozyme using molecular dynamics with a polarizable force field

Lysozyme is a well‐studied enzyme that hydrolyzes the β‐(1,4)‐glycosidic linkage of N‐acetyl‐β‐glucosamine (NAG)_n oligomers. The active site of hen egg‐white lysozyme (HEWL) is believed to consist of six subsites, A‐F that can accommodate six sugar residues. We present studies exploring the use of polarizable force fields in conjunction with all‐atom molecular dynamics (MD) simulations to analyze binding structures of complexes of lysozyme and NAG trisaccharide, (NAG)₃. MD trajectories are applied to analyze structures and conformation of the complex as well as protein–ligand interactions, including the hydrogen‐bonding network in the binding pocket. Two binding modes (ABC and BCD) of (NAG)₃ are investigated independently based on a fixed‐charge model and a polarizable model. We also apply molecular mechanics with generalized born and surface area (MM‐GBSA) methods based on MD using both nonpolarizable and polarizable force fields to compute binding free energies. We also study the correlation between root‐mean‐squared deviation and binding free energies of the wildtype and W62Y mutant; we find that for this prototypical system, approaches using the MD trajectories coupled with implicit solvent models are equivalent for polarizable and fixed‐charge models. © 2012 Wiley Periodicals, Inc.     

DIRECT‐ID: An automated method to identify and quantify conformational variations—application to β2‐adrenergic GPCR

The conformational dynamics of a macromolecule can be modulated by a number of factors, including changes in environment, ligand binding, and interactions with other macromolecules, among others. We present a method that quantifies the differences in macromolecular conformational dynamics and automatically extracts the structural features responsible for these changes. Given a set of molecular dynamics (MD) simulations of a macromolecule, the norms of the differences in covariance matrices are calculated for each pair of trajectories. A matrix of these norms thus quantifies the differences in conformational dynamics across the set of simulations. For each pair of trajectories, covariance difference matrices are parsed to extract structural elements that undergo changes in conformational properties. As a demonstration of its applicability to biomacromolecular systems, the method, referred to as DIRECT‐ID, was used to identify relevant ligand‐modulated structural variations in the β₂‐adrenergic (β₂AR) G‐protein coupled receptor. Micro‐second MD simulations of the β₂AR in an explicit lipid bilayer were run in the apo state and complexed with the ligands: BI‐167107 (agonist), epinephrine (agonist), salbutamol (long‐acting partial agonist), or carazolol (inverse agonist). Each ligand modulated the conformational dynamics of β₂AR differently and DIRECT‐ID analysis of the inverse‐agonist vs. agonist‐modulated β₂AR identified residues known through previous studies to selectively propagate deactivation/activation information, along with some previously unidentified ligand‐specific microswitches across the GPCR. This study demonstrates the utility of DIRECT‐ID to rapidly extract functionally relevant conformational dynamics information from extended MD simulations of large and complex macromolecular systems. © 2015 Wiley Periodicals, Inc.     

STOCK: Structure mapper and online coarse‐graining kit for molecular simulations

We present a web toolkit STructure mapper and Online Coarse‐graining Kit for setting up coarse‐grained molecular simulations. The kit consists of two tools: structure mapping and Boltzmann inversion tools. The aim of the first tool is to define a molecular mapping from high, for example, all‐atom, to low, that is, coarse‐grained, resolution. Using a graphical user interface it generates input files, which are compatible with standard coarse‐graining packages, for example, Versatile Object‐oriented Toolkit for Coarse‐graining Applications and DL_CGMAP. Our second tool generates effective potentials for coarse‐grained simulations preserving the structural properties, for example, radial distribution functions, of the underlying higher resolution model. The required distribution functions can be provided by any simulation package. Simulations are performed on a local machine and only the distributions are uploaded to the server. The applicability of the toolkit is validated by mapping atomistic pentane and polyalanine molecules to a coarse‐grained representation. Effective potentials are derived for systems of TIP3P (transferable intermolecular potential 3 point) water molecules and salt solution. The presented coarse‐graining web toolkit is available at http://stock.cmm.ki.si . © 2014 Wiley Periodicals, Inc.     

The dynamo library for molecular simulations using hybrid quantum mechanical and molecular mechanical potentials

The Dynamo module library has been developed for the simulation of molecular systems using hybrid quantum mechanical (QM) and molecular mechanical (MM) potentials. Dynamo is not a program package but is a library of Fortran 90 modules that can be employed by those interested in writing their own programs for performing molecular simulations. The library supports a range of different types of molecular calculation including geometry optimizations, reaction‐path determinations and molecular dynamics and Monte Carlo simulations. This article outlines the general structure and capabilities of the library and describes in detail Dynamo's semiempirical QM/MM hybrid potential. Results are presented to indicate three particular aspects of this implementation—the handling of long‐range nonbonding interactions, the nature of the boundary between the quantum mechanical and molecular mechanical atoms and how to perform path‐integral hybrid‐potential molecular dynamics simulations. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 1088–1100, 2000     

Large‐scale asynchronous and distributed multidimensional replica exchange molecular simulations and efficiency analysis

We describe methods to perform replica exchange molecular dynamics (REMD) simulations asynchronously (ASyncRE). The methods are designed to facilitate large scale REMD simulations on grid computing networks consisting of heterogeneous and distributed computing environments as well as on homogeneous high‐performance clusters. We have implemented these methods on NSF (National Science Foundation) XSEDE (Extreme Science and Engineering Discovery Environment) clusters and BOINC (Berkeley Open Infrastructure for Network Computing) distributed computing networks at Temple University and Brooklyn College at CUNY (the City University of New York). They are also being implemented on the IBM World Community Grid. To illustrate the methods, we have performed extensive (more than 60 ms in aggregate) simulations for the beta‐cyclodextrin‐heptanoate host‐guest system in the context of one‐ and two‐dimensional ASyncRE, and we used the results to estimate absolute binding free energies using the binding energy distribution analysis method. We propose ways to improve the efficiency of REMD simulations: these include increasing the number of exchanges attempted after a specified molecular dynamics (MD) period up to the fast exchange limit and/or adjusting the MD period to allow sufficient internal relaxation within each thermodynamic state. Although ASyncRE simulations generally require long MD periods (>picoseconds) per replica exchange cycle to minimize the overhead imposed by heterogeneous computing networks, we found that it is possible to reach an efficiency similar to conventional synchronous REMD, by optimizing the combination of the MD period and the number of exchanges attempted per cycle. © 2015 Wiley Periodicals, Inc.     

Structural properties of D‐mannopyranosyl rings containing O‐acetyl side‐chains

The crystal structures of 1,2,3,4,6‐penta‐O‐acetyl‐α‐d ‐mannopyranose, C₁₆H₂₂O₁₁, and 2,3,4,6‐tetra‐O‐acetyl‐α‐d ‐mannopyranosyl‐(1→2)‐3,4,6‐tri‐O‐acetyl‐α‐d ‐mannopyranosyl‐(1→3)‐1,2,4,6‐tetra‐O‐acetyl‐α‐d ‐mannopyranose, C₄₀H₅₄O₂₇, were determined and compared to those of methyl 2,3,4,6‐tetra‐O‐acetyl‐α‐d ‐mannopyranoside, methyl α‐d ‐mannopyranoside and methyl α‐d ‐mannopyranosyl‐(1→2)‐α‐d ‐mannopyranoside to evaluate the effects of O‐acetylation on bond lengths, bond angles and torsion angles. In general, O‐acetylation exerts little effect on the exo‐ and endocyclic C—C and endocyclic C—O bond lengths, but the exocyclic C—O bonds involved in O‐acetylation are lengthened by ～0.02 Å. The conformation of the O‐acetyl side‐chains is highly conserved, with the carbonyl O atom either eclipsing the H atom attached to a 2°‐alcoholic C atom or bisecting the H—C—H bond angle of a 1°‐alcoholic C atom. Of the two C—O bonds that determine O‐acetyl side‐chain conformation, that involving the alcoholic C atom exhibits greater rotational variability than that involving the carbonyl C atom. These findings are in good agreement with recent solution NMR studies of O‐acetyl side‐chain conformations in saccharides. Experimental evidence was also obtained to confirm density functional theory (DFT) predictions of C—O and O—H bond‐length behavior in a C—O—H fragment involved in hydrogen bonding.