Computer simulations of membrane proteins

Computing the Structural Dynamics of RVFV L Protein Domain in Aqueous Glycerol Solutions

Many biological and biotechnological processes are controlled by protein–protein and protein–solvent interactions. In order to understand, predict, and optimize such processes, it is important to understand how solvents affect protein structure during protein–solvent interactions. In this study, all-atom molecular dynamics are used to investigate the structural dynamics and energetic properties of a C-terminal domain of the Rift Valley Fever Virus L protein solvated in glycerol and aqueous glycerol solutions in different concentrations by molecular weight. The Generalized Amber Force Field is modified by including restrained electrostatic potential atomic charges for the glycerol molecules. The peptide is considered in detail by monitoring properties like the root-mean-squared deviation, root-mean-squared fluctuation, radius of gyration, hydrodynamic radius, end-to-end distance, solvent-accessible surface area, intra-potential energy, and solvent–peptide interaction energies for hundreds of nanoseconds. Secondary structure analysis is also performed to examine the extent of conformational drift for the individual helices and sheets. We predict that the peptide helices and sheets are maintained only when the modeling strategy considers the solvent with lower glycerol concentration. We also find that the solvent-peptide becomes more cohesive with decreasing glycerol concentrations. The density and radial distribution function of glycerol solvent calculated when modeled with the modified atomic charges show a very good agreement with experimental results and other simulations at 298.15K.

Theoretical study of the interactions between peptide tyrosine tyrosine [PYY (1-36)], a newly identified modulator in type 2 diabetes pathophysiology, with receptors NPY1R and NPY4R

PURPOSE

Diabetes mellitus is a common condition in the clinically obese. Bariatric surgery is one of the ways to put type 2 diabetes in remission. Recent findings propose the appetite-regulator peptide tyrosine tyrosine (PYY) as a therapeutic option for patients with type 2 diabetes. This novel gut hormone restores impaired insulin and glucagon secretion in pancreatic islets and is implicated in type 2 diabetes reversal after bariatric surgery. The current study elucidates the interactions between PYY and the NPY1R and NPY4R receptors using computational methods.

METHODS

Protein structure prediction, molecular docking simulation, and molecular dynamics (MD) simulation were performed to elucidate the interactions of PYY with NPY1R and NPY4R.

RESULTS

The predicted binding models of PYY-NPY receptors are in agreement with those described in the literature, although different interaction partners are presented for the C-terminal tail of PYY. Non-polar interactions are predicted to drive the formation of the protein complex. The calculated binding energies show that PYY has higher affinity for NPY4R (ΔG_GBSA = -65.08 and ΔG_PBSA = -87.62 kcal/mol) than for NPY1R (ΔG_GBSA = -23.11 and ΔG_PBSA = -50.56 kcal/mol).

CONCLUSIONS

Based on the constructed models, the binding conformations obtained from docking and MD simulation for both the PYY-NPY1R and PYY-NPY4R complexes provide a detailed map of possible interactions. The calculated binding energies show a higher affinity of PYY for NPY4R. These findings may help to understand the mechanisms behind the improvement of diabetes following bariatric surgery.

Influence of Lipid Composition on the Insertion Process of Glyphosate into Membranes: A Thermodynamic Study

In this work, molecular dynamics simulations were applied to investigate the influence of lipid composition of the model membrane on the insertion of glyphosate (in its charged state, GLYP2-). The profiles of free energy, entropy and enthalpy were obtained through umbrella sampling calculations, for lipid bilayers composed by only 1,2-dipalmitoyl-sn-glycerol-3-phosphocholine (DPPC), only 1,2-dipalmitoyl-sn-glycerol-3-phosphoserine (DPPS) or a symmetric binary mixture of DPPC and DPPS. In general, the location, the values of minima and maxima of the free energy, and the trend of free energy profiles are influenced by the lipid composition of the lipid bilayer. The driving force in the glyphosate insertion process depends on the lipid composition of the membrane model. If the lipid bilayer is composed solely of DPPS or DPPC, GLYP2- insertion is driven by a favorable enthalpic change. However, if the membrane is composed of a mixture of both lipids, this process is driven by a favorable entropic change. In the lipid bilayer containing DPPS, the glyphosate was found to penetrate hydrated and coordinated with Na+ ions, in contrast to the pure zwitterionic lipid bilayer which penetrated only hydrated. This effect is independent of the concentration of sodium ions present in the bulk solution.

Molecular Dynamics Simulations in Drug Discovery and Pharmaceutical Development

Molecular dynamics (MD) simulations have become increasingly useful in the modern drug development process. In this review, we give a broad overview of the current application possibilities of MD in drug discovery and pharmaceutical development. Starting from the target validation step of the drug development process, we give several examples of how MD studies can give important insights into the dynamics and function of identified drug targets such as sirtuins, RAS proteins, or intrinsically disordered proteins. The role of MD in antibody design is also reviewed. In the lead discovery and lead optimization phases, MD facilitates the evaluation of the binding energetics and kinetics of the ligand-receptor interactions, therefore guiding the choice of the best candidate molecules for further development. The importance of considering the biological lipid bilayer environment in the MD simulations of membrane proteins is also discussed, using G-protein coupled receptors and ion channels as well as the drug-metabolizing cytochrome P450 enzymes as relevant examples. Lastly, we discuss the emerging role of MD simulations in facilitating the pharmaceutical formulation development of drugs and candidate drugs. Specifically, we look at how MD can be used in studying the crystalline and amorphous solids, the stability of amorphous drug or drug-polymer formulations, and drug solubility. Moreover, since nanoparticle drug formulations are of great interest in the field of drug delivery research, different applications of nano-particle simulations are also briefly summarized using multiple recent studies as examples. In the future, the role of MD simulations in facilitating the drug development process is likely to grow substantially with the increasing computer power and advancements in the development of force fields and enhanced MD methodologies.

Whole-Cell Models and Simulations in Molecular Detail

Comprehensive data about the composition and structure of cellular components have enabled the construction of quantitative whole-cell models. While kinetic network-type models have been established, it is also becoming possible to build physical, molecular-level models of cellular environments. This review outlines challenges in constructing and simulating such models and discusses near- and long-term opportunities for developing physical whole-cell models that can connect molecular structure with biological function.

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.

The Control of Colloidal Grain Boundaries through Evaporative Vertical Self-Assembly

Self-assembly continuously gains attention as an excellent method to create novel nanoscale structures with a wide range of applications in photonics, optoelectronics, biomedical engineering, and heat transfer applications. However, self-assembly is governed by a diversity of complex interparticle forces that cause fabricating defectless large scale (>1 cm) colloidal crystals, or opals, to be a daunting challenge. Despite numerous efforts to find an optimal method that offers the perfect colloidal crystal by minimizing defects, it has been difficult to provide physical interpretations that govern the development of defects such as grain boundaries. This study reports the control over grain domains and intentional defect characteristics that develop during evaporative vertical deposition. The degree of particle crystallinity and evaporation conditions is shown to govern the grain domain characteristics, such as shapes and sizes. In particular, the grains fabricated with 300 and 600 nm sphere diameters can be tuned into single-column structures exceeding ≈1 mm by elevating heating temperature up to 93 °C. The understanding of self-assembly physics presented in this work will enable the fabrication of novel self-assembled structures with high periodicity and offer fundamental groundworks for developing large-scale crack-free structures.

Multiscale Modeling and Simulation Approaches to Lipid-Protein Interactions

Lipid membranes play a crucial role in living systems by compartmentalizing biological processes and forming a barrier between these processes and the environment. Naturally, a large apparatus of biomolecules is responsible for construction, maintenance, transport, and degradation of these lipid barriers. Additional classes of biomolecules are tasked with transport of specific substances or transduction of signals from the environment across lipid membranes. In this article, we intend to describe a set of techniques that enable one to build accurate models of lipid systems and their associated proteins, and to simulate their dynamics over a variety of time and length scales. We discuss the methods and challenges that allow us to derive structural, mechanistic, and thermodynamic information from these models. We also show how these models have recently been applied in research to study some of the most complex lipid-protein systems to date, including bacterial and viral envelopes, neuronal membranes, and mammalian signaling systems.

Dynamic Interactions between Lipid-Tethered DNA and Phospholipid Membranes

Lipid-anchored DNA can attach functional cargo to bilayer membranes in DNA nanotechnology, synthetic biology, and cell biology research. To optimize DNA anchoring, an understanding of DNA-membrane interactions in terms of binding strength, extent, and structural dynamics is required. Here we use experiments and molecular dynamics (MD) simulations to determine how the membrane binding of cholesterol-modified DNA depends on electrostatic and steric factors involving the lipid headgroup charge, duplexed or single-stranded DNA, and the buffer composition. The experiments distinguish between free and membrane vesicle-bound DNA and thereby reveal the surface density of anchored DNA and its binding affinity, something which had previously not been known. The K_d values range from 8.5 ± 4.9 to 466 ± 134 μM whereby negatively charged headgroups led to weak binding due to the electrostatic repulsion with respect to the negatively charged DNA. Atomistic MD simulations explain the findings and elucidate the dynamic nature of anchored DNA such as the mushroom-like conformation of single-stranded DNA hovering over the bilayer surface in contrast to a straight-up conformation of double-stranded DNA. The biophysical insight into the binding strength to membranes as well as the molecular accessibility of DNA for hybridization to molecular cargo is expected to facilitate the creation of biomimetic DNA versions of natural membrane nanopores and cytoskeletons for research and nanobiotechnology.

Molecular dynamics simulations of glyphosate in a DPPC lipid bilayer

Extensive molecular dynamics simulations have been performed to study the effect of glyphosate (in their neutral and charged forms, GLYP and GLYP2-, respectively) on fully hydrated DiPalmitoylPhosphatidylCholine (DPPC) lipid bilayer. First, we calculated the free energy profile (using the Umbrella Sampling technique) for both states of charge of glyphosate. The minimum value for the free energy for GLYP is ∼-60 kJ mol-1 located at z = ±1.7 nm (from the lipid bilayer center), and there is almost no maximum at the center of the lipid bilayer. By contrast, the minimum for GLYP2- is ∼-35 kJ mol-1 located at z = ± 1.4 nm (from the lipid bilayer center), and the maximum reaches ∼35 kJ mol-1 at the center of the lipid bilayer. Then, different lipid bilayer properties were analyzed for different glyphosate:lipid (G:L) ratios. The mean area per lipid was slightly affected, increasing only 5% (in the presence of glyphosate at high concentrations), which is in agreement with the slight decrease in deuterium order parameters. As for the thickness of the bilayer, it is observed that the state of charge produces opposite effects. On one hand, the neutral state produces an increase in the thickness of the lipid bilayer; on the other, the charged form produces a decrease in the thickness, which not depend linearly on the G:L ratios, either. The orientation of the DPPC head groups is practically unaffected throughout the range of the G:L ratios studied. Finally, the mobility of the lipids of the bilayer is strongly affected by the presence of glyphosate, considerably increasing its lateral diffusion coefficient noteworthy (one order of magnitude), with increasing G:L ratio.

Folding a viral peptide in different membrane environments: pathway and sampling analyses

Flock House virus (FHV) is a well-characterized model system to study infection mechanisms in non-enveloped viruses. A key stage of the infection cycle is the disruption of the endosomal membrane by a component of the FHV capsid, the membrane active γ peptide. In this study, we perform all-atom molecular dynamics simulations of the 21 N-terminal residues of the γ peptide interacting with membranes of differing compositions. We carry out umbrella sampling calculations to study the folding of the peptide to a helical state in homogenous and heterogeneous membranes consisting of neutral and anionic lipids. From the trajectory data, we evaluate folding energetics and dissect the mechanism of folding in the different membrane environments. We conclude the study by analyzing the extent of configurational sampling by performing time-lagged independent component analysis.

Developing force fields when experimental data is sparse: AMBER/GAFF-compatible parameters for inorganic and alkyl oxoanions

We present a set of Lennard-Jones parameters for classical, all-atom models of acetate and various alkylated and non-alkylated forms of sulfate, sulfonate and phosphate ions, optimized to reproduce their interactions with water and with the physiologically relevant sodium, ammonium and methylammonium cations. The parameters are internally consistent and are fully compatible with the Generalized Amber Force Field (GAFF), the AMBER force field for proteins, the accompanying TIP3P water model and the sodium model of Joung and Cheatham. The parameters were developed primarily relying on experimental information - hydration free energies and solution activity derivatives at 0.5 m concentration - with ab initio, gas phase calculations being used for the cases where experimental information is missing. The ab initio parameterization scheme presented here is distinct from other approaches because it explicitly connects gas phase binding energies to intermolecular interactions in solution. We demonstrate that the original GAFF/AMBER parameters often overestimate anion-cation interactions, leading to an excessive number of contact ion pairs in solutions of carboxylate ions, and to aggregation in solutions of divalent ions. GAFF/AMBER parameters lead to excessive numbers of salt bridges in proteins and of contact ion pairs between sodium and acidic protein groups, issues that are resolved by using the optimized parameters presented here.

H. R, Thamleena A. H, Joseph A, Ben A, V. U. K. Roles of different amino-acid residues towards binding and selective transport of K+ through KcsA K+-ion channel

Each amino acid in the selectivity filter plays a distinct role towards binding and transport of K⁺ ion through KcsA.

Molecular modelling studies of kdr mutations in voltage gated sodium channel revealed significant conformational variations contributing to insecticide resistance

Voltage gated sodium channels (VGSC) of mosquito vectors are the primary targets of dichlorodiphenyltrichloroethane (DDT) and other synthetic pyrethroids used in public health programmes. The knockdown resistant (kdr) mutations in VGSC are associated with the insecticide resistance especially in Anophelines. The present study is aimed to emphasize and demarcate the impact of three kdr-mutations such as L1014S, L1014F and L1014H on insecticide resistance. The membrane model of sodium transport domain of VGSC (STD-VGSC) was constructed using de novo approach based on domain and trans-membrane predictions. The comparative molecular modelling studies of wild type and mutant models of STD-VGSC revealed that L1014F mutant was observed to be near native to the wild type model in all the respects, but, L1014S and L1014H mutations showed drastic variations in the energy levels, root mean square fluctuations (RMSF) that resulted in conformational variations. The predicted binding sites also showed variable cavity volumes and RMSF in L1014S and L1014H mutants. Further, DDT also found be bound in near native manner to wild type in L1014F mutant and with variable orientation and affinities in L1014S and L1014H mutants. The variations and fluctuations observed in mutant structures explained that each mutation has its specific impact on the conformation of VGSC and its binding with DDT. The study provides new insights into the structure-function-correlations of mutant STD-VGSC structures and demonstrates the role and effects of kdr mutations on insecticide resistance in mosquito vectors.

Molecular Modelling, Dynamics, and Docking of Membrane Proteins

Computational tools and techniques are now most popular and promising to progress the research at rapid rate. Molecular modelling studies contribute their maximum role in wide variety of disciplines especially in proteomics and drug discovery strategies. Molecular dynamics and molecular docking algorithms are now became an essential part in daily research activities of every laboratory throughout the world. These strategies are now well established and standardised to study any specific protein of interest and drug molecule. But still there exist considerable drawbacks in a special concern with membrane proteins as the presently available tools and methods cannot be applied directly to them. Modelling, dynamics and docking studies of membrane proteins need a special care and attention as several challenges are to be crossed with an intensive care to produce a reliable result. This chapter is aimed to discuss such challenges and solutions to handle membrane proteins.

Colored spectrum characteristics of thermal noise on the molecular scale

Thermal noise is of fundamental importance to many processes. Traditionally, thermal noise has been treated as white noise on the macroscopic scale. Using molecular dynamics simulations and power spectrum analysis, we show that the thermal noise of solute molecules in water is non-white on the molecular scale, which is in contrast to the conventional theory. In the frequency domain from 2 × 10¹¹ Hz to 10¹³ Hz, the power spectrum of thermal noise for polar solute molecules resembles the spectrum of 1/f noise. The power spectrum of thermal noise for non-polar solute molecules deviates only slightly from the spectrum of white noise. The key to this phenomenon is the existence of hydrogen bonds between polar solute molecules and solvent water molecules. Furthermore, for polar solute molecules, the degree of power spectrum deviation from that of white noise is associated with the average lifetime of the hydrogen bonds between the solute and the solvent molecules.

All-Atom Continuous Constant pH Molecular Dynamics With Particle Mesh Ewald and Titratable Water

Development of a pH stat to properly control solution pH in biomolecular simulations has been a long-standing goal in the community. Toward this goal recent years have witnessed the emergence of the so-called constant pH molecular dynamics methods. However, the accuracy and generality of these methods have been hampered by the use of implicit-solvent models or truncation-based electrostatic schemes. Here we report the implementation of the particle mesh Ewald (PME) scheme into the all-atom continuous constant pH molecular dynamics (CpHMD) method, enabling CpHMD to be performed with a standard MD engine at a fractional added computational cost. We demonstrate the performance using pH replica-exchange CpHMD simulations with titratable water for a stringent test set of proteins, HP36, BBL, HEWL, and SNase. With the sampling time of 10 ns per replica, most pK_a's are converged, yielding the average absolute and root-mean-square deviations of 0.61 and 0.77, respectively, from experiment. Linear regression of the calculated vs experimental pK_a shifts gives a correlation coefficient of 0.79, a slope of 1, and an intercept near 0. Analysis reveals inadequate sampling of structure relaxation accompanying a protonation-state switch as a major source of the remaining errors, which are reduced as simulation prolongs. These data suggest PME-based CpHMD can be used as a general tool for pH-controlled simulations of macromolecular systems in various environments, enabling atomic insights into pH-dependent phenomena involving not only soluble proteins but also transmembrane proteins, nucleic acids, surfactants, and polysaccharides.

Conserved methionine dictates substrate preference in Nramp-family divalent metal transporters

Natural resistance-associated macrophage protein (Nramp) family transporters catalyze uptake of essential divalent transition metals like iron and manganese. To discriminate against abundant competitors, the Nramp metal-binding site should favor softer transition metals, which interact either covalently or ionically with coordinating molecules, over hard calcium and magnesium, which interact mainly ionically. The metal-binding site contains an unusual, but conserved, methionine, and its sulfur coordinates transition metal substrates, suggesting a vital role in their transport. Using a bacterial Nramp model system, we show that, surprisingly, this conserved methionine is dispensable for transport of the physiological manganese substrate and similar divalents iron and cobalt, with several small amino acid replacements still enabling robust uptake. Moreover, the methionine sulfur's presence makes the toxic metal cadmium a preferred substrate. However, a methionine-to-alanine substitution enables transport of calcium and magnesium. Thus, the putative evolutionary pressure to maintain the Nramp metal-binding methionine likely exists because it-more effectively than any other amino acid-increases selectivity for low-abundance transition metal transport in the presence of high-abundance divalents like calcium and magnesium.

Molecular Simulations of Gram-Negative Bacterial Membranes: A Vignette of Some Recent Successes

In the following review we use recent examples from the literature to discuss progress in the area of atomistic and coarse-grained molecular dynamics simulations of selected bacterial membranes and proteins, with a particular focus on Gram-negative bacteria. As structural biology continues to provide increasingly high-resolution data on the proteins that reside within these membranes, simulations have an important role to play in linking these data with the dynamical behavior and function of these proteins. In particular, in the last few years there has been significant progress in addressing the issue of biochemical complexity of bacterial membranes such that the heterogeneity of the lipid and protein components of these membranes are now being incorporated into molecular-level models. Thus, in future we can look forward to complementary data from structural biology and molecular simulations combining to provide key details of structure-dynamics-function relationships in bacterial membranes.

Glu106 targeted inhibitors of ORAI1 as potential Ca2+ release-activated Ca2+ (CRAC) channel blockers - molecular modeling and docking studies

Calcium release-activated calcium modulator 1(ORAI1) is an integral component of the calcium release-activated calcium channel (CRAC) channel complex and plays a central role in regulating Ca^2 +concentrations in T-lymphocytes. It is critical for many physiological processes, including cell-proliferation, cytokine production and activation of the immune system. Loss of ORAI1 function is linked with rheumatoid arthritis (RA) and hence pharmacological blockers of ORAI1 could be potential therapeutic agents for the treatment of RA. In this study, we have used a high-throughput screening approach to inhibit the binding of Ca²⁺toward ORAI1 and the interactions are verified through induced fit docking. The results hint that these compounds act by possibly binding with, and thereby blocking Ca²⁺-binding with ORAI1 (E106). The molecular dynamics (MD) simulations shows strong support toward the hit compounds by showing the ligand potency throughout the simulation timescale of 30 ns. We have thus identified a novel class of highly stable, potential lead compounds that directly bind with the selectivity filter region E106 and block Ca²⁺ binding on ORAI1. This resulting alteration in the pore geometry of ORAI1 due to the strong blocking mechanism of lead compounds will greatly diminish its function and the downstream activities that result from the same including decreased production of cytokines in autoimmune disorders. This study may lay the foundation for finding novel lead compounds for clinical trials that could positively modulate the course of autoimmune disorders with ORAI1 as its specific target.

Living on the edge: Simulations of bacterial outer-membrane proteins

Gram-negative bacteria are distinguished in part by a second, outer membrane surrounding them. This membrane is distinct from others, possessing an outer leaflet composed not of typical phospholipids but rather large, highly charged molecules known as lipopolysaccharides. Therefore, modeling the structure and dynamics of proteins embedded in the outer membrane requires careful consideration of their native environment. In this review, we examine how simulations of such outer-membrane proteins have evolved over the last two decades, culminating most recently in detailed, highly accurate atomistic models of the outer membrane. We also draw attention to how the simulations have coupled with experiments to produce novel insights unattainable through a single approach. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Protein-fluctuation-induced water-pore formation in ion channel voltage-sensor translocation across a lipid bilayer membrane

We have applied a combined fluorescence microscopy and single-ion-channel electric current recording approach, correlating with molecular dynamics (MD) simulations, to study the mechanism of voltage-sensor domain translocation across a lipid bilayer. We use the colicin Ia ion channel as a model system, and our experimental and simulation results show the following: (1) The open-close activity of an activated colicin Ia is not necessarily sensitive to the amplitude of the applied cross-membrane voltage when the cross-membrane voltage is around the resting potential of excitable membranes; and (2) there is a significant probability that the activation of colicin Ia occurs by forming a transient and fluctuating water pore of ∼15 Å diameter in the lipid bilayer membrane. The location of the water-pore formation is nonrandom and highly specific, right at the insertion site of colicin Ia charged residues in the lipid bilayer membrane, and the formation is intrinsically associated with the polypeptide conformational fluctuations and solvation dynamics. Our results suggest an interesting mechanistic pathway for voltage-sensitive ion channel activation, and specifically for translocation of charged polypeptide chains across the lipid membrane under a transmembrane electric field: the charged polypeptide domain facilitates the formation of hydrophilic water pore in the membrane and diffuses through the hydrophilic pathway across the membrane; i.e., the charged polypeptide chain can cross a lipid membrane without entering into the hydrophobic core of the lipid membrane but entirely through the aqueous and hydrophilic environment to achieve a cross-membrane translocation. This mechanism sheds light on the intensive and fundamental debate on how a hydrophilic and charged peptide domain diffuses across the biologically inaccessible high-energy barrier of the hydrophobic core of a lipid bilayer: The peptide domain does not need to cross the hydrophobic core to move across a lipid bilayer.

Fluctuating hydrodynamics of multicomponent membranes with embedded proteins

A simulation method for the dynamics of inhomogeneous lipid bilayer membranes is presented. The membrane is treated using stochastic Saffman-Delbrück hydrodynamics, coupled to a phase-field description of lipid composition and discrete membrane proteins. Multiple applications are considered to validate and parameterize the model. The dynamics of membrane composition fluctuations above the critical point and phase separation dynamics below the critical point are studied in some detail, including the effects of adding proteins to the mixture.

Biomembranes in atomistic and coarse-grained simulations

The architecture of biological membranes is tightly coupled to the localization, organization, and function of membrane proteins. The organelle-specific distribution of lipids allows for the formation of functional microdomains (also called rafts) that facilitate the segregation and aggregation of membrane proteins and thus shape their function. Molecular dynamics simulations enable to directly access the formation, structure, and dynamics of membrane microdomains at the molecular scale and the specific interactions among lipids and proteins on timescales from picoseconds to microseconds. This review focuses on the latest developments of biomembrane force fields for both atomistic and coarse-grained molecular dynamics (MD) simulations, and the different levels of coarsening of biomolecular structures. It also briefly introduces scale-bridging methods applicable to biomembrane studies, and highlights selected recent applications.

Calculating hydrodynamic interactions for membrane-embedded objects

A recently introduced numerical scheme for calculating self-diffusion coefficients of solid objects embedded in lipid bilayer membranes is extended to enable calculation of hydrodynamic interactions between multiple objects. The method is used to validate recent analytical predictions by Oppenheimer and Diamant [Biophys. J. 96, 3041 2009] related to the coupled diffusion of membrane embedded proteins and is shown to converge to known near-field lubrication results as objects closely approach one another; however, the present methodology also applies outside of the limiting regimes where analytical results are available. Multiple different examples involving pairs of disk-like objects with various constraints imposed on their relative motions demonstrate the importance of hydrodynamic interactions in the dynamics of proteins and lipid domains on membrane surfaces. It is demonstrated that the relative change in self-diffusion of a membrane embedded object upon perturbation by a similar proximal solid object displays a maximum for object sizes comparable to the Saffman-Delbrück length of the membrane.

Molecular dynamics simulations of membrane proteins

Membrane protein structures are underrepresented in the Protein Data Bank (PDB) due to difficulties associated with expression and crystallization. As such, it is one area where computational studies, particularly Molecular Dynamics (MD) simulations, can provide useful additional information. Recently, there has been substantial progress in the simulation of lipid bilayers and membrane proteins embedded within them. Initial efforts at simulating membrane proteins embedded within a lipid bilayer were relatively slow and interactive processes, but recent advances now mean that the setup and running of membrane protein simulations is somewhat more straightforward, though not without its problems. In this chapter, we outline practical methods for setting up and running MD simulations of a membrane protein embedded within a lipid bilayer and discuss methodologies that are likely to contribute future improvements.

Instability of carbon nanoparticles interacting with lipid bilayers

For nanoparticles penetrating biological tissue, modelling indicates that without external forces, carbon nanoparticles will remain trapped in lipid bilayers.

Lipopolysaccharide membrane building and simulation

While membrane simulations are widely employed to study the structure and dynamics of various lipid bilayers and membrane proteins in the bilayers, simulations of lipopolysaccharides (LPS) in membrane environments have been limited due to their structural complexity, difficulties in building LPS-membrane systems, and lack of the appropriate molecular force fields. In this work, as a first step to extend CHARMM-GUI Membrane Builder to incorporate LPS molecules and to explore their structures and dynamics in membrane environments using molecular dynamics simulations, we describe step-by-step procedures to build LPS bilayer systems using CHARMM and the recently developed CHARMM carbohydrate and lipid force fields. Such procedures are illustrated by building various bilayers of Escherichia coli R1.O6 LPS and the presentation of preliminary simulation results in terms of per-LPS area and density distributions of various components along the membrane normal.

Exploring the sequence space for (tri-)peptide self-assembly to design and discover new hydrogels

Peptides that self-assemble into nanostructures are of tremendous interest for biological, medical, photonic and nanotechnological applications. The enormous sequence space that is available from 20 amino acids probably harbours many interesting candidates, but it is currently not possible to predict supramolecular behaviour from sequence alone. Here, we demonstrate computational tools to screen for the aqueous self-assembly propensity in all of the 8,000 possible tripeptides and evaluate these by comparison with known examples. We applied filters to select for candidates that simultaneously optimize the apparently contradicting requirements of aggregation propensity and hydrophilicity, which resulted in a set of design rules for self-assembling sequences. A number of peptides were subsequently synthesized and characterized, including the first reported tripeptides that are able to form a hydrogel at neutral pH. These tools, which enable the peptide sequence space to be searched for supramolecular properties, enable minimalistic peptide nanotechnology to deliver on its promise.

Assembly and stability of Salmonella enterica ser. Typhi TolC protein in POPE and DMPE

In this work we assessed the suitability of two different lipid membranes for the simulation of a TolC protein from Salmonella enterica serovar Typhi. The TolC protein family is found in many pathogenic Gram-negative bacteria including Vibrio cholera and Pseudomonas aeruginosa and acts as an outer membrane channel for expulsion of drug and toxin from the cell. In S. typhi, the causative agent for typhoid fever, the TolC outer membrane protein is an antigen for the pathogen. The lipid environment is an important modulator of membrane protein structure and function. We evaluated the conformation of the TolC protein in the presence of DMPE and POPE bilayers using molecular dynamics simulation. The S. typhi TolC protein exhibited similar conformational dynamics to TolC and its homologues. Conformational flexibility of the protein is seen in the C-terminal, extracellular loops, and α-helical region. Despite differences in the two lipids, significant similarities in the motion of the protein in POPE and DMPE were observed, including the rotational motion of the C-terminal residues and the partially open extracellular loops. However, analysis of the trajectories demonstrated effects of hydrophobic matching of the TolC protein in the membrane, particularly in the lengthening of the lipids and subtle movements of the protein's β-barrel towards the lower leaflet in DMPE. The study exhibited the use of molecular dynamics simulation in revealing the differential effect of membrane proteins and lipids on each other. In this study, POPE is potentially a more suitable model for future simulation of the S. typhi TolC protein.

A structural model of the active ribosome-bound membrane protein insertase YidC

The integration of most membrane proteins into the cytoplasmic membrane of bacteria occurs co-translationally. The universally conserved YidC protein mediates this process either individually as a membrane protein insertase, or in concert with the SecY complex. Here, we present a structural model of YidC based on evolutionary co-variation analysis, lipid-versus-protein-exposure and molecular dynamics simulations. The model suggests a distinctive arrangement of the conserved five transmembrane domains and a helical hairpin between transmembrane segment 2 (TM2) and TM3 on the cytoplasmic membrane surface. The model was used for docking into a cryo-electron microscopy reconstruction of a translating YidC-ribosome complex carrying the YidC substrate FOc. This structure reveals how a single copy of YidC interacts with the ribosome at the ribosomal tunnel exit and identifies a site for membrane protein insertion at the YidC protein-lipid interface. Together, these data suggest a mechanism for the co-translational mode of YidC-mediated membrane protein insertion.

Synthesis of coumarin derivatives as fluorescent probes for membrane and cell dynamics studies

Three coumarin-derived fluorescent probes, 3-acetyl-7-[(6-bromohexyl)oxy]-2H-chromen-2-one (FM1), 7-[(6-bromohexyl)oxy]-4-methyl-2H-chromen-2-one (FM2) and ethyl 2-{7-[(6-bromohexyl)oxy]-2-oxo-2H-chromen-4-yl}acetate (FM3), are described, with their photophysical constants. The compounds were tested in preliminary studies employing epifluorescence microscopy demonstrating that they allow the imaging of human neuroblastoma SH-SY5Y cell membranes. The structure of FM3 was confirmed by X-ray crystallographic analysis. Molecular dynamics (MD) simulations were used to characterize the localization and interactions of the studied compounds with a lipid bilayer model of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC).

Unveiling the gating mechanism of ECF transporter RibU

Energy-coupling factor (ECF) transporters are responsible for uptake of micronutrients in prokaryotes. The recently reported crystal structure of an ECF transporter RibU provided a foundation for understanding the structure and transport mechanism of ECF transporters. In the present study, molecular dynamics (MD) was carried out to study the conformational changes of the S component RibU upon binding by riboflavin. Our result and analysis revealed a critically important gating mechanism, in which part of loop5 (L5′) (eleven residues, missing in the crystal structure) between TM5 and TM6 is dynamically flexible and serves as a gate. Specifically, the L5′ opens a large cavity accessible to riboflavin from the extracellular space in Apo-RibU and closes the cavity upon riboflavin binding through hydrophobic packing with riboflavin. Thus, L5′is proposed to be the gate for riboflavin binding. In addition, steered molecular dynamics (SMD) simulation is employed to investigate the translocation dynamics of RibU during riboflavin transport. The simulation result does not show evidence that the S component alone can carry out the transport function. Since loop regions are very flexible and therefore could not be resolved by crystallography, their dynamics are hard to predict based on crystal structure alone.

Effects of PEGylation on the binding interaction of magainin 2 and tachyplesin I with lipid bilayer surface

Poly(ethylene glycol) (PEG)-grafted magainin 2 and tachyplesin I were simulated with lipid bilayers. In the simulations of PEGylated magainin 2 and tachyplesin I in water, both peptides are wrapped by PEG chains. The α-helical structure of PEGylated magainin 2 is broken in water, while the β-sheet of PEGylated tachyplesin I keeps stable, similar to the structural behavior of unPEGylated peptides, in agreement with experiments. Simulations of PEGylated peptides with lipid bilayers show that PEG chains block the electrostatic interaction between cationic residues of peptides and anionic phosphates of lipids, leading to the less binding of the peptide to the bilayer surface, which is observed more significantly for magainin 2 than for tachyplesin I. Since the random-coiled magainin 2 can be more completely covered by PEGs than does the β-sheet tachyplesin I, the PEGylation effect on the decreased binding is larger for magainin 2, showing the dependence of PEGylation on the peptide structure. These simulation findings qualitatively support the experimental observation of the different extents of the reduced membrane-permeabilizing activity for PEGylated magainin 2 and tachyplesin I.

Simulative and experimental investigation on the cleavage site that generates the soluble human LOX-1

Lectin-like oxidized low-density lipoprotein receptor-1 (LOX-1) is a scavenger receptor that mediates the recognition, the binding and internalization of ox-LDL. A truncated soluble form of LOX-1 (sLOX-1) has been identified that, at elevated levels, has been associated to acute coronary syndrome. Human sLOX-1 is the extracellular part of membrane LOX-1 which is cleaved in the NECK domain with a mechanism that has not yet been identified. Purification of human sLOX-1 has been carried out to experimentally identify the cleavage site region within the NECK domain. Molecular modelling and classical molecular dynamics simulation techniques have been used to characterize the structural and dynamical properties of the LOX-1 NECK domain in the presence and absence of the CTLD recognition region, taking into account the obtained proteolysis results. The simulative data indicate that the NECK domain is stabilized by the coiled-coil heptad repeat motif along the simulations, shows a definite flexibility pattern and is characterized by specific electrostatic potentials. The detection of a mobile inter-helix space suggests an explanation for the in vivo susceptibility of the NECK domain to the proteolytic cleavage, validating the assumption that the NECK domain sequence is composed of a coiled-coil motif destabilized in specific regions of functional significance.

Membrane perturbation action mode and structure-activity relationships of Protonectin, a novel antimicrobial peptide from the venom of the neotropical social wasp Agelaia pallipes pallipes

With the extensive use of antibiotics, multidrug-resistant bacteria emerge frequently. New antimicrobial agents with novel modes of action are urgently needed. It is now widely accepted that antimicrobial peptides (AMPs) could be promising alternatives to conventional antibiotics. In this study, we aimed to study the antimicrobial activity and mechanism of action of protonectin, a cationic peptide from the venom of the neotropical social wasp Agelaia pallipes pallipes. We demonstrated that protonectin exhibits potent antimicrobial activity against a spectrum of bacteria, including multidrug-resistant strains. To further understand this mechanism, the structural features of protonectin and its analogs were studied by circular dichroism (CD). The CD spectra demonstrated that protonectin and its natural analog polybia-CP formed a typical α-helical conformation in the membrane-mimicking environment, while its proline-substituted analog had much lower or even no α-helix conformation. Molecular dynamics simulations indicated that the α-helical conformation in the membrane is required for the exhibition of antibacterial activity. In conclusion, protonectin exhibits potent antimicrobial activity by disruption of the integrity of the bacterial membrane, and its α-helical confirmation in the membrane is essential for this action.

Membrane undulation induced by NS4A of Dengue virus: a molecular dynamics simulation study

Nonstructural protein 4A (NS4A) of Dengue virus (DENV) is a membrane protein involved in rearrangements of the endoplasmic reticulum membrane that are required for formation of replication vesicles. NS4A is composed most likely of three membrane domains. The N- and C-terminal domains are supposed to traverse the lipid membrane whereas the central one is thought to reside on the membrane surface, thus forming a u-shaped protein. All three membrane domains are proposed to be helical by secondary structure prediction programs. After performing multi nanosecond molecular dynamics (MD) simulations at various temperatures (300, 310, and 315.15 K) with each of the individual domains, they are used in a docking approach to define putative association motifs of the transmembrane domains (TMDs). Two structures of the u-shaped protein are generated by separating two assembled TMDs linking them with the membrane-attached domain. Lipid undulation is monitored with the structures embedded in a fully hydrated lipid bilayer applying multiple 200 ns MD simulations at 310 K. An intact structure of the protein supports membrane undulation. The strong unwinding of the helices in the domain-linking section of one of the structures lowers its capability to induce membrane curvature. Unwinding of the link region is due to interactions of two tryptophan residues, Trp-96 and 104. These results provide first insights into the membrane-altering properties of DENV NS4A.