1
|
Woods H, Leman JK, Meiler J. Modeling membrane geometries implicitly in Rosetta. Protein Sci 2024; 33:e4908. [PMID: 38358133 PMCID: PMC10868433 DOI: 10.1002/pro.4908] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2023] [Revised: 01/05/2024] [Accepted: 01/08/2024] [Indexed: 02/16/2024]
Abstract
Interactions between membrane proteins (MPs) and lipid bilayers are critical for many cellular functions. In the Rosetta molecular modeling suite, the implicit membrane energy function is based on a "slab" model, which represent the membrane as a flat bilayer. However, in nature membranes often have a curvature that is important for function and/or stability. Even more prevalent, in structural biology research MPs are reconstituted in model membrane systems such as micelles, bicelles, nanodiscs, or liposomes. Thus, we have modified the existing membrane energy potentials within the RosettaMP framework to allow users to model MPs in different membrane geometries. We show that these modifications can be utilized in core applications within Rosetta such as structure refinement, protein-protein docking, and protein design. For MP structures found in curved membranes, refining these structures in curved, implicit membranes produces higher quality models with structures closer to experimentally determined structures. For MP systems embedded in multiple membranes, representing both membranes results in more favorable scores compared to only representing one of the membranes. Modeling MPs in geometries mimicking the membrane model system used in structure determination can improve model quality and model discrimination.
Collapse
Affiliation(s)
- Hope Woods
- Center of Structural Biology, Vanderbilt UniversityNashvilleTennesseeUSA
- Chemical and Physical Biology ProgramVanderbilt UniversityNashvilleTennesseeUSA
| | | | - Jens Meiler
- Center of Structural Biology, Vanderbilt UniversityNashvilleTennesseeUSA
- Department of ChemistryVanderbilt UniversityNashvilleTennesseeUSA
- Institute for Drug Discovery, Leipzig University Medical SchoolLeipzigGermany
| |
Collapse
|
2
|
Maurer M, Lazaridis T. Transmembrane β-Barrel Models of α-Synuclein Oligomers. J Chem Inf Model 2023; 63:7171-7179. [PMID: 37963823 DOI: 10.1021/acs.jcim.3c00997] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2023]
Abstract
The aggregation of α-synuclein is implicated in a number of neurodegenerative diseases, such as Parkinson's and Multiple System Atrophy, but the role of these aggregates in disease development is not clear. One possible mechanism of cytotoxicity is the disturbance or permeabilization of cell membranes by certain types of oligomers. However, no high-resolution structure of such membrane-embedded complexes has ever been determined. Here we construct and evaluate putative transmembrane β-barrels formed by this protein. Examination of the α-synuclein sequence reveals two regions that could form membrane-embedded β-hairpins: 64-92 (the NAC), and 35-56, which harbors many familial Parkinson's mutations. The stability of β-barrels formed by these hairpins is examined first in implicit membrane pores and then by multimicrosecond all-atom simulations. We find that a NAC region barrel remains stably inserted and hydrated for at least 10 μs. A 35-56 barrel remains stably inserted in the membrane but dehydrates and collapses if all His50 are neutral or if His50 is replaced by Q. If half of the His50 are doubly protonated, the barrel takes an oval shape but remains hydrated for at least 10 μs. Possible implications of these findings for α-synuclein pathology are discussed.
Collapse
Affiliation(s)
- Manuela Maurer
- Department of Chemistry & Biochemistry, City College of New York/CUNY, 160 Convent Ave, New York, New York 10031, United States
| | - Themis Lazaridis
- Department of Chemistry & Biochemistry, City College of New York/CUNY, 160 Convent Ave, New York, New York 10031, United States
| |
Collapse
|
3
|
Vasquez Rodriguez SY, Lazaridis T. Simulations suggest a scaffolding mechanism of membrane deformation by the caveolin 8S complex. Biophys J 2023; 122:4082-4090. [PMID: 37742070 PMCID: PMC10598286 DOI: 10.1016/j.bpj.2023.09.008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2023] [Revised: 09/11/2023] [Accepted: 09/14/2023] [Indexed: 09/25/2023] Open
Abstract
Caveolins form complexes of various sizes that deform membranes into polyhedral shapes. However, the recent structure of the 8S complex was disk-like with a flat membrane-binding surface. How can a flat complex deform membranes into nonplanar structures? Molecular dynamics simulations revealed that the 8S complex rapidly takes the form of a suction cup. Simulations on implicit membrane vesicles determined that binding is stronger when E140 gets protonated. In that case, the complex binds much more strongly to 5- and 10-nm-radius vesicles. A concave membrane-binding surface readily explains the membrane-deforming ability of caveolins by direct scaffolding. We propose that the 8S complex sits at the vertices of the caveolar polyhedra, rather than at the center of the polyhedral faces.
Collapse
Affiliation(s)
| | - Themis Lazaridis
- Department of Chemistry, City College of New York/CUNY, New York, New York; Graduate Programs in Chemistry, Biochemistry, and Physics, The Graduate Center, City University of New York, New York, New York.
| |
Collapse
|
4
|
Jung JH, Seo PJ, Oh E, Kim J. Temperature perception by plants. TRENDS IN PLANT SCIENCE 2023; 28:924-940. [PMID: 37045740 DOI: 10.1016/j.tplants.2023.03.006] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/18/2022] [Revised: 02/16/2023] [Accepted: 03/09/2023] [Indexed: 06/19/2023]
Abstract
Plants constantly face fluctuating ambient temperatures and must adapt to survive under stressful conditions. Temperature affects many aspects of plant growth and development through a complex network of transcriptional responses. Although temperature sensing is a crucial primary step in initiating transcriptional responses via Ca2+ and/or reactive oxygen species signaling, an understanding of how plants perceive temperature has remained elusive. However, recent studies have yielded breakthroughs in our understanding of temperature sensors and thermosensation mechanisms. We review recent findings on potential temperature sensors and emerging thermosensation mechanisms, including biomolecular condensate formation through phase separation in plants. We also compare the temperature perception mechanisms of plants with those of other organisms to provide insights into understanding temperature sensing by plants.
Collapse
Affiliation(s)
- Jae-Hoon Jung
- Department of Biological Sciences, Sungkyunkwan University, Suwon 16419, Korea
| | - Pil Joon Seo
- Department of Chemistry, Seoul National University, Seoul 08826, Korea
| | - Eunkyoo Oh
- Department of Life Sciences, Korea University, Seoul 02841, Korea
| | - Jungmook Kim
- Department of Bioenergy Science and Technology, Chonnam National University, Gwangju 61186, Korea; Department of Integrative Food, Bioscience, and Technology, Chonnam National University, Gwangju 61186, Korea.
| |
Collapse
|
5
|
Li R, Pang L. Comparing the effects of proteins with IDRs on membrane system in yeast, mammalian cells, and the model plant Arabidopsis. CURRENT OPINION IN PLANT BIOLOGY 2023; 74:102375. [PMID: 37172364 DOI: 10.1016/j.pbi.2023.102375] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/21/2022] [Revised: 03/26/2023] [Accepted: 04/07/2023] [Indexed: 05/14/2023]
Abstract
Membrane vesiculation is an energy-costing process. Previous studies paid much attention to proteins with curvature-inducing motifs. Recent publications reveal that the liquid-like protein assembly on membrane surfaces provides an efficient yet structure-independent mechanism for increasing the membrane curvature, which plays important roles in vesicle transport in many aspects. Intrinsically disordered regions (IDRs) within the proteins are highly potent drivers of membrane curvature by providing large hydrodynamic radii to generate steric pressure. Biomolecular condensates formed by phase separation can provide a reaction platform for sequential processes or generate a wetting surface to sequestrate cargos and trigger membrane remodeling. We review the latest progress in yeast and mammalian cells, focus on the mechanism of clathrin-mediated endocytosis (CME) and autophagy initiation, and compare with what we know in model plant Arabidopsis. The comparison may give important insights into the understanding of basic membrane trafficking mechanisms in plant cells.
Collapse
Affiliation(s)
- Ruixi Li
- Key Laboratory of Molecular Design for Plant Cell Factory of Guangdong Higher Education Institutes, Institute of Plant and Food Science, School of Life Sciences, Southern University of Science and Technology, Shenzhen, 518055, China.
| | - Lei Pang
- Key Laboratory of Molecular Design for Plant Cell Factory of Guangdong Higher Education Institutes, Institute of Plant and Food Science, School of Life Sciences, Southern University of Science and Technology, Shenzhen, 518055, China
| |
Collapse
|
6
|
Cyanobacterial membrane dynamics in the light of eukaryotic principles. Biosci Rep 2023; 43:232406. [PMID: 36602300 PMCID: PMC9950537 DOI: 10.1042/bsr20221269] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2022] [Revised: 12/23/2022] [Accepted: 01/03/2023] [Indexed: 01/06/2023] Open
Abstract
Intracellular compartmentalization is a hallmark of eukaryotic cells. Dynamic membrane remodeling, involving membrane fission/fusion events, clearly is crucial for cell viability and function, as well as membrane stabilization and/or repair, e.g., during or after injury. In recent decades, several proteins involved in membrane stabilization and/or dynamic membrane remodeling have been identified and described in eukaryotes. Yet, while typically not having a cellular organization as complex as eukaryotes, also bacteria can contain extra internal membrane systems besides the cytoplasmic membranes (CMs). Thus, also in bacteria mechanisms must have evolved to stabilize membranes and/or trigger dynamic membrane remodeling processes. In fact, in recent years proteins, which were initially defined being eukaryotic inventions, have been recognized also in bacteria, and likely these proteins shape membranes also in these organisms. One example of a complex prokaryotic inner membrane system is the thylakoid membrane (TM) of cyanobacteria, which contains the complexes of the photosynthesis light reaction. Cyanobacteria are evolutionary closely related to chloroplasts, and extensive remodeling of the internal membrane systems has been observed in chloroplasts and cyanobacteria during membrane biogenesis and/or at changing light conditions. We here discuss common principles guiding eukaryotic and prokaryotic membrane dynamics and the proteins involved, with a special focus on the dynamics of the cyanobacterial TMs and CMs.
Collapse
|
7
|
Xue M, Cao Y, Shen C, Guo W. Computational Advances of Protein/Neurotransmitter-membrane Interactions Involved in Vesicle Fusion and Neurotransmitter Release. J Mol Biol 2023; 435:167818. [PMID: 36089056 DOI: 10.1016/j.jmb.2022.167818] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2022] [Revised: 08/22/2022] [Accepted: 09/04/2022] [Indexed: 02/04/2023]
Abstract
Vesicle fusion is of crucial importance to neuronal communication at neuron terminals. The exquisite but complex fusion machinery for neurotransmitter release is tightly controlled and regulated by protein/neurotransmitter-membrane interactions. Computational 'microscopies', in particular molecular dynamics simulations and related techniques, have provided notable insight into the physiological process over the past decades, and have made enormous contributions to fields such as neurology, pharmacology and pathophysiology. Here we review the computational advances of protein/neurotransmitter-membrane interactions related to presynaptic vesicle-membrane fusion and neurotransmitter release, and outline the in silico challenges ahead for understanding this important physiological process.
Collapse
Affiliation(s)
- Minmin Xue
- Key Laboratory for Intelligent Nano Materials and Devices of Ministry of Education, State Key Laboratory of Mechanics and Control of Mechanical Structures, and Institute for Frontier Science, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
| | - Yuwei Cao
- State Key Laboratory of Tribology, Department of Mechanical Engineering, Tsinghua University, Beijing, China
| | - Chun Shen
- Key Laboratory for Intelligent Nano Materials and Devices of Ministry of Education, State Key Laboratory of Mechanics and Control of Mechanical Structures, and Institute for Frontier Science, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China.
| | - Wanlin Guo
- Key Laboratory for Intelligent Nano Materials and Devices of Ministry of Education, State Key Laboratory of Mechanics and Control of Mechanical Structures, and Institute for Frontier Science, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China; State Key Laboratory of Tribology, Department of Mechanical Engineering, Tsinghua University, Beijing, China.
| |
Collapse
|
8
|
Nepal B, Das S, Reith ME, Kortagere S. Overview of the structure and function of the dopamine transporter and its protein interactions. Front Physiol 2023; 14:1150355. [PMID: 36935752 PMCID: PMC10020207 DOI: 10.3389/fphys.2023.1150355] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2023] [Accepted: 02/21/2023] [Indexed: 03/06/2023] Open
Abstract
The dopamine transporter (DAT) plays an integral role in dopamine neurotransmission through the clearance of dopamine from the extracellular space. Dysregulation of DAT is central to the pathophysiology of numerous neuropsychiatric disorders and as such is an attractive therapeutic target. DAT belongs to the solute carrier family 6 (SLC6) class of Na+/Cl- dependent transporters that move various cargo into neurons against their concentration gradient. This review focuses on DAT (SCL6A3 protein) while extending the narrative to the closely related transporters for serotonin and norepinephrine where needed for comparison or functional relevance. Cloning and site-directed mutagenesis experiments provided early structural knowledge of DAT but our contemporary understanding was achieved through a combination of crystallization of the related bacterial transporter LeuT, homology modeling, and subsequently the crystallization of drosophila DAT. These seminal findings enabled a better understanding of the conformational states involved in the transport of substrate, subsequently aiding state-specific drug design. Post-translational modifications to DAT such as phosphorylation, palmitoylation, ubiquitination also influence the plasma membrane localization and kinetics. Substrates and drugs can interact with multiple sites within DAT including the primary S1 and S2 sites involved in dopamine binding and novel allosteric sites. Major research has centered around the question what determines the substrate and inhibitor selectivity of DAT in comparison to serotonin and norepinephrine transporters. DAT has been implicated in many neurological disorders and may play a role in the pathology of HIV and Parkinson's disease via direct physical interaction with HIV-1 Tat and α-synuclein proteins respectively.
Collapse
Affiliation(s)
- Binod Nepal
- Department of Microbiology and Immunology, Drexel University College of Medicine, Philadelphia, PA, United States
| | - Sanjay Das
- Department of Microbiology and Immunology, Drexel University College of Medicine, Philadelphia, PA, United States
| | - Maarten E. Reith
- Department of Psychiatry, New York University School of Medicine, New York City, NY, United States
| | - Sandhya Kortagere
- Department of Microbiology and Immunology, Drexel University College of Medicine, Philadelphia, PA, United States
- *Correspondence: Sandhya Kortagere,
| |
Collapse
|
9
|
Ouyang J, Sheng Y, Wang W. Recent Advances of Studies on Cell-Penetrating Peptides Based on Molecular Dynamics Simulations. Cells 2022; 11:cells11244016. [PMID: 36552778 PMCID: PMC9776715 DOI: 10.3390/cells11244016] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2022] [Revised: 12/09/2022] [Accepted: 12/10/2022] [Indexed: 12/14/2022] Open
Abstract
With the ability to transport cargo molecules across cell membranes with low toxicity, cell-penetrating peptides (CPPs) have become promising candidates for next generation peptide-based drug delivery vectors. Over the past three decades since the first CPP was discovered, a great deal of work has been done on the cellular uptake mechanisms and the applications for the delivery of therapeutic molecules, and significant advances have been made. But so far, we still do not have a precise and unified understanding of the structure-activity relationship of the CPPs. Molecular dynamics (MD) simulations provide a method to reveal peptide-membrane interactions at the atomistic level and have become an effective complement to experiments. In this paper, we review the progress of the MD simulations on CPP-membrane interactions, including the computational methods and technical improvements in the MD simulations, the research achievements in the CPP internalization mechanism, CPP decoration and coupling, and the peptide-induced membrane reactions during the penetration process, as well as the comparison of simulated and experimental results.
Collapse
Affiliation(s)
- Jun Ouyang
- School of Public Courses, Bengbu Medical College, Bengbu 233030, China
- Collaborative Innovation Center of Advanced Microstructures, National Laboratory of Solid State Microstructure, Department of Physics, Nanjing University, Nanjing 210093, China
| | - Yuebiao Sheng
- Collaborative Innovation Center of Advanced Microstructures, National Laboratory of Solid State Microstructure, Department of Physics, Nanjing University, Nanjing 210093, China
- High Performance Computing Center, Nanjing University, Nanjing 210093, China
- Correspondence: (Y.S.); (W.W.)
| | - Wei Wang
- Collaborative Innovation Center of Advanced Microstructures, National Laboratory of Solid State Microstructure, Department of Physics, Nanjing University, Nanjing 210093, China
- Correspondence: (Y.S.); (W.W.)
| |
Collapse
|
10
|
van Hilten N, Stroh KS, Risselada HJ. Efficient Quantification of Lipid Packing Defect Sensing by Amphipathic Peptides: Comparing Martini 2 and 3 with CHARMM36. J Chem Theory Comput 2022; 18:4503-4514. [PMID: 35709386 PMCID: PMC9281404 DOI: 10.1021/acs.jctc.2c00222] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In biological systems, proteins can be attracted to curved or stretched regions of lipid bilayers by sensing hydrophobic defects in the lipid packing on the membrane surface. Here, we present an efficient end-state free energy calculation method to quantify such sensing in molecular dynamics simulations. We illustrate that lipid packing defect sensing can be defined as the difference in mechanical work required to stretch a membrane with and without a peptide bound to the surface. We also demonstrate that a peptide's ability to concurrently induce excess leaflet area (tension) and elastic softening─a property we call the "characteristic area of sensing" (CHAOS)─and lipid packing sensing behavior are in fact two sides of the same coin. In essence, defect sensing displays a peptide's propensity to generate tension. The here-proposed mechanical pathway is equally accurate yet, computationally, about 40 times less costly than the commonly used alchemical pathway (thermodynamic integration), allowing for more feasible free energy calculations in atomistic simulations. This enabled us to directly compare the Martini 2 and 3 coarse-grained and the CHARMM36 atomistic force fields in terms of relative binding free energies for six representative peptides including the curvature sensor ALPS and two antiviral amphipathic helices (AH). We observed that Martini 3 qualitatively reproduces experimental trends while producing substantially lower (relative) binding free energies and shallower membrane insertion depths compared to atomistic simulations. In contrast, Martini 2 tends to overestimate (relative) binding free energies. Finally, we offer a glimpse into how our end-state-based free energy method can enable the inverse design of optimal lipid packing defect sensing peptides when used in conjunction with our recently developed evolutionary molecular dynamics (Evo-MD) method. We argue that these optimized defect sensors─aside from their biomedical and biophysical relevance─can provide valuable targets for the development of lipid force fields.
Collapse
Affiliation(s)
- Niek van Hilten
- Leiden Institute of Chemistry, Leiden University, Leiden 2300 RA, The Netherlands
| | - Kai Steffen Stroh
- Department of Physics, Technical University Dortmund, Dortmund 44221, Germany.,Institute for Theoretical Physics, Georg-August-University Göttingen, Göttingen 37077, Germany
| | - Herre Jelger Risselada
- Leiden Institute of Chemistry, Leiden University, Leiden 2300 RA, The Netherlands.,Department of Physics, Technical University Dortmund, Dortmund 44221, Germany.,Institute for Theoretical Physics, Georg-August-University Göttingen, Göttingen 37077, Germany
| |
Collapse
|
11
|
Fu Y, Zeno WF, Stachowiak JC, Johnson ME. A continuum membrane model can predict curvature sensing by helix insertion. SOFT MATTER 2021; 17:10649-10663. [PMID: 34792524 PMCID: PMC8877990 DOI: 10.1039/d1sm01333e] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Protein domains, such as ENTH (epsin N-terminal homology) and BAR (bin/amphiphysin/rvs), contain amphipathic helices that drive preferential binding to curved membranes. However, predicting how the physical parameters of these domains control this 'curvature sensing' behavior is challenging due to the local membrane deformations generated by the nanoscopic helix on the surface of a large sphere. We here use a deformable continuum model that accounts for the physical properties of the membrane and the helix insertion to predict curvature sensing behavior, with direct validation against multiple experimental datasets. We show that the insertion can be modeled as a local change to the membrane's spontaneous curvature, cins0, producing excellent agreement with the energetics extracted from experiments on ENTH binding to vesicles and cylinders, and of ArfGAP helices to vesicles. For small vesicles with high curvature, the insertion lowers the membrane energy by relieving strain on a membrane that is far from its preferred curvature of zero. For larger vesicles, however, the insertion has the inverse effect, de-stabilizing the membrane by introducing more strain. We formulate here an empirical expression that accurately captures numerically calculated membrane energies as a function of both basic membrane properties (bending modulus κ and radius R) as well as stresses applied by the inserted helix (cins0 and area Ains). We therefore predict how these physical parameters will alter the energetics of helix binding to curved vesicles, which is an essential step in understanding their localization dynamics during membrane remodeling processes.
Collapse
Affiliation(s)
- Yiben Fu
- T. C. Jenkins Department of Biophysics, The Johns Hopkins University, 3400 N. Charles St., Baltimore, Maryland 21218, USA.
| | - Wade F Zeno
- Mork Family Department of Chemical Engineering and Materials Science, The University of Southern California, Los Angeles, California, 90089, USA
| | - Jeanne C Stachowiak
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA
| | - Margaret E Johnson
- T. C. Jenkins Department of Biophysics, The Johns Hopkins University, 3400 N. Charles St., Baltimore, Maryland 21218, USA.
| |
Collapse
|
12
|
Mechanism of negative membrane curvature generation by I-BAR domains. Structure 2021; 29:1440-1452.e4. [PMID: 34520736 DOI: 10.1016/j.str.2021.07.010] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2021] [Revised: 04/16/2021] [Accepted: 07/22/2021] [Indexed: 11/23/2022]
Abstract
The membrane sculpting ability of BAR domains has been attributed to the intrinsic curvature of their banana-shaped dimeric structure. However, there is often a mismatch between this intrinsic curvature and the diameter of the membrane tubules generated. I-BAR domains are especially mysterious since they are almost flat but generate high negative membrane curvature. Here, we use atomistic implicit-solvent computer modeling to show that the membrane bending of the IRSp53 I-BAR domain is dictated by its higher oligomeric structure, whose curvature is completely unrelated to the intrinsic curvature of the dimer. Two other I-BARs give similar results, whereas a flat F-BAR sheet develops a concave membrane-binding interface, consistent with its observed positive membrane curvature generation. Laterally interacting helical spirals of I-BAR dimers on tube interiors are stable and have an enhanced binding energy that is sufficient for membrane bending to experimentally observed tubule diameters at a reasonable surface density.
Collapse
|
13
|
Stroh KS, Risselada HJ. Quantifying Membrane Curvature Sensing of Peripheral Proteins by Simulated Buckling and Umbrella Sampling. J Chem Theory Comput 2021; 17:5276-5286. [PMID: 34261315 DOI: 10.1021/acs.jctc.1c00021] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Membrane curvature plays an essential role in the organization and trafficking of membrane associated proteins. Comparison or prediction of the experimentally resolved protein concentrations adopted at different membrane curvatures requires direct quantification of the relative partitioning free energy. Here, we present a highly efficient and simple to implement a free-energy calculation method which is able to directly resolve the relative partitioning free energy of proteins as a direct function of membrane curvature, i.e., a curvature sensing profile, within (coarse-grained) molecular dynamics simulations. We demonstrate its utility by resolving these profiles for two known curvature sensing peptides, namely ALPS and α-synuclein, for a membrane curvature ranging from -1/6.5 to +1/6.5 nm-1. We illustrate that the difference in relative partitioning (binding) free energy between these two extrema is only about 13 kBT for both peptides, illustrating that the driving force of curvature sensing is subtle. Furthermore, we illustrate that ALPS and α-synuclein sense curvature via a contrasting mechanism, which is differentially affected by membrane composition. In addition, we demonstrate that the intrinsic spontaneous curvature of both of these peptides lies beyond the range of membrane curvature accessible in micropipette aspiration experiments, being about 1/7 nm -1. Our approach offers an efficient and simple to implement in silico tool for exploring and screening the membrane curvature sensing mechanisms of proteins.
Collapse
Affiliation(s)
- Kai Steffen Stroh
- Institute for Theoretical Physics, Georg-August University Göttingen, Göttingen, Germany
| | - Herre Jelger Risselada
- Institute for Theoretical Physics, Georg-August University Göttingen, Göttingen, Germany.,Leiden Institute of Chemistry, Leiden University, Leiden, The Netherlands.,Leibniz Institute for Surface Engineering, Leipzig, Germany
| |
Collapse
|
14
|
Amos SBA, Schwarz TC, Shi J, Cossins BP, Baker TS, Taylor RJ, Konrat R, Sansom MSP. Membrane Interactions of α-Synuclein Revealed by Multiscale Molecular Dynamics Simulations, Markov State Models, and NMR. J Phys Chem B 2021; 125:2929-2941. [PMID: 33719460 PMCID: PMC8006134 DOI: 10.1021/acs.jpcb.1c01281] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2021] [Revised: 03/01/2021] [Indexed: 01/30/2023]
Abstract
α-Synuclein (αS) is a presynaptic protein that binds to cell membranes and is linked to Parkinson's disease (PD). Binding of αS to membranes is a likely first step in the molecular pathophysiology of PD. The αS molecule can adopt multiple conformations, being largely disordered in water, adopting a β-sheet conformation when present in amyloid fibrils, and forming a dynamic multiplicity of α-helical conformations when bound to lipid bilayers and related membrane-mimetic surfaces. Multiscale molecular dynamics simulations in conjunction with nuclear magnetic resonance (NMR) and cross-linking mass spectrometry (XLMS) measurements are used to explore the interactions of αS with an anionic lipid bilayer. The simulations and NMR measurements together reveal a break in the helical structure of the central non-amyloid-β component (NAC) region of αS in the vicinity of residues 65-70, which may facilitate subsequent oligomer formation. Coarse-grained simulations of αS starting from the structure of αS when bound to a detergent micelle reveal the overall pattern of protein contacts to anionic lipid bilayers, while subsequent all-atom simulations provide details of conformational changes upon membrane binding. In particular, simulations and NMR data for liposome-bound αS indicate incipient β-strand formation in the NAC region, which is supported by intramolecular contacts seen via XLMS and simulations. Markov state models based on the all-atom simulations suggest a mechanism of conformational change of membrane-bound αS via a dynamic helix break in the region of residue 65 in the NAC region. The emergent dynamic model of membrane-interacting αS advances our understanding of the mechanism of PD, potentially aiding the design of novel therapeutic approaches.
Collapse
Affiliation(s)
- Sarah-Beth
T. A. Amos
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K.
| | - Thomas C. Schwarz
- Department
of Structural and Computational
Biology, Max Perutz Laboratories, University
of Vienna, Campus Vienna
Biocenter 5, Vienna A-1030, Austria
| | - Jiye Shi
- UCB
Pharma, 208 Bath Road, Slough SL1 3WE, U.K.
| | | | | | | | - Robert Konrat
- Department
of Structural and Computational
Biology, Max Perutz Laboratories, University
of Vienna, Campus Vienna
Biocenter 5, Vienna A-1030, Austria
| | - Mark S. P. Sansom
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K.
| |
Collapse
|
15
|
Moqadam M, Tubiana T, Moutoussamy EE, Reuter N. Membrane models for molecular simulations of peripheral membrane proteins. ADVANCES IN PHYSICS: X 2021. [DOI: 10.1080/23746149.2021.1932589] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022] Open
Affiliation(s)
- Mahmoud Moqadam
- Department of Chemistry, University of Bergen, Bergen, Norway
- Computational Biology Unit, Department of Informatics, University of Bergen, Bergen, Norway
| | - Thibault Tubiana
- Department of Chemistry, University of Bergen, Bergen, Norway
- Computational Biology Unit, Department of Informatics, University of Bergen, Bergen, Norway
| | - Emmanuel E. Moutoussamy
- Computational Biology Unit, Department of Informatics, University of Bergen, Bergen, Norway
- Department of Biological Sciences, University of Bergen, Bergen, Norway
| | - Nathalie Reuter
- Department of Chemistry, University of Bergen, Bergen, Norway
- Computational Biology Unit, Department of Informatics, University of Bergen, Bergen, Norway
| |
Collapse
|
16
|
van Hilten N, Stroh KS, Risselada HJ. Membrane Thinning Induces Sorting of Lipids and the Amphipathic Lipid Packing Sensor (ALPS) Protein Motif. Front Physiol 2020; 11:250. [PMID: 32372966 PMCID: PMC7177014 DOI: 10.3389/fphys.2020.00250] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2020] [Accepted: 03/05/2020] [Indexed: 12/16/2022] Open
Abstract
Heterogeneities (e.g., membrane proteins and lipid domains) and deformations (e.g., highly curved membrane regions) in biological lipid membranes cause lipid packing defects that may trigger functional sorting of lipids and membrane-associated proteins. To study these phenomena in a controlled and efficient way within molecular simulations, we developed an external field protocol that artificially enhances packing defects in lipid membranes by enforcing local thinning of a flat membrane region. For varying lipid compositions, we observed strong thinning-induced depletion or enrichment, depending on the lipid's intrinsic shape and its effect on a membrane's elastic modulus. In particular, polyunsaturated and lysolipids are strongly attracted to regions high in packing defects, whereas phosphatidylethanolamine (PE) lipids and cholesterol are strongly repelled from it. Our results indicate that externally imposed changes in membrane thickness, area, and curvature are underpinned by shared membrane elastic principles. The observed sorting toward the thinner membrane region is in line with the sorting expected for a positively curved membrane region. Furthermore, we have demonstrated that the amphipathic lipid packing sensor (ALPS) protein motif, a known curvature and packing defect sensor, is effectively attracted to thinner membrane regions. By extracting the force that drives amphipathic molecules toward the thinner region, our thinning protocol can directly quantify and score the lipid packing sensing of different amphipathic molecules. In this way, our protocol paves the way toward high-throughput exploration of potential defect- and curvature-sensing motifs, making it a valuable addition to the molecular simulation toolbox.
Collapse
Affiliation(s)
- Niek van Hilten
- Leiden Institute of Chemistry, Leiden University, Leiden, Netherlands
| | - Kai Steffen Stroh
- Institute for Theoretical Physics, Georg August University Göttingen, Göttingen, Germany
| | - Herre Jelger Risselada
- Leiden Institute of Chemistry, Leiden University, Leiden, Netherlands.,Institute for Theoretical Physics, Georg August University Göttingen, Göttingen, Germany
| |
Collapse
|
17
|
Nepal B, Sepehri A, Lazaridis T. Mechanisms of negative membrane curvature sensing and generation by ESCRT III subunit Snf7. Protein Sci 2020; 29:1473-1485. [PMID: 32142182 DOI: 10.1002/pro.3851] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2019] [Revised: 03/03/2020] [Accepted: 03/04/2020] [Indexed: 12/13/2022]
Abstract
Certain proteins have the propensity to bind to negatively curved membranes and generate negative membrane curvature. The mechanism of action of these proteins is much less studied and understood than those that sense and generate positive curvature. In this work, we use implicit membrane modeling to explore the mechanism of an important negative curvature sensing and generating protein: the main ESCRT III subunit Snf7. We find that Snf7 monomers alone can sense negative curvature and that curvature sensitivity increases for dimers and trimers. We have observed spontaneous bending of Snf7 oligomers into circular structures with preferred radius of ~20 nm. The preferred curvature of Snf7 filaments is further confirmed by the simulations of preformed spirals on a cylindrical membrane surface. Snf7 filaments cannot bind with the same interface to flat and curved membranes. We find that even when a filament has the preferred radius, it is always less stable on the flat membrane surface than on the interior cylindrical membrane surface. This provides an additional energy for membrane bending which has not been considered in the spiral spring model. Furthermore, the rings on the cylindrical spirals are bridged together by helix 4 and hence are extra stabilized compared to the spirals on the flat membrane surface.
Collapse
Affiliation(s)
- Binod Nepal
- Department of Chemistry, City College of New York, New York, New York, USA
| | - Aliasghar Sepehri
- Department of Chemistry, City College of New York, New York, New York, USA
| | - Themis Lazaridis
- Department of Chemistry, City College of New York, New York, New York, USA.,Graduate Programs in Chemistry, Biochemistry, and Physics, The Graduate Center, City University of New York, New York, New York, USA
| |
Collapse
|
18
|
Protein Structure Prediction and Design in a Biologically Realistic Implicit Membrane. Biophys J 2020; 118:2042-2055. [PMID: 32224301 DOI: 10.1016/j.bpj.2020.03.006] [Citation(s) in RCA: 45] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2019] [Revised: 02/04/2020] [Accepted: 03/09/2020] [Indexed: 11/19/2022] Open
Abstract
Protein design is a powerful tool for elucidating mechanisms of function and engineering new therapeutics and nanotechnologies. Although soluble protein design has advanced, membrane protein design remains challenging because of difficulties in modeling the lipid bilayer. In this work, we developed an implicit approach that captures the anisotropic structure, shape of water-filled pores, and nanoscale dimensions of membranes with different lipid compositions. The model improves performance in computational benchmarks against experimental targets, including prediction of protein orientations in the bilayer, ΔΔG calculations, native structure discrimination, and native sequence recovery. When applied to de novo protein design, this approach designs sequences with an amino acid distribution near the native amino acid distribution in membrane proteins, overcoming a critical flaw in previous membrane models that were prone to generating leucine-rich designs. Furthermore, the proteins designed in the new membrane model exhibit native-like features including interfacial aromatic side chains, hydrophobic lengths compatible with bilayer thickness, and polar pores. Our method advances high-resolution membrane protein structure prediction and design toward tackling key biological questions and engineering challenges.
Collapse
|
19
|
Amphipathic Helices of Cellular Proteins Can Replace the Helix in M2 of Influenza A Virus with Only Small Effects on Virus Replication. J Virol 2020; 94:JVI.01605-19. [PMID: 31694941 DOI: 10.1128/jvi.01605-19] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2019] [Accepted: 11/04/2019] [Indexed: 11/20/2022] Open
Abstract
M2 of influenza virus functions as a proton channel during virus entry. In addition, an amphipathic helix in its cytoplasmic tail plays a role during budding. It targets M2 to the assembly site where it inserts into the inner membrane leaflet to induce curvature that causes virus scission. Since vesicularization of membranes can be performed by a variety of amphiphilic peptides, we used reverse genetics to investigate whether the peptides can substitute for M2's helix. Virus could not be generated if M2's helix was deleted or replaced by a peptide predicted not to form an amphiphilic helix. In contrast, viruses could be rescued if the M2 helix was exchanged by helices known to induce membrane curvature. Infectious virus titers were marginally reduced if M2 contains the helix of the amphipathic lipid packing sensor from the Epsin N-terminal homology domain or the nonnatural membrane inducer RW16. Transmission electron microscopy of infected cells did not reveal unequivocal evidence that virus budding or membrane scission was disturbed in any of the mutants. Instead, individual virus mutants exhibit other defects in M2, such as reduced surface expression, incorporation into virus particles, and ion channel activity. The protein composition and specific infectivity were also altered for mutant virions. We conclude that the presence of an amphiphilic helix in M2 is essential for virus replication but that other helices can replace its basic (curvature-inducing) function.IMPORTANCE Influenza virus is unique among enveloped viruses since it does not rely on the cellular ESCRT machinery for budding. Instead, viruses encode their own scission machine, the M2 protein. M2 is targeted to the edge of the viral assembly site, where it inserts an amphiphilic helix into the membrane to induce curvature. Cellular proteins utilize a similar mechanism for scission of vesicles. We show that the helix of M2 can be replaced by helices from cellular proteins with only small effects on virus replication. No evidence was obtained that budding is disturbed, but individual mutants exhibit other defects in M2 that explain the reduced virus titers. In contrast, no virus could be generated if the helix of M2 is deleted or replaced by irrelevant sequences. These experiments support the concept that M2 requires an amphiphilic helix to induce membrane curvature, but its biophysical properties are more important than the amino acid sequence.
Collapse
|
20
|
Sato Y, Kuwahara K, Mogami K, Takahashi K, Nishizawa S. Amphipathic helical peptide-based fluorogenic probes for a marker-free analysis of exosomes based on membrane-curvature sensing. RSC Adv 2020; 10:38323-38327. [PMID: 35517518 PMCID: PMC9057301 DOI: 10.1039/d0ra07763a] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2020] [Accepted: 10/11/2020] [Indexed: 01/08/2023] Open
Abstract
With increasing knowledge about the diverse roles of exosomes in the biological process, much attention has been paid to develop analytical methods for detection and quantification of exosomes. Immunoassays based on the recognition of exosomal protein markers by antibodies were widely used. However, considering that exosomal protein composition varies with the cell type, the protein markers should be carefully selected for a sensitive and selective analysis of target exosomes. Herein, we developed a new class of exosome-binding fluorogenic probes based on membrane curvature (MC) sensing of amphipathic helical (AH) peptides for exosome analysis without the need to use protein markers on the exosomal membranes. The C-terminal region of apolipoprotein A-I labeled with Nile red (ApoC-NR) exhibited a significant fluorescence enhancement upon selective binding to the highly curved membranes of synthetic vesicles. Circular dichroism (CD) measurements involving 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC)/1-2-dioleoyl-sn-glycerol (DOG) vesicles suggested that ApoC-NR recognizes the lipid packing defects in the surface of highly curved membranes via the hydrophobic insertion of the α-helix structure of the ApoC unit. ApoC-NR exhibited a stronger binding affinity for exosome-sized vesicles and a higher MC selectivity compared to all other previously reported peptide probes. ApoC-NR can be used in a simple and rapid “mix and read” analysis of various kinds of exosomes derived from different cell types (limit of detection: –105 particles/μL) without being influenced by the variation in the expression of the surface proteins of the exosomes, which stands in sharp contrast to immunoassays. Fluorogenic probes based on membrane curvature sensing-amphipathic helical peptides have been developed for a marker-free exosome analysis.![]()
Collapse
Affiliation(s)
- Yusuke Sato
- Department of Chemistry
- Graduate School of Science
- Tohoku University
- Sendai 980-8578
- Japan
| | - Kazuki Kuwahara
- Department of Chemistry
- Graduate School of Science
- Tohoku University
- Sendai 980-8578
- Japan
| | - Kenta Mogami
- Department of Chemistry
- Graduate School of Science
- Tohoku University
- Sendai 980-8578
- Japan
| | - Kenta Takahashi
- Department of Chemistry
- Graduate School of Science
- Tohoku University
- Sendai 980-8578
- Japan
| | - Seiichi Nishizawa
- Department of Chemistry
- Graduate School of Science
- Tohoku University
- Sendai 980-8578
- Japan
| |
Collapse
|
21
|
Mori T, Sugita Y. Implicit Micelle Model for Membrane Proteins Using Superellipsoid Approximation. J Chem Theory Comput 2019; 16:711-724. [DOI: 10.1021/acs.jctc.9b00783] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Takaharu Mori
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako-shi, Saitama, 351-0198, Japan
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako-shi, Saitama, 351-0198, Japan
- RIKEN Center for Computational Science, 7-1-26 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
- RIKEN Center for Biosystems Dynamics Research, 7-1-26 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
22
|
Marrink SJ, Corradi V, Souza PC, Ingólfsson HI, Tieleman DP, Sansom MS. Computational Modeling of Realistic Cell Membranes. Chem Rev 2019; 119:6184-6226. [PMID: 30623647 PMCID: PMC6509646 DOI: 10.1021/acs.chemrev.8b00460] [Citation(s) in RCA: 410] [Impact Index Per Article: 82.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2018] [Indexed: 12/15/2022]
Abstract
Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.
Collapse
Affiliation(s)
- Siewert J. Marrink
- Groningen
Biomolecular Sciences and Biotechnology Institute & Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Valentina Corradi
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Paulo C.T. Souza
- Groningen
Biomolecular Sciences and Biotechnology Institute & Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Helgi I. Ingólfsson
- Biosciences
and Biotechnology Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94550, United States
| | - D. Peter Tieleman
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Mark S.P. Sansom
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K.
| |
Collapse
|
23
|
Avci FG, Akbulut BS, Ozkirimli E. Membrane Active Peptides and Their Biophysical Characterization. Biomolecules 2018; 8:biom8030077. [PMID: 30135402 PMCID: PMC6164437 DOI: 10.3390/biom8030077] [Citation(s) in RCA: 92] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2018] [Revised: 08/08/2018] [Accepted: 08/13/2018] [Indexed: 12/12/2022] Open
Abstract
In the last 20 years, an increasing number of studies have been reported on membrane active peptides. These peptides exert their biological activity by interacting with the cell membrane, either to disrupt it and lead to cell lysis or to translocate through it to deliver cargos into the cell and reach their target. Membrane active peptides are attractive alternatives to currently used pharmaceuticals and the number of antimicrobial peptides (AMPs) and peptides designed for drug and gene delivery in the drug pipeline is increasing. Here, we focus on two most prominent classes of membrane active peptides; AMPs and cell-penetrating peptides (CPPs). Antimicrobial peptides are a group of membrane active peptides that disrupt the membrane integrity or inhibit the cellular functions of bacteria, virus, and fungi. Cell penetrating peptides are another group of membrane active peptides that mainly function as cargo-carriers even though they may also show antimicrobial activity. Biophysical techniques shed light on peptide–membrane interactions at higher resolution due to the advances in optics, image processing, and computational resources. Structural investigation of membrane active peptides in the presence of the membrane provides important clues on the effect of the membrane environment on peptide conformations. Live imaging techniques allow examination of peptide action at a single cell or single molecule level. In addition to these experimental biophysical techniques, molecular dynamics simulations provide clues on the peptide–lipid interactions and dynamics of the cell entry process at atomic detail. In this review, we summarize the recent advances in experimental and computational investigation of membrane active peptides with particular emphasis on two amphipathic membrane active peptides, the AMP melittin and the CPP pVEC.
Collapse
Affiliation(s)
- Fatma Gizem Avci
- Bioengineering Department, Marmara University, Kadikoy, 34722 Istanbul, Turkey.
| | | | - Elif Ozkirimli
- Chemical Engineering Department, Bogazici University, Bebek, 34342 Istanbul, Turkey.
| |
Collapse
|