Search

The fusion of the viral and target cell membranes is a key step in the life cycle of all enveloped viruses. Here, a range of structural data is used to generate an evidence-based model of the active conformation of an archetypical type-I fusion protein, the Ebola glycoprotein 2 (GP2). The stability of the trimeric complex is demonstrated using molecular dynamics and validated by simulating the interaction of the complex with a lipid bilayer. In this model, the fusion peptides project away from the central helix bundle parallel to the target membrane. This maximizes contact with the host membrane, enhances lateral stability, and would explain why, when activated, viral fusion proteins are trimeric.

α-Helical membrane-active antimicrobial peptides (AMPs) are known to act via a range of mechanisms, including the formation of barrel-stave and toroidal pores and the micellisation of the membrane (carpet mechanism). Different mechanisms imply that the peptides adopt different 3D structures when bound at the water-membrane interface, a highly amphipathic environment. Here, an evolutionary algorithm is used to predict the 3D structure of a range of α-helical membrane-active AMPs at the water-membrane interface by optimising amphipathicity. This amphipathic structure prediction (ASP) is capable of distinguishing between curved and linear peptides solved experimentally, potentially allowing the activity and mechanism of action of different membrane-active AMPs to be predicted. The ASP algorithm is accessible via a web interface at http://atb.uq.edu.au/asp/.

E-cadherin is a transmembrane glycoprotein that facilitates inter-cellular adhesion in the epithelium. The ectodomain of the native structure is comprised of five repeated immunoglobulin-like domains. All E-cadherin crystal structures show the protein in one of three alternative conformations: a monomer, a strand-swapped trans homodimer and the so-called X-dimer, which is proposed to be a kinetic intermediate to forming the strand-swapped trans homodimer. However, previous studies have indicated that even once the trans strand-swapped dimer is formed, the complex is highly dynamic and the E-cadherin monomers may reorient relative to each other. Here, molecular dynamics simulations have been used to investigate the stability and conformational flexibility of the human E-cadherin trans strand-swapped dimer. In four independent, 100 ns simulations, the dimer moved away from the starting structure and converged to a previously unreported structure, which we call the Y-dimer. The Y-dimer was present for over 90% of the combined simulation time, suggesting that it represents a stable conformation of the E-cadherin dimer in solution. The Y-dimer conformation is stabilised by interactions present in both the trans strand-swapped dimer and X-dimer crystal structures, as well as additional interactions not found in any E-cadherin dimer crystal structures. The Y-dimer represents a previously unreported, stable conformation of the human E-cadherin trans strand-swapped dimer and suggests that the available crystal structures do not fully capture the conformations that the human E-cadherin trans homodimer adopts in solution.

Based on differences between the x-ray crystal structures of ligand-bound and unbound forms, the activation of the erythropoietin receptor (EPOR) was initially proposed to involve a cross-action scissorlike motion. However, the validity of the motions involved in the scissorlike model has been recently challenged. Here, atomistic molecular dynamics simulations are used to examine the structure of the extracellular domain of the EPOR dimer in the presence and absence of erythropoietin and a series of agonistic or antagonistic mimetic peptides free in solution. The simulations suggest that in the absence of crystal packing effects, the EPOR chains in the different dimers adopt very similar conformations with no clear distinction between the agonist and antagonist-bound complexes. This questions whether the available x-ray crystal structures of EPOR truly represent active or inactive conformations. The study demonstrates the difficulty in using such structures to infer a mechanism of action, especially in the case of membrane receptors where just part of the structure has been considered in addition to potential confounding effects that arise from the comparison of structures in a crystal as opposed to a membrane environment. The work highlights the danger of assigning functional significance to small differences between structures of proteins bound to different ligands in a crystal environment without consideration of the effects of the crystal lattice and thermal motion.

Biobutanol production by fermentation is potentially a sustainable alternative to butanol production from fossil fuels. However, the toxicity of butanol to fermentative bacteria, resulting largely from cell membrane fluidization, limits production titers and is a major factor limiting the uptake of the technology. Here, studies were undertaken, in vitro and in silico, on the butanol effects on a representative bacterial (i.e. Escherichia coli) inner cell membrane. A critical butanol : lipid ratio for stability of 2 : 1 was observed, computationally, consistent with complete interdigitation. However, at this ratio the bilayer was ∼20% thicker than for full interdigitation. Furthermore, butanol intercalation induced acyl chain bending and increased disorder, measured as a 27% lateral diffusivity increase experimentally in a supported lipid bilayer. There was also a monophasic Tm reduction in butanol-treated large unilamellar vesicles. Both behaviours are inconsistent with an interdigitated gel. Butanol thus causes only partial interdigitation at physiological temperatures, due to butanol accumulating at the phospholipid headgroups. Acyl tail disordering (i.e. splaying and bending) fills the subsequent voids. Finally, butanol short-circuits the bilayer and creates a coupled system where interdigitated and splayed phospholipids coexist. These findings will inform the design of strategies targeting bilayer stability for increasing biobutanol production titers. Ministry of Education (MOE) Nanyang Technological University National Research Foundation (NRF) Accepted version The computational work for this article was performed on resources of the National Supercomputing Centre, Singapore (https://www.nscc.sg). This work was supported by the Ministry of Education, Singapore (MOE) Grants M4360005 and Tier 1 RG146/17. SCELSE is funded by Singapore’s Ministry of Education, National Research Federation, Nanyang Technological University (NTU), and National University of Singapore (NUS) and hosted by NTU in partnership with NUS.

Biological membranes display a great diversity in lipid composition and lateral structure that is crucial in a variety of cellular functions. Simulations of membranes have contributed significantly to the understanding of the properties, functions and behaviour of membranes and membrane-protein assemblies. This success relies on the ability of the force field used to describe lipid-lipid and lipid-environment interactions accurately, reproducibly and realistically. In this review, we present some recent progress in lipid force-field development and validation strategies. In particular, we highlight how a range of properties obtained from various experimental techniques on lipid bilayers and membranes, can be used to assess the quality of a force field. We discuss the limitations and assumptions that are inherent to both computational and experimental approaches and how these can influence the comparison between simulations and experimental data. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

The apparent activity of the multidrug transporter P-glycoprotein (P-gp) is enhanced by the presence of cholesterol. Whether this is due to the direct effect of cholesterol on the activity of P-gp, its effect on the local concentration of substrate in the membrane, or its effect on the rate of entry of the drug into the cell, is unknown. In this study, molecular dynamics simulation techniques coupled with potential of mean force calculations have been used to investigate the role of cholesterol in the movement of four P-gp substrates across a POPC bilayer in the presence or absence of 10% cholesterol. The simulations suggest that the presence of cholesterol lowers the free energy associated with entering the middle of the bilayer in a substrate-specific manner. These findings suggest that P-gp substrates may preferentially accumulate in cholesterol-rich regions of the membrane, which may explain its enhanced transport activity.

Free to read at publisher's site. Many venom peptides are potent and selective inhibitors of voltage-gated ion channels, including channels that are validated therapeutic targets for treatment of a wide range of human diseases. However, the development of novel venom-peptide-based therapeutics requires an understanding of their mechanism of action. In the case of voltage-gated ion channels, venom peptides act either as pore blockers that bind to the extracellular side of the channel pore or gating modifiers that bind to one or more of the membrane-embedded voltage sensor domains. In the case of gating modifiers, it has been debated whether the peptide must partition into the membrane to reach its binding site. In this study, we used surface plasmon resonance, fluorescence spectroscopy and molecular dynamics to directly compare the lipid-binding properties of two gating modifiers (μ-TRTX-Hd1a and ProTx-I) and two pore blockers (ShK and KIIIA). Only ProTx-I was found to bind to model membranes. Our results provide further evidence that the ability to insert into the lipid bilayer is not a requirement to be a gating modifier. In addition, we characterised the surface of ProTx-I that mediates its interaction with neutral and anionic phospholipid membranes and show that it preferentially interacts with anionic lipids.

The human multidrug transporter P-glycoprotein (P-gp) transports over 200 chemically diverse substrates, influencing their bioavailability and tissue distribution. Pharmacological studies have identified both competitive and noncompetitive P-gp substrates, but neither the precise location of the substrate binding sites, nor the basis of competitive and noncompetitive interactions has been fully characterized. Here, potential of mean force (PMF) calculations are used to identify the transport-competent minimum free energy binding locations of five compounds, Hoechst 33342, Rhodamine 123, paclitaxel, tariquidar, and verapamil to P-gp. Unrestrained molecular dynamics simulations were also performed to confirm the substrates were stable in the energy wells determined using the PMF calculations. All compounds had energy minima within the P-gp transmembrane (TM) pore. For Hoechst 33342 and Rhodamine 123, a second minimum outside the TM pore was also identified. Based on this and previous studies of nicardipine and morphine [ Subramanian et al. J. Chem. Inf. Model. 2015 , 55 , 1202 ], a general scheme that accounts for the observed noncompetitive and competitive substrate interactions with P-gp is proposed.

Background and Purpose The spider‐venom peptide PcTx1 is the most potent and selective inhibitor of acid‐sensing ion channel (ASIC) 1a. It has centrally acting analgesic activity and is neuroprotective in rodent models of ischaemic stroke. Understanding the molecular details of the PcTx1 : ASIC1a interaction should facilitate development of therapeutically useful ASIC1a modulators. Previously, we showed that several key pharmacophore residues of PcTx1 reside in a dynamic β‐hairpin loop; conclusions confirmed by recent crystal structures of the complex formed between PcTx1 and chicken ASIC1 (cASIC1). Numerous peptide : channel contacts were observed in these crystal structures, but it remains unclear which of these are functionally important.

The multidrug transporter P-glycoprotein (P-gp) is central to the development of multidrug resistance in cancer. While residues essential for transport and binding have been identified, the location, composition, and specificity of potential drug binding sites are uncertain. Here molecular dynamics simulations are used to calculate the free energy profile for the binding of morphine and nicardipine to P-gp. We show that morphine and nicardipine primarily interact with key residues implicated in binding and transport from mutational studies, binding at different but overlapping sites within the transmembrane pore. Their permeation pathways were distinct but involved overlapping sets of residues. The results indicate that the binding location and permeation pathways of morphine and nicardipine are not well separated and cannot be considered as unique. This has important implications for our understanding of substrate uptake and transport by P-gp. Our results are independent of the choice of starting structure and consistent with a range of experimental studies.

Phox-homology (PX) domains target proteins to the organelles of the secretary and endocytic systems by binding to phosphatidylinositol phospholipids (PIPs). Among all the structures of PX domains that have been solved, only three have been solved in a complex with the main physiological ligand: PtdIns3P. In this work, molecular dynamic simulations have been used to explore the structure and dynamics of the p40(phox)-PX domain and the SNX17-PX domain and their interaction with membrane-bound PtdIns3P. In the simulations, both PX domains associated spontaneously with the membrane-bound PtdIns3P and formed stable complexes. The interaction between the p40(phox)-PX domain and PtdIns3P in the membrane was found to be similar to the crystal structure of the p40(phox)-PX-PtdIns3P complex that is available. The interaction between the SNX17-PX domain and PtdIns3P was similar to that observed in the p40(phox)-PX-PtdIns3P complex; however, some residues adopted different orientations. The simulations also showed that nonspecific interactions between the beta 1-beta 2 loop and the membrane play an important role in the interaction of membrane bound PtdIns3P and different PX domains. The behaviour of unbound PtdIns3P within a 2-oleoyl-1-palmitoyl-sn-glycero-3-phosphocholine (POPC) membrane environment was also examined and compared to the available experimental data and simulation studies of related molecules. (C) 2014 Wiley Periodicals, Inc.

Crystallographic studies of ligands bound to biological macromolecules (proteins and nucleic acids) represent\ud an important source of information concerning drug-target interactions, providing atomic level insights\ud into the physical chemistry of complex formation between macromolecules and ligands. Of the\ud more than 115,000 entries extant in the Protein Data Bank (PDB) archive, ~75% include at least one non-polymeric\ud ligand. Ligand geometrical and stereochemical quality, the suitability of ligand models for in silico drug\ud discovery and design, and the goodness-of-fit of ligand models to electron-density maps vary widely across\ud the archive. We describe the proceedings and conclusions from the first Worldwide PDB/Cambridge Crystallographic\ud Data Center/Drug Design Data Resource (wwPDB/CCDC/D3R) Ligand Validation Workshop\ud held at the Research Collaboratory for Structural Bioinformatics at Rutgers University on July 30–31, 2015.\ud Experts in protein crystallography from academe and industry came together with non-profit and for-profit\ud software providers for crystallography and with experts in computational chemistry and data archiving to\ud discuss and make recommendations on best practices, as framed by a series of questions central to structural\ud studies of macromolecule-ligand complexes. What data concerning bound ligands should be archived\ud in the PDB? How should the ligands be best represented? How should structural models of macromoleculeligand\ud complexes be validated? What supplementary information should accompany publications of structural\ud studies of biological macromolecules? Consensus recommendations on best practices developed in\ud response to each of these questions are provided, together with some details regarding implementation.\ud Important issues addressed but not resolved at the workshop are also enumerated.

The Type-I Cytokine Receptors regulate diverse cellular functions that encompass immune responses, growth, metabolism, lactation and red blood cell production. A sub-group of receptors that are models for receptor activation includes the growth hormone receptor (GHR), prolactin receptor (PRLR) and erythropoietin receptor (EPOR). While the cellular effects of receptor activation have been well characterised, the atomic details of how these receptors mechanically couple cytokine binding to the extracellular domains through the plasma membrane has remained elusive. Using fully atomistic molecular dynamics simulations, we have shown the extracellular domains of the GHR and PRLR homodimers undergo a relative rigid body rotation upon activation. Recently, we have tested if the rotation-based mechanism could be extended to erythropoietin (EPO) induced activation of the EPOR homodimer along with agonistic and antagonistic EPO mimetic peptides. In the simulations, the EPOR protomers rotate as rigid bodies similarly to the GHR and PRLR upon cytokine activation. However, analysis of the mimetic-bound EPOR could not distinguish between agonised or antagonised conformations of the dimer. Simulations of the GHR, PRLR and EPOR embedded in membranes exhibited different behavior in cholesterol rich membranes, suggesting a key role of cholesterol in receptor activation. The simulations demonstrate a shared cytokine induced rotational-based mechanism of activation of Type-I Cytokine Receptors. The atomic resolution of the mechanical coupling of the receptors through the plasma membrane provides clearer understanding of how the signal of cytokine binding is transmitted from the extracellular to the cytosolic domains.

The growth hormone receptor (GHR) has been the archetype for the class I cytokine receptor family. The original paradigm for how growth hormone (GH) binding to the GHR led to activation of the associated JAK2 was based on biophysical and structural studies over 20 years ago which showed that growth hormone binding to the extracellular domain of the receptor occurred in a sequential fashion, first binding to a high affinity site on a single receptor (site 1) and then binding to a lower affinity site (site 2) on a second receptor. We propose a new model for JAK2 activation by the growth hormone receptor. We utilised multiple approaches to define this model including FRET to monitor positioning of the JAK2 binding motif, leucine-zipper receptor constructs to control receptor transmembrane helix position, atomistic modelling of TM helix interactions and docking of JAK2 kinase and inhibitory pseudokinase crystal structures with an opposing pair in trans. Surprisingly, we identified that receptor activation induces a separation of its JAK2 binding motifs which leads to removal of the pseudokinase domain from the kinase domain of the partner JAK2 and pairing of the two kinase domains, facilitating trans-activation. In addition, we have recently shown that the receptor intracellular domain is intrinsically disordered and has specific interactions with the lipids of the inner membrane. This model of receptor activation may represent a common mechanism to other class I cytokine receptors.

Free to read at publisher's site. ProTx-II is a disulfide-rich peptide toxin from tarantula venom able to inhibit the human voltage-gated sodium channel 1.7 (hNaV1.7), a channel reported to be involved in nociception, and thus it might have potential as a pain therapeutic. ProTx-II acts by binding to the membrane-embedded voltage sensor domain of hNaV1.7, but the precise peptide channel-binding site and the importance of membrane binding on the inhibitory activity of ProTx-II remain unknown. In this study, we examined the structure and membrane-binding properties of ProTx-II and several analogues using NMR spectroscopy, surface plasmon resonance, fluorescence spectroscopy, and molecular dynamics simulations. Our results show a direct correlation between ProTx-II membrane binding affinity and its potency as an hNaV1.7 channel inhibitor. The data support a model whereby a hydrophobic patch on the ProTx-II surface anchors the molecule at the cell surface in a position that optimizes interaction of the peptide with the binding site on the voltage sensor domain. This is the first study to demonstrate that binding of ProTx-II to the lipid membrane is directly linked to its potency as an hNaV1.7 channel inhibitor.

Acid-sensing ion channel 1a (ASIC1a) is a primary acid sensor in the peripheral and central nervous system. It has been implicated as a novel therapeutic target for a broad range of pathophysiological conditions including pain, ischemic stroke, depression, and autoimmune diseases such as multiple sclerosis. The only known selective blocker of ASIC1a is π-TRTX-Pc1a (PcTx1), a disulfide-rich 40-residue peptide isolated from spider venom. π-TRTX-Pc1a is an effective analgesic in rodent models of acute pain and it provides neuroprotection in a mouse model of ischemic stroke. Thus, understanding the molecular basis of the π-TRTX-Pc1a-ASIC1a interaction should facilitate development of therapeutically useful ASIC1a blockers. We therefore developed an efficient bacterial expression system to produce a panel of π-TRTX-Pc1a mutants for probing structure-activity relationships as well as isotopically labeled toxin for determination of its solution structure and dynamics. We demonstrate that the toxin pharmacophore resides in a β-hairpin loop that was revealed to be mobile over a wide range of time scales using molecular dynamics simulations in combination with NMR spin relaxation and relaxation dispersion measurements. The toxin-receptor interaction was modeled by in silico docking of the toxin structure onto a homology model of rat ASIC1a in a restraints-driven approach that was designed to take account of the dynamics of the toxin pharmacophore and the consequent remodeling of side-chain conformations upon receptor binding. The resulting model reveals new insights into the mechanism of action of π-TRTX-Pc1a and provides an experimentally validated template for the rational design of therapeutically useful π-TRTX-Pc1a mimetics.

Short cationic antimicrobial peptides (AMPs) are believed to act either by inducing transmembrane pores or disrupting membranes in a detergent-like manner. For example, the antimicrobial peptides aurein 1.2, citropin 1.1, maculatin 1.1 and caerin 1.1, despite being closely related, appear to act by fundamentally different mechanisms depending on their length. Using molecular dynamics simulations, the structural properties of these four peptides have been examined in solution as well as in a variety of membrane environments. It is shown that each of the peptides has a strong preference for binding to regions of high membrane curvature and that the structure of the peptides is dependent on the degree of local curvature. This suggests that the shorter peptides aurein 1.2 and citropin 1.1 act via a detergent-like mechanism because they can induce high local, but not long-range curvature, whereas the longer peptides maculatin 1.1 and caerin 1.1 require longer range curvature to fold and thus bind to and stabilize transmembrane pores.

Atomistic molecular dynamics simulations have been used to investigate the conformational changes associated with the binding of human growth hormone (hGH) to the extracellular domains (ECD) of the human growth hormone receptor (hGHR), thereby shedding light on the mechanism of activation. It is shown that the removal of hGH from the hormone-bound receptor complex results in a counter-clockwise rotation of the twosubunits relative to each other by 30°–64° (average 45° ± 14°), in close agreement in terms of both the magnitude and direction of the rotation with that proposed based on mutagenesis experiments. In addition to providing evidence to support a rotational activation mechanism, the simulations have enabled the nature of the interaction interfaces in both the cytokine-bound and unliganded hGHR states to be analyzed in detail. Proteins 2010. © 2009 Wiley-Liss, Inc.

A large variety of antimicrobial peptides have been shown to act, at least in vitro, by potation of the lipid membrane. The nanometre size of these pores, however, complicates their structural characterization by experimental techniques. Here we use molecular dynamics simulations, to study the interaction of a specific class of antimicrobial peptides, melittin, with a dipalmitoylphosphatidylcholine bilayer in atomic detail. We show that transmembrane pores spontaneously form above a critical peptide to lipid ratio. The lipid molecules bend inwards to form a toroidally shaped pore but with only one or two peptides lining the pore. This is in strong contrast to the traditional models of toroidal pores in which the peptides are assumed to adopt a transmembrane orientation. We find that peptide aggregation, either prior or after binding to the membrane surface. is a prerequisite to pore formation. The presence of a stable helical secondary structure of the peptide, however is not. Furthermore, results obtained with modified peptides point to the importance of electrostatic interactions in the potation process. Removing the charges of the basic amino-acid residues of melittin prevents pore formation. It was also found that in the absence of counter ions pores not only form more rapidly but lead to membrane rupture. The rupture process occurs via a novel recursive potation pathway, which we coin the Droste mechanism. (c) 2008 Elsevier B.V. All rights reserved.

Class I hydrophobins are fungal proteins that self-assemble into robust amphipathic rodlet monolayers on the surface of aerial structures such as spores and fruiting bodies. These layers share many structural characteristics with amyloid fibrils and belong to the growing family of functional amyloid-like materials produced by microorganisms. Although the three-dimensional structure of the soluble monomeric form of a class I hydrophobin has been determined, little is known about the molecular structure of the rodlets or their assembly mechanism. Several models have been proposed, some of which suggest that the Cys3–Cys4 loop has a critical role in the initiation of assembly or in the polymeric structure. In order to provide insight into the relationship between hydrophobin sequence and rodlet assembly, we investigated the role of the Cys3–Cys4 loop in EAS, a class I hydrophobin from Neurospora crassa. Remarkably, deletion of up to 15 residues from this 25-residue loop does not impair rodlet formation or reduce the surface activity of the protein, and the physicochemical properties of rodlets formed by this mutant are indistinguishable from those of its full-length counterpart. In addition, the core structure of the truncation mutant is essentially unchanged. Molecular dynamics simulations carried out on the full-length protein and this truncation mutant binding to an air–water interface show that, although it is hydrophobic, the loop does not play a role in positioning the protein at the surface. These results demonstrate that the Cys3–Cys4 loop does not have an integral role in the formation or structure of the rodlets and that the major determinant of the unique properties of these proteins is the amphipathic core structure, which is likely to be preserved in all hydrophobins despite the high degree of sequence variation across the family.

The Death Receptor 5 (DR5), a member of tumor necrosis factor receptor (TNFR) superfamily of receptors, triggers apoptosis (programmed cell death) when stimulated by its tridentate ligand TRAIL. Until recently it was generally assumed that the activation of DR5 resulted from the recruitment of three independent receptor units, leading to the trimerization of intracellular domains. However, there is mounting evidence to suggest that, in the absence of ligand, such cytokine receptors primarily reside as preformed complexes. In this work, molecular dynamics simulations of the TRALL-DR5 complex, the unbound receptor trimer and individual receptor monomers are compared to gain insight in the mechanism of activation. The results suggest that, in the absence of TRAIL, DR5 has a strong propensity to self-associate and that this is primarily mediated through interactions of the membrane proximal domains. The association of the free receptors leads to a loss of the threefold symmetry found within the receptor-ligand complex. The simulations suggest that the primary role of TRAIL is to induce threefold-symmetry within the DR5 complex and to constrain the receptor to a specific conformation. The implications of this in terms of the mechanism by which the receptor switches from an inactive to an active state are discussed.

The question of why plants release isoprene when heat stressed has been continuously debated for more than half a century. In this work we use molecular dynamics simulation techniques to directly investigate the interaction between isoprene and a model phospholipid membrane in atomic detail. It is found that isoprene partitions preferentially in the center of the membrane and in a dose dependent manner enhances the order within the membrane without significantly changing the dynamical properties of the system. At a concentration of 20 mol% isoprene (16 isoprene molecules per 64 lipid molecules) the effect of the addition of isoprene on the membrane order is equivalent to a reduction in temperature of 10 K, rising to a reduction of 30 K at 43 mol% isoprene. The significance of the work is that it provides for the first time direct evidence that isoprene stabilizes lipid membranes and reduces the likelihood of a phospholipid membrane undergoing a heat induced phase transition. Furthermore it provides a clear mechanistic picture as to why plants specifically utilize isoprene for this purpose.

The stability of the crystal structure of the multidrug transporter P-glycoprotein proposed by Aller et al. (PDBid 3G5U ) has been examined under different environmental conditions using molecular dynamics. We show that in the presence of the detergent cholate, the structure of P-glycoprotein solved at pH 7.5 is stable. However, when incorporated into a cholesterol-enriched POPC membrane in the presence of 150 mM NaCl, the structure rapidly deforms. Only when the simulation conditions closely matched the experimental conditions under which P-glycoprotein is transport active was a stable conformation obtained. Specifically, the presence of Mg(2+), which bound to distinct sites in the nucleotide binding domains (NBDs), and the double protonation of the catalytic histidines (His583 and His1228) and His149 were required. While the structure obtained in a membrane environment under these conditions is very similar to the crystal structure of Aller et al., there are several key differences. The NBDs are in direct contact, reminiscent of the open state of MalK. The angle between the transmembrane domains is also increased, resulting in an outward motion of the intracellular loops. Notably, the structures obtained from the simulations provide a better match to a range of experimental cross-linking data than does the original 3G5U-a crystal structure. This work highlights the effect small changes in environmental conditions can have of the conformation of a membrane protein and the importance of representing the experimental conditions appropriately in modeling studies.

Molecular dynamics simulations of the magainin MG-H2 peptide interacting with a model phospholipid membrane have been used to investigate the mechanism by which antimicrobial peptides act. Multiple copies of the peptide were randomly placed in solution close to the membrane. The peptide readily bound to the membrane, and above a certain concentration, the peptide was observed to cooperatively induce the formation of a nanometer- sized, toroidally shaped pore in the bilayer. In sharp contrast with the commonly accepted model of a toroidal pore, only one peptide was typically found near the center of the pore. The remaining peptides lay close to the edge of the pore, maintaining a predominantly parallel orientation with respect to the membrane.

Bacterial chaperonin, GroEL, together with its co-chaperonin, GroES, facilitates the folding of a variety of polypeptides. Experiments suggest that GroEL stimulates protein folding by multiple cycles of binding and release. Misfolded proteins first bind to an exposed hydrophobic surface on GroEL. GroES then encapsulates the substrate and triggers its release into the central cavity of the GroEL/ES complex for folding. In this work, we investigate the possibility to facilitate protein folding in molecular dynamics simulations by mimicking the effects of GroEL/ES namely, repeated binding and release, together with spatial confinement. During the binding stage, the (metastable) partially folded proteins are allowed to attach spontaneously to a hydrophobic surface within the simulation box. This destabilizes the structures, which are then transferred into a spatially confined cavity for folding. The approach has been tested by attempting to refine protein structural models generated using the ROSETTA procedure for ab initio structure prediction. Dramatic improvements in regard to the deviation of protein models from the corresponding experimental structures were observed. The results suggest that the primary effects of the GroEL/ES system can be mimicked in a simple coarse-grained manner and be used to facilitate protein folding in molecular dynamics simulations. Furthermore, the results Sur port the assumption that the spatial confinement in GroEL/ES assists the folding of encapsulated proteins.

The formation, relative stability, and possible stoichiometries of two (self-)complementary peptide sequences (B and E) designed to form either a parallel homodimeric (B + B) or an antiparallel heterodimeric (B + E) coiled coil have been investigated. Peptide B shows a characteristic coiled coil pattern in circular dichroism spectra at pH 7.4, whereas peptide E is apparently random coiled under these conditions. The peptides are complementary to each other, with peptide E forming a coiled coil when mixed with peptide B. Molecular dynamics simulations show that combinations of B + B and B + E readily form a dimeric coiled coil, whereas E + E does not fall in line with the experimental data. However, the simulations strongly suggest the preferred orientation of the helices in the homodimeric coiled coil is antiparallel, with interactions at the interface quite different to that of the idealized model. In addition, molecular dynamics simulations suggest equilibrium between dimers, trimers, and tetramers of alpha-helices for peptide B.

HIV-1 protease is most active under weakly acidic conditions (pH 3.5-6.5), when the catalytic Asp25 and Asp25' residues share 1 proton. At neutral pH, this proton is lost and the stability of the structure is reduced. Here we present an investigation of the effect of pH on the dynamics of HIV-1 protease using MD simulation techniques. MD simulations of the solvated HIV-1 protease with the Asp25/25' residues monoprotonated and deprotonated have been performed. In addition we investigated the effect of the inclusion of Na+ and Cl- ions to mimic physiological salt conditions. The simulations of the monoprotonated form and deprotonated form including Na+ show very similar behavior. In both cases the protein remained stable in the compact, "self-blocked" conformation in which the active site is blocked by the tips of the flaps. In the deprotonated system a Na+ ion binds tightly to the catalytic dyad shielding the repulsion between the COO- groups. Ab initio calculations also suggest the geometry of the active site with the Na+ bound closely resembles that of the monoprotonated case. In the simulations of the deprotonated form (without Na+ ions), a water molecule bound between the Asp25 Asp25' side-chains. This disrupted the dimerization interface and eventually led to a fully open conformation. (C) 2004 Wiley-Liss, Inc.

A novel method for the refinement of misfolded protein structures is proposed in which the properties of the solvent environment are oscillated in order to mimic some aspects of the role of molecular chaperones play in protein folding in vivo. Specifically, the hydrophobicity of the solvent is cycled by repetitively altering the partial charges on solvent molecules (water) during a molecular dynamics simulation. During periods when the hydrophobicity of the solvent is increased, intramolecular hydrogen bonding and secondary structure formation are promoted. During periods of increased solvent polarity, poorly packed regions of secondary structures are destabilized, promoting structural rearrangement. By cycling between these two extremes, the aim is to minimize the formation of long-lived intermediates. The approach has been applied to the refinement of structural models of three proteins generated by using the ROSETTA procedure for ab initio structure prediction. A significant improvement in the deviation of the model structures from the corresponding experimental structures was observed. Although preliminary, the results indicate computationally mimicking some functions of molecular chaperones in molecular dynamics simulations can promote the correct formation of secondary structure and thus be of general use in protein folding simulations and in the refinement of structural models of small- to medium-size proteins.

Here, we use coarse grained molecular dynamics (MD) simulations to study the spontaneous aggregation of dipalmitoylphosphatidylcholine (DPPC) lipids into small unilamellar vesicles. We show that the aggregation process occurs on a nanosecond time scale, with bicelles and cuplike vesicles formed at intermediate stages. Formation of hemifused vesicles is also observed at higher lipid concentration. With either 25% dipalmitoylphosphatidylethanolamine (DPPE) or lysoPC mixed into the system, the final stages of the aggregation process occur significantly faster. The structure of the spontaneously formed vesicles is analyzed in detail. Microsecond simulations of isolated vesicles reveal significant differences in the packing of the lipids between the inner and outer monolayers, and between PC, PE, and lysoPC. Due to the small size of the vesicles they remain almost perfectly spherical, undergoing very limited shape fluctuations or bilayer undulations. The lipid lateral diffusion rate is found to be faster in the outer than in the inner monolayer. The water permeability coefficient of the pure DPPC vesicles is of the order of 10(-3) cm s(-1), in agreement with experimental measurements.

The initial stages of the gating of the mechanoselective channel of large conductance from Mycobacterium tuberculosis have been studied in atomic detail using molecular dynamics simulation techniques. A truncated form of the protein complex embedded in a palmitoyloleoylphosphatidylcholine lipid bilayer and surrounded by explicit water was simulated under different pressure conditions to mimic the effects of tension and compression within the membrane on the protein. As a direct result of lateral tension being applied to the membrane, an increase in the tilt of a subset of the transmembrane helices was observed. This in turn led to the enlargement of the pore and the disruption of the hydrophobic gate consisting of residues Ile-14 and Val-21. The simulations suggest that opening occurs in a sequential staged process. Such a mechanism could explain the partial opening or staged conductance observed in patch-clamp experiments using related large conductance mechanosensitive channel complexes.

Molecular dynamics simulation techniques have been used to investigate the effect of 2,2,2-trifluoroethanol (TFE) as a cosolvent on the stability of three different secondary structure-forming peptides: the α-helix from Melittin, the three-stranded β-sheet peptide Betanova, and the β-hairpin 41–56 from the B1 domain of protein G. The peptides were studied in pure water and 30% (vol/vol) TFE/water mixtures at 300 K. The simulations suggest that the stabilizing effect of TFE is induced by the preferential aggregation of TFE molecules around the peptides. This coating displaces water, thereby removing alternative hydrogen-bonding partners and providing a low dielectric environment that favors the formation of intrapeptide hydrogen bonds. Because TFE interacts only weakly with nonpolar residues, hydrophobic interactions within the peptides are not disrupted. As a consequence, TFE promotes stability rather than inducing denaturation.

The cell surface epidermal growth factor receptor (EGFR) plays a critical role in cell development and oncogenesis. The binding of growth factors to the EGFR results in a mechanical signal being transmitted through the plasma membrane. In this study, atomistic molecular dynamics simulations have been used to investigate the conformational changes associated with the binding of the epidermal growth factor (EGF) and transforming growth factor α (TGFα) to the EGFR. In the simulations, the removal of the EGF and TGFα from the extracellular domain of the EGFR homodimer led to a relative rotation of the protomers of 16-35° about the dimerization axis. The three N-terminal domains that make up the extracellular region of the receptor undergo essentially rigid-body motion. The dimerization interface itself was found to be largely unaffected by the removal of the ligand. In most simulations, the rotation within the dimer was associated with an opening of the cytokine-binding sites. On the basis of these simulations, a simple mechanical model that explains the coupling between the binding of ligand and the motions in the extracellular domains is proposed.

Modeling the mechanisms of protofibril twisting: Molecular-dynamics simulations of simian viral peptide aggregates show that β sheets of 10 to 30 chains form left-handed helical ribbons with saddlelike curvature (see picture). These structures are highly dynamic, with oscillations around an average twist angle of 9-10°, and a pitch of 15-20 nm, depending on β-sheet length. The peptides studied are key to viral entry into host cells.

ATP Binding Cassette (ABC) transporters couple the binding and hydrolysis of ATP to the transport of substrate molecules across the membrane. The mechanism by which ATP binding and/or hydrolysis drives the conformational changes associated with substrate transport has not yet been characterized fully. Here, changes in the conformation of the ABC export protein P-glycoprotein on ATP binding are examined in a series of molecular dynamics simulations. When one molecule of ATP is placed at the ATP binding site associated with each of the two nucleotide binding domains (NBDs), the membrane-embedded P-glycoprotein crystal structure adopts two distinct metastable conformations. In one, each ATP molecule interacts primarily with the Walker A motif of the corresponding NBD. In the other, the ATP molecules interacts with both Walker A motif of one NBD and the Signature motif of the opposite NBD inducing the partial dimerization of the NBDs. This interaction is more extensive in one of the two ATP binding site, leading to an asymmetric structure. The overall conformation of the transmembrane domains is not altered in either of these metastable states, indicating that the conformational changes associated with ATP binding observed in the simulations in the absence of substrate do not lead to the outward-facing conformation and thus would be insufficient in themselves to drive transport. Nevertheless, the metastable intermediate ATP-bound conformations observed are compatible with a wide range of experimental cross-linking data demonstrating the simulations do capture physiologically important conformations. Analysis of the interaction between ATP and its cofactor Mg(2+) with each NBD indicates that the coordination of ATP and Mg(2+) differs between the two NBDs. The role structural asymmetry may play in ATP binding and hydrolysis is discussed. Furthermore, we demonstrate that our results are not heavily influenced by the crystal structure chosen for initiation of the simulations.

This chapter discusses the study of the interaction of peptides and proteins at interfaces using molecular dynamics (MD) simulation techniques. First, the chapter discusses how computational methods, in particular MD simulation techniques, complement experimental methods. Then it focuses on three main areas of research: the interaction of peptides and proteins with (i) biological membranes and lipid bilayers, (ii) airâwater and oilâwater interfaces, and (iii) organic and inorganic sorbents.

Molecular dynamics simulations have been used to investigate the effect of DMSO on 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC) and 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) phospholipid bilayers. The concentration of DMSO was varied between 0 and 25.0 mol %. For both lipids, DMSO causes the membrane to expand in the plane of the membrane while thinning normal to that plane. Above a critical concentration, pores in the membrane form spontaneously, and if the concentration is increased further, then the bilayer structure is destroyed. Even at concentrations below those required to induce pores, DMSO readily diffuses across the bilayers. The free-energy profile associated with the diffusion of a DMSO molecules across the membrane has been calculated. The simulations suggest that the DOPC bilayer is more resistant to the deleterious effects of DMSO, both increasing the stability of the membranes and decreasing the rate at which DMSO diffuses across the membrane. In this way, the work highlights the importance of investigating the lipid composition of cell membranes when characterizing the effects of cryosolvents.

Biofilm formation, in which bacteria are embedded within an extracellular matrix, is the default form of microbial life in most natural and engineered habitats. In this work, atomistic molecular dynamics simulations were employed to examine the self-assembly of the polysaccharide Granulan to provide insight into the molecular interactions that lead to biofilm formation. Granulan is a major gel forming matrix component of granular microbial biofilms found in used-water treatment systems. Molecular dynamics simulations showed that Granulan forms an antiparallel double helix stabilized by complementary hydrogen bonds between the β-glucosamine of one strand and the N-acetyl-β-galactosamine-2-acetoamido-2-deoxy-α-galactopyranuronic pair of the other in both the presence and absence of Ca(2+). It is shown that Ca(2+) binds primarily to the carboxyl group of the terminal hexuronic acid of the sugar branch and that interactions between branches mediated by Ca(2+) suggest a possible mechanism for strengthening gels by facilitating interhelical bridging.

European Biophysics Journal, 40 (7), ISSN:0175-7571, ISSN:1432-1017

Resorcinolic lipids, or resorcinols, are commonly found in plant membranes. They consist of a substituted benzene ring forming the hydrophilic lipid head, attached to an alkyl chain forming the hydrophobic tail. Experimental results show alternative effects of resorcinols on lipid membranes. Depending on whether they are added to lipid solutions before or after the formation of the liposomes, they either stabilize or destabilize these liposomes. Here we use atomistic molecular dynamics simulations to elucidate the molecular nature of this dual effect. Systems composed of either one of three resorcinol homologs, differing in the alkyl tail length, interacting with dimyristoylphosphatidylcholine lipid bilayers were studied. It is shown that resorcinols preincorporated into bilayers induce order within the lipid acyl chains, decrease the hydration of the lipid headgroups, and make the bilayers less permeable to water. In contrast, simulations in which the resorcinols are incorporated from the aqueous solution into a preformed phospholipid bilayer induce local disruption, leading to either transient pore formation or even complete rupture of the membrane. In line with the experimental data, our simulations thus demonstrate that resorcinols can either disturb or stabilize the membrane structure, and offer a detailed view of the underlying molecular mechanism.

Computer-based molecular simulation techniques are increasingly used to interpret experimental data on biomolecular systems at an atomic level. Direct comparison between experiment and simulation is, however, seldom straightforward. The available experimental data are limited in scope and generally correspond to averages over both time and space. A critical analysis of the various factors that may influence the apparent degree of agreement between the results of simulations and experimentally measured quantities is presented and illustrated using examples from recent literature.

Molecular dynamics simulations have been used to study the spontaneous aggregation of a concentrated solution of dipalmitoylphosphatidylcholine (DPPC) molecules in water into a small vesicle. The molecules were represented in atomistic detail. Starting from a DPPC solution in water, an oblong vesicle with a long axis of 15 nm and short axes of 10 nm was formed spontaneously. After 90 ns of simulation, the vesicle contained a number of water pores. Water pores were shown to facilitate exchange of lipids between inner and outer leaflets. Lipid tails were shown to be less ordered in the inner leaflet of the vesicle, as compared to those in the outer leaflet of the vesicle and an equilibrated lamellar bilayer.

Hydrophilic pores are formed in peptide free lipid bilayers under mechanical stress. It has been proposed that the transport of ionic species across such membranes is largely determined by the existence of such meta-stable hydrophilic pores. To study the properties of these structures and understand the mechanism by which pore expansion leads to membrane rupture, a series of molecular dynamics simulations of a dipalmitoylphosphatidylcholine (DPPC) bilayer have been conducted. The system was simulated in two different states; first, as a bilayer containing a meta-stable pore and second, as an equilibrated bilayer without a pore. Surface tension in both cases was applied to study the formation and stability of hydrophilic pores inside the bilayers. It is observed that below a critical threshold tension of approximately 38 mN/m the pores are stabilized. The minimum radius at which a pore can be stabilized is 0.7 nm. Based on the critical threshold tension the line tension of the bilayer was estimated to be approximately 3 x 10(-11) N, in good agreement with experimental measurements. The flux of water molecules through these stabilized pores was analyzed, and the structure and size of the pores characterized. When the lateral pressure exceeds the threshold tension, the pores become unstable and start to expand causing the rupture of the membrane. In the simulations the mechanical threshold tension necessary to cause rupture of the membrane on a nanosecond timescale is much higher in the case of the equilibrated bilayers, as compared with membranes containing preexisting pores.

Molecular dynamics simulations of pore formation and mem-brane rupture in phospholipid bilayers under mechanical andelectrical stress at an atomic level are presented. Pore formationcan be induced on a nanosecond time scale in simulations wherethe lateral pressure exceeds -200 bar or where an electric field of0.5 V/nm is applied across the membrane.Lipid bilayer membranes are remarkable structures consistingof two leaflets of phospholipids. Their mechanical properties arecentral to understanding the behavior of cell membranes. Forexample, the formation of transient water pores is believed tounderlie the passive transport of protons and hydrophilic compoundsthrough the bilayer. Pore formation is also relevant during processessuch as cell fusion and is important for drug release from liposomes.Experimentally, pores can be induced in membranes by applyingmechanical stress (i.e., pipet aspiration experiments) or an electricfield (electroporation). The precise mechanism by which pores formand their size, structure, and stability are, however, poorlyunderstood.

Hydrophobins are fungal proteins that self-assemble at hydrophilic/hydrophobic interfaces into amphipathic membranes. These assemblages are extremely stable and posses the remarkable ability to invert the polarity of the surface on which they are adsorbed. Neither the three-dimensional structure of a hydrophobin nor the mechanism by which they function is known. Nevertheless, there are experimental indications that the self-assembled form of the hydrophobins SC3 and EAS at a water/air interface is rich with beta-sheet secondary structure. In this paper we report results from molecular dynamics simulations, showing that fully extended SC3 undergoes fast (approximately 100 ns) folding at a water/hexane interface to an elongated planar structure with extensive beta-sheet secondary elements. Simulations in each of the bulk solvents result in a mainly unstructured globular protein. The dramatic enhancement in secondary structure, whether kinetic or thermodynamic in origin, highlights the role interfaces between phases with large differences in polarity can have on folding. The partitioning of the residue side-chains to one of the two phases can serve as a strong driving force to initiate secondary structure formation. The interactions of the side-chains with the environment at an interface can also stabilize configurations that otherwise would not occur in a homogenous solution.

Molecular dynamics (MD) simulations and polarized subnanosecond time-resolved flavin fluorescence spectroscopy have been used to study the conformational dynamics of the flavin adenine dinucleotide (FAD) cofactor in aqueous solution. FAD displays a highly heterogeneous fluorescence intensity decay, resulting in lifetime spectra with two major components: a dominant 7-ps contribution that is characteristic of ultrafast fluorescence quenching and a 2.7-ns contribution resulting from moderate quenching. MD simulations were performed in both the ground state and first excited state. The simulations showed transitions from "open" conformations to "closed" conformations in which the flavin and adenine ring systems stack coplanarly. Stacking generally occurred within the lifetime of the flavin excited state (4.7 ns in water), and yielded a simulated fluorescence lifetime on the order of the nanosecond lifetime that was observed experimentally. Hydrogen bonds in the ribityl-pyrophosphate-ribofuranosyl chain connecting both ring systems form highly stable cooperative networks and dominate the conformational transitions of the molecule. Fluorescence quenching in FAD is mainly determined by the coplanar stacking of the flavin and adenine ring systems, most likely through a mechanism of photoinduced electron transfer. Whereas in stacked conformations fluorescence is quenched nearly instantaneously, open fluorescent conformations can stack during the lifetime of the flavin excited state, resulting in immediate fluorescence quenching upon stacking.

Molecular dynamics simulations of protein folding and unfolding are often carried out at temperatures (400–600K) that are much higher than physiological or room temperature to speed up the (un)folding process. Use of such high temperatures changes both the protein and solvent properties considerably, compared to physiological or room temperature. Water models designed for use in conjunction with biomolecules, such as the simple point charge (SPC) model, have generally been calibrated at room temperature and pressure. To determine the distortive effect of high simulation temperatures on the behavior of such “room temperature” water models, the structural, dynamic, and thermodynamic properties of the much-used SPC water model are investigated in the temperature range from 300 to 500K. Both constant pressure and constant volume conditions, as used in protein simulations, were analyzed. We found that all properties analyzed change markedly with increasing temperature, but no phase transition in this temperature range was observed.

Antimicrobial peptides of the innate immune system represent our first line of defense against infection. They exhibit a broad range of microbicidal activity encompassing Gram-positive and -negative bacteria, mycobacteria and spirochetes as well as some fungi and enveloped viruses. Common antimicrobial peptides act by forming transmembrane channels or pores.1,2 Despite these seemingly generic properties their activity toward certain cell types is often highly sequence dependent.3 This raises critical questions as to how these peptides discriminate between different cell types based on the membrane composition. Using molecular dynamics simulations, the structural properties of a range of pore forming peptides has been examined in solution and a variety of membrane environments. It is shown that pore forming peptides have a strong preference for regions of high membrane curvature and that the structure adopted is dependent on the degree of local curvature. The simulations further suggest that it is a correspondence between curvature of the peptide and the intrinsic curvature of the membrane that could account for the effect of variation in head group and tail length on antimicrobial activity.[1] Leontiadou, H, Mark, AE and Marrink, SJ Antimicrobial peptides in action. J. Am. Chem. Soc. 2006, 128, 12156-12161.[2] Sengupta, D., Leontiadou, H., Mark, A. E. and Marrink, S. J. Toroidal Pores Formed by Antimicrobial Peptides Show Significant Disorder. Biochim. Biophys. Acta - Biomembranes. 2008, 1778, 2308-2317.[3] Yesylevskyy, S., Marrink, S. J. and Mark, A. E. Alternative mechanisms for the interaction of the cell-penetrating peptides Penetratin and the TAT peptide with lipid bilayers. Biophys. J. 2009, 97, 40-49.[4] Chen, R., Mark, A. E. The effect of membrane curvature on the conformation of antimicrobial peptides: implications for binding and the mechanism of action. Eur. Biophys. J., 2011, 40, 545-553.