1
|
Archaeal Lipids Regulating the Trimeric Structure Dynamics of Bacteriorhodopsin for Efficient Proton Release and Uptake. Int J Mol Sci 2022; 23:ijms23136913. [PMID: 35805918 PMCID: PMC9278134 DOI: 10.3390/ijms23136913] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2022] [Revised: 06/19/2022] [Accepted: 06/20/2022] [Indexed: 11/30/2022] Open
Abstract
S-TGA-1 and PGP-Me are native archaeal lipids associated with the bacteriorhodopsin (bR) trimer and contribute to protein stabilization and native dynamics for proton transfer. However, little is known about the underlying molecular mechanism of how these lipids regulate bR trimerization and efficient photocycling. Here, we explored the specific binding of S-TGA-1 and PGP-Me with the bR trimer and elucidated how specific interactions modulate the bR trimeric structure and proton release and uptake using long-term atomistic molecular dynamic simulations. Our results showed that S-TGA-1 and PGP-Me are essential for stabilizing the bR trimer and maintaining the coherent conformational dynamics necessary for proton transfer. The specific binding of S-TGA-1 with W80 and K129 regulates proton release on the extracellular surface by forming a “Glu-shared” model. The interaction of PGP-Me with K40 ensures proton uptake by accommodating the conformation of the helices to recruit enough water molecules on the cytoplasmic side. The present study results could fill in the theoretical gaps of studies on the functional role of archaeal lipids and could provide a reference for other membrane proteins containing similar archaeal lipids.
Collapse
|
2
|
Manukyan AK. Structural aspects and activation mechanism of human secretory group IIA phospholipase. EUROPEAN BIOPHYSICS JOURNAL: EBJ 2020; 49:511-531. [DOI: 10.1007/s00249-020-01458-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/21/2020] [Revised: 08/12/2020] [Accepted: 08/14/2020] [Indexed: 11/30/2022]
|
3
|
Good D, Pham C, Jagas J, Lewandowski JR, Ladizhansky V. Solid-State NMR Provides Evidence for Small-Amplitude Slow Domain Motions in a Multispanning Transmembrane α-Helical Protein. J Am Chem Soc 2017; 139:9246-9258. [PMID: 28613900 PMCID: PMC5510093 DOI: 10.1021/jacs.7b03974] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2017] [Indexed: 02/06/2023]
Abstract
Proteins are dynamic entities and populate ensembles of conformations. Transitions between states within a conformational ensemble occur over a broad spectrum of amplitude and time scales, and are often related to biological function. Whereas solid-state NMR (SSNMR) spectroscopy has recently been used to characterize conformational ensembles of proteins in the microcrystalline states, its applications to membrane proteins remain limited. Here we use SSNMR to study conformational dynamics of a seven-helical transmembrane (TM) protein, Anabaena Sensory Rhodopsin (ASR) reconstituted in lipids. We report on site-specific measurements of the 15N longitudinal R1 and rotating frame R1ρ relaxation rates at two fields of 600 and 800 MHz and at two temperatures of 7 and 30 °C. Quantitative analysis of the R1 and R1ρ values and of their field and temperature dependencies provides evidence of motions on at least two time scales. We modeled these motions as fast local motions and slower collective motions of TM helices and of structured loops, and used the simple model-free and extended model-free analyses to fit the data and estimate the amplitudes, time scales and activation energies. Faster picosecond (tens to hundreds of picoseconds) local motions occur throughout the protein and are dominant in the middle portions of the TM helices. In contrast, the amplitudes of the slower collective motions occurring on the nanosecond (tens to hundreds of nanoseconds) time scales, are smaller in the central parts of helices, but increase toward their cytoplasmic sides as well as in the interhelical loops. ASR interacts with a soluble transducer protein on its cytoplasmic surface, and its binding affinity is modulated by light. The larger amplitude of motions on the cytoplasmic side of the TM helices correlates with the ability of ASR to undergo large conformational changes in the process of binding/unbinding the transducer.
Collapse
Affiliation(s)
- Daryl Good
- Department
of Physics and Biophysics Interdepartmental Group, University
of Guelph, Guelph, Ontario N1G 2W1, Canada
| | - Charlie Pham
- Department
of Physics and Biophysics Interdepartmental Group, University
of Guelph, Guelph, Ontario N1G 2W1, Canada
| | - Jacob Jagas
- Department
of Physics and Biophysics Interdepartmental Group, University
of Guelph, Guelph, Ontario N1G 2W1, Canada
| | - Józef R. Lewandowski
- Department
of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Vladimir Ladizhansky
- Department
of Physics and Biophysics Interdepartmental Group, University
of Guelph, Guelph, Ontario N1G 2W1, Canada
| |
Collapse
|
4
|
Jang H, Muratcioglu S, Gursoy A, Keskin O, Nussinov R. Membrane-associated Ras dimers are isoform-specific: K-Ras dimers differ from H-Ras dimers. Biochem J 2016; 473:1719-32. [PMID: 27057007 PMCID: PMC7830773 DOI: 10.1042/bcj20160031] [Citation(s) in RCA: 80] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2016] [Accepted: 04/07/2016] [Indexed: 02/07/2023]
Abstract
Are the dimer structures of active Ras isoforms similar? This question is significant since Ras can activate its effectors as a monomer; however, as a dimer, it promotes Raf's activation and MAPK (mitogen-activated protein kinase) cell signalling. In the present study, we model possible catalytic domain dimer interfaces of membrane-anchored GTP-bound K-Ras4B and H-Ras, and compare their conformations. The active helical dimers formed by the allosteric lobe are isoform-specific: K-Ras4B-GTP favours the α3 and α4 interface; H-Ras-GTP favours α4 and α5. Both isoforms also populate a stable β-sheet dimer interface formed by the effector lobe; a less stable β-sandwich interface is sustained by salt bridges of the β-sheet side chains. Raf's high-affinity β-sheet interaction is promoted by the active helical interface. Collectively, Ras isoforms' dimer conformations are not uniform; instead, the isoform-specific dimers reflect the favoured interactions of the HVRs (hypervariable regions) with cell membrane microdomains, biasing the effector-binding site orientations, thus isoform binding selectivity.
Collapse
Affiliation(s)
- Hyunbum Jang
- Cancer and Inflammation Program, National Cancer Institute at Frederick, Basic Science Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research, Frederick, MD 21702, U.S.A
| | - Serena Muratcioglu
- Department of Chemical and Biological Engineering, Koc University, Istanbul, Turkey
| | - Attila Gursoy
- Department of Computer Engineering, Koc University, Istanbul, Turkey
| | - Ozlem Keskin
- Department of Chemical and Biological Engineering, Koc University, Istanbul, Turkey
| | - Ruth Nussinov
- Cancer and Inflammation Program, National Cancer Institute at Frederick, Basic Science Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research, Frederick, MD 21702, U.S.A. Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine, Tel Aviv University, Tel Aviv 69978, Israel
| |
Collapse
|
5
|
Jang H, Banerjee A, Chavan TS, Lu S, Zhang J, Gaponenko V, Nussinov R. The higher level of complexity of K-Ras4B activation at the membrane. FASEB J 2016; 30:1643-55. [PMID: 26718888 PMCID: PMC4799498 DOI: 10.1096/fj.15-279091] [Citation(s) in RCA: 69] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2015] [Accepted: 12/14/2015] [Indexed: 12/19/2022]
Abstract
Is nucleotide exchange sufficient to activate K-Ras4B? To signal, oncogenic rat sarcoma (Ras) anchors in the membrane and recruits effectors by exposing its effector lobe. With the use of NMR and molecular dynamics (MD) simulations, we observed that in solution, farnesylated guanosine 5'-diphosphate (GDP)-bound K-Ras4B is predominantly autoinhibited by its hypervariable region (HVR), whereas the GTP-bound state favors an activated, HVR-released state. On the anionic membrane, the catalytic domain adopts multiple orientations, including parallel (∼180°) and perpendicular (∼90°) alignments of the allosteric helices, with respect to the membrane surface direction. In the autoinhibited state, the HVR is sandwiched between the effector lobe and the membrane; in the active state, with membrane-anchored farnesyl and unrestrained HVR, the catalytic domain fluctuates reinlessly, exposing its effector-binding site. Dimerization and clustering can reduce the fluctuations. This achieves preorganized, productive conformations. Notably, we also observe HVR-autoinhibited K-Ras4B-GTP states, with GDP-bound-like orientations of the helices. Thus, we propose that the GDP/GTP exchange may not be sufficient for activation; instead, our results suggest that the GDP/GTP exchange, HVR sequestration, farnesyl insertion, and orientation/localization of the catalytic domain at the membrane conjointly determine the active or inactive state of K-Ras4B. Importantly, K-Ras4B-GTP can exist in active and inactive states; on its own, GTP binding may not compel K-Ras4B activation.-Jang, H., Banerjee, A., Chavan, T. S, Lu, S., Zhang, J., Gaponenko, V., Nussinov, R. The higher level of complexity of K-Ras4B activation at the membrane.
Collapse
Affiliation(s)
- Hyunbum Jang
- *Basic Science Program, Leidos Biomedical Research, Frederick National Laboratory for Cancer Research, Cancer and Inflammation Program, National Cancer Institute at Frederick, Frederick, Maryland, USA; Department of Chemistry, Department of Medicinal Chemistry, and Department of Biochemistry and Molecular Genetics, University of Illinois at Chicago, Chicago, Illinois, USA; Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis, Chinese Ministry of Education, Shanghai JiaoTong University, School of Medicine, Shanghai, China; and Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Avik Banerjee
- *Basic Science Program, Leidos Biomedical Research, Frederick National Laboratory for Cancer Research, Cancer and Inflammation Program, National Cancer Institute at Frederick, Frederick, Maryland, USA; Department of Chemistry, Department of Medicinal Chemistry, and Department of Biochemistry and Molecular Genetics, University of Illinois at Chicago, Chicago, Illinois, USA; Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis, Chinese Ministry of Education, Shanghai JiaoTong University, School of Medicine, Shanghai, China; and Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Tanmay S Chavan
- *Basic Science Program, Leidos Biomedical Research, Frederick National Laboratory for Cancer Research, Cancer and Inflammation Program, National Cancer Institute at Frederick, Frederick, Maryland, USA; Department of Chemistry, Department of Medicinal Chemistry, and Department of Biochemistry and Molecular Genetics, University of Illinois at Chicago, Chicago, Illinois, USA; Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis, Chinese Ministry of Education, Shanghai JiaoTong University, School of Medicine, Shanghai, China; and Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Shaoyong Lu
- *Basic Science Program, Leidos Biomedical Research, Frederick National Laboratory for Cancer Research, Cancer and Inflammation Program, National Cancer Institute at Frederick, Frederick, Maryland, USA; Department of Chemistry, Department of Medicinal Chemistry, and Department of Biochemistry and Molecular Genetics, University of Illinois at Chicago, Chicago, Illinois, USA; Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis, Chinese Ministry of Education, Shanghai JiaoTong University, School of Medicine, Shanghai, China; and Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Jian Zhang
- *Basic Science Program, Leidos Biomedical Research, Frederick National Laboratory for Cancer Research, Cancer and Inflammation Program, National Cancer Institute at Frederick, Frederick, Maryland, USA; Department of Chemistry, Department of Medicinal Chemistry, and Department of Biochemistry and Molecular Genetics, University of Illinois at Chicago, Chicago, Illinois, USA; Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis, Chinese Ministry of Education, Shanghai JiaoTong University, School of Medicine, Shanghai, China; and Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Vadim Gaponenko
- *Basic Science Program, Leidos Biomedical Research, Frederick National Laboratory for Cancer Research, Cancer and Inflammation Program, National Cancer Institute at Frederick, Frederick, Maryland, USA; Department of Chemistry, Department of Medicinal Chemistry, and Department of Biochemistry and Molecular Genetics, University of Illinois at Chicago, Chicago, Illinois, USA; Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis, Chinese Ministry of Education, Shanghai JiaoTong University, School of Medicine, Shanghai, China; and Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Ruth Nussinov
- *Basic Science Program, Leidos Biomedical Research, Frederick National Laboratory for Cancer Research, Cancer and Inflammation Program, National Cancer Institute at Frederick, Frederick, Maryland, USA; Department of Chemistry, Department of Medicinal Chemistry, and Department of Biochemistry and Molecular Genetics, University of Illinois at Chicago, Chicago, Illinois, USA; Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis, Chinese Ministry of Education, Shanghai JiaoTong University, School of Medicine, Shanghai, China; and Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
| |
Collapse
|
6
|
Jang H, Abraham SJ, Chavan TS, Hitchinson B, Khavrutskii L, Tarasova NI, Nussinov R, Gaponenko V. Mechanisms of membrane binding of small GTPase K-Ras4B farnesylated hypervariable region. J Biol Chem 2015; 290:9465-77. [PMID: 25713064 PMCID: PMC4392252 DOI: 10.1074/jbc.m114.620724] [Citation(s) in RCA: 90] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2014] [Revised: 02/19/2015] [Indexed: 01/08/2023] Open
Abstract
K-Ras4B belongs to a family of small GTPases that regulates cell growth, differentiation and survival. K-ras is frequently mutated in cancer. K-Ras4B association with the plasma membrane through its farnesylated and positively charged C-terminal hypervariable region (HVR) is critical to its oncogenic function. However, the structural mechanisms of membrane association are not fully understood. Here, using confocal microscopy, surface plasmon resonance, and molecular dynamics simulations, we observed that K-Ras4B can be distributed in rigid and loosely packed membrane domains. Its membrane binding domain interaction with phospholipids is driven by membrane fluidity. The farnesyl group spontaneously inserts into the disordered lipid microdomains, whereas the rigid microdomains restrict the farnesyl group penetration. We speculate that the resulting farnesyl protrusion toward the cell interior allows oligomerization of the K-Ras4B membrane binding domain in rigid microdomains. Unlike other Ras isoforms, K-Ras4B HVR contains a single farnesyl modification and positively charged polylysine sequence. The high positive charge not only modulates specific HVR binding to anionic phospholipids but farnesyl membrane orientation. Phosphorylation of Ser-181 prohibits spontaneous farnesyl membrane insertion. The mechanism illuminates the roles of HVR modifications in K-Ras4B targeting microdomains of the plasma membrane and suggests an additional function for HVR in regulation of Ras signaling.
Collapse
Affiliation(s)
- Hyunbum Jang
- From the Basic Science Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research and Cancer and Inflammation Program, NCI-Frederick, National Institutes of Health, Frederick, Maryland 21702
| | - Sherwin J Abraham
- Department of Molecular and Cellular Physiology, Stanford University, Stanford, California 94305, Departments of Biochemistry and Molecular Genetics and
| | - Tanmay S Chavan
- Department of Molecular and Cellular Physiology, Stanford University, Stanford, California 94305, Medicinal Chemistry, University of Illinois, Chicago, Illinois 60607, and
| | | | - Lyuba Khavrutskii
- From the Basic Science Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research and Cancer and Inflammation Program, NCI-Frederick, National Institutes of Health, Frederick, Maryland 21702
| | - Nadya I Tarasova
- Cancer and Inflammation Program, NCI-Frederick, National Institutes of Health, Frederick, Maryland 21702,
| | - Ruth Nussinov
- From the Basic Science Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research and Cancer and Inflammation Program, NCI-Frederick, National Institutes of Health, Frederick, Maryland 21702, Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine, Tel Aviv University, Tel Aviv 69978, Israel
| | - Vadim Gaponenko
- Departments of Biochemistry and Molecular Genetics and Medicinal Chemistry, University of Illinois, Chicago, Illinois 60607, and
| |
Collapse
|
7
|
Wang T, Sessions AO, Lunde CS, Rouhani S, Glaeser RM, Duan Y, Facciotti MT. Deprotonation of D96 in bacteriorhodopsin opens the proton uptake pathway. Structure 2013; 21:290-7. [PMID: 23394942 DOI: 10.1016/j.str.2012.12.018] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2012] [Revised: 12/08/2012] [Accepted: 12/12/2012] [Indexed: 11/18/2022]
Abstract
Despite extensive investigation, the precise mechanism controlling the opening of the cytoplasmic proton uptake pathway in bacteriorhodopsin (bR) has remained a mystery. From an analysis of the X-ray structure of the D96G/F171C/F219L triple mutant of bR and 60 independent molecular dynamics simulations of bR photointermediates, we report that the deprotonation of D96, a key residue in proton transfer reactions, serves two roles that occur sequentially. First, D96 donates a proton to the Schiff base. Subsequently, the deprotonation of D96 serves to "unlatch" the cytoplasmic side. The latching function of D96 appears to be remarkably robust, functioning to open hydration channels in all photointermediate structures. These results suggest that the protonation state of D96 may be the critical biophysical cue controlling the opening and closing of the cytoplasmic half-channel in bR. We suspect that this protonation-switch mechanism could also be utilized in other proton pumps to minimize backflow and reinforce directionality.
Collapse
Affiliation(s)
- Ting Wang
- Genome Center and Department of Biomedical Engineering, University of California at Davis, One Shields Avenue, Davis, CA 95616, USA
| | | | | | | | | | | | | |
Collapse
|
8
|
Yang J, Aslimovska L, Glaubitz C. Molecular Dynamics of Proteorhodopsin in Lipid Bilayers by Solid-State NMR. J Am Chem Soc 2011; 133:4874-81. [DOI: 10.1021/ja109766n] [Citation(s) in RCA: 73] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Jun Yang
- Institute for Biophysical Chemistry & Centre for Biomolecular Magnetic Resonance, Goethe University Frankfurt/M., Max-von-Laue Strasse 9, 60438 Frankfurt/M., Germany
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Center for Magnetic Resonance, Wuhan Institute of Physics and Mathematics, The Chinese Academy of Sciences, Wuhan 430071, China
| | - Lubica Aslimovska
- Institute for Biophysical Chemistry & Centre for Biomolecular Magnetic Resonance, Goethe University Frankfurt/M., Max-von-Laue Strasse 9, 60438 Frankfurt/M., Germany
| | - Clemens Glaubitz
- Institute for Biophysical Chemistry & Centre for Biomolecular Magnetic Resonance, Goethe University Frankfurt/M., Max-von-Laue Strasse 9, 60438 Frankfurt/M., Germany
| |
Collapse
|
9
|
Khelashvili G, Grossfield A, Feller SE, Pitman MC, Weinstein H. Structural and dynamic effects of cholesterol at preferred sites of interaction with rhodopsin identified from microsecond length molecular dynamics simulations. Proteins 2009; 76:403-17. [PMID: 19173312 PMCID: PMC4101808 DOI: 10.1002/prot.22355] [Citation(s) in RCA: 100] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
An unresolved question about GPCR function is the role of membrane components in receptor stability and activation. In particular, cholesterol is known to affect the function of membrane proteins, but the details of its effect on GPCRs are still elusive. Here, we describe how cholesterol modulates the behavior of the TM1-TM2-TM7-helix 8(H8) functional network that comprises the highly conserved NPxxY(x)(5,6)F motif, through specific interactions with the receptor. The inferences are based on the analysis of microsecond length molecular dynamics (MD) simulations of rhodopsin in an explicit membrane environment. Three regions on the rhodopsin exhibit the highest cholesterol density throughout the trajectory: the extracellular end of TM7, a location resembling the high-density sterol area from the electron microscopy data; the intracellular parts of TM1, TM2, and TM4, a region suggested as the cholesterol binding site in the recent X-ray crystallography data on beta(2)-adrenergic GPCR; and the intracellular ends of TM2-TM3, a location that was categorized as the high cholesterol density area in multiple independent 100 ns MD simulations of the same system. We found that cholesterol primarily affects specific local perturbations of the helical TM domains such as the kinks in TM1, TM2, and TM7. These local distortions, in turn, relate to rigid-body motions of the TMs in the TM1-TM2-TM7-H8 bundle. The specificity of the effects stems from the nonuniform distribution of cholesterol around the protein. Through correlation analysis we connect local effects of cholesterol on structural perturbations with a regulatory role of cholesterol in the structural rearrangements involved in GPCR function.
Collapse
Affiliation(s)
- George Khelashvili
- Department of Physiology and Biophysics, Weill Medical College of Cornell University, New York, New York 10021, USA.
| | | | | | | | | |
Collapse
|
10
|
Wassenaar TA, Daura X, Padrós E, Mark AE. Calcium binding to the purple membrane: A molecular dynamics study. Proteins 2009; 74:669-81. [PMID: 18704943 DOI: 10.1002/prot.22182] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The purple membrane (PM) is a specialized membrane patch found in halophilic archaea, containing the photoreceptor bacteriorhodopsin (bR). It is long known that calcium ions bind to the PM, but their position and role remain elusive to date. Molecular dynamics simulations in conjunction with a highly detailed model of the PM have been used to investigate the stability of calcium ions placed at three proposed cation binding sites within bR, one near the Schiff base, one in the region of the proton release group, and one near Glu9. The simulations suggest that, of the sites investigated, the binding of calcium ions was most likely at the proton release group. Binding in the region of the Schiff base, while possible, was associated with significant changes in local geometry. Calcium ions placed near Glu9 in the interior of bR (simultaneously to a Ca(2+) near the Schiff base and another one near the Glu194-Glu204 site) were not stable. The results obtained are discussed in relation to recent experimental observations and theoretical considerations.
Collapse
Affiliation(s)
- Tsjerk A Wassenaar
- Department of Biophysical Chemistry, Groningen Biomolecular Sciences and Biotechnology Institute (GBB), University of Groningen, 9747 AG Groningen, The Netherlands
| | | | | | | |
Collapse
|
11
|
Abstract
The protein bacteriorhodopsin pumps protons across a bacterial membrane; its pumping cycle is triggered by the photoisomerization of a retinal cofactor and involves multiple proton-transfer reactions between intermittent protonation sites. These transfers are either direct or mediated by hydrogen-bonded networks, which may include internal water molecules. The terminal step of the proton-transfer sequence is the proton release from a pocket near Glu194 and Glu204 to the extracellular bulk during the transition from the L to the M photointermediate states. The polar and charged side chains connecting these two regions in the crystal structures show no structural changes between the initial bR state and the L/M states, and no intermittent protonation changes have been detected so far in this region. Based on biomolecular simulations, we propose two potential proton-release channels, which connect the release pocket to the extracellular medium. In simulations of the L photointermediate we observe bulk water entering these channels and forming transient hydrogen-bonded networks, which could serve as fast deprotonation pathways from the release pocket to the bulk via a Grotthuss mechanism. For the first channel, we find that the triple Arg7, Glu9, and Tyr79 acts as a valve, thereby gating water uptake and release. The second channel has two release paths, which split at the position Asn76/Pro77 underneath the release group. Here, water molecules either exchange directly with the bulk or diffuse within the protein towards Arg 134/Lys129, where the exchange with the bulk occurs.
Collapse
Affiliation(s)
- Alain Chaumont
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | | | | | | |
Collapse
|
12
|
Jang H, Michaud-Agrawal N, Johnston JM, Woolf TB. How to lose a kink and gain a helix: pH independent conformational changes of the fusion domains from influenza hemagglutinin in heterogeneous lipid bilayers. Proteins 2008; 72:299-312. [PMID: 18214961 DOI: 10.1002/prot.21925] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
We have simulated two conformations of the fusion domain of influenza hemagglutinin (HA) within explicit water, salt, and heterogeneous lipid bilayers composed of POPC:POPG (4:1). Each conformation has seven different starting points in which the initial peptide structure is the same for each conformation, but the location across the membrane normal and lipid arrangement around the peptide are varied, giving a combined total simulation time of 140 ns. For the HA5 conformation (primary structure from recent NMR spectroscopy at pH = 5), the peptide exhibits a stable and less kinked structure in the lipid bilayer compared to that from the NMR studies. The relative fusogenic behavior of the different conformations has been investigated by calculation of the relative free energy of insertion into the hydrophobic region of lipid bilayer as a function of the depth of immersion. For the HA7 conformations (primary structure from recent NMR spectroscopy at pH = 7.4), while the N-terminal helix preserves its initial structure, the flexible C-terminal chain produces a transient helical motif inside the lipid bilayer. This conformational change is pH-independent, and is closely related to the peptide insertion into the lipid bilayer.
Collapse
Affiliation(s)
- Hyunbum Jang
- Department of Physiology, Johns Hopkins University School of Medicine, Baltimore, Maryland 21205, USA
| | | | | | | |
Collapse
|
13
|
Konidala P, Niemeyer B. Molecular dynamics simulations of pea (Pisum sativum) lectin structure with octyl glucoside detergents: the ligand interactions and dynamics. Biophys Chem 2007; 128:215-30. [PMID: 17532552 DOI: 10.1016/j.bpc.2007.04.006] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2007] [Revised: 04/18/2007] [Accepted: 04/18/2007] [Indexed: 11/23/2022]
Abstract
The mitogenic pea (Pisum sativum) lectin is a legume protein of non-immunoglobulin nature capable of specific recognition of glucose derivatives without altering its structure. Molecular dynamics simulations were performed in a realistic environment to investigate the structure and interaction properties of pea lectin with various concentrations of n-octyl-beta-d-glucopyranoside (OG) detergent monomers distributed inside explicit solvent cell. In addition, the diffusion coefficients of the ligands (OG, Ca2+, Mn2+, and Cl-) and the water molecules were also reported. The structural flexibility of the lectin was conserved in all simulations. The self-assembly of OG monomers into a small micelle at the hydrophobic site of the lectin was noticed in the simulation with 20 OG monomers. The interaction energy analysis concludes that the lectin was appropriately termed an adaptive structure. One or rarely two binding sites were observed at an instant in each simulation that were electrostatically favoured for the OG to interact with the surface amino acid residues. Enhanced binding of OG to the pea lectin was quantified in the system containing only Ca2+ divalent ions. Interestingly, no binding was observed in the simulation without divalent ions. Furthermore, the lectin-ligand complex was stabilized by multiple hydrogen bonds and at least one water bridge. Finally, the work was also in accordance with the published work elsewhere that the simulations performed with different initial conditions and using higher nonbonded cutoffs for the van der Waals and electrostatic interactions provide more accurate information and clues than the single large simulation of the biomolecular system of interest.
Collapse
Affiliation(s)
- Praveen Konidala
- Institute of Thermodynamics, Helmut-Schmidt-University / University of the Federal Armed Forces Hamburg, Holstenhofweg 85, D-22043 Hamburg, Germany
| | | |
Collapse
|
14
|
Jang H, Zheng J, Nussinov R. Models of beta-amyloid ion channels in the membrane suggest that channel formation in the bilayer is a dynamic process. Biophys J 2007; 93:1938-49. [PMID: 17526580 PMCID: PMC1959551 DOI: 10.1529/biophysj.107.110148] [Citation(s) in RCA: 156] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Here we model the Alzheimer beta-peptide ion channel with the goal of obtaining insight into the mechanism of amyloid toxicity. The models are built based on NMR data of the oligomers, with the universal U-shaped (strand-turn-strand) motif. After 30-ns simulations in the bilayer, the channel dimensions, shapes and subunit organization are in good agreement with atomic force microscopy (AFM). The models use the Abeta(17-42) pentamer NMR-based coordinates. Extension and bending of the straight oligomers lead to two channel topologies, depending on the direction of the curvature: 1), the polar/charged N-terminal beta-strand of Abeta(17-42) faces the water-filled pore, and the hydrophobic C-terminal beta-strand faces the bilayer (CNpNC; p for pore); and 2), the C-terminal beta-strand faces the solvated pore (NCpCN). In the atomistic simulations in a fully solvated DOPC lipid bilayer, the first (CNpNC) channel preserves the pore and conducts solvent; by contrast, hydrophobic collapse blocks the NCpCN channel. AFM demonstrated open pores and collapsed complexes. The final averaged CNpNC pore dimensions (outer diameter 8 nm; inner diameter approximately 2.5 nm) are in the AFM range (8-12 nm; approximately 2 nm, respectively). Further, in agreement with high-resolution AFM images, during the simulations, the channels spontaneously break into ordered subunits in the bilayer; however, we also observe that the subunits are loosely connected by partially disordered inner beta-sheet, suggesting subunit mobility in the bilayer. The cationic channel has strong selective affinity for Ca(2+), supporting experimental calcium-selective beta-amyloid channels. Membrane permeability and consequent disruption of calcium homeostasis were implicated in cellular degeneration. Consequently, the CNpNC channel topology can sign cell death, offering insight into amyloid toxicity via an ion "trap-release" transport mechanism. The observed loosely connected subunit organization suggests that amyloid channel formation in the bilayer is a dynamic, fluid process involving subunit association, dissociation, and channel rearrangements.
Collapse
Affiliation(s)
- Hyunbum Jang
- Center for Cancer Research Nanobiology Program, SAIC-Frederick, National Cancer Institute-Frederick, Frederick, MD 21702, USA
| | | | | |
Collapse
|
15
|
Jang H, Ma B, Woolf TB, Nussinov R. Interaction of protegrin-1 with lipid bilayers: membrane thinning effect. Biophys J 2006; 91:2848-59. [PMID: 16861271 PMCID: PMC1578484 DOI: 10.1529/biophysj.106.084046] [Citation(s) in RCA: 61] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2006] [Accepted: 07/12/2006] [Indexed: 11/18/2022] Open
Abstract
Protegrins (PG) are important in defending host tissues, preventing infection via an attack on the membrane surface of invading microorganisms. Protegrins have powerful antibiotic abilities, but the molecular-level mechanisms underlying the interactions of their beta-sheet motifs with the membrane are not known. Protegrin-1 (PG-1) is composed of 18 amino acids with a high content of basic residues and two disulfide bonds. Here we focused on the stability of PG-1 at the amphipathic interface in lipid bilayers and on the details of the peptide-membrane interactions. We simulated all-atom models of the PG-1 monomer with explicit water and lipid bilayers composed of both homogeneous POPC (palmitoyl-oleyl-phosphatidylcholine) lipids and a mixture of POPC/POPG (palmitoyl-oleyl-phosphatidylglycerol) (4:1) lipids. We observed that local thinning of the lipid bilayers mediated by the peptide is enhanced in the lipid bilayer containing POPG, consistent with experimental results of selective membrane targeting. The beta-hairpin motif of PG-1 is conserved in both lipid settings, whereas it is highly bent in aqueous solution. The conformational dynamics of PG-1, especially the highly charged beta-hairpin turn region, are found to be mostly responsible for disturbing the membrane. Even though the eventual membrane disruption requires PG-1 oligomers, our simulations clearly show the first step of the monomeric effects. The thinning effects in the bilayer should relate to pore/channel formation in the lipid bilayer and thus be responsible for further defects in the membrane caused by oligomer.
Collapse
Affiliation(s)
- Hyunbum Jang
- Center for Cancer Research Nanobiology Program, SAIC-Frederick, Inc., NCI-Frederick, Frederick, Maryland 21702, USA
| | | | | | | |
Collapse
|
16
|
Tanizaki S, Feig M. Molecular dynamics simulations of large integral membrane proteins with an implicit membrane model. J Phys Chem B 2006; 110:548-56. [PMID: 16471567 DOI: 10.1021/jp054694f] [Citation(s) in RCA: 61] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The heterogeneous dielectric generalized Born (HDGB) methodology is an the extension of the GBMV model for the simulation of integral membrane proteins with an implicit membrane environment. Three large integral membrane proteins, the bacteriorhodopsin monomer and trimer and the BtuCD protein, were simulated with the HDGB model in order to evaluate how well thermodynamic and dynamic properties are reproduced. Effects of the truncation of electrostatic interactions were examined. For all proteins, the HDGB model was able to generate stable trajectories that remained close to the starting experimental structures, in excellent agreement with explicit membrane simulations. Dynamic properties evaluated through a comparison of B-factors are also in good agreement with experiment and explicit membrane simulations. However, overall flexibility was slightly underestimated with the HDGB model unless a very large electrostatic cutoff is employed. Results with the HDGB model are further compared with equivalent simulations in implicit aqueous solvent, demonstrating that the membrane environment leads to more realistic simulations.
Collapse
Affiliation(s)
- Seiichiro Tanizaki
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824-1319, USA
| | | |
Collapse
|
17
|
Tanizaki S, Feig M. A generalized Born formalism for heterogeneous dielectric environments: application to the implicit modeling of biological membranes. J Chem Phys 2006; 122:124706. [PMID: 15836408 DOI: 10.1063/1.1865992] [Citation(s) in RCA: 153] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Reliable computer simulations of complex biological environments such as integral membrane proteins with explicit water and lipid molecules remain a challenging task. We propose a modification of the standard generalized Born theory of homogeneous solvent for modeling the heterogeneous dielectric environments such as lipid/water interfaces. Our model allows the representation of biological membranes in the form of multiple layered dielectric regions with dielectric constants that are different from the solute cavity. The proposed new formalism is shown to predict the electrostatic component of solvation free energy with a relative error of 0.17% compared to exact finite-difference solutions of the Poisson equation for a transmembrane helix test system. Molecular dynamics simulations of melittin and bacteriorhodopsin are carried out and performed over 10 ns and 7 ns of simulation time, respectively. The center of melittin along the membrane normal in these stable simulations is in excellent agreement with the relevant experimental data. Simulations of bacteriorhodopsin started from the experimental structure remained stable and in close agreement with experiment. We also examined the free energy profiles of water and amino acid side chain analogs upon membrane insertion. The results with our implicit membrane model agree well with the experimental transfer free energy data from cyclohexane to water as well as explicit solvent simulations of water and selected side chain analogs.
Collapse
Affiliation(s)
- Seiichiro Tanizaki
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824-1319, USA
| | | |
Collapse
|
18
|
Gumbart J, Wang Y, Aksimentiev A, Tajkhorshid E, Schulten K. Molecular dynamics simulations of proteins in lipid bilayers. Curr Opin Struct Biol 2005; 15:423-31. [PMID: 16043343 PMCID: PMC2474857 DOI: 10.1016/j.sbi.2005.07.007] [Citation(s) in RCA: 129] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2005] [Revised: 06/14/2005] [Accepted: 07/13/2005] [Indexed: 11/19/2022]
Abstract
With recent advances in X-ray crystallography of membrane proteins promising many new high-resolution structures, molecular dynamics simulations will become increasingly valuable for understanding membrane protein function, as they can reveal the dynamic behavior concealed in the static structures. Dramatic increases in computational power, in synergy with more efficient computational methodologies, now allow us to carry out molecular dynamics simulations of any structurally known membrane protein in its native environment, covering timescales of up to 0.1 micros. At the frontiers of membrane protein simulations are ion channels, aquaporins, passive and active transporters, and bioenergetic proteins.
Collapse
Affiliation(s)
- James Gumbart
- Theoretical and Computational Biophysics Group, Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | | | | | | | | |
Collapse
|
19
|
Kaledin AL, Kaledin M, Bowman JM. All-Atom Calculation of the Normal Modes of Bacteriorhodopsin Using a Sliding Block Iterative Diagonalization Method. J Chem Theory Comput 2005; 2:166-74. [DOI: 10.1021/ct050161z] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Alexey L. Kaledin
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computing, Emory University, Atlanta, Georgia 30322
| | - Martina Kaledin
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computing, Emory University, Atlanta, Georgia 30322
| | - Joel M. Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computing, Emory University, Atlanta, Georgia 30322
| |
Collapse
|
20
|
Lee AG. How lipids and proteins interact in a membrane: a molecular approach. MOLECULAR BIOSYSTEMS 2005; 1:203-12. [PMID: 16880984 DOI: 10.1039/b504527d] [Citation(s) in RCA: 101] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Membrane proteins in a biological membrane are surrounded by a shell or annulus of 'solvent' lipid molecules. These lipid molecules in general interact rather non-specifically with the protein molecules, although a few 'hot-spots' may be present on the protein where anionic lipids bind with high affinity. Because of the low structural specificity of most of the annular sites, the composition of the lipid annulus will be rather similar to the bulk lipid composition of the membrane. The structures of the solvent lipid molecules are important in determining the conformational state of a membrane protein, and hence its activity, through charge and hydrogen bonding interactions between the lipid headgroups and residues in the protein, and through hydrophobic matching between the protein and the surrounding lipid bilayer. Evidence is also accumulating for the presence of 'co-factor' lipid molecules binding with high specificity to membrane proteins, often between transmembrane alpha-helices, and often being essential for activity.
Collapse
Affiliation(s)
- Anthony G Lee
- School of Biological Sciences, University of Southampton, Southampton, UK.
| |
Collapse
|
21
|
Karnaukhova E, Vasileiou C, Wang A, Berova N, Nakanishi K, Borhan B. Circular dichroism of heterochromophoric and partially regenerated purple membrane: Search for exciton coupling. Chirality 2005; 18:72-83. [PMID: 16385624 DOI: 10.1002/chir.20222] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
In order to determine the origin of the bisignate CD spectra of native purple membrane, heterochromophoric analogues containing bacteriorhodopsin regenerated with native all-trans-retinal and retinal analogues were investigated. The data collected for the purple membrane samples containing two different chromophores suggest the additive character of the CD spectra. This conclusion was supported by a series of spectra using 5,6-dihydroretinal and 3-dehydroretinal and by using 33% regenerated PM in buffer and in presence of osmolytes. Our results support the idea of conformational heterogeneity of the chromophores in the bR in the trimer, suggesting that the three bR subunits in the trimer are not conformationally equal, and therefore, the bisignate CD spectrum of bR in the purple membrane occurs rather due to a superposition of the CD spectra from variously distorted bR subunits in the trimer than interchromophoric exciton-coupling interactions.
Collapse
Affiliation(s)
- Elena Karnaukhova
- Center for Biologics Evaluation and Research, Food and Drug Administration, Bethesda, MD 20892, USA
| | | | | | | | | | | |
Collapse
|