Wonpil Im

CHARMM: The biomolecular simulation program

CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecular simulation program. It has been developed over the last three decades with a primary focus on molecules of biological interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estimators, molecular minimization, dynamics, and analysis techniques, and model‐building capabilities. The CHARMM program is applicable to problems involving a much broader class of many‐particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical‐molecular mechanical force fields, to all‐atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numerous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.

CHARMM‐GUI: A web‐based graphical user interface for CHARMM

CHARMM is an academic research program used widely for macromolecular mechanics and dynamics with versatile analysis and manipulation tools of atomic coordinates and dynamics trajectories. CHARMM‐GUI, http://www.charmm‐gui.org, has been developed to provide a web‐based graphical user interface to generate various input files and molecular systems to facilitate and standardize the usage of common and advanced simulation techniques in CHARMM. The web environment provides an ideal platform to build and validate a molecular model system in an interactive fashion such that, if a problem is found through visual inspection, one can go back to the previous setup and regenerate the whole system again. In this article, we describe the currently available functional modules of CHARMM‐GUI Input Generator that form a basis for the advanced simulation techniques. Future directions of the CHARMM‐GUI development project are also discussed briefly together with other features in the CHARMM‐GUI website, such as Archive and Movie Gallery.

Generalized born model with a simple smoothing function

Based on recent developments in generalized Born (GB) theory that employ rapid volume integration schemes (M. S. Lee, F. R. Salabury, Jr., and C. L. Brooks III, J Chem Phys 2002, 116, 10606) we have recast the calculation of the self‐electrostatic solvation energy to utilize a simple smoothing function at the dielectric boundary. The present GB model is formulated in this manner to provide consistency with the Poisson–Boltzmann (PB) theory previously developed to yield numerically stable electrostatic solvation forces based on finite‐difference methods (W. Im, D. Beglov, and B. Roux, Comp Phys Commun 1998, 111, 59). Our comparisons show that the present GB model is indeed an efficient and accurate approach to reproduce corresponding PB solvation energies and forces. With only two adjustable parameters—a0 to modulate the Coulomb field term, and a1 to include a correction term beyond Coulomb field—the PB solvation energies are reproduced within 1% error on average for a variety of proteins. Detailed analysis shows that the PB energy can be reproduced within 2% absolute error with a confidence of about 95%. In addition, the solvent‐exposed surface area of a biomolecule, as commonly used in calculations of the nonpolar solvation energy, can be calculated accurately and efficiently using the simple smoothing function and the volume integration method. Our implicit solvent GB calculations are about 4.5 times slower than the corresponding vacuum calculations. Using the simple smoothing function makes the present GB model roughly three times faster than GB models, which attempt to mimic the Lee–Richards molecular volume.

CHARMM-GUI Membrane Builder for Mixed Bilayers and Its Application to Yeast Membranes

The CHARMM-GUI Membrane Builder (http://www.charmm-gui.org/input/membrane), an intuitive, straightforward, web-based graphical user interface, was expanded to automate the building process of heterogeneous lipid bilayers, with or without a protein and with support for up to 32 different lipid types. The efficacy of these new features was tested by building and simulating lipid bilayers that resemble yeast membranes, composed of cholesterol, dipalmitoylphosphatidylcholine, dioleoylphosphatidylcholine, palmitoyloleoylphosphatidylethanolamine, palmitoyloleoylphosphatidylamine, and palmitoyloleoylphosphatidylserine. Four membranes with varying concentrations of cholesterol and phospholipids were simulated, for a total of 170 ns at 303.15 K. Unsaturated phospholipid chain concentration had the largest influence on membrane properties, such as average lipid surface area, density profiles, deuterium order parameters, and cholesterol tilt angle. Simulations with a high concentration of unsaturated chains (73%, membrane(unsat)) resulted in a significant increase in lipid surface area and a decrease in deuterium order parameters, compared with membranes with a high concentration of saturated chains (60-63%, membrane(sat)). The average tilt angle of cholesterol with respect to bilayer normal was largest, and the distribution was significantly broader for membrane(unsat). Moreover, short-lived cholesterol orientations parallel to the membrane surface existed only for membrane(unsat). The membrane(sat) simulations were in a liquid-ordered state, and agree with similar experimental cholesterol-containing membranes.

Performance comparison of generalized born and Poisson methods in the calculation of electrostatic solvation energies for protein structures

This study compares generalized Born (GB) and Poisson (PB) methods for calculating electrostatic solvation energies of proteins. A large set of GB and PB implementations from our own laboratories as well as others is applied to a series of protein structure test sets for evaluating the performance of these methods. The test sets cover a significant range of native protein structures of varying size, fold topology, and amino acid composition as well as nonnative extended and misfolded structures that may be found during structure prediction and folding/unfolding studies. We find that the methods tested here span a wide range from highly accurate and computationally demanding PB‐based methods to somewhat less accurate but more affordable GB‐based approaches and a few fast, approximate PB solvers. Compared with PB solvation energies, the latest, most accurate GB implementations were found to achieve errors of 1% for relative solvation energies between different proteins and 0.4% between different conformations of the same protein. This compares to accurate PB solvers that produce results with deviations of less than 0.25% between each other for both native and nonnative structures. The performance of the best GB methods is discussed in more detail for the application for force field‐based minimizations or molecular dynamics simulations.

Continuum Solvation Model: computation of electrostatic forces from numerical solutions to the Poisson-Boltzmann equation

A rigorous formulation of the solvation forces (first derivatives) associated with the electrostatic free energy calculated from numerical solutions of the linearized Poisson-Boltzmann equation on a discrete grid is described. The solvation forces are obtained from the formal solution of the linearized Poisson-Boltzmann equation written in terms of the Green function. An intermediate region for the solute-solvent dielectric boundary is introduced to yield a continuous solvation free energy and accurate solvation forces. A series of numerical tests show that the calculated forces agree extremely well with finite-difference derivatives of the solvation free energy. To gain a maximum efficiency, the nonpolar contribution to the free energy is expressed in terms of the discretized grid used for the electrostatic problem. The current treatment of solvation forces can be used to introduce the influence of a continuum solvation model in molecular mechanics calculations of large biological systems.

Automated Builder and Database of Protein/Membrane Complexes for Molecular Dynamics Simulations

Molecular dynamics simulations of membrane proteins have provided deeper insights into their functions and interactions with surrounding environments at the atomic level. However, compared to solvation of globular proteins, building a realistic protein/membrane complex is still challenging and requires considerable experience with simulation software. Membrane Builder in the CHARMM-GUI website (http://www.charmm-gui.org) helps users to build such a complex system using a web browser with a graphical user interface. Through a generalized and automated building process including system size determination as well as generation of lipid bilayer, pore water, bulk water, and ions, a realistic membrane system with virtually any kinds and shapes of membrane proteins can be generated in 5 minutes to 2 hours depending on the system size. Default values that were elaborated and tested extensively are given in each step to provide reasonable options and starting points for both non-expert and expert users. The efficacy of Membrane Builder is illustrated by its applications to 12 transmembrane and 3 interfacial membrane proteins, whose fully equilibrated systems with three different types of lipid molecules (DMPC, DPPC, and POPC) and two types of system shapes (rectangular and hexagonal) are freely available on the CHARMM-GUI website. One of the most significant advantages of using the web environment is that, if a problem is found, users can go back and re-generate the whole system again before quitting the browser. Therefore, Membrane Builder provides the intuitive and easy way to build and simulate the biologically important membrane system.

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 P21 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.

An Implicit Membrane Generalized Born Theory for the Study of Structure, Stability, and Interactions of Membrane Proteins

Exploiting recent developments in generalized Born (GB) electrostatics theory, we have reformulated the calculation of the self-electrostatic solvation energy to account for the influence of biological membranes. Consistent with continuum Poisson-Boltzmann (PB) electrostatics, the membrane is approximated as an solvent-inaccessible infinite planar low-dielectric slab. The present membrane GB model closely reproduces the PB electrostatic solvation energy profile across the membrane. The nonpolar contribution to the solvation energy is taken to be proportional to the solvent-exposed surface area (SA) with a phenomenological surface tension coefficient. The proposed membrane GB/SA model requires minor modifications of the pre-existing GB model and appears to be quite efficient. By combining this implicit model for the solvent/bilayer environment with advanced computational sampling methods, like replica-exchange molecular dynamics, we are able to fold and assemble helical membrane peptides. We examine the reliability of this model and approach by applications to three membrane peptides: melittin from bee venom, the transmembrane domain of the M2 protein from Influenza A (M2-TMP), and the transmembrane domain of glycophorin A (GpA). In the context of these proteins, we explore the role of biological membranes (represented as a low-dielectric medium) in affecting the conformational changes in melittin, the tilt of transmembrane peptides with respect to the membrane normal (M2-TMP), helix-to-helix interactions in membranes (GpA), and the prediction of the configuration of transmembrane helical bundles (GpA). The present method is found to perform well in each of these cases and is anticipated to be useful in the study of folding and assembly of membrane proteins as well as in structure refinement and modeling of membrane proteins where a limited number of experimental observables are available.

Ion Permeation and Selectivity of OmpF Porin: A Theoretical Study Based on Molecular Dynamics, Brownian Dynamics, and Continuum Electrodiffusion Theory

Three different theoretical approaches are used and compared to refine our understanding of ion permeation through the channel formed by OmpF porin from Escherichia coli. Those approaches are all-atom molecular dynamics (MD) in which ions, solvent, and lipids are represented explicitly, Brownian dynamics (BD) in which ions are represented explicitly, while solvent and lipids are represented as featureless dielectrics, and Poisson-Nernst-Planck (PNP) electrodiffusion theory in which both solvent and local ion concentrations are represented as a continuum. First, the ability of the different theoretical approaches in reproducing the equilibrium average ion density distribution in OmpF porin bathed by a 1M KCl symmetric salt solution is examined. Under those conditions the PNP theory is equivalent to the non-linear Poisson-Boltzmann (PB) theory. Analysis shows that all the three approaches are able to capture the important electrostatic interactions between ions and the charge distribution of the channel that govern ion permeation and selectivity in OmpF. The K(+) and Cl(-) density distributions obtained from the three approaches are very consistent with one another, which suggests that a treatment on the basis of a rigid protein and continuum dielectric solvent is valid in the case of OmpF. Interestingly, both BD and continuum electrostatics reproduce the distinct left-handed twisted ion pathways for K(+) and Cl(-) extending over the length of the pore which were observed previously in MD. Equilibrium BD simulations in the grand canonical ensemble indicate that the channel is very attractive for cations, particularly at low salt concentration. On an average there is 1.55 K(+) inside the pore in 10mM KCl. Remarkably, there is still 0.17 K(+) on average inside the pore even at a concentration as low as 1microM KCl. Secondly, non-equilibrium ion flow through OmpF is calculated using BD and PNP and compared with experimental data. The channel conductance in 0.2M and 1M KCl calculated using BD is in excellent accord with the experimental data. The calculations reproduce the experimentally well-known conductance-concentration relation and also reveal an asymmetry in the channel conductance (a larger conductance is observed under a positive transmembrane potential). Calculations of the channel conductance for three mutants (R168A, R132A, and K16A) in 1M KCl suggest that the asymmetry in the channel conductance arises mostly from the permanent charge distribution of the channel rather than the shape of the pore itself. Lastly, the calculated reversal potential in a tenfold salt gradient (0.1:1M KCl) is 27.4(+/-1.3)mV (BD) and 22.1(+/-0.6)mV (PNP), in excellent accord with the experimental value of 24.3mV. Although most of the results from PNP are qualitatively reasonable, the calculated channel conductance is about 50% higher than that calculated from BD probably because of a lack of some dynamical ion-ion correlations.