CHARMM-GUI HMMM Builder for Membrane Simulations with the Highly Mobile Membrane-Mimetic Model

CHARMM-GUI DEER facilitator for spin-pair distance distribution calculations and preparation of restrained-ensemble molecular dynamics simulations

The double electron-electron resonance (DEER) is a powerful structural biology technique to obtain distance information in the range of 18 to 80 å by measuring the dipolar coupling between two unpaired electron spins. The distance distributions obtained from the experiment provide valuable structural information about the protein in its native environment that can be exploited using restrained ensemble molecular dynamics (reMD) simulations. We present a new tool DEER Facilitator in CHARMM-GUI that consists of two modules Spin-Pair Distributor and reMD Prepper to setup simulations that utilize information from DEER experiments. Spin-Pair Distributor provides a web-based interface to calculate the spin-pair distance distribution of labeled sites in a protein using MD simulations. The calculated distribution can be used to guide the selection of the labeling sites in experiments as well as validate different protein structure models. reMD Prepper facilities the setup of reMD simulations using different types of spin labels in four different environments including vacuum, solution, micelle, and bilayer. The applications of these two modules are demonstrated with several test cases. Spin-Pair Distributor and reMD Prepper are available at http://www.charmm-gui.org/input/deer and http://www.charmm-gui.org/input/deerre. DEER Facilitator is expected to facilitate advanced biomolecular modeling and simulation, thereby leading to an improved understanding of the structure and dynamics of complex biomolecular systems based on experimental DEER data. © 2019 Wiley Periodicals, Inc.

Interfacial Aromatics Mediating Cation-π Interactions with Choline-Containing Lipids Can Contribute as Much to Peripheral Protein Affinity for Membranes as Aromatics Inserted below the Phosphates

Membrane-binding interfaces of peripheral proteins are restricted to a small part of their exposed surface, so the ability to engage in strong selective interactions with membrane lipids at various depths in the interface, both below and above the phosphates, is an advantage. Driven by their hydrophobicity, aromatic amino acids preferentially partition into membrane interfaces, often below the phosphates, yet enthalpically favorable interactions with the lipid headgroups, above the phosphate plane, are likely to further stabilize high interfacial positions. Using free-energy perturbation, we calculate the energetic cost of alanine substitution for 11 interfacial aromatic amino acids from 3 peripheral proteins. We show that the involvement in cation-π interactions with the headgroups (i) increases the ΔΔG_transfer as compared with insertion at the same depth without cation-π stabilization and (ii) can contribute at least as much as deeper insertion below the phosphates, highlighting the multiple roles of aromatics in peripheral membrane protein affinity.

CHARMM-GUI Glycan Modeler for modeling and simulation of carbohydrates and glycoconjugates

Characterizing glycans and glycoconjugates in the context of three-dimensional structures is important in understanding their biological roles and developing efficient therapeutic agents. Computational modeling and molecular simulation have become an essential tool complementary to experimental methods. Here, we present a computational tool, Glycan Modeler for in silico N-/O-glycosylation of the target protein and generation of carbohydrate-only systems. In our previous study, we developed Glycan Reader, a web-based tool for detecting carbohydrate molecules from a PDB structure and generation of simulation system and input files. As integrated into Glycan Reader in CHARMM-GUI, Glycan Modeler (Glycan Reader & Modeler) enables to generate the structures of glycans and glycoconjugates for given glycan sequences and glycosylation sites using PDB glycan template structures from Glycan Fragment Database (http://glycanstructure.org/fragment-db). Our benchmark tests demonstrate the universal applicability of Glycan Reader & Modeler to various glycan sequences and target proteins. We also investigated the structural properties of modeled glycan structures by running 2-μs molecular dynamics simulations of HIV envelope protein. The simulations show that the modeled glycan structures built by Glycan Reader & Modeler have the similar structural features compared to the ones solved by X-ray crystallography. We also describe the representative examples of glycoconjugate modeling with video demos to illustrate the practical applications of Glycan Reader & Modeler. Glycan Reader & Modeler is freely available at http://charmm-gui.org/input/glycan.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

CHARMM-GUI Nanodisc Builder for modeling and simulation of various nanodisc systems

Nanodiscs are discoidal protein-lipid complexes that have wide applications in membrane protein studies. Modeling and simulation of nanodiscs are challenging due to the absence of structures of many membrane scaffold proteins (MSPs) that wrap around the membrane bilayer. We have developed CHARMM-GUI Nanodisc Builder (http://www.charmm-gui.org/input/nanodisc) to facilitate the setup of nanodisc simulation systems by modeling the MSPs with defined size and known structural features. A total of 11 different nanodiscs with a diameter from 80 to 180 Å are made available in both the all-atom CHARMM and two coarse-grained (PACE and Martini) force fields. The usage of the Nanodisc Builder is demonstrated with various simulation systems. The structures and dynamics of proteins and lipids in these systems were analyzed, showing similar behaviors to those from previous all-atom and coarse-grained nanodisc simulations. We expect the Nanodisc Builder to be a convenient and reliable tool for modeling and simulation of nanodisc systems. © 2019 Wiley Periodicals, Inc.

Phosphatidic acid induces conformational changes in Sec18 protomers that prevent SNARE priming

Eukaryotic cell homeostasis requires transfer of cellular components among organelles and relies on membrane fusion catalyzed by SNARE proteins. Inactive SNARE bundles are reactivated by hexameric N-ethylmaleimide-sensitive factor, vesicle-fusing ATPase (Sec18/NSF)-driven disassembly that enables a new round of membrane fusion. We previously found that phosphatidic acid (PA) binds Sec18 and thereby sequesters it from SNAREs and that PA dephosphorylation dissociates Sec18 from the membrane, allowing it to engage SNARE complexes. We now report that PA also induces conformational changes in Sec18 protomers and that hexameric Sec18 cannot bind PA membranes. Molecular dynamics (MD) analyses revealed that the D1 and D2 domains of Sec18 contain PA-binding sites and that the residues needed for PA binding are masked in hexameric Sec18. Importantly, these simulations also disclosed that a major conformational change occurs in the linker region between the D1 and D2 domains, which is distinct from the conformational changes that occur in hexameric Sec18 during SNARE priming. Together, these findings indicate that PA regulates Sec18 function by altering its architecture and stabilizing membrane-bound Sec18 protomers.

Molecular mechanisms of pore formation and membrane disruption by the antimicrobial lantibiotic peptide Mutacin 1140

The emergence of antibiotic-resistance is a major concern to global human health and identification of novel antibiotics is critical to mitigate the threat.

Setting Up All-Atom Molecular Dynamics Simulations to Study the Interactions of Peripheral Membrane Proteins with Model Lipid Bilayers

All-atom molecular dynamics (MD) simulations enable the study of biological systems at atomic detail, complement the understanding gained from experiment, and can also motivate experimental techniques to further examine a given biological process. This method is based on statistical mechanics; it predicts the trajectory of atoms over time by solving Newton's Laws of motion taking into account all forces. Here, we describe the use of this methodology to study the interaction between peripheral membrane proteins and a lipid bilayer. Specifically, we provide step-by-step instructions to set up MD simulations to study the binding and interaction of the amphipathic helix of Osh4, a lipid transport protein, and Thanatin, an antimicrobial peptide (AMP), with model lipid bilayers using both fully detailed lipid tails and the highly mobile membrane-mimetic (HMMM) method to enhance conformational sampling.

Calcium-Induced Lipid Nanocluster Structures: Sculpturing of the Plasma Membrane

The plasma membrane of the cell is a complex, tightly regulated, heterogeneous environment shaped by proteins, lipids, and small molecules. Ca²⁺ ions are important cellular messengers, spatially separated from anionic lipids. After cell injury, disease, or apoptotic events, anionic lipids are externalized to the outer leaflet of the plasma membrane and encounter Ca²⁺, resulting in dramatic changes in the plasma membrane structure and initiation of signaling cascades. Despite the high chemical and biological significance, the structures of lipid-Ca²⁺ nanoclusters are still not known. Previously, we demonstrated by solid-state nuclear magnetic resonance (NMR) spectroscopy that upon binding to Ca²⁺, individual phosphatidylserine lipids populate two distinct yet-to-be-characterized structural environments. Here, we concurrently employ extensive all-atom molecular dynamics (MD) simulations with our accelerated membrane mimetic and detailed NMR measurements to identify lipid-Ca²⁺ nanocluster conformations. We find that major structural characteristics of these nanoclusters, including interlipid pair distances and chemical shifts, agree with observable NMR parameters. Simulations reveal that lipid-ion nanoclusters are shaped by two characteristic, long-lived lipid structures induced by divalent Ca²⁺. Using ab initio quantum mechanical calculations of chemical shifts on MD-captured lipid-ion complexes, we show that computationally observed conformations are validated by experimental NMR data. Both NMR measurements of diluted specifically labeled lipids and MD simulations reveal that the basic structural unit that reshapes the membrane is a Ca²⁺-coordinated phosphatidylserine tetramer. Our combined computational and experimental approach presented here can be applied to other complex systems in which charged membrane-active molecular agents leave structural signatures on lipids.

Influence of membrane lipid composition on the structure and activity of γ-secretase

γ-Secretase (GS) is a multi-subunit membrane-embedded aspartyl protease that cleaves more than 80 integral membrane proteins, including the amyloid precursor protein (APP) to produce the amyloid-β (Aβ) peptide. Oligomerization and aggregation of the 42-amino acid length Aβ isoform in the brain has been associated with the development and progression of Alzheimer's disease (AD). Based on recent experimental structural studies and using multiscale computational modeling approaches, the conformational states and protein-membrane interactions of the GS complex embedded in six homogeneous and six heterogeneous lipid bilayers were characterized. In order to identify potential lipid and cholesterol binding sites, GS regions with high lipid/cholesterol occupancy values were analyzed using atomistic and coarse-grained simulations. Long lipid residence times were observed to be correlated with a large number of hydrogen bonds between the charged headgroups and key GS amino acids. This observation provides a plausible explanation for the inhibition of GS by charged lipids observed in previous experimental studies. Computed lateral pressure profiles suggest that higher transmembrane pressures favor active state conformations of the catalytic subunit. A probable mechanism for the regulation of the local stress response in cholesterol-rich multicomponent lipid bilayers is identified. Finally, it is demonstrated that interactions between the nicastrin extracellular domain and lipid headgroups leads to a compact structural conformation of the GS complex. Overall, this study provides valuable insight into the effect of bilayer lipid composition on the GS structural ensemble and its function.

Sensitivity of peripheral membrane proteins to the membrane context: A case study of phosphatidylserine and the TIM proteins

There is a diverse class of peripheral membrane-binding proteins that specifically bind phosphatidylserine (PS), a lipid that signals apoptosis or cell fusion depending on the membrane context of its presentation. PS-receptors are specialized for particular PS-presenting pathways, indicating that they might be sensitive to the membrane context. In this review, we describe a combination of thermodynamic, structural, and computational techniques that can be used to investigate the mechanisms underlying this sensitivity. As an example, we focus on three PS-receptors of the T-cell Immunoglobulin and Mucin containing (TIM) protein family, which we have previously shown to differ in their sensitivity to PS surface density.

Molecular details of spontaneous insertion and interaction of HCV non-structure 3 protease protein domain with PIP2-containing membrane

Hepatitis C virus (HCV), known as the leading cause of liver cirrhosis, viral hepatitis, and hepatocellular carcinoma, has been affecting more than 150 million people globally. The HCV non-structure 3 (NS3) protease protein domain plays a key role in HCV replication and pathogenesis; and is currently a primary target for HCV antiviral therapy. Through unbiased molecular dynamics simulations which take advantage of the novel highly mobile membrane mimetic model, we constructed the membrane-bound state of the protein domain at the atomic level. Our results indicated that protease domain of HCV NS3 protein can spontaneously bind and penetrate to an endoplasmic reticulum complex membrane containing phosphatidylinositol 4,5-bisphosphate (PIP2). An amphipathic helix α₀ and loop S1 show their anchoring role to keep the protein on the membrane surface. Proper orientation of the protein domain at membrane surface was identified through measuring tilt angles of two specific vectors, wherein residue R161 plays a crucial role in its final orientation. Remarkably, PIP2 molecules were observed to bind to three main sites of the protease domain via specific electrostatic contacts and hydrogen bonds. PIP2-interaction determines the protein orientation at the membrane while both hydrophobic interplay and PIP2-interaction can stabilize the NS3 - membrane complex. Simulated results provide us with a detailed characterization of insertion, orientation and PIP2-interaction of NS3 protease domain at membrane environment, thus enhancing our understanding of structural functions and mechanism for the association of HCV non-structure 3 protein with respect to ER membranes.

Coupling X-Ray Reflectivity and In Silico Binding to Yield Dynamics of Membrane Recognition by Tim1

The dynamic nature of lipid membranes presents significant challenges with respect to understanding the molecular basis of protein/membrane interactions. Consequently, there is relatively little known about the structural mechanisms by which membrane-binding proteins might distinguish subtle variations in lipid membrane composition and/or structure. We have previously developed a multidisciplinary approach that combines molecular dynamics simulation with interfacial x-ray scattering experiments to produce an atomistic model for phosphatidylserine recognition by the immune receptor Tim4. However, this approach requires a previously determined protein crystal structure in a membrane-bound conformation. Tim1, a Tim4 homolog with distinct differences in both immunological function and sensitivity to membrane composition, was crystalized in a closed-loop conformation that is unlikely to support membrane binding. Here we have used a previously described highly mobile membrane mimetic membrane in combination with a conventional lipid bilayer model to generate a membrane-bound configuration of Tim1 in silico. This refined structure provided a significantly improved fit of experimental x-ray reflectivity data. Moreover, the coupling of the x-ray reflectivity analysis with both highly mobile membrane mimetic membranes and conventional lipid bilayer molecular dynamics simulations yielded a dynamic model of phosphatidylserine membrane recognition by Tim1 with atomic-level detail. In addition to providing, to our knowledge, new insights into the molecular mechanisms that distinguish the various Tim receptors, these results demonstrate that in silico membrane-binding simulations can remove the requirement that the existing crystal structure be in the membrane-bound conformation for effective x-ray reflectivity analysis. Consequently, this refined methodology has the potential for much broader applicability with respect to defining the atomistic details of membrane-binding proteins.

Lipid specificity of the membrane binding domain of coagulation factor X

Essentials Membrane-binding GLA domains of coagulation factors are essential for proper clot formation. Factor X (FX) is specific to phosphatidylserine (PS) lipids through unknown atomic-level interactions. Molecular dynamics simulations were used to develop the first membrane-bound model of FX-GLA. PS binding modes of FX-GLA were described, and potential PS-specific binding sites identified.

SUMMARY

Background Factor X (FX) binds to cell membranes in a highly phospholipid-dependent manner and, in complex with tissue factor and factor VIIa (FVIIa), initiates the clotting cascade. Experimental information concerning the membrane-bound structure of FX with atomic resolution has remained elusive because of the fluid nature of cellular membranes. FX is known to bind preferentially to phosphatidylserine (PS). Objectives To develop the first membrane-bound model of the FX-GLA domain to PS at atomic level, and to identify PS-specific binding sites of the FX-GLA domain. Methods Molecular dynamics (MD) simulations were performed to develop an atomic-level model for the FX-GLA domain bound to PS bilayers. We utilized a membrane representation with enhanced lipid mobility, termed the highly mobile membrane mimetic (HMMM), permitting spontaneous membrane binding and insertion by FX-GLA in multiple 100-ns simulations. In 14 independent simulations, FX-GLA bound spontaneously to the membrane. The resulting membrane-bound models were converted from HMMM to conventional membrane and simulated for an additional 100 ns. Results The final membrane-bound FX-GLA model allowed for detailed characterization of the orientation, insertion depth and lipid interactions of the domain, providing insight into the molecular basis of its PS specificity. All binding simulations converged to the same configuration despite differing initial orientations. Conclusions Analysis of interactions between residues in FX-GLA and lipid-charged groups allowed for potential PS-specific binding sites to be identified. This new structural and dynamic information provides an additional step towards a full understanding of the role of atomic-level lipid-protein interactions in regulating the critical and complex clotting cascade.

Spontaneous membrane insertion of a dengue virus NS2A peptide

Non-structural NS2A protein of Dengue virus is essential for viral replication but poorly characterized because of its high hydrophobicity. We have previously shown experimentally that NS2A possess a segment, peptide dens25, known to insert into membranes and interact specifically with negatively-charged phospholipids. To characterize its membrane interaction we have used two types of molecular dynamics membrane model systems, a highly mobile membrane mimetic (HMMM) and an endoplasmic reticulum (ER) membrane-like model. Using the HMMM system, we have been able of demonstrating the spontaneous binding of dens25 to the negatively-charged phospholipid 1,2-divaleryl-sn-glycero-3-phosphate containing membrane whereas no binding was observed for the membrane containing the zwitterionic one 1,2-divaleryl-sn-glycero-3-phosphocholine. Using the ER-like membrane model system, we demonstrate the spontaneous insertion of dens25 into the middle of the membrane, it maintained its three-dimensional structure and presented a nearly parallel orientation with respect to the membrane surface. Both charged and hydrophobic amino acids, presenting an interfacial/hydrophobic pattern characteristic of a membrane-proximal segment, are responsible for membrane binding and insertion. Dens25 might control protein/membrane interaction and be involved in membrane rearrangements critical for the viral cycle. These data should help us in the development of inhibitor molecules that target NS2A segments involved in membrane reorganisation.

CHARMM-GUI 10 years for biomolecular modeling and simulation

CHARMM-GUI, http://www.charmm-gui.org, is a web-based graphical user interface that prepares complex biomolecular systems for molecular simulations. CHARMM-GUI creates input files for a number of programs including CHARMM, NAMD, GROMACS, AMBER, GENESIS, LAMMPS, Desmond, OpenMM, and CHARMM/OpenMM. Since its original development in 2006, CHARMM-GUI has been widely adopted for various purposes and now contains a number of different modules designed to set up a broad range of simulations: (1) PDB Reader & Manipulator, Glycan Reader, and Ligand Reader & Modeler for reading and modifying molecules; (2) Quick MD Simulator, Membrane Builder, Nanodisc Builder, HMMM Builder, Monolayer Builder, Micelle Builder, and Hex Phase Builder for building all-atom simulation systems in various environments; (3) PACE CG Builder and Martini Maker for building coarse-grained simulation systems; (4) DEER Facilitator and MDFF/xMDFF Utilizer for experimentally guided simulations; (5) Implicit Solvent Modeler, PBEQ-Solver, and GCMC/BD Ion Simulator for implicit solvent related calculations; (6) Ligand Binder for ligand solvation and binding free energy simulations; and (7) Drude Prepper for preparation of simulations with the CHARMM Drude polarizable force field. Recently, new modules have been integrated into CHARMM-GUI, such as Glycolipid Modeler for generation of various glycolipid structures, and LPS Modeler for generation of lipopolysaccharide structures from various Gram-negative bacteria. These new features together with existing modules are expected to facilitate advanced molecular modeling and simulation thereby leading to an improved understanding of the structure and dynamics of complex biomolecular systems. Here, we briefly review these capabilities and discuss potential future directions in the CHARMM-GUI development project. © 2016 Wiley Periodicals, Inc.

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.

Extension of the Highly Mobile Membrane Mimetic to Transmembrane Systems through Customized in Silico Solvents

The mechanics of the protein-lipid interactions of transmembrane proteins are difficult to capture with conventional atomic molecular dynamics, due to the slow lateral diffusion of lipids restricting sampling to states near the initial membrane configuration. The highly mobile membrane mimetic (HMMM) model accelerates lipid dynamics by modeling the acyl tails nearest to the membrane center as a fluid organic solvent while maintaining an atomic description of the lipid headgroups and short acyl tails. The HMMM has been applied to many peripheral protein systems; however, the organic solvent used to date caused deformations in transmembrane proteins by intercalating into the protein and disrupting interactions between individual side chains. We ameliorate the effect of the solvent on transmembrane protein structure through the development of two new in silico Lennard-Jones solvents. The parameters for the new solvents were determined through an extensive parameter search in order to match the bulk properties of alkanes in a highly simplified model. Using these new solvents, we substantially improve the insertion free energy profiles of 10 protein side chain analogues across the entire bilayer. In addition, we reduce the intercalation of solvent into transmembrane systems, resulting in native-like transmembrane protein structures from five different topological classes within a HMMM bilayer. The parametrization of the solvents, in addition to their computed physical properties, is discussed. By combining high lipid lateral diffusion with intact transmembrane proteins, we foresee the developed solvents being useful to efficiently identify membrane composition inhomogeneities and lipid binding caused by the presence of membrane proteins.

Interaction of KRas4b with anionic membranes: A special role for PIP2

KRas4b is a small G-protein whose constitutively active oncogenic mutants are present in 90% of pancreatic cancers. Using fully post-translationally modified KRAS4b, we investigated the role of lipid identity in the recruitment of KRas4b to a membrane surface of defined composition. Application of a newly developed single frequency fluorescence anisotropy decay experiment to this system revealed that KRas4b has a significant binding preference for Nanodisc bilayers containing PIP₂. We conducted molecular dynamics simulations to look for an origin of this specificity. In the case of membranes containing PIP₂ the protein formed long-lived salt bridges with PIP₂ head groups but not the monovalent DMPS, explaining the experimentally observed lipid specificity. Additionally, we report that PIP₂ forms key contacts with Helix-4 on the catalytic domain of KRas4b that orient the protein in a manner expected to facilitate association with upstream and downstream signaling partners.

Brain grants permission of access to Zika virus but denies entry to drugs: a molecular modeling perspective to infiltrate the boundary

Thein silicodesign of targeted Zika virus inhibitors.

Effects of Spin-Labels on Membrane Burial Depth of MARCKS-ED Residues

Site-directed spin-labeling electron paramagnetic resonance spectroscopy is a useful tool to obtain information about the environment of specific residues. One of its applications is to investigate membrane protein topology based on the accessibility of the spin label, with the assumption that the position of the spin label in the membrane is close to that of the native residue. This assumption is valid in proteins with well-ordered structures, but could be problematic in small peptides because the labeling may cause a perturbation that is large enough to change local interactions between the peptide and the membrane. To quantitatively characterize such effects, we have simulated the association of a 25-amino-acid peptide, MARCKS-ED, to membranes with and without spin labels. Our simulations show that the depths of spin labels are ∼6-17 Å deeper than the unlabeled charged and polar residues in the wild-type. When the hydrophobic residue Phe is labeled, however, the spin-label depth is close to that of the native residue as well as the experimental value. Our study suggests that one should be cautious in interpretation of spin label data when charged and polar residues in small peptides are labeled.

Challenges in structural approaches to cell modeling

Computational modeling is essential for structural characterization of biomolecular mechanisms across the broad spectrum of scales. Adequate understanding of biomolecular mechanisms inherently involves our ability to model them. Structural modeling of individual biomolecules and their interactions has been rapidly progressing. However, in terms of the broader picture, the focus is shifting toward larger systems, up to the level of a cell. Such modeling involves a more dynamic and realistic representation of the interactomes in vivo, in a crowded cellular environment, as well as membranes and membrane proteins, and other cellular components. Structural modeling of a cell complements computational approaches to cellular mechanisms based on differential equations, graph models, and other techniques to model biological networks, imaging data, etc. Structural modeling along with other computational and experimental approaches will provide a fundamental understanding of life at the molecular level and lead to important applications to biology and medicine. A cross section of diverse approaches presented in this review illustrates the developing shift from the structural modeling of individual molecules to that of cell biology. Studies in several related areas are covered: biological networks; automated construction of three-dimensional cell models using experimental data; modeling of protein complexes; prediction of non-specific and transient protein interactions; thermodynamic and kinetic effects of crowding; cellular membrane modeling; and modeling of chromosomes. The review presents an expert opinion on the current state-of-the-art in these various aspects of structural modeling in cellular biology, and the prospects of future developments in this emerging field.

Microscopic Characterization of Membrane Transporter Function by In Silico Modeling and Simulation

Membrane transporters mediate one of the most fundamental processes in biology. They are the main gatekeepers controlling active traffic of materials in a highly selective and regulated manner between different cellular compartments demarcated by biological membranes. At the heart of the mechanism of membrane transporters lie protein conformational changes of diverse forms and magnitudes, which closely mediate critical aspects of the transport process, most importantly the coordinated motions of remotely located gating elements and their tight coupling to chemical processes such as binding, unbinding and translocation of transported substrate and cotransported ions, ATP binding and hydrolysis, and other molecular events fueling uphill transport of the cargo. An increasing number of functional studies have established the active participation of lipids and other components of biological membranes in the function of transporters and other membrane proteins, often acting as major signaling and regulating elements. Understanding the mechanistic details of these molecular processes require methods that offer high spatial and temporal resolutions. Computational modeling and simulations technologies empowered by advanced sampling and free energy calculations have reached a sufficiently mature state to become an indispensable component of mechanistic studies of membrane transporters in their natural environment of the membrane. In this article, we provide an overview of a number of major computational protocols and techniques commonly used in membrane transporter modeling and simulation studies. The article also includes practical hints on effective use of these methods, critical perspectives on their strengths and weak points, and examples of their successful applications to membrane transporters, selected from the research performed in our own laboratory.

Atomic-level description of protein-lipid interactions using an accelerated membrane model

Peripheral membrane proteins are structurally diverse proteins that are involved in fundamental cellular processes. Their activity of these proteins is frequently modulated through their interaction with cellular membranes, and as a result techniques to study the interfacial interaction between peripheral proteins and the membrane are in high demand. Due to the fluid nature of the membrane and the reversibility of protein-membrane interactions, the experimental study of these systems remains a challenging task. Molecular dynamics simulations offer a suitable approach to study protein-lipid interactions; however, the slow dynamics of the lipids often prevents sufficient sampling of specific membrane-protein interactions in atomistic simulations. To increase lipid dynamics while preserving the atomistic detail of protein-lipid interactions, in the highly mobile membrane-mimetic (HMMM) model the membrane core is replaced by an organic solvent, while short-tailed lipids provide a nearly complete representation of natural lipids at the organic solvent/water interface. Here, we present a brief introduction and a summary of recent applications of the HMMM to study different membrane proteins, complementing the experimental characterization of the presented systems, and we offer a perspective of future applications of the HMMM to study other classes of membrane proteins. This article is part of a Special Issue entitled: Membrane proteins edited by J.C. Gumbart and Sergei Noskov.

Ionic Hydrogen Bonds and Lipid Packing Defects Determine the Binding Orientation and Insertion Depth of RecA on Multicomponent Lipid Bilayers

We describe a computational and experimental approach for probing the binding properties of the RecA protein at the surface of anionic membranes. Fluorescence measurements indicate that RecA behaves differently when bound to phosphatidylglycerol (PG)- and cardiolipin (CL)-containing liposomes. We use a multistage computational protocol that integrates an implicit membrane/solvent model, the highly mobile mimetic membrane model, and the full atomistic membrane model to study how different anionic lipids perturb RecA binding to the membrane. With anionic lipids studied here, the binding interface involves three key regions: the N-terminal helix, the DNA binding loop L2, and the M-M7 region. The nature of binding involves both electrostatic interactions between cationic protein residues and lipid polar/charged groups and insertion of hydrophobic residues. The L2 loop contributes more to membrane insertion than the N-terminal helix. More subtle aspects of RecA-membrane interaction are influenced by specific properties of anionic lipids. Ionic hydrogen bonds between the carboxylate group in phosphatidylserine and several lysine residues in the C-terminal region of RecA stabilize the parallel (∥) binding orientation, which is not locally stable on PG- and CL-containing membranes despite similarity in the overall charge density. Lipid packing defects, which are more prevalent in the presence of conical lipids, are observed to enhance the insertion depth of hydrophobic motifs. The computational finding that RecA binds in a similar orientation to PG- and CL-containing membranes is consistent with the fact that PG alone is sufficient to induce RecA polar localization, although CL might be more effective because of its tighter binding to RecA. The different fluorescence behaviors of RecA upon binding to PG- and CL-containing liposomes is likely due to the different structures and flexibility of the C-terminal region of RecA when it binds to different anionic phospholipids.

Efficient preparation and analysis of membrane and membrane protein systems

Molecular dynamics (MD) simulations have become a highly important technique to consider lipid membrane systems, and quite often they provide considerable added value to laboratory experiments. Rapid development of both software and hardware has enabled the increase of time and size scales reachable by MD simulations to match those attainable by several accurate experimental techniques. However, until recently, the quality and maturity of software tools available for building membrane models for simulations as well as analyzing the results of these simulations have seriously lagged behind. Here, we discuss the recent developments of such tools from the end-users' point of view. In particular, we review the software that can be employed to build lipid bilayers and other related structures with or without embedded membrane proteins to be employed in MD simulations. Additionally, we provide a brief critical insight into force fields and MD packages commonly used for membrane and membrane protein simulations. Finally, we list analysis tools that can be used to study the properties of membrane and membrane protein systems. In all these points we comment on the respective compatibility of the covered tools. We also share our opinion on the current state of the available software. We briefly discuss the most commonly employed tools and platforms on which new software can be built. We conclude the review by providing a few ideas and guidelines on how the development of tools can be further boosted to catch up with the rapid pace at which the field of membrane simulation progresses. This includes improving the compatibility between software tools and promoting the openness of the codes on which these applications rely. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field

Proper treatment of nonbonded interactions is essential for the accuracy of molecular dynamics (MD) simulations, especially in studies of lipid bilayers. The use of the CHARMM36 force field (C36 FF) in different MD simulation programs can result in disagreements with published simulations performed with CHARMM due to differences in the protocols used to treat the long-range and 1-4 nonbonded interactions. In this study, we systematically test the use of the C36 lipid FF in NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM. A wide range of Lennard-Jones (LJ) cutoff schemes and integrator algorithms were tested to find the optimal simulation protocol to best match bilayer properties of six lipids with varying acyl chain saturation and head groups. MD simulations of a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer were used to obtain the optimal protocol for each program. MD simulations with all programs were found to reasonably match the DPPC bilayer properties (surface area per lipid, chain order parameters, and area compressibility modulus) obtained using the standard protocol used in CHARMM as well as from experiments. The optimal simulation protocol was then applied to the other five lipid simulations and resulted in excellent agreement between results from most simulation programs as well as with experimental data. AMBER compared least favorably with the expected membrane properties, which appears to be due to its use of the hard-truncation in the LJ potential versus a force-based switching function used to smooth the LJ potential as it approaches the cutoff distance. The optimal simulation protocol for each program has been implemented in CHARMM-GUI. This protocol is expected to be applicable to the remainder of the additive C36 FF including the proteins, nucleic acids, carbohydrates, and small molecules.

CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field

Proper treatment of nonbonded interactions is essential for the accuracy of molecular dynamics (MD) simulations, especially in studies of lipid bilayers. The use of the CHARMM36 force field (C36 FF) in different MD simulation programs can result in disagreements with published simulations performed with CHARMM due to differences in the protocols used to treat the long-range and 1-4 nonbonded interactions. In this study, we systematically test the use of the C36 lipid FF in NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM. A wide range of Lennard-Jones (LJ) cutoff schemes and integrator algorithms were tested to find the optimal simulation protocol to best match bilayer properties of six lipids with varying acyl chain saturation and head groups. MD simulations of a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer were used to obtain the optimal protocol for each program. MD simulations with all programs were found to reasonably match the DPPC bilayer properties (surface area per lipid, chain order parameters, and area compressibility modulus) obtained using the standard protocol used in CHARMM as well as from experiments. The optimal simulation protocol was then applied to the other five lipid simulations and resulted in excellent agreement between results from most simulation programs as well as with experimental data. AMBER compared least favorably with the expected membrane properties, which appears to be due to its use of the hard-truncation in the LJ potential versus a force-based switching function used to smooth the LJ potential as it approaches the cutoff distance. The optimal simulation protocol for each program has been implemented in CHARMM-GUI. This protocol is expected to be applicable to the remainder of the additive C36 FF including the proteins, nucleic acids, carbohydrates, and small molecules.