1
|
Sahoo AR, Souza PCT, Meng Z, Buck M. Transmembrane dimers of type 1 receptors sample alternate configurations: MD simulations using coarse grain Martini 3 versus AlphaFold2 Multimer. Structure 2023; 31:735-745.e2. [PMID: 37075749 PMCID: PMC10833135 DOI: 10.1016/j.str.2023.03.014] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2022] [Revised: 02/07/2023] [Accepted: 03/23/2023] [Indexed: 04/21/2023]
Abstract
Structures and dynamics of transmembrane (TM) receptor regions are key to understanding their signaling mechanism across membranes. Here we examine configurations of TM region dimers, assembled using the recent Martini 3 force field for coarse-grain (CG) molecular dynamics simulations. At first glance, our results show only a reasonable agreement with ab initio predictions using PREDDIMER and AlphaFold2 Multimer and with nuclear magnetic resonance (NMR)-derived structures. 5 of 11 CG TM structures are similar to the NMR structures (within <3.5 Å root-mean-square deviation [RMSD]) compared with 10 and 9 using PREDDIMER and AlphaFold2, respectively (with 8 structures of the later within 1.5 Å). Surprisingly, AlphaFold2 predictions are closer to NMR structures when the 2001 instead of 2020 database is used for training. The CG simulations reveal that alternative configurations of TM dimers readily interconvert with a predominant population. The implications for transmembrane signaling are discussed, including for the development of peptide-based pharmaceuticals.
Collapse
Affiliation(s)
- Amita R Sahoo
- Department of Physiology and Biophysics, Case Western Reserve University, School of Medicine, 10900 Euclid Avenue, Cleveland, OH 44106, USA
| | - Paulo C T Souza
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS & University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| | - Zhiyuan Meng
- Department of Physiology and Biophysics, Case Western Reserve University, School of Medicine, 10900 Euclid Avenue, Cleveland, OH 44106, USA
| | - Matthias Buck
- Department of Physiology and Biophysics, Case Western Reserve University, School of Medicine, 10900 Euclid Avenue, Cleveland, OH 44106, USA.
| |
Collapse
|
2
|
Padhi S, Ramakrishna S, Priyakumar UD. Prediction of the structures of helical membrane proteins based on a minimum unfavorable contacts approach. J Comput Chem 2015; 36:539-52. [PMID: 25565454 DOI: 10.1002/jcc.23828] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2014] [Revised: 11/20/2014] [Accepted: 11/21/2014] [Indexed: 11/12/2022]
Abstract
An understanding of structure-function relationships of membrane proteins continues to be a challenging problem, owing to the difficulty in obtaining their structures experimentally. This study suggests a method for modeling membrane protein structures that can be used to generate a reliable initial conformation prior to the use of other approaches for sampling conformations. It involves optimizing the orientation of hydrophilic residues so as to minimize unfavorable contacts with the hydrophobic tails of the lipid bilayer. Starting with the optimized initial conformation for three different proteins modeled based on this method, two independent approaches have been used for sampling the conformational space of the proteins. Both approaches are able to predict structures reasonably close to experimental structures, indicating that the initial structure enables the sampling of conformations that are close to the native structure. Possible improvements in the method for making it broadly applicable to helical membrane proteins are discussed.
Collapse
Affiliation(s)
- Siladitya Padhi
- Centre for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad, 500032, India
| | | | | |
Collapse
|
3
|
Critical glycosylated residues in exon three of erythrocyte glycophorin A engage Plasmodium falciparum EBA-175 and define receptor specificity. mBio 2014; 5:e01606-14. [PMID: 25205096 PMCID: PMC4173783 DOI: 10.1128/mbio.01606-14] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023] Open
Abstract
Erythrocyte invasion is an essential step in the pathogenesis of malaria. The erythrocyte binding-like (EBL) family of Plasmodium falciparum proteins recognizes glycophorins (Gp) on erythrocytes and plays a critical role in attachment during invasion. However, the molecular basis for specific receptor recognition by each parasite ligand has remained elusive, as is the case with the ligand/receptor pair P. falciparum EBA-175 (PfEBA-175)/GpA. This is due largely to difficulties in producing properly glycosylated and functional receptors. Here, we developed an expression system to produce recombinant glycosylated and functional GpA, as well as mutations and truncations. We identified the essential binding region and determinants for PfEBA-175 engagement, demonstrated that these determinants are required for the inhibition of parasite growth, and identified the glycans important in mediating the PfEBA-175–GpA interaction. The results suggest that PfEBA-175 engages multiple glycans of GpA encoded by exon 3 and that the presentation of glycans is likely required for high-avidity binding. The absence of exon 3 in GpB and GpE due to a splice site mutation confers specific recognition of GpA by PfEBA-175. We speculate that GpB and GpE may have arisen due to selective pressure to lose the PfEBA-175 binding site in GpA. The expression system described here has wider application for examining other EBL members important in parasite invasion, as well as additional pathogens that recognize glycophorins. The ability to define critical binding determinants in receptor-ligand interactions, as well as a system to genetically manipulate glycosylated receptors, opens new avenues for the design of interventions that disrupt parasite invasion. Plasmodium falciparum uses distinct ligands that bind host cell receptors for invasion of red blood cells (RBCs) during malaria infection. A key entry pathway involves P. falciparum EBA-175 (PfEBA-175) recognizing glycophorin A (GpA) on RBCs. Despite knowledge of this protein-protein interaction, the complete mechanism for specific receptor engagement is not known. PfEBA-175 recognizes GpA but is unable to engage the related RBC receptor GpB or GpE. Understanding the necessary elements that enable PfEBA-175 to specifically recognize GpA is critical in developing specific and potent inhibitors of PfEBA-175 that disrupt host cell invasion and aid in malaria control. Here, we describe a novel system to produce and manipulate the host receptor GpA. Using this system, we probed the elements in GpA necessary for engagement and thus for host cell invasion. These studies have important implications for understanding how ligands and receptors interact and for the future development of malaria interventions.
Collapse
|
4
|
Yuzlenko O, Lazaridis T. Interactions between ionizable amino acid side chains at a lipid bilayer-water interface. J Phys Chem B 2011; 115:13674-84. [PMID: 21985663 DOI: 10.1021/jp2052213] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
Potentials of mean force (PMF) between ionizable amino acid side chains (Arg, Lys, His, Glu) in the headgroup area of a palmitoyl oleoyl phosphatidylcholine lipid bilayer were obtained from all-atom molecular dynamics simulations and the adaptive biasing force method. Simulations in bulk water were also performed for comparison. Side chains were constrained in collinear, stacking, and orthogonal (T-shaped) orientations. The most structured and attractive PMFs were observed for hydrogen-bonded side chains. Contact minima occurred at a distance of 2.6-3.1 Å between selected atoms or centers of mass with the most attractive interaction (-9.6 kcal/mol) observed between Arg(+) and Glu(-). Hydrogen bonds play a significant role in stabilizing these interactions. Interactions between like charged side chains can also be very attractive if the charges are screened by surrounding molecules or groups (e.g., the PMF value at the contact minimum for Arg(+)···Arg(+) is -7.6 kcal/mol). Like charged side chains can have contact minima as close as 3.6 Å. The PMFs depend strongly on the relative orientation of the side chains. In agreement with experimental studies and other simulations, we found the stacking arrangement of like charged side chains to be the most favorable orientation. Interaction energies and Lennard-Jones energies between side chains, headgroups, and water molecules were analyzed in order to rationalize the observed PMFs and their dependence on orientation. In general, the results cannot be explained by simple dielectric arguments.
Collapse
Affiliation(s)
- Olga Yuzlenko
- Department of Chemistry, City College of the City University of New York, 160 Convent Avenue, New York, New York 10031, USA
| | | |
Collapse
|
5
|
Frank A, Andricioaei I. A comparative study on the ability of two implicit solvent lipid models to predict transmembrane helix tilt angles. J Membr Biol 2010; 239:57-62. [PMID: 21152910 PMCID: PMC3030950 DOI: 10.1007/s00232-010-9325-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2010] [Accepted: 11/05/2010] [Indexed: 11/30/2022]
Abstract
Free-energy profiles describing the relative orientation of membrane proteins along predefined coordinates can be efficiently calculated by means of umbrella simulations. Such simulations generate reliable orientational distributions but are difficult to converge because of the very long equilibration times of the solvent and the lipid bilayer in explicit representation. Two implicit lipid membrane models are here applied in combination with the umbrella sampling strategy to the simulation of the transmembrane (TM) helical segment from virus protein U (Vpu). The models are used to study both orientation and energetics of this α-helical peptide as a function of hydrophobic mismatch. We observe that increasing the degree of positive hydrophobic mismatch increased the tilt angle of Vpu. These findings agree well with experimental data and as such validate the solvation models used in this study.
Collapse
Affiliation(s)
- Aaron Frank
- Department of Chemistry, University of California, Irvine, CA 92697, USA
| | | |
Collapse
|
6
|
Janosi L, Prakash A, Doxastakis M. Lipid-modulated sequence-specific association of glycophorin A in membranes. Biophys J 2010; 99:284-92. [PMID: 20655857 DOI: 10.1016/j.bpj.2010.04.005] [Citation(s) in RCA: 49] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2010] [Revised: 04/01/2010] [Accepted: 04/05/2010] [Indexed: 11/26/2022] Open
Abstract
Protein association in lipid membranes is a complex process with thermodynamics directed by a multitude of different factors. Amino-acid sequence is a molecular parameter that affects dimerization as shown by limited directed mutations along the transmembrane domains. Membrane-mediated interactions are also important although details of such contributions remain largely unclear. In this study, we probe directly the free energy of association of Glycophorin A by means of extensive parallel Monte Carlo simulations with recently developed methods and a model that accounts for sequence-specificity while representing lipid membranes faithfully. We find that lipid-induced interactions are significant both at short and intermediate separations. The ability of molecules to tilt in a specific hydrophobic environment extends their accessible interfaces, leading to intermittent contacts during protein recognition. The dimer with the lowest free energy is largely determined by the favorable lipid-induced attractive interactions at the closest distance. Finally, the coarse-grained model employed herein, together with the extensive sampling performed, provides estimates of the free energy of association that are in excellent agreement with existing data.
Collapse
Affiliation(s)
- Lorant Janosi
- Department of Chemical and Biomolecular Engineering, University of Houston, Houston, Texas, USA
| | | | | |
Collapse
|
7
|
Varriale S, Merlino A, Coscia MR, Mazzarella L, Oreste U. An evolutionary conserved motif is responsible for immunoglobulin heavy chain packing in the B cell membrane. Mol Phylogenet Evol 2010; 57:1238-44. [PMID: 20937398 DOI: 10.1016/j.ympev.2010.09.022] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2010] [Revised: 09/27/2010] [Accepted: 09/28/2010] [Indexed: 11/18/2022]
Abstract
All species of vertebrates synthesize immunoglobulin molecules, which differ in an number of aspects but also share a few common features responsible for their function, such as the presence of a transmembrane domain in the membrane bound form of the immunoglobulin heavy chain (IgTMD) that ensures communication with the signal transducing Igα-Igβ peptides. We have analyzed the gene sequence encoding the IgTMD of different heavy chain isotypes of very distant species, from shark to mammals. The IgTMD sequences show a high degree of sequence identity and their encoding nucleotide sequences were shown to be subject to purifying selection at most sites. We have built molecular models of seven IgTMDs from different vertebrate species and have investigated the formation of homodimer in a palmitoyl oleoyl phosphatidylcholine (POPC) lipid bilayer by molecular dynamics simulations. We found that the conserved FXXXFXXS/TXXXS motif, never observed to date in protein transmembrane chains, is responsible for the two heavy chains association through two pairs of Phe-Phe hydrophobic interactions and two pairs of Ser/Thr-Ser/Ser hydrogen bonds. This interaction pattern, which stabilizes the dimer conformation in the lipid bilayer, was unique, being different from any other pattern identified in transmembrane helices to date.
Collapse
Affiliation(s)
- Sonia Varriale
- Institute of Protein Biochemistry, CNR, via P. Castellino 111, 80131 Napoli, Italy
| | | | | | | | | |
Collapse
|
8
|
Chugunov AO, Efremov RG. [Prediction of the spatial structure of proteins: emphasis on membrane targets]. RUSSIAN JOURNAL OF BIOORGANIC CHEMISTRY 2010; 35:744-60. [PMID: 20208575 DOI: 10.1134/s106816200906003x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Knowledge of the spatial structure of proteins is a prerequisite for both awareness of their functional mechanisms and the framework for rational drug discovery and design. Meanwhile, direct structural determination is often hampered or impractical due to the complexity, expensiveness, and limited capabilities of experimental techniques. These issues are especially pronounced for integral membrane proteins. On numerous occasions, the theoretical prediction of protein structures may facilitate the process by exploiting physical or empirical principles. This paper surveys modern techniques for the prediction of the spatial structure of proteins using computer algorithms, and the main emphasis is placed on the most "complex" targets - membrane proteins (MPs). The first part of the review describes de novo methods based on empirical physical principles; in the second part, a comparative modeling philosophy, which accounts for the structure of related proteins, is described. Special focus is made regarding pharmacologically relevant classes of G-coupled receptors, receptor tyrosine ki-nases, and other MPs. Algorithms for the assessment of the models quality and potential fields of application of computer models are discussed.
Collapse
|
9
|
Jusoh SA, Welsch C, Siu SWI, Böckmann RA, Helms V. Contribution of charged and polar residues for the formation of the E1-E2 heterodimer from Hepatitis C Virus. J Mol Model 2010; 16:1625-37. [PMID: 20195665 DOI: 10.1007/s00894-010-0672-1] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2009] [Accepted: 01/25/2010] [Indexed: 02/05/2023]
Abstract
The transmembrane domains of the envelope glycoprotein E1 and E2 have crucial multifunctional roles in the biogenesis of hepatitis C virus. We have performed molecular dynamics simulations to investigate a structural model of the transmembrane segments of the E1-E2 heterodimer. The simulations support the key role of the Lys370-Asp728 ion pair for mediating the E1-E2 heterodimerization. In comparison to these two residues, the simulation results also reveal the differential effect of the conserved Arg730 residue that has been observed in experimental studies. Furthermore, we discovered the formation of inter-helical hydrogen bonds via Asn367 that stabilize dimer formation. Simulations of single and double mutants further demonstrate the importance of the ion-pair and polar interactions between the interacting helix monomers. The conformation of the E1 fragment in the simulation of the E1-E2 heterodimer is in close agreement with an NMR structure of the E1 transmembrane segment. The proposed model of the E1-E2 heterodimer supports the postulated cooperative insertion of both helices by the translocon complex into the bilayer.
Collapse
Affiliation(s)
- Siti Azma Jusoh
- Center for Bioinformatics, Saarland University, 66041, Saarbruecken, Germany
| | | | | | | | | |
Collapse
|
10
|
Sengupta D, Marrink SJ. Lipid-mediated interactions tune the association of glycophorin A helix and its disruptive mutants in membranes. Phys Chem Chem Phys 2010; 12:12987-96. [DOI: 10.1039/c0cp00101e] [Citation(s) in RCA: 97] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
11
|
Zhang J, Lazaridis T. Transmembrane helix association affinity can be modulated by flanking and noninterfacial residues. Biophys J 2009; 96:4418-27. [PMID: 19486666 DOI: 10.1016/j.bpj.2009.03.008] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2008] [Revised: 02/13/2009] [Accepted: 03/03/2009] [Indexed: 11/24/2022] Open
Abstract
The GxxxG sequence motif mediates the association of transmembrane (TM) helices by providing a site of close contact between them. However, it is not sufficient for strong association. For example, both bacteriophage M13 major coat protein (MCP) and human erythrocyte protein glycophorin A (GpA) contain a GxxxG motif in their TM domains and form a homodimer, but the association affinity of MCP, measured by the ToxCAT in vivo assay, is dramatically weaker than that of GpA. Even when all interfacial residues of MCP were substituted for those of GpA (MCP-GpA), association remained significantly weaker than in GpA. Here we provide an explanation for these experimental observations using molecular dynamics simulations in an implicit membrane (IMM1-GC). The association free energies of GpA29 (GpA with 29 residues all from the wild-type sequence), GpA15p11 (GpA with 15 residues from the wild-type sequence plus 11 flanking residues from the ToxCAT construct), MCP, and MCP-GpA TM helices were calculated and compared. MCP and MCP-GpA have the same flanking residues used in the ToxCAT assay as those in GpA15p11, but the position of the flanking residues relative to the GxxxG motif is different. The calculated association free energies follow experimental observations: the association affinity of MCP-GpA falls between those of GpA15p11 and MCP wild-type. MCP exhibits an equally strong interhelical interaction in the TM domain. A major reason for the weaker association of MCP in the calculations was the noninterfacial residue Lys-40, which in the dimer structure is forced to be buried in the membrane interior. To alleviate the desolvation cost, in MCP and MCP-GpA dimers, Lys-40 gets deprotonated. A second factor that modulates association affinity is the flanking residues. Thanks to them, GpA15p11 exhibits a much stronger association affinity than GpA29. The positioning of the flanking residues is also important, as evidenced by the difference in association affinity between MCP and MCP-GpA on one hand and GpA15p11 on the other. Thus, residues outside the contact interface can exert a significant influence on transmembrane helix association affinity.
Collapse
Affiliation(s)
- Jinming Zhang
- Department of Chemistry, City College of New York/CUNY, New York, NY 10031, USA
| | | |
Collapse
|
12
|
Sayyed-Ahmad A, Khandelia H, Kaznessis YN. Relative free energy of binding between antimicrobial peptides and SDS or DPC micelles. MOLECULAR SIMULATION 2009; 35:986-997. [PMID: 21113423 DOI: 10.1080/08927020902902742] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Abstract
We present relative binding free energy calculations for six antimicrobial peptide-micelle systems, three peptides interacting with two types of micelles. The peptides are the scorpion derived antimicrobial peptide (AMP), IsCT and two of its analogues. The micelles are dodecylphosphatidylcholine (DPC) and sodium dodecylsulphate (SDS) micelles. The interfacial electrostatic properties of DPC and SDS micelles are assumed to be similar to those of zwitterionic mammalian and anionic bacterial membrane interfaces, respectively. We test the hypothesis that the binding strength between peptides and the anionic micelle SDS can provide information on peptide antimicrobial activity, since it is widely accepted that AMPs function by binding to and disrupting the predominantly anionic lipid bilayer of the bacterial cytoplasmic membrane. We also test the hypothesis that the binding strength between peptides and the zwitterionic micelle DPC can provide information on peptide haemolytic activities, since it is accepted that they also bind to and disrupt the zwitterionic membrane of mammalian cells. Equilibrium structures of the peptides, micelles and peptide-micelle complexes are obtained from more than 300 ns of molecular dynamics simulations. A thermodynamic cycle is introduced to compute the binding free energy from electrostatic, non-electrostatic and entropic contributions. We find relative binding free energy strengths between peptides and SDS to correlate with the experimentally measured rankings for peptide antimicrobial activities, and relative free energy binding strengths between peptides and DPC to correlate with the observed rankings for peptide haemolytic toxicities. These findings point to the importance of peptide-membrane binding strength for antimicrobial activity and haemolytic activity.
Collapse
Affiliation(s)
- Abdallah Sayyed-Ahmad
- Department of Chemical Engineering and Materials Science, and the Digital Technology Center, University of Minnesota, 421 Washington Ave. SE, Minneapolis, MN 55455, USA
| | | | | |
Collapse
|
13
|
Affiliation(s)
- J. Krüger
- Institute of Biophotonics, School of Medical Science and Engineering, National Yang-Ming University, 155, Section 2, Li-Nong Street, Taipei 112, Taiwan
| | - W. B. Fischer
- Institute of Biophotonics, School of Medical Science and Engineering, National Yang-Ming University, 155, Section 2, Li-Nong Street, Taipei 112, Taiwan
| |
Collapse
|
14
|
Brooks B, Brooks C, MacKerell A, Nilsson L, Petrella R, Roux B, Won Y, Archontis G, Bartels C, Boresch S, Caflisch A, Caves L, Cui Q, Dinner A, Feig M, Fischer S, Gao J, Hodoscek M, Im W, Kuczera K, Lazaridis T, Ma J, Ovchinnikov V, Paci E, Pastor R, Post C, Pu J, Schaefer M, Tidor B, Venable RM, Woodcock HL, Wu X, Yang W, York D, Karplus M. CHARMM: the biomolecular simulation program. J Comput Chem 2009; 30:1545-614. [PMID: 19444816 PMCID: PMC2810661 DOI: 10.1002/jcc.21287] [Citation(s) in RCA: 6075] [Impact Index Per Article: 405.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecular simulation program. It has been developed over the last three decades with a primary focus on molecules of biological interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estimators, molecular minimization, dynamics, and analysis techniques, and model-building capabilities. The CHARMM program is applicable to problems involving a much broader class of many-particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical-molecular mechanical force fields, to all-atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numerous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.
Collapse
Affiliation(s)
- B.R. Brooks
- Laboratory of Computational Biology, National Heart, Lung, and
Blood Institute, National Institutes of Health, Bethesda, MD 20892
| | - C.L. Brooks
- Departments of Chemistry & Biophysics, University of
Michigan, Ann Arbor, MI 48109
| | - A.D. MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy,
University of Maryland, Baltimore, MD, 21201
| | - L. Nilsson
- Karolinska Institutet, Department of Biosciences and Nutrition,
SE-141 57, Huddinge, Sweden
| | - R.J. Petrella
- Department of Chemistry and Chemical Biology, Harvard University,
Cambridge, MA 02138
- Department of Medicine, Harvard Medical School, Boston, MA
02115
| | - B. Roux
- Department of Biochemistry and Molecular Biology, University of
Chicago, Gordon Center for Integrative Science, Chicago, IL 60637
| | - Y. Won
- Department of Chemistry, Hanyang University, Seoul
133–792 Korea
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | - M. Karplus
- Department of Chemistry and Chemical Biology, Harvard University,
Cambridge, MA 02138
- Laboratoire de Chimie Biophysique, ISIS, Université de
Strasbourg, 67000 Strasbourg France
| |
Collapse
|
15
|
Vereshaga YA, Volynsky PE, Pustovalova JE, Nolde DE, Arseniev AS, Efremov RG. Specificity of helix packing in transmembrane dimer of the cell death factor BNIP3: a molecular modeling study. Proteins 2009; 69:309-25. [PMID: 17600828 DOI: 10.1002/prot.21555] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
BNIP3 is a mitochondrial 19-kDa proapoptotic protein, a member of the Bcl-2 family. It has a single COOH-terminal transmembrane (TM) alpha-helical domain, which is required for membrane targeting, proapoptotic activity, hetero- and homo-dimerization in membrane. The role and the molecular details of association of TM helices of BNIP3 are yet to be established. Here, we present a molecular modeling study of helix interactions in its membrane domain. The approach combines Monte Carlo conformational search in an implicit hydrophobic slab followed by molecular dynamics simulations in a hydrated full-atom lipid bilayer. The former technique was used for exhaustive sampling of the peptides' conformational space and for generation of putative "native-like" structures of the dimer. The latter ones were taken as realistic starting points to assess stability and dynamic behavior of the complex in explicit lipid-water surrounding. As a result, several groups of tightly packed right-handed structures of the dimer were proposed. They have almost similar helix-helix interface, which includes the motif A(176)xxxG(180)xxxG(184) and agrees well with previous mutagenesis data and preliminary NMR analysis. Molecular dynamics simulations of these structures reveal perfect adaptation of most of them to heterogeneous membrane environment. A remarkable feature of the predicted dimeric structures is the occurrence of a cluster of H-bonded histidine 173 and serines 168 and 172 on the helix interface, near the N-terminus. Because of specific polar interactions between the monomers, this part of the dimer has no such dense packing as the C-terminal one, thus allowing penetration of water from the extramembrane side into the membrane interior. We propose that the ionization state of His(173) can mediate structural and dynamic properties of the dimer. This, in turn, may be related to pH-dependent proapoptotic activity of BNIP3, which is triggering on by acidosis appearing under hypoxic conditions.
Collapse
Affiliation(s)
- Yana A Vereshaga
- M.M. Shemyakin and Yu.A. Ovchinnikov Institute of Bioorganic Chemistry, Russian Academy of Sciences, Moscow V-437, 117997 GSP, Russia
| | | | | | | | | | | |
Collapse
|
16
|
Gervais C, Wüst T, Landau DP, Xu Y. Application of the Wang-Landau algorithm to the dimerization of glycophorin A. J Chem Phys 2009; 130:215106. [PMID: 19508105 PMCID: PMC2719476 DOI: 10.1063/1.3148186] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2009] [Accepted: 05/13/2009] [Indexed: 11/15/2022] Open
Abstract
A two-step Monte Carlo procedure is developed to investigate the dimerization process of the homodimer glycophorin A. In the first step, the energy density of states of the system is estimated by the Wang-Landau algorithm. In the second step, a production run is performed during which various energetical and structural observables are sampled to provide insight into the thermodynamics of the system. All seven residues LIxxGVxxGVxxT constituting the contact interface play a dominating role in the dimerization, however at different stages of the process. The leucine motif and to some extent the GxxxG motif are involved at the very beginning of the dimerization when the two helices come into contact, ensuring an interface already similar to the native one. At a lower temperature, the threonine motif stabilizes by hydrogen bonding the dimer, which finally converges toward its native state at around 300 K. The power and flexibility of the procedure employed here makes it an interesting alternative to other Monte Carlo methods for the study of similar protein systems.
Collapse
Affiliation(s)
- Claire Gervais
- Department of Biochemistry and Molecular Biology, Computational Systems Biology Laboratory and Institute of Bioinformatics, The University of Georgia, Athens, Georgia 30602, USA.
| | | | | | | |
Collapse
|
17
|
Pendley SS, Yu YB, Cheatham TE. Molecular dynamics guided study of salt bridge length dependence in both fluorinated and non-fluorinated parallel dimeric coiled-coils. Proteins 2009; 74:612-29. [PMID: 18704948 PMCID: PMC2692595 DOI: 10.1002/prot.22177] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
The alpha-helical coiled-coil is one of the most common oligomerization motifs found in both native and engineered proteins. To better understand the stability and dynamics of the coiled-coil motifs, including those modified by fluorination, several fluorinated and nonfluorinated parallel dimeric coiled-coil protein structures were designed and modeled. We also attempt to investigate how changing the length and geometry of the important stabilizing salt bridges influences the coiled-coil protein structure. Molecular dynamics (MD) and free energy simulations with AMBER used a particle mesh Ewald treatment of the electrostatics in explicit TIP3P solvent with balanced force field treatments. Preliminary studies with legacy force fields (ff94, ff96, and ff99) show a profound instability of the coiled-coil structures in short MD simulation. Significantly, better behavior is evident with the more balanced ff99SB and ff03 protein force fields. Overall, the results suggest that the coiled-coil structures can readily accommodate the larger acidic arginine or S-2,7-diaminoheptanedoic acid mutants in the salt bridge, whereas substitution of the smaller L-ornithine residue leads to rapid disruption of the coiled-coil structure on the MD simulation time scale. This structural distortion of the secondary structure allows both the formation of large hydration pockets proximal to the charged groups and within the hydrophobic core. Moreover, the increased structural fluctuations and movement lead to a decrease in the water occupancy lifetimes in the hydration pockets. In contrast, analysis of the hydration in the stable dimeric coiled-coils shows high occupancy water sites along the backbone residues with no water occupancy in the hydrophobic core, although transitory water interactions with the salt bridge residues are evident. The simulations of the fluorinated coiled-coils suggest that in some cases fluorination electrostatically stabilizes the intermolecular coiled-coil salt bridges. Structural analyses also reveal different side chain rotamer preferences for leucine when compared with 5,5,5,5',5',5'-hexafluoroleucine mutants. These observed differences in the side chain rotamer populations suggest differential changes in the side chain conformational entropy upon coiled-coil formation when the protein is fluorinated. The free energy of hydration of the isolated 5,5,5,5',5',5'-hexafluoroleucine amino acid is calculated to be 1.1 kcal/mol less stable than leucine; this hydrophobic penalty in the monomer may provide a driving force for coiled-coil dimer formation. Estimation of the ellipticity at 222 nm from a series of snapshots from the MD simulations with DicroCalc shows distinct increases in the ellipticity when the coiled-coil is fluorinated, which suggests that the helicity in the folded coiled-coils is greater when fluorinated.
Collapse
Affiliation(s)
- Scott S. Pendley
- Departments of Pharmaceutics and Pharmaceutical Chemistry, University of Utah, 2000 South 30 East, Skaggs Hall 201, Salt Lake City, UT 84112
| | - Yihua B. Yu
- Departments of Pharmaceutics and Pharmaceutical Chemistry, University of Utah, 2000 South 30 East, Skaggs Hall 201, Salt Lake City, UT 84112
- Departments of Pharmaceutical Sciences and Bioengineering, University of Maryland, University of Maryland, 20 Penn Street, Rm. 635, Baltimore, MD 21201
| | - Thomas E. Cheatham
- Departments of Pharmaceutics and Pharmaceutical Chemistry, University of Utah, 2000 South 30 East, Skaggs Hall 201, Salt Lake City, UT 84112
- Department of Medicinal Chemistry, University of Utah, 2000 South 30 East, Skaggs Hall 201, Salt Lake City, UT 84112
- Department of Bioengineering, University of Utah, 2000 South 30 East, Skaggs Hall 201, Salt Lake City, UT 84112
| |
Collapse
|
18
|
Structure and dimerization of the teleost transmembrane immunoglobulin region. J Mol Graph Model 2008; 27:401-7. [DOI: 10.1016/j.jmgm.2008.07.001] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2008] [Revised: 07/08/2008] [Accepted: 07/10/2008] [Indexed: 12/20/2022]
|
19
|
Kao L, Sassani P, Azimov R, Pushkin A, Abuladze N, Peti-Peterdi J, Liu W, Newman D, Kurtz I. Oligomeric structure and minimal functional unit of the electrogenic sodium bicarbonate cotransporter NBCe1-A. J Biol Chem 2008; 283:26782-94. [PMID: 18658147 DOI: 10.1074/jbc.m804006200] [Citation(s) in RCA: 52] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
The electrogenic sodium bicarbonate cotransporter NBCe1-A mediates the basolateral absorption of sodium and bicarbonate in the proximal tubule. In this study the oligomeric state and minimal functional unit of NBCe1-A were investigated. Wild-type (wt) NBCe1-A isolated from mouse kidney or heterologously expressed in HEK293 cells was predominantly in a dimeric state as was shown using fluorescence energy transfer, pulldown, immunoprecipitation, cross-linking experiments, and nondenaturing perfluorooctanoate-PAGE. NBCe1-A monomers were found to be covalently linked by S-S bonds. When each of the 15 native cysteine residues were individually removed on a wt-NBCe1-A backbone, dimerization of the cotransporter was not affected. In experiments involving multiple native cysteine residue removal, both Cys(630) and Cys(642) in extracellular loop 3 were shown to mediate S-S bond formation between NBCe1-A monomers. When native NBCe1-A cysteine residues were individually reintroduced into a cysteineless NBCe1-A mutant backbone, the finding that a Cys(992) construct that lacked S-S bonds functioned normally indicated that stable covalent linkage of NBCe1-A monomers was not a necessary requirement for functional activity of the cotransporter. Studies using concatameric constructs of wt-NBCe1-A, whose activity is resistant to methanesulfonate reagents, and an NBCe1-A(T442C) mutant, whose activity is completely inhibited by methanesulfonate reagents, confirmed that NBCe1-A monomers are functional. Our results demonstrate that wt-NBCe1-A is predominantly a homodimer, dependent on S-S bond formation that is composed of functionally active monomers.
Collapse
Affiliation(s)
- Liyo Kao
- Division of Nephrology, David Geffen School Medicine, UCLA, Los Angeles, California 90095-1689, USA
| | | | | | | | | | | | | | | | | |
Collapse
|
20
|
Dell'Orco D, De Benedetti PG, Fanelli F. In Silico Screening of Mutational Effects on Transmembrane Helix Dimerization: Insights from Rigid-Body Docking and Molecular Dynamics Simulations. J Phys Chem B 2007; 111:9114-24. [PMID: 17602582 DOI: 10.1021/jp071383r] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In this study, a docking-based protocol has been probed for its ability to predict the effects of 32 single and double mutations on glycophorin A (GpA) homodimerization. Rigid-body docking simulations have been paralleled by molecular dynamics (MD) simulations in implicit membrane. The rigid-body docking-based approach proved effective in reconstituting the native architecture of the GpA dimer for the wild type and the wild-type-like mutants. The good correlative models between the intermolecular interaction descriptors derived both by rigid-body docking simulations and by MD simulations and experimental relative free energies support the assumption that the mutation-induced changes in the association free energy of GpA helices are essentially ascribed to differences in the packing interactions, whereas almost all the variations in the entropic contributions to the association free energy are constant and/or negligible. The MD-based models achieved provide insights into the structural determinants for disruptive and restoring mutational effects. The computational approaches presented in this study are fast and effective, and constitute simple and promising tools for in silico screening of mutational effects on the association properties of integral membrane proteins.
Collapse
Affiliation(s)
- Daniele Dell'Orco
- Department of Chemistry and Dulbecco Telethon Institute, University of Modena and Reggio Emilia, via Campi 183, 41100 Modena, Italy
| | | | | |
Collapse
|
21
|
Johansson ACV, Lindahl E. Amino-acid solvation structure in transmembrane helices from molecular dynamics simulations. Biophys J 2006; 91:4450-63. [PMID: 17012325 PMCID: PMC1779925 DOI: 10.1529/biophysj.106.092767] [Citation(s) in RCA: 67] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Understanding the solvation of amino acids in biomembranes is an important step to better explain membrane protein folding. Several experimental studies have shown that polar residues are both common and important in transmembrane segments, which means they have to be solvated in the hydrophobic membrane, at least until helices have aggregated to form integral proteins. In this work, we have used computer simulations to unravel these interactions on the atomic level, and classify intramembrane solvation properties of amino acids. Simulations have been performed for systematic mutations in poly-Leu helices, including not only each amino acid type, but also every z-position in a model helix. Interestingly, many polar or charged residues do not desolvate completely, but rather retain hydration by snorkeling or pulling in water/headgroups--even to the extent where many of them exist in a microscopic polar environment, with hydration levels corresponding well to experimental hydrophobicity scales. This suggests that even for polar/charged residues a large part of solvation cost is due to entropy, not enthalpy loss. Both hydration level and hydrogen bonding exhibit clear position-dependence. Basic side chains cause much less membrane distortion than acidic, since they are able to form hydrogen bonds with carbonyl groups instead of water or headgroups. This preference is supported by sequence statistics, where basic residues have increased relative occurrence at carbonyl z-coordinates. Snorkeling effects and N-/C-terminal orientation bias are directly observed, which significantly reduces the effective thickness of the hydrophobic core. Aromatic side chains intercalate efficiently with lipid chains (improving Trp/Tyr anchoring to the interface) and Ser/Thr residues are stabilized by hydroxyl groups sharing hydrogen bonds to backbone oxygens.
Collapse
Affiliation(s)
- Anna C V Johansson
- Stockholm Bioinformatics Center, Stockholm University, Stockholm, Sweden
| | | |
Collapse
|
22
|
Fleishman SJ, Ben-Tal N. Progress in structure prediction of α-helical membrane proteins. Curr Opin Struct Biol 2006; 16:496-504. [PMID: 16822664 DOI: 10.1016/j.sbi.2006.06.003] [Citation(s) in RCA: 44] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2006] [Revised: 06/05/2006] [Accepted: 06/16/2006] [Indexed: 10/24/2022]
Abstract
Transmembrane (TM) proteins comprise 20-30% of the genome but, because of experimental difficulties, they represent less than 1% of the Protein Data Bank. The dearth of membrane protein structures makes computational prediction a potentially important means of obtaining novel structures. Recent advances in computational methods have been combined with experimental data to constrain the modeling of three-dimensional structures. Furthermore, threading and ab initio modeling approaches that were effective for soluble proteins have been applied to TM domains. Surprisingly, experimental structures, proteomic analyses and bioinformatics have revealed unexpected architectures that counter long-held views on TM protein structure and stability. Future computational and experimental studies aimed at understanding the thermodynamic and evolutionary bases of these architectural details will greatly enhance predictive capabilities.
Collapse
Affiliation(s)
- Sarel J Fleishman
- Department of Biochemistry, George S. Wise Faculty of Life Sciences, Tel-Aviv University Ramat Aviv 69978, Israel
| | | |
Collapse
|