1
|
Kiani YS, Jabeen I. Challenges of Protein-Protein Docking of the Membrane Proteins. Methods Mol Biol 2024; 2780:203-255. [PMID: 38987471 DOI: 10.1007/978-1-0716-3985-6_12] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/12/2024]
Abstract
Despite the recent advances in the determination of high-resolution membrane protein (MP) structures, the structural and functional characterization of MPs remains extremely challenging, mainly due to the hydrophobic nature, low abundance, poor expression, purification, and crystallization difficulties associated with MPs. Whereby the major challenges/hurdles for MP structure determination are associated with the expression, purification, and crystallization procedures. Although there have been significant advances in the experimental determination of MP structures, only a limited number of MP structures (approximately less than 1% of all) are available in the Protein Data Bank (PDB). Therefore, the structures of a large number of MPs still remain unresolved, which leads to the availability of widely unplumbed structural and functional information related to MPs. As a result, recent developments in the drug discovery realm and the significant biological contemplation have led to the development of several novel, low-cost, and time-efficient computational methods that overcome the limitations of experimental approaches, supplement experiments, and provide alternatives for the characterization of MPs. Whereby the fine tuning and optimizations of these computational approaches remains an ongoing endeavor.Computational methods offer a potential way for the elucidation of structural features and the augmentation of currently available MP information. However, the use of computational modeling can be extremely challenging for MPs mainly due to insufficient knowledge of (or gaps in) atomic structures of MPs. Despite the availability of numerous in silico methods for 3D structure determination the applicability of these methods to MPs remains relatively low since all methods are not well-suited or adequate for MPs. However, sophisticated methods for MP structure predictions are constantly being developed and updated to integrate the modifications required for MPs. Currently, different computational methods for (1) MP structure prediction, (2) stability analysis of MPs through molecular dynamics simulations, (3) modeling of MP complexes through docking, (4) prediction of interactions between MPs, and (5) MP interactions with its soluble partner are extensively used. Towards this end, MP docking is widely used. It is notable that the MP docking methods yet few in number might show greater potential in terms of filling the knowledge gap. In this chapter, MP docking methods and associated challenges have been reviewed to improve the applicability, accuracy, and the ability to model macromolecular complexes.
Collapse
Affiliation(s)
- Yusra Sajid Kiani
- School of Interdisciplinary Engineering and Sciences (SINES), National University of Sciences and Technology (NUST), Islamabad, Pakistan
| | - Ishrat Jabeen
- School of Interdisciplinary Engineering and Sciences (SINES), National University of Sciences and Technology (NUST), Islamabad, Pakistan.
| |
Collapse
|
2
|
Molecular Insights into Substrate Binding of the Outer Membrane Enzyme OmpT. Catalysts 2023. [DOI: 10.3390/catal13020214] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
The enzyme OmpT of the outer membrane of Escherichia coli shows proteolytic activity and cleaves peptides and proteins. Using molecular dynamics simulations in a fully hydrated lipid bilayer on a time scale of hundreds of nanoseconds, we draw a detailed atomic picture of substrate recognition in the OmpT-holo enzyme complex. Hydrogen bonds and salt bridges are essential for maintaining the integrity of the active site and play a central role for OmpT in recognizing its substrate. Electrostatic interactions are critical at all stages from approaching the substrate to docking at the active site. Computational alanine scanning based on the Molecular Mechanics Generalized Born Surface Area (MM-GBSA) approach confirms the importance of multiple residues in the active site that form salt bridges. The substrate fluctuates along the axis of the β-barrel, which is associated with oscillations of the binding cleft formed by the residue pairs D210-H212 and D83-D85. Principal component analysis suggests that substrate and protein movements are correlated. We observe the transient presence of putative catalytic water molecules near the active site, which may be involved in the nucleophilic attack on the cleavable peptide bond of the substrate.
Collapse
|
3
|
Koehler Leman J, Bonneau R. Specificities of Modeling of Membrane Proteins Using Multi-Template Homology Modeling. Methods Mol Biol 2023; 2627:141-166. [PMID: 36959446 DOI: 10.1007/978-1-0716-2974-1_8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/25/2023]
Abstract
Structures of membrane proteins are challenging to determine experimentally and currently represent only about 2% of the structures in the Protein Data Bank. Because of this disparity, methods for modeling membrane proteins are fewer and of lower quality than those for modeling soluble proteins. However, better expression, crystallization, and cryo-EM techniques have prompted a recent increase in experimental structures of membrane proteins, which can act as templates to predict the structure of closely related proteins through homology modeling. Because homology modeling relies on a structural template, it is easier and more accurate than fold recognition methods or de novo modeling, which are used when the sequence similarity between the query sequence and the sequence of related proteins in structural databases is below 25%. In homology modeling, a query sequence is mapped onto the coordinates of a single template and refined. With the increase in available templates, several templates often cover overlapping segments of the query sequence. Multi-template modeling can be used to identify the best template for local segments and join them into a single model. Here we provide a protocol for modeling membrane proteins from multiple templates in the Rosetta software suite. This approach takes advantage of several integrated frameworks, namely, RosettaScripts, RosettaCM, and RosettaMP with the membrane scoring function.
Collapse
Affiliation(s)
- Julia Koehler Leman
- Center for Computational Biology, Flatiron Institute, Simons Foundation, New York, NY, USA
| | - Richard Bonneau
- Center for Computational Biology, Flatiron Institute, Simons Foundation, New York, NY, USA
- Department of Biology, New York University, New York, NY, USA
- Department of Computer Science, New York University, New York, NY, USA
- Center for Data Science, New York University, New York, NY, USA
- Prescient Design, a Genentech Accelerator, New York, NY, USA
| |
Collapse
|
4
|
Hasanzadeh A, Hamblin MR, Kiani J, Noori H, Hardie JM, Karimi M, Shafiee H. Could artificial intelligence revolutionize the development of nanovectors for gene therapy and mRNA vaccines? NANO TODAY 2022; 47:101665. [PMID: 37034382 PMCID: PMC10081506 DOI: 10.1016/j.nantod.2022.101665] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
Gene therapy enables the introduction of nucleic acids like DNA and RNA into host cells, and is expected to revolutionize the treatment of a wide range of diseases. This growth has been further accelerated by the discovery of CRISPR/Cas technology, which allows accurate genomic editing in a broad range of cells and organisms in vitro and in vivo. Despite many advances in gene delivery and the development of various viral and non-viral gene delivery vectors, the lack of highly efficient non-viral systems with low cellular toxicity remains a challenge. The application of cutting-edge technologies such as artificial intelligence (AI) has great potential to find new paradigms to solve this issue. Herein, we review AI and its major subfields including machine learning (ML), neural networks (NNs), expert systems, deep learning (DL), computer vision and robotics. We discuss the potential of AI-based models and algorithms in the design of targeted gene delivery vehicles capable of crossing extracellular and intracellular barriers by viral mimicry strategies. We finally discuss the role of AI in improving the function of CRISPR/Cas systems, developing novel nanobots, and mRNA vaccine carriers.
Collapse
Affiliation(s)
- Akbar Hasanzadeh
- Cellular and Molecular Research Center, Iran University of Medical Sciences, Tehran 1449614535, Iran
- Department of Medical Nanotechnology, Faculty of Advanced Technologies in Medicine, Iran University of Medical Sciences, Tehran 1449614535, Iran
| | - Michael R Hamblin
- Laser Research Centre, Faculty of Health Science, University of Johannesburg, Doornfontein 2028, South Africa
- Radiation Biology Research Center, Iran University of Medical Sciences, Tehran, Iran
| | - Jafar Kiani
- Oncopathology Research Center, Iran University of Medical Sciences, Tehran 1449614535, Iran
- Department of Molecular Medicine, Faculty of Advanced Technologies in Medicine, Iran University of Medical Sciences, Tehran, Iran
| | - Hamid Noori
- Cellular and Molecular Research Center, Iran University of Medical Sciences, Tehran 1449614535, Iran
- Department of Medical Nanotechnology, Faculty of Advanced Technologies in Medicine, Iran University of Medical Sciences, Tehran 1449614535, Iran
| | - Joseph M. Hardie
- Division of Engineering in Medicine, Department of Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, MA, 02139 USA
| | - Mahdi Karimi
- Cellular and Molecular Research Center, Iran University of Medical Sciences, Tehran 1449614535, Iran
- Department of Medical Nanotechnology, Faculty of Advanced Technologies in Medicine, Iran University of Medical Sciences, Tehran 1449614535, Iran
- Oncopathology Research Center, Iran University of Medical Sciences, Tehran 1449614535, Iran
- Research Center for Science and Technology in Medicine, Tehran University of Medical Sciences, Tehran 141556559, Iran
- Applied Biotechnology Research Centre, Tehran Medical Science, Islamic Azad University, Tehran 1584743311, Iran
| | - Hadi Shafiee
- Division of Engineering in Medicine, Department of Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, MA, 02139 USA
| |
Collapse
|
5
|
Zhang J, Yuan H, Yao X, Chen S. Endogenous ion channels expressed in human embryonic kidney (HEK-293) cells. Pflugers Arch 2022; 474:665-680. [PMID: 35567642 DOI: 10.1007/s00424-022-02700-z] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2022] [Revised: 04/25/2022] [Accepted: 04/30/2022] [Indexed: 12/21/2022]
Abstract
Mammalian expression systems, particularly the human embryonic kidney (HEK-293) cells, combined with electrophysiological studies, have greatly benefited our understanding of the function, characteristic, and regulation of various ion channels. It was previously assumed that the existence of endogenous ion channels in native HEK-293 cells could be negligible. Still, more and more ion channels are gradually reported in native HEK-293 cells, which should draw our attention. In this regard, we summarize the different ion channels that are endogenously expressed in HEK-293 cells, including voltage-gated Na+ channels, Ca2+ channels, K+ channels, Cl- channels, nonselective cation channels, TRP channels, acid-sensitive ion channels, and Piezo channels, which may complicate the recording of the heterogeneously expressed ion channels to a certain degree. We noted that the expression patterns and channel profiles varied with different studies, which may be due to the distinct originality of the cells, cell culture conditions, passage numbers, and different recording protocols. Therefore, a better knowledge of endogenous ion channels may help minimize potential problems in characterizing heterologously expressed ion channels. Based on this, it is recommended that HEK-293 cells from unknown sources should be examined before transfection for the characterization of their functional profile, especially when the expression level of exogenous ion channels does not overwhelm the endogenous ion channels largely, or the current amplitude is not significantly higher than the native currents.
Collapse
Affiliation(s)
- Jun Zhang
- School of Biomedical Sciences and Li Ka Shing Institute of Health Science, The Chinese University of Hong Kong, Hong Kong, China
| | - Huikai Yuan
- Department of General Surgery, The First Affiliated Hospital of Harbin Medical University, Harbin, China
| | - Xiaoqiang Yao
- School of Biomedical Sciences and Li Ka Shing Institute of Health Science, The Chinese University of Hong Kong, Hong Kong, China
| | - Shuo Chen
- Department of Biopharmaceutical Sciences, School of Pharmacy, Harbin Medical University at Daqing, No. 39 Xinyang Rd, High-tech District, Daqing, 163319, Heilongjiang Province, China.
| |
Collapse
|
6
|
Wang XR, Cao TT, Jia CM, Tian XM, Wang Y. Quantitative prediction model for affinity of drug-target interactions based on molecular vibrations and overall system of ligand-receptor. BMC Bioinformatics 2021; 22:497. [PMID: 34649499 PMCID: PMC8515642 DOI: 10.1186/s12859-021-04389-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2021] [Accepted: 09/20/2021] [Indexed: 12/27/2022] Open
Abstract
Background The study of drug–target interactions (DTIs) affinity plays an important role in safety assessment and pharmacology. Currently, quantitative structure–activity relationship (QSAR) and molecular docking (MD) are most common methods in research of DTIs affinity. However, they often built for a specific target or several targets, and most QSAR and MD methods were based either on structure of drug molecules or on structure of receptors with low accuracy and small scope of application. How to construct quantitative prediction models with high accuracy and wide applicability remains a challenge. To this end, this paper screened molecular descriptors based on molecular vibrations and took molecule-target as a whole system to construct prediction models with high accuracy-wide applicability based on dissociation constant (Kd) and concentration for 50% of maximal effect (EC50), and to provide reference for quantifying affinity of DTIs. Results After comprehensive comparison, the results showed that RF models are optimal models to analyze and predict DTIs affinity with coefficients of determination (R2) are all greater than 0.94. Compared to the quantitative models reported in literatures, the RF models developed in this paper have higher accuracy and wide applicability. In addition, E-state molecular descriptors associated with molecular vibrations and normalized Moreau-Broto autocorrelation (G3), Moran autocorrelation (G4), transition-distribution (G7) protein descriptors are of higher importance in the quantification of DTIs. Conclusion Through screening molecular descriptors based on molecular vibrations and taking molecule-target as whole system, we obtained optimal models based on RF with more accurate-widely applicable, which indicated that selection of molecular descriptors associated with molecular vibrations and the use of molecular-target as whole system are reliable methods for improving performance of models. It can provide reference for quantifying affinity of DTIs. Supplementary Information The online version contains supplementary material available at 10.1186/s12859-021-04389-w.
Collapse
Affiliation(s)
- Xian-Rui Wang
- Key Laboratory of TCM-Information Engineer of State Administration of TCM, School of Chinese Pharmacy, Beijing University of Chinese Medicine, Beijing, 100102, China
| | - Ting-Ting Cao
- Key Laboratory of TCM-Information Engineer of State Administration of TCM, School of Chinese Pharmacy, Beijing University of Chinese Medicine, Beijing, 100102, China
| | - Cong Min Jia
- Key Laboratory of TCM-Information Engineer of State Administration of TCM, School of Chinese Pharmacy, Beijing University of Chinese Medicine, Beijing, 100102, China
| | - Xue-Mei Tian
- Key Laboratory of TCM-Information Engineer of State Administration of TCM, School of Chinese Pharmacy, Beijing University of Chinese Medicine, Beijing, 100102, China
| | - Yun Wang
- Key Laboratory of TCM-Information Engineer of State Administration of TCM, School of Chinese Pharmacy, Beijing University of Chinese Medicine, Beijing, 100102, China.
| |
Collapse
|
7
|
Vickery ON, Stansfeld PJ. CG2AT2: an Enhanced Fragment-Based Approach for Serial Multi-scale Molecular Dynamics Simulations. J Chem Theory Comput 2021; 17:6472-6482. [PMID: 34492188 PMCID: PMC8515810 DOI: 10.1021/acs.jctc.1c00295] [Citation(s) in RCA: 46] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
![]()
Coarse-grained molecular
dynamics provides a means for simulating
the assembly and interactions of macromolecular complexes at a reduced
level of representation, thereby allowing both longer timescale and
larger sized simulations. Here, we describe an enhanced fragment-based
protocol for converting macromolecular complexes from coarse-grained
to atomistic resolution, for further refinement and analysis. While
the focus is upon systems that comprise an integral membrane protein
embedded in a phospholipid bilayer, the technique is also suitable
for membrane-anchored and soluble protein/nucleotide complexes. Overall,
this provides a method for generating an accurate and well-equilibrated
atomic-level description of a macromolecular complex. The approach
is evaluated using a diverse test set of 11 system configurations
of varying size and complexity. Simulations are assessed in terms
of protein stereochemistry, conformational drift, lipid/protein interactions,
and lipid dynamics.
Collapse
Affiliation(s)
- Owen N Vickery
- School of Life Sciences & Department of Chemistry, University of Warwick, Gibbet Hill Campus, Coventry CV4 7AL, U.K
| | - Phillip J Stansfeld
- School of Life Sciences & Department of Chemistry, University of Warwick, Gibbet Hill Campus, Coventry CV4 7AL, U.K
| |
Collapse
|
8
|
Li B, Mendenhall J, Capra JA, Meiler J. A Multitask Deep-Learning Method for Predicting Membrane Associations and Secondary Structures of Proteins. J Proteome Res 2021; 20:4089-4100. [PMID: 34236204 PMCID: PMC8650144 DOI: 10.1021/acs.jproteome.1c00410] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Prediction of residue-level structural attributes and protein-level structural classes helps model protein tertiary structures and understand protein functions. Existing methods are either specialized on only one class of proteins or developed to predict only a specific type of residue-level attribute. In this work, we develop a new deep-learning method, named Membrane Association and Secondary Structure Predictor (MASSP), for accurately predicting both residue-level structural attributes (secondary structure, location, orientation, and topology) and protein-level structural classes (bitopic, α-helical, β-barrel, and soluble). MASSP integrates a multilayer two-dimensional convolutional neural network (2D-CNN) with a long short-term memory (LSTM) neural network into a multitasking framework. Our comparison shows that MASSP performs equally well or better than the state-of-the-art methods in predicting residue-level secondary structures, boundaries of transmembrane segments, and topology. Furthermore, it achieves outstanding accuracy in predicting protein-level structural classes. MASSP automatically distinguishes the structural classes of input sequences and identifies transmembrane segments and topologies if present, making it broadly applicable to different classes of proteins. In summary, MASSP's good performance and broad applicability make it well suited for annotating residue-level attributes and protein-level structural classes at the proteome scale.
Collapse
Affiliation(s)
- Bian Li
- Department of Biological Sciences, Vanderbilt University, Nashville, Tennessee 37203, United States.,Center for Structural Biology, Vanderbilt University, Nashville, Tennessee 37203, United States
| | - Jeffrey Mendenhall
- Center for Structural Biology, Vanderbilt University, Nashville, Tennessee 37203, United States.,Department of Chemistry, Vanderbilt University, Nashville, Tennessee 37203, United States
| | - John A Capra
- Bakar Computational Health Sciences Institute and Department of Epidemiology and Biostatistics, University of California, San Francisco, California 94143, United States
| | - Jens Meiler
- Center for Structural Biology, Vanderbilt University, Nashville, Tennessee 37203, United States.,Department of Chemistry, Vanderbilt University, Nashville, Tennessee 37203, United States.,Institute for Drug Discovery, University Leipzig Medical School, Leipzig 04109, Germany
| |
Collapse
|
9
|
Lachance J, Matteau D, Brodeur J, Lloyd CJ, Mih N, King ZA, Knight TF, Feist AM, Monk JM, Palsson BO, Jacques P, Rodrigue S. Genome-scale metabolic modeling reveals key features of a minimal gene set. Mol Syst Biol 2021; 17:e10099. [PMID: 34288418 PMCID: PMC8290834 DOI: 10.15252/msb.202010099] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2020] [Revised: 06/18/2021] [Accepted: 06/22/2021] [Indexed: 12/19/2022] Open
Abstract
Mesoplasma florum, a fast-growing near-minimal organism, is a compelling model to explore rational genome designs. Using sequence and structural homology, the set of metabolic functions its genome encodes was identified, allowing the reconstruction of a metabolic network representing ˜ 30% of its protein-coding genes. Growth medium simplification enabled substrate uptake and product secretion rate quantification which, along with experimental biomass composition, were integrated as species-specific constraints to produce the functional iJL208 genome-scale model (GEM) of metabolism. Genome-wide expression and essentiality datasets as well as growth data on various carbohydrates were used to validate and refine iJL208. Discrepancies between model predictions and observations were mechanistically explained using protein structures and network analysis. iJL208 was also used to propose an in silico reduced genome. Comparing this prediction to the minimal cell JCVI-syn3.0 and its parent JCVI-syn1.0 revealed key features of a minimal gene set. iJL208 is a stepping-stone toward model-driven whole-genome engineering.
Collapse
Affiliation(s)
| | - Dominick Matteau
- Département de BiologieUniversité de SherbrookeSherbrookeQCCanada
| | - Joëlle Brodeur
- Département de BiologieUniversité de SherbrookeSherbrookeQCCanada
| | - Colton J Lloyd
- Department of BioengineeringUniversity of CaliforniaSan Diego, La JollaCAUSA
| | - Nathan Mih
- Department of BioengineeringUniversity of CaliforniaSan Diego, La JollaCAUSA
| | - Zachary A King
- Department of BioengineeringUniversity of CaliforniaSan Diego, La JollaCAUSA
| | | | - Adam M Feist
- Department of BioengineeringUniversity of CaliforniaSan Diego, La JollaCAUSA
- Department of PediatricsUniversity of CaliforniaSan Diego, La JollaCAUSA
| | - Jonathan M Monk
- Department of BioengineeringUniversity of CaliforniaSan Diego, La JollaCAUSA
| | - Bernhard O Palsson
- Department of BioengineeringUniversity of CaliforniaSan Diego, La JollaCAUSA
- Department of PediatricsUniversity of CaliforniaSan Diego, La JollaCAUSA
- Bioinformatics and Systems Biology ProgramUniversity of CaliforniaSan Diego, La JollaCAUSA
- Novo Nordisk Foundation Center for BiosustainabilityTechnical University of DenmarkLyngbyDenmark
| | | | | |
Collapse
|
10
|
Jiang V, Khare SD, Banta S. Computational structure prediction provides a plausible mechanism for electron transfer by the outer membrane protein Cyc2 from Acidithiobacillus ferrooxidans. Protein Sci 2021; 30:1640-1652. [PMID: 33969560 DOI: 10.1002/pro.4106] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2021] [Revised: 04/30/2021] [Accepted: 05/03/2021] [Indexed: 12/14/2022]
Abstract
Cyc2 is the key protein in the outer membrane of Acidithiobacillus ferrooxidans that mediates electron transfer between extracellular inorganic iron and the intracellular central metabolism. This cytochrome c is specific for iron and interacts with periplasmic proteins to complete a reversible electron transport chain. A structure of Cyc2 has not yet been characterized experimentally. Here we describe a structural model of Cyc2, and associated proteins, to highlight a plausible mechanism for the ferrous iron electron transfer chain. A comparative modeling protocol specific for trans membrane beta barrel (TMBB) proteins in acidophilic conditions (pH ~ 2) was applied to the primary sequence of Cyc2. The proposed structure has three main regimes: Extracellular loops exposed to low-pH conditions, a TMBB, and an N-terminal cytochrome-like region within the periplasmic space. The Cyc2 model was further refined by identifying likely iron and heme docking sites. This represents the first computational model of Cyc2 that accounts for the membrane microenvironment and the acidity in the extracellular matrix. This approach can be used to model other TMBBs which can be critical for chemolithotrophic microbial growth.
Collapse
Affiliation(s)
- Virginia Jiang
- Department of Chemical Engineering, Columbia University in the City of New York, New York, New York, USA
| | - Sagar D Khare
- Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, Piscataway, New Jersey, USA
| | - Scott Banta
- Department of Chemical Engineering, Columbia University in the City of New York, New York, New York, USA
| |
Collapse
|
11
|
Mesdaghi S, Murphy DL, Sánchez Rodríguez F, Burgos-Mármol JJ, Rigden DJ. In silico prediction of structure and function for a large family of transmembrane proteins that includes human Tmem41b. F1000Res 2021; 9:1395. [PMID: 33520197 PMCID: PMC7818093 DOI: 10.12688/f1000research.27676.1] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Accepted: 11/23/2020] [Indexed: 01/07/2023] Open
Abstract
Background: Recent strides in computational structural biology have opened up an opportunity to understand previously uncharacterised proteins. The under-representation of transmembrane proteins in the Protein Data Bank highlights the need to apply new and advanced bioinformatics methods to shed light on their structure and function. This study focuses on a family of transmembrane proteins containing the Pfam domain PF09335 ('SNARE_ASSOC'/ 'VTT '/'Tvp38'). One prominent member, Tmem41b, has been shown to be involved in early stages of autophagosome formation and is vital in mouse embryonic development as well as being identified as a viral host factor of SARS-CoV-2. Methods: We used evolutionary covariance-derived information to construct and validate ab initio models, make domain boundary predictions and infer local structural features. Results: The results from the structural bioinformatics analysis of Tmem41b and its homologues showed that they contain a tandem repeat that is clearly visible in evolutionary covariance data but much less so by sequence analysis. Furthermore, cross-referencing of other prediction data with covariance analysis showed that the internal repeat features two-fold rotational symmetry. Ab initio modelling of Tmem41b and homologues reinforces these structural predictions. Local structural features predicted to be present in Tmem41b were also present in Cl -/H + antiporters. Conclusions: The results of this study strongly point to Tmem41b and its homologues being transporters for an as-yet uncharacterised substrate and possibly using H + antiporter activity as its mechanism for transport.
Collapse
Affiliation(s)
- Shahram Mesdaghi
- Institute of Systems, Molecular and Integrative Biology, University of Liverpool, Liverpool, L69 7ZB, UK
| | - David L. Murphy
- Institute of Systems, Molecular and Integrative Biology, University of Liverpool, Liverpool, L69 7ZB, UK
| | - Filomeno Sánchez Rodríguez
- Institute of Systems, Molecular and Integrative Biology, University of Liverpool, Liverpool, L69 7ZB, UK
| | - J. Javier Burgos-Mármol
- Institute of Systems, Molecular and Integrative Biology, University of Liverpool, Liverpool, L69 7ZB, UK
| | - Daniel J. Rigden
- Institute of Systems, Molecular and Integrative Biology, University of Liverpool, Liverpool, L69 7ZB, UK,
| |
Collapse
|
12
|
Gulsevin A, Meiler J. Prediction of amphipathic helix-membrane interactions with Rosetta. PLoS Comput Biol 2021; 17:e1008818. [PMID: 33730029 PMCID: PMC8007005 DOI: 10.1371/journal.pcbi.1008818] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2020] [Revised: 03/29/2021] [Accepted: 02/18/2021] [Indexed: 01/17/2023] Open
Abstract
Amphipathic helices have hydrophobic and hydrophilic/charged residues situated on opposite faces of the helix. They can anchor peripheral membrane proteins to the membrane, be attached to integral membrane proteins, or exist as independent peptides. Despite the widespread presence of membrane-interacting amphipathic helices, there is no computational tool within Rosetta to model their interactions with membranes. In order to address this need, we developed the AmphiScan protocol with PyRosetta, which runs a grid search to find the most favorable position of an amphipathic helix with respect to the membrane. The performance of the algorithm was tested in benchmarks with the RosettaMembrane, ref2015_memb, and franklin2019 score functions on six engineered and 44 naturally-occurring amphipathic helices using membrane coordinates from the OPM and PDBTM databases, OREMPRO server, and MD simulations for comparison. The AmphiScan protocol predicted the coordinates of amphipathic helices within less than 3Å of the reference structures and identified membrane-embedded residues with a Matthews Correlation Constant (MCC) of up to 0.57. Overall, AmphiScan stands as fast, accurate, and highly-customizable protocol that can be pipelined with other Rosetta and Python applications. Amphipathic helices are important targets as antibacterial peptides and as domains of membrane proteins that play a role in sensing the membrane environment. Understanding how amphipathic helices interact with membrane enables us to design better peptides and understand how membrane proteins use them to interact with their environment. However, there is a limited number of tools available for the modeling of amphipathic helices in membranes. Implicit membrane models can be used for this purpose as simplistic representations of the membrane environment. In this work, we developed the AmphiScan protocol that can be used to predict membrane coordinates of amphipathic helices starting with a helix structure in an implicit membrane environment. We benchmarked the performance of AmphiScan on engineered LK peptides, naturally-occurring amphipathic helices, and hydrophobic and hydrophilic peptides. Our approach provides a reliable and customizable tool to model amphipathic helix–membrane interactions, and pose a platform for the screening of amphipathic helix properties in silico.
Collapse
Affiliation(s)
- Alican Gulsevin
- Department of Chemistry, Center for Structural Biology, Vanderbilt University, Nashville, Tennessee, United States of America
| | - Jens Meiler
- Department of Chemistry, Center for Structural Biology, Vanderbilt University, Nashville, Tennessee, United States of America
- Institute for Drug Discovery, Leipzig University Medical School, 04103 Leipzig, Germany
- * E-mail:
| |
Collapse
|
13
|
Wada A, Prates ÉT, Hirano R, Werner AZ, Kamimura N, Jacobson DA, Beckham GT, Masai E. Characterization of aromatic acid/proton symporters in Pseudomonas putida KT2440 toward efficient microbial conversion of lignin-related aromatics. Metab Eng 2021; 64:167-179. [PMID: 33549838 DOI: 10.1016/j.ymben.2021.01.013] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2020] [Revised: 12/10/2020] [Accepted: 01/30/2021] [Indexed: 11/18/2022]
Abstract
Pseudomonas putida KT2440 (hereafter KT2440) is a well-studied platform bacterium for the production of industrially valuable chemicals from heterogeneous mixtures of aromatic compounds obtained from lignin depolymerization. KT2440 can grow on lignin-related monomers, such as ferulate (FA), 4-coumarate (4CA), vanillate (VA), 4-hydroxybenzoate (4HBA), and protocatechuate (PCA). Genes associated with their catabolism are known, but knowledge about the uptake systems remains limited. In this work, we studied the KT2440 transporters of lignin-related monomers and their substrate selectivity. Based on the inhibition by protonophores, we focused on five genes encoding aromatic acid/H+ symporter family transporters categorized into major facilitator superfamily that uses the proton motive force. The mutants of PP_1376 (pcaK) and PP_3349 (hcnK) exhibited significantly reduced growth on PCA/4HBA and FA/4CA, respectively, while no change was observed on VA for any of the five gene mutants. At pH 9.0, the conversion of these compounds by hcnK mutant (FA/4CA) and vanK mutant (VA) was dramatically reduced, revealing that these transporters are crucial for the uptake of the anionic substrates at high pH. Uptake assays using 14C-labeled substrates in Escherichia coli and biosensor-based assays confirmed that PcaK, HcnK, and VanK have ability to take up PCA, FA/4CA, and VA/PCA, respectively. Additionally, analyses of the predicted protein structures suggest that the size and hydropathic properties of the substrate-binding sites of these transporters determine their substrate preferences. Overall, this study reveals that at physiological pH, PcaK and HcnK have a major role in the uptake of PCA/4HBA and FA/4CA, respectively, and VanK is a VA/PCA transporter. This information can contribute to the engineering of strains for the efficient conversion of lignin-related monomers to value-added chemicals.
Collapse
Affiliation(s)
- Ayumu Wada
- Department of Bioengineering, Nagaoka University of Technology, Nagaoka, Niigata, Japan
| | - Érica T Prates
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, USA; Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, USA
| | - Ryo Hirano
- Department of Bioengineering, Nagaoka University of Technology, Nagaoka, Niigata, Japan
| | - Allison Z Werner
- Renewable Resources and Enabling Sciences Center, National Renewable Energy Laboratory, Golden, CO, 80401, USA; Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, USA
| | - Naofumi Kamimura
- Department of Bioengineering, Nagaoka University of Technology, Nagaoka, Niigata, Japan
| | - Daniel A Jacobson
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, USA; Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, USA
| | - Gregg T Beckham
- Renewable Resources and Enabling Sciences Center, National Renewable Energy Laboratory, Golden, CO, 80401, USA; Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, USA
| | - Eiji Masai
- Department of Bioengineering, Nagaoka University of Technology, Nagaoka, Niigata, Japan.
| |
Collapse
|
14
|
Lazareva EA, Lezzhov AA, Chergintsev DA, Golyshev SA, Dolja VV, Morozov SY, Heinlein M, Solovyev AG. Reticulon-like properties of a plant virus-encoded movement protein. THE NEW PHYTOLOGIST 2021; 229:1052-1066. [PMID: 32866987 DOI: 10.1111/nph.16905] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/26/2020] [Accepted: 08/20/2020] [Indexed: 06/11/2023]
Abstract
Plant viruses encode movement proteins (MPs) that ensure the transport of viral genomes through plasmodesmata (PD) and use cell endomembranes, mostly the endoplasmic reticulum (ER), for delivery of viral genomes to PD and formation of PD-anchored virus replication compartments. Here, we demonstrate that the Hibiscus green spot virus BMB2 MP, an integral ER protein, induces constrictions of ER tubules, decreases the mobility of ER luminal content, and exhibits an affinity to highly curved membranes. These properties are similar to those described for reticulons, cellular proteins that induce membrane curvature to shape the ER tubules. Similar to reticulons, BMB2 adopts a W-like topology within the ER membrane. BMB2 targets PD and increases their size exclusion limit, and these BMB2 activities correlate with the ability to induce constrictions of ER tubules. We propose that the induction of ER constrictions contributes to the BMB2-dependent increase in PD permeability and formation of the PD-associated replication compartments, therefore facilitating the virus intercellular spread. Furthermore, we show that the ER tubule constrictions also occur in cells expressing TGB2, one of the three MPs of Potato virus X (PVX), and in PVX-infected cells, suggesting that reticulon-like MPs are employed by diverse RNA viruses.
Collapse
Affiliation(s)
- Ekaterina A Lazareva
- Department of Virology, Biological Faculty, Moscow State University, Moscow, 119234, Russia
| | - Alexander A Lezzhov
- Faculty of Bioengineering and Bioinformatics, Moscow State University, Moscow, 119991, Russia
| | - Denis A Chergintsev
- Department of Plant Physiology, Biological Faculty, Moscow State University, Moscow, 119234, Russia
| | - Sergei A Golyshev
- A.N. Belozersky Institute of Physico-Chemical Biology, Moscow State University, Moscow, 119992, Russia
| | - Valerian V Dolja
- Department of Botany and Plant Pathology, Oregon State University, Corvallis, OR, 97331, USA
| | - Sergey Y Morozov
- Department of Virology, Biological Faculty, Moscow State University, Moscow, 119234, Russia
- A.N. Belozersky Institute of Physico-Chemical Biology, Moscow State University, Moscow, 119992, Russia
| | - Manfred Heinlein
- Institute for Plant Molecular Biology (IBMP-CNRS), University of Strasbourg, Strasbourg, 67000, France
| | - Andrey G Solovyev
- Department of Virology, Biological Faculty, Moscow State University, Moscow, 119234, Russia
- A.N. Belozersky Institute of Physico-Chemical Biology, Moscow State University, Moscow, 119992, Russia
- Institute of Molecular Medicine, Sechenov First Moscow State Medical University, Moscow, 119991, Russia
- Institute of Agricultural Biotechnology, Russian Academy of Agricultural Sciences, Moscow, 127550, Russia
| |
Collapse
|
15
|
Rosário-Ferreira N, Marques-Pereira C, Gouveia RP, Mourão J, Moreira IS. Guardians of the Cell: State-of-the-Art of Membrane Proteins from a Computational Point-of-View. Methods Mol Biol 2021; 2315:3-28. [PMID: 34302667 DOI: 10.1007/978-1-0716-1468-6_1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Membrane proteins (MPs) encompass a large family of proteins with distinct cellular functions, and although representing over 50% of existing pharmaceutical drug targets, their structural and functional information is still very scarce. Over the last years, in silico analysis and algorithm development were essential to characterize MPs and overcome some limitations of experimental approaches. The optimization and improvement of these methods remain an ongoing process, with key advances in MPs' structure, folding, and interface prediction being continuously tackled. Herein, we discuss the latest trends in computational methods toward a deeper understanding of the atomistic and mechanistic details of MPs.
Collapse
Affiliation(s)
- Nícia Rosário-Ferreira
- Coimbra Chemistry Center, Department of Chemistry, University of Coimbra, Coimbra, Portugal.,Center for Neuroscience and Cell Biology, Center for Innovative Biomedicine and Biotechnology, University of Coimbra, Coimbra, Portugal
| | - Catarina Marques-Pereira
- Center for Neuroscience and Cell Biology, Center for Innovative Biomedicine and Biotechnology, University of Coimbra, Coimbra, Portugal.,PhD Programme in Experimental Biology and Biomedicine, Institute for Interdisciplinary Research (IIIUC), University of Coimbra, Coimbra, Portugal
| | - Raquel P Gouveia
- Center for Neuroscience and Cell Biology, Center for Innovative Biomedicine and Biotechnology, University of Coimbra, Coimbra, Portugal
| | - Joana Mourão
- Center for Neuroscience and Cell Biology, University of Coimbra, Coimbra, Portugal
| | - Irina S Moreira
- Department of Life Sciences, University of Coimbra, Coimbra, Portugal.
| |
Collapse
|
16
|
Abstract
Protein engineering can yield new molecular tools for nanotechnology and therapeutic applications through modulating physiochemical and biological properties. Engineering membrane proteins is especially attractive because they perform key cellular processes including transport, nutrient uptake, removal of toxins, respiration, motility, and signaling. In this chapter, we describe two protocols for membrane protein engineering with the Rosetta software: (1) ΔΔG calculations for single point mutations and (2) sequence optimization in different membrane lipid compositions. These modular protocols are easily adaptable for more complex problems and serve as a foundation for efficient membrane protein engineering calculations.
Collapse
Affiliation(s)
- Rebecca F Alford
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Jeffrey J Gray
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, MD, USA. .,Program in Molecular Biophysics, Johns Hopkins University, Baltimore, MD, USA.
| |
Collapse
|
17
|
Integrative modeling of membrane-associated protein assemblies. Nat Commun 2020; 11:6210. [PMID: 33277503 PMCID: PMC7718903 DOI: 10.1038/s41467-020-20076-5] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2020] [Accepted: 11/13/2020] [Indexed: 01/03/2023] Open
Abstract
Membrane proteins are among the most challenging systems to study with experimental structural biology techniques. The increased number of deposited structures of membrane proteins has opened the route to modeling their complexes by methods such as docking. Here, we present an integrative computational protocol for the modeling of membrane-associated protein assemblies. The information encoded by the membrane is represented by artificial beads, which allow targeting of the docking toward the binding-competent regions. It combines efficient, artificial intelligence-based rigid-body docking by LightDock with a flexible final refinement with HADDOCK to remove potential clashes at the interface. We demonstrate the performance of this protocol on eighteen membrane-associated complexes, whose interface lies between the membrane and either the cytosolic or periplasmic regions. In addition, we provide a comparison to another state-of-the-art docking software, ZDOCK. This protocol should shed light on the still dark fraction of the interactome consisting of membrane proteins.
Collapse
|
18
|
Vázquez-Martínez JA, Gómez-Lim MA, Morales-Ríos E, Gonzalez-y-Merchand JA, Ortiz-Navarrete V. Short Disordered Epitope of CRTAM Ig-Like V Domain as a Potential Target for Blocking Antibodies. Int J Mol Sci 2020; 21:ijms21228798. [PMID: 33233764 PMCID: PMC7699905 DOI: 10.3390/ijms21228798] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2020] [Revised: 11/16/2020] [Accepted: 11/16/2020] [Indexed: 11/25/2022] Open
Abstract
Class-I Restricted T Cell-Associated Molecule (CRTAM) is a protein that is expressed after T cell activation. The interaction of CRTAM with its ligand, nectin-like 2 (Necl2), is required for the efficient production of IL-17, IL-22, and IFNγ by murine CD4 T cells, and it plays a role in optimal CD8 T and NK cell cytotoxicity. CRTAM promotes the pro-inflammatory cytokine profile; therefore, it may take part in the immunopathology of autoimmune diseases such as diabetes type 1 or colitis. Thus, antibodies that block the interaction between CRTAM and Necl2 would be useful for controlling the production of these inflammatory cytokines. In this work, using bioinformatics predictions, we identified three short disordered epitopes (sDE1-3) that are located in the Ig-like domains of murine CRTAM and are conserved in mammalian species. We performed a structural analysis by molecular dynamics simulations of sDE1 (QHPALKSSKY, Ig-like V), sDE2 (QRNGEKSVVK, Ig-like C1), and sDE3 (CSTERSKKPPPQI, Ig-like C1). sDE1, which is located within a loop of the contact interface of the heterotypic interaction with Nectl2, undergoes an order–disorder transition. On the contrary, even though sDE2 and sDE3 are flexible and also located within loops, they do not undergo order–disorder transitions. We evaluated the immunogenicity of sDE1 and sDE3 through the expression of these epitopes in chimeric L1 virus-like particles. We confirmed that sDE1 induces polyclonal antibodies that recognize the native folding of CRTAM expressed in activated murine CD4 T cells. In contrast, sDE3 induces polyclonal antibodies that recognize the recombinant protein hCRTAM-Fc, but not the native CRTAM. Thus, in this study, an exposed disordered epitope in the Ig-like V domain of CRTAM was identified as a potential site for therapeutic antibodies.
Collapse
Affiliation(s)
- Julio Angel Vázquez-Martínez
- Departamento de Microbiología, Escuela Nacional de Ciencias Biológicas, Instituto Politécnico Nacional, 11340 Ciudad de Mexico, Mexico; (J.A.V.-M.); (J.A.G.-y.-M.)
- Departamento de Biomedicina Molecular, Centro de Investigación y Estudios Avanzados del Instituto Politécnico Nacional (CINVESTAV), 07360 Ciudad de Mexico, Mexico
- Departamento de Ingeniería Genética, Centro de Investigación y Estudios Avanzados del Instituto Politécnico Nacional (CINVESTAV), 36821 Irapuato, Guanajuato, Mexico;
| | - Miguel Angel Gómez-Lim
- Departamento de Ingeniería Genética, Centro de Investigación y Estudios Avanzados del Instituto Politécnico Nacional (CINVESTAV), 36821 Irapuato, Guanajuato, Mexico;
| | - Edgar Morales-Ríos
- Departamento de Bioquímica, Centro de Investigación y Estudios Avanzados del Instituto Politécnico Nacional (CINVESTAV), 07360 Ciudad de Mexico, Mexico;
| | - Jorge Alberto Gonzalez-y-Merchand
- Departamento de Microbiología, Escuela Nacional de Ciencias Biológicas, Instituto Politécnico Nacional, 11340 Ciudad de Mexico, Mexico; (J.A.V.-M.); (J.A.G.-y.-M.)
| | - Vianney Ortiz-Navarrete
- Departamento de Biomedicina Molecular, Centro de Investigación y Estudios Avanzados del Instituto Politécnico Nacional (CINVESTAV), 07360 Ciudad de Mexico, Mexico
- Correspondence:
| |
Collapse
|
19
|
Determinants of Endoplasmic Reticulum-to-Lipid Droplet Protein Targeting. Dev Cell 2020; 54:471-487.e7. [PMID: 32730754 DOI: 10.1016/j.devcel.2020.07.001] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2019] [Revised: 05/11/2020] [Accepted: 07/07/2020] [Indexed: 02/06/2023]
Abstract
Lipid droplet (LD) formation from the endoplasmic reticulum (ER) is accompanied by the targeting and accumulation of specific hydrophobic, membrane-embedded proteins on LDs. The determinants of this process are unknown. Here, we study the hydrophobic membrane motifs of two Drosophila melanogaster proteins, GPAT4 and ALG14, that utilize this pathway, and we identify crucial sequence features that mediate LD accumulation. Molecular dynamics simulations and studies in cells reveal that LD targeting of these motifs requires deeply inserted tryptophans that have lower free energy in the LD oil phase and positively charged residues near predicted hairpin hinges that become less constrained in the LD environment. Analyzing hydrophobic motifs from similar LD-targeting proteins, it appears that the distribution of tryptophan and positively charged residues distinguishes them from non-LD-targeting membrane motifs. Our studies identify specific sequence features and principles of hydrophobic membrane motifs that mediate their accumulation on LDs.
Collapse
|
20
|
Leman JK, Weitzner BD, Lewis SM, Adolf-Bryfogle J, Alam N, Alford RF, Aprahamian M, Baker D, Barlow KA, Barth P, Basanta B, Bender BJ, Blacklock K, Bonet J, Boyken SE, Bradley P, Bystroff C, Conway P, Cooper S, Correia BE, Coventry B, Das R, De Jong RM, DiMaio F, Dsilva L, Dunbrack R, Ford AS, Frenz B, Fu DY, Geniesse C, Goldschmidt L, Gowthaman R, Gray JJ, Gront D, Guffy S, Horowitz S, Huang PS, Huber T, Jacobs TM, Jeliazkov JR, Johnson DK, Kappel K, Karanicolas J, Khakzad H, Khar KR, Khare SD, Khatib F, Khramushin A, King IC, Kleffner R, Koepnick B, Kortemme T, Kuenze G, Kuhlman B, Kuroda D, Labonte JW, Lai JK, Lapidoth G, Leaver-Fay A, Lindert S, Linsky T, London N, Lubin JH, Lyskov S, Maguire J, Malmström L, Marcos E, Marcu O, Marze NA, Meiler J, Moretti R, Mulligan VK, Nerli S, Norn C, Ó'Conchúir S, Ollikainen N, Ovchinnikov S, Pacella MS, Pan X, Park H, Pavlovicz RE, Pethe M, Pierce BG, Pilla KB, Raveh B, Renfrew PD, Burman SSR, Rubenstein A, Sauer MF, Scheck A, Schief W, Schueler-Furman O, Sedan Y, Sevy AM, Sgourakis NG, Shi L, Siegel JB, Silva DA, Smith S, Song Y, Stein A, Szegedy M, Teets FD, Thyme SB, Wang RYR, Watkins A, Zimmerman L, Bonneau R. Macromolecular modeling and design in Rosetta: recent methods and frameworks. Nat Methods 2020; 17:665-680. [PMID: 32483333 PMCID: PMC7603796 DOI: 10.1038/s41592-020-0848-2] [Citation(s) in RCA: 409] [Impact Index Per Article: 102.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2019] [Accepted: 04/22/2020] [Indexed: 12/12/2022]
Abstract
The Rosetta software for macromolecular modeling, docking and design is extensively used in laboratories worldwide. During two decades of development by a community of laboratories at more than 60 institutions, Rosetta has been continuously refactored and extended. Its advantages are its performance and interoperability between broad modeling capabilities. Here we review tools developed in the last 5 years, including over 80 methods. We discuss improvements to the score function, user interfaces and usability. Rosetta is available at http://www.rosettacommons.org.
Collapse
Affiliation(s)
- Julia Koehler Leman
- Center for Computational Biology, Flatiron Institute, Simons Foundation, New York, NY, USA.
- Department of Biology, New York University, New York, New York, USA.
| | - Brian D Weitzner
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, MD, USA
- Department of Biochemistry, University of Washington, Seattle, WA, USA
- Institute for Protein Design, University of Washington, Seattle, WA, USA
- Lyell Immunopharma Inc., Seattle, WA, USA
| | - Steven M Lewis
- Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
- Department of Biochemistry, Duke University, Durham, NC, USA
- Cyrus Biotechnology, Seattle, WA, USA
| | - Jared Adolf-Bryfogle
- Department of Immunology and Microbiology, The Scripps Research Institute, La Jolla, CA, USA
| | - Nawsad Alam
- Department of Microbiology and Molecular Genetics, IMRIC, Ein Kerem Faculty of Medicine, Hebrew University of Jerusalem, Jerusalem, Israel
| | - Rebecca F Alford
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Melanie Aprahamian
- Department of Chemistry and Biochemistry, Ohio State University, Columbus, OH, USA
| | - David Baker
- Department of Biochemistry, University of Washington, Seattle, WA, USA
- Institute for Protein Design, University of Washington, Seattle, WA, USA
| | - Kyle A Barlow
- Graduate Program in Bioinformatics, University of California San Francisco, San Francisco, CA, USA
| | - Patrick Barth
- Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
- Baylor College of Medicine, Department of Pharmacology, Houston, TX, USA
| | - Benjamin Basanta
- Department of Biochemistry, University of Washington, Seattle, WA, USA
- Biological Physics Structure and Design PhD Program, University of Washington, Seattle, WA, USA
| | - Brian J Bender
- Department of Pharmacology, Vanderbilt University, Nashville, TN, USA
| | - Kristin Blacklock
- Institute of Quantitative Biomedicine, Rutgers, The State University of New Jersey, Piscataway, NJ, USA
| | - Jaume Bonet
- Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Scott E Boyken
- Institute for Protein Design, University of Washington, Seattle, WA, USA
- Lyell Immunopharma Inc., Seattle, WA, USA
| | - Phil Bradley
- Fred Hutchinson Cancer Research Center, Seattle, WA, USA
| | - Chris Bystroff
- Department of Biological Sciences, Rensselaer Polytechnic Institute, Troy, NY, USA
| | - Patrick Conway
- Department of Biochemistry, University of Washington, Seattle, WA, USA
| | - Seth Cooper
- Khoury College of Computer Sciences, Northeastern University, Boston, MA, USA
| | - Bruno E Correia
- Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Brian Coventry
- Department of Biochemistry, University of Washington, Seattle, WA, USA
| | - Rhiju Das
- Department of Biochemistry, Stanford University School of Medicine, Stanford, CA, USA
| | | | - Frank DiMaio
- Department of Biochemistry, University of Washington, Seattle, WA, USA
- Institute for Protein Design, University of Washington, Seattle, WA, USA
| | - Lorna Dsilva
- Khoury College of Computer Sciences, Northeastern University, Boston, MA, USA
| | - Roland Dunbrack
- Institute for Cancer Research, Fox Chase Cancer Center, Philadelphia, PA, USA
| | - Alexander S Ford
- Department of Biochemistry, University of Washington, Seattle, WA, USA
| | - Brandon Frenz
- Institute for Protein Design, University of Washington, Seattle, WA, USA
- Cyrus Biotechnology, Seattle, WA, USA
| | - Darwin Y Fu
- Department of Chemistry, Vanderbilt University, Nashville, TN, USA
| | - Caleb Geniesse
- Department of Biochemistry, Stanford University School of Medicine, Stanford, CA, USA
| | | | - Ragul Gowthaman
- University of Maryland Institute for Bioscience and Biotechnology Research, Rockville, MD, USA
- Department of Cell Biology and Molecular Genetics, University of Maryland, College Park, MD, USA
| | - Jeffrey J Gray
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, MD, USA
- Program in Molecular Biophysics, Johns Hopkins University, Baltimore, MD, USA
| | - Dominik Gront
- Faculty of Chemistry, Biological and Chemical Research Centre, University of Warsaw, Warsaw, Poland
| | - Sharon Guffy
- Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
| | - Scott Horowitz
- Department of Chemistry & Biochemistry, University of Denver, Denver, CO, USA
- The Knoebel Institute for Healthy Aging, University of Denver, Denver, CO, USA
| | - Po-Ssu Huang
- Department of Biochemistry, University of Washington, Seattle, WA, USA
| | - Thomas Huber
- Research School of Chemistry, Australian National University, Canberra, Australian Capital Territory, Australia
| | - Tim M Jacobs
- Program in Bioinformatics and Computational Biology, Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
| | | | - David K Johnson
- Center for Computational Biology, University of Kansas, Lawrence, KS, USA
| | - Kalli Kappel
- Biophysics Program, Stanford University, Stanford, CA, USA
| | - John Karanicolas
- Institute for Cancer Research, Fox Chase Cancer Center, Philadelphia, PA, USA
| | - Hamed Khakzad
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
- Institute for Computational Science, University of Zurich, Zurich, Switzerland
- S3IT, University of Zurich, Zurich, Switzerland
| | - Karen R Khar
- Cyrus Biotechnology, Seattle, WA, USA
- Center for Computational Biology, University of Kansas, Lawrence, KS, USA
| | - Sagar D Khare
- Institute of Quantitative Biomedicine, Rutgers, The State University of New Jersey, Piscataway, NJ, USA
- Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, Piscataway, NJ, USA
- Department of Chemistry and Chemical Biology, The State University of New Jersey, Piscataway, NJ, USA
- Center for Integrative Proteomics Research, Rutgers, The State University of New Jersey, Piscataway, NJ, USA
- Computational Biology and Molecular Biophysics Program, Rutgers, The State University of New Jersey, Piscataway, NJ, USA
| | - Firas Khatib
- Department of Computer and Information Science, University of Massachusetts Dartmouth, Dartmouth, MA, USA
| | - Alisa Khramushin
- Department of Microbiology and Molecular Genetics, IMRIC, Ein Kerem Faculty of Medicine, Hebrew University of Jerusalem, Jerusalem, Israel
| | - Indigo C King
- Department of Biochemistry, University of Washington, Seattle, WA, USA
- Cyrus Biotechnology, Seattle, WA, USA
| | - Robert Kleffner
- Khoury College of Computer Sciences, Northeastern University, Boston, MA, USA
| | - Brian Koepnick
- Department of Biochemistry, University of Washington, Seattle, WA, USA
| | - Tanja Kortemme
- Department of Bioengineering and Therapeutic Sciences, University of California San Francisco, San Francisco, CA, USA
| | - Georg Kuenze
- Department of Chemistry, Vanderbilt University, Nashville, TN, USA
- Center for Structural Biology, Vanderbilt University, Nashville, TN, USA
| | - Brian Kuhlman
- Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
| | - Daisuke Kuroda
- Medical Device Development and Regulation Research Center, School of Engineering, University of Tokyo, Tokyo, Japan
- Department of Bioengineering, School of Engineering, University of Tokyo, Tokyo, Japan
| | - Jason W Labonte
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, MD, USA
- Department of Chemistry, Franklin & Marshall College, Lancaster, PA, USA
| | - Jason K Lai
- Baylor College of Medicine, Department of Pharmacology, Houston, TX, USA
| | - Gideon Lapidoth
- Department of Biomolecular Sciences, Weizmann Institute of Science, Rehovot, Israel
| | - Andrew Leaver-Fay
- Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
| | - Steffen Lindert
- Department of Chemistry and Biochemistry, Ohio State University, Columbus, OH, USA
| | - Thomas Linsky
- Department of Biochemistry, University of Washington, Seattle, WA, USA
- Institute for Protein Design, University of Washington, Seattle, WA, USA
| | - Nir London
- Department of Microbiology and Molecular Genetics, IMRIC, Ein Kerem Faculty of Medicine, Hebrew University of Jerusalem, Jerusalem, Israel
| | - Joseph H Lubin
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Sergey Lyskov
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Jack Maguire
- Program in Bioinformatics and Computational Biology, Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
| | - Lars Malmström
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
- Institute for Computational Science, University of Zurich, Zurich, Switzerland
- S3IT, University of Zurich, Zurich, Switzerland
- Division of Infection Medicine, Department of Clinical Sciences Lund, Faculty of Medicine, Lund University, Lund, Sweden
| | - Enrique Marcos
- Department of Biochemistry, University of Washington, Seattle, WA, USA
- Institute for Research in Biomedicine Barcelona, The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Orly Marcu
- Department of Microbiology and Molecular Genetics, IMRIC, Ein Kerem Faculty of Medicine, Hebrew University of Jerusalem, Jerusalem, Israel
| | - Nicholas A Marze
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Jens Meiler
- Center for Structural Biology, Vanderbilt University, Nashville, TN, USA
- Departments of Chemistry, Pharmacology and Biomedical Informatics, Vanderbilt University, Nashville, TN, USA
- Institute for Chemical Biology, Vanderbilt University, Nashville, TN, USA
| | - Rocco Moretti
- Department of Chemistry, Vanderbilt University, Nashville, TN, USA
| | - Vikram Khipple Mulligan
- Center for Computational Biology, Flatiron Institute, Simons Foundation, New York, NY, USA
- Department of Biochemistry, University of Washington, Seattle, WA, USA
- Institute for Protein Design, University of Washington, Seattle, WA, USA
| | - Santrupti Nerli
- Department of Computer Science, University of California Santa Cruz, Santa Cruz, CA, USA
| | - Christoffer Norn
- Department of Biomolecular Sciences, Weizmann Institute of Science, Rehovot, Israel
| | - Shane Ó'Conchúir
- Department of Bioengineering and Therapeutic Sciences, University of California San Francisco, San Francisco, CA, USA
| | - Noah Ollikainen
- Department of Bioengineering and Therapeutic Sciences, University of California San Francisco, San Francisco, CA, USA
| | - Sergey Ovchinnikov
- Department of Biochemistry, University of Washington, Seattle, WA, USA
- Institute for Protein Design, University of Washington, Seattle, WA, USA
- Molecular and Cellular Biology Program, University of Washington, Seattle, WA, USA
| | - Michael S Pacella
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Xingjie Pan
- Department of Bioengineering and Therapeutic Sciences, University of California San Francisco, San Francisco, CA, USA
| | - Hahnbeom Park
- Department of Biochemistry, University of Washington, Seattle, WA, USA
| | - Ryan E Pavlovicz
- Department of Biochemistry, University of Washington, Seattle, WA, USA
- Institute for Protein Design, University of Washington, Seattle, WA, USA
- Cyrus Biotechnology, Seattle, WA, USA
| | - Manasi Pethe
- Department of Chemistry and Chemical Biology, The State University of New Jersey, Piscataway, NJ, USA
- Center for Integrative Proteomics Research, Rutgers, The State University of New Jersey, Piscataway, NJ, USA
| | - Brian G Pierce
- University of Maryland Institute for Bioscience and Biotechnology Research, Rockville, MD, USA
- Department of Cell Biology and Molecular Genetics, University of Maryland, College Park, MD, USA
| | - Kala Bharath Pilla
- Research School of Chemistry, Australian National University, Canberra, Australian Capital Territory, Australia
| | - Barak Raveh
- Department of Microbiology and Molecular Genetics, IMRIC, Ein Kerem Faculty of Medicine, Hebrew University of Jerusalem, Jerusalem, Israel
| | - P Douglas Renfrew
- Center for Computational Biology, Flatiron Institute, Simons Foundation, New York, NY, USA
| | - Shourya S Roy Burman
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Aliza Rubenstein
- Institute of Quantitative Biomedicine, Rutgers, The State University of New Jersey, Piscataway, NJ, USA
- Computational Biology and Molecular Biophysics Program, Rutgers, The State University of New Jersey, Piscataway, NJ, USA
| | - Marion F Sauer
- Chemical and Physical Biology Program, Vanderbilt Vaccine Center, Vanderbilt University, Nashville, TN, USA
| | - Andreas Scheck
- Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - William Schief
- Department of Immunology and Microbiology, The Scripps Research Institute, La Jolla, CA, USA
| | - Ora Schueler-Furman
- Department of Microbiology and Molecular Genetics, IMRIC, Ein Kerem Faculty of Medicine, Hebrew University of Jerusalem, Jerusalem, Israel
| | - Yuval Sedan
- Department of Microbiology and Molecular Genetics, IMRIC, Ein Kerem Faculty of Medicine, Hebrew University of Jerusalem, Jerusalem, Israel
| | - Alexander M Sevy
- Chemical and Physical Biology Program, Vanderbilt Vaccine Center, Vanderbilt University, Nashville, TN, USA
| | - Nikolaos G Sgourakis
- Department of Chemistry and Biochemistry, University of California Santa Cruz, Santa Cruz, CA, USA
| | - Lei Shi
- Department of Biochemistry, University of Washington, Seattle, WA, USA
- Institute for Protein Design, University of Washington, Seattle, WA, USA
| | - Justin B Siegel
- Department of Chemistry, University of California, Davis, Davis, CA, USA
- Department of Biochemistry and Molecular Medicine, University of California, Davis, Davis, California, USA
- Genome Center, University of California, Davis, Davis, CA, USA
| | | | - Shannon Smith
- Department of Chemistry, Vanderbilt University, Nashville, TN, USA
| | - Yifan Song
- Department of Biochemistry, University of Washington, Seattle, WA, USA
- Institute for Protein Design, University of Washington, Seattle, WA, USA
- Cyrus Biotechnology, Seattle, WA, USA
| | - Amelie Stein
- Department of Bioengineering and Therapeutic Sciences, University of California San Francisco, San Francisco, CA, USA
| | - Maria Szegedy
- Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, Piscataway, NJ, USA
| | - Frank D Teets
- Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
| | - Summer B Thyme
- Department of Biochemistry, University of Washington, Seattle, WA, USA
| | - Ray Yu-Ruei Wang
- Department of Biochemistry, University of Washington, Seattle, WA, USA
| | - Andrew Watkins
- Department of Biochemistry, Stanford University School of Medicine, Stanford, CA, USA
| | - Lior Zimmerman
- Department of Microbiology and Molecular Genetics, IMRIC, Ein Kerem Faculty of Medicine, Hebrew University of Jerusalem, Jerusalem, Israel
| | - Richard Bonneau
- Center for Computational Biology, Flatiron Institute, Simons Foundation, New York, NY, USA.
- Department of Biology, New York University, New York, New York, USA.
- Department of Computer Science, New York University, New York, NY, USA.
- Center for Data Science, New York University, New York, NY, USA.
| |
Collapse
|
21
|
Finding the Keys to the CAR: Identifying Novel Target Antigens for T Cell Redirection Immunotherapies. Int J Mol Sci 2020; 21:ijms21020515. [PMID: 31947597 PMCID: PMC7014258 DOI: 10.3390/ijms21020515] [Citation(s) in RCA: 48] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2019] [Revised: 01/08/2020] [Accepted: 01/09/2020] [Indexed: 02/06/2023] Open
Abstract
Oncology immunotherapy has been a significant advancement in cancer treatment and involves harnessing and redirecting a patient’s immune response towards their own tumour. Specific recognition and elimination of tumour cells was first proposed over a century ago with Paul Erlich’s ‘magic bullet’ theory of therapy. In the past decades, targeting cancer antigens by redirecting T cells with antibodies using either bispecific T cell engagers (BiTEs) or chimeric antigen receptor (CAR) T cell therapy has achieved impressive clinical responses. Despite recent successes in haematological cancers, linked to a high and uniformly expressed CD19 antigen, the efficacy of T cell therapies in solid cancers has been disappointing, in part due to antigen escape. Targeting heterogeneous solid tumours with T cell therapies will require the identification of novel tumour specific targets. These targets can be found among a range of cell-surface expressed antigens, including proteins, glycolipids or carbohydrates. In this review, we will introduce the current tumour target antigen classification, outline existing approaches to discover novel tumour target antigens and discuss considerations for future design of antibodies with a focus on their use in CAR T cells.
Collapse
|
22
|
Altan I, Khan AR, James S, Quinn MK, McManus JJ, Charbonneau P. Using Schematic Models to Understand the Microscopic Basis for Inverted Solubility in γD-Crystallin. J Phys Chem B 2019; 123:10061-10072. [PMID: 31557434 DOI: 10.1021/acs.jpcb.9b07774] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Inverted solubility-melting a crystal by cooling-is observed in a handful of proteins, such as carbomonoxy hemoglobin C and γD-crystallin. In human γD-crystallin, the phenomenon is associated with the mutation of the 23rd residue, a proline, to a threonine, serine, or valine. One proposed microscopic mechanism entails an increase in surface hydrophobicity upon mutagenesis. Recent crystal structures of a double mutant that includes the P23T mutation allow for a more careful investigation of this proposal. Here, we first measure the surface hydrophobicity of various mutant structures of γD-crystallin and discern no notable increase in hydrophobicity upon mutating the 23rd residue. We then investigate the solubility inversion regime with a schematic patchy particle model that includes one of three variants of temperature-dependent patch energies: two of the hydrophobic effect, and one of a more generic nature. We conclude that, while solubility inversion due to the hydrophobic effect may be possible, microscopic evidence to support it in γD-crystallin is weak. More generally, we find that solubility inversion requires a fine balance between patch strengths and their temperature-dependent component, which may explain why inverted solubility is not commonly observed in proteins. We also find that the temperature-dependent interaction has only a negligible impact on liquid-liquid phase boundaries of γD-crystallin, in line with previous experimental observations.
Collapse
Affiliation(s)
| | - Amir R Khan
- School of Biochemistry and Immunology , Trinity College Dublin , Dublin , Ireland
| | - Susan James
- Department of Chemistry , Maynooth University , Maynooth , Ireland
| | - Michelle K Quinn
- Department of Chemistry , Maynooth University , Maynooth , Ireland
| | | | | |
Collapse
|
23
|
Differences in protein structural regions that impact functional specificity in GT2 family β-glucan synthases. PLoS One 2019; 14:e0224442. [PMID: 31665152 PMCID: PMC6821405 DOI: 10.1371/journal.pone.0224442] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2019] [Accepted: 10/14/2019] [Indexed: 12/16/2022] Open
Abstract
Most cell wall and secreted β-glucans are synthesised by the CAZy Glycosyltransferase 2 family (www.cazy.org), with different members catalysing the formation of (1,4)-β-, (1,3)-β-, or both (1,4)- and (1,3)-β-glucosidic linkages. Given the distinct physicochemical properties of each of the resultant β-glucans (cellulose, curdlan, and mixed linkage glucan, respectively) are crucial to their biological and biotechnological functions, there is a desire to understand the molecular evolution of synthesis and how linkage specificity is determined. With structural studies hamstrung by the instability of these proteins to solubilisation, we have utilised in silico techniques and the crystal structure for a bacterial cellulose synthase to further understand how these enzymes have evolved distinct functions. Sequence and phylogenetic analyses were performed to determine amino acid conservation, both family-wide and within each sub-family. Further structural analysis centred on comparison of a bacterial curdlan synthase homology model with the bacterial cellulose synthase crystal structure, with molecular dynamics simulations performed with their respective β-glucan products bound in the trans-membrane channel. Key residues that differentially interact with the different β-glucan chains and have sub-family-specific conservation were found to reside at the entrance of the trans-membrane channel. The linkage-specific catalytic activity of these enzymes and hence the type of β-glucan chain built is thus likely determined by the different interactions between the proteins and the first few glucose residues in the channel, which in turn dictates the position of the acceptor glucose. The sequence-function relationships for the bacterial β-glucan synthases pave the way for extending this understanding to other kingdoms, such as plants.
Collapse
|
24
|
Weinstein JY, Elazar A, Fleishman SJ. A lipophilicity-based energy function for membrane-protein modelling and design. PLoS Comput Biol 2019; 15:e1007318. [PMID: 31461441 PMCID: PMC6736313 DOI: 10.1371/journal.pcbi.1007318] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2019] [Revised: 09/10/2019] [Accepted: 08/01/2019] [Indexed: 01/14/2023] Open
Abstract
Membrane-protein design is an exciting and increasingly successful research area which has led to landmarks including the design of stable and accurate membrane-integral proteins based on coiled-coil motifs. Design of topologically more complex proteins, such as most receptors, channels, and transporters, however, demands an energy function that balances contributions from intra-protein contacts and protein-membrane interactions. Recent advances in water-soluble all-atom energy functions have increased the accuracy in structure-prediction benchmarks. The plasma membrane, however, imposes different physical constraints on protein solvation. To understand these constraints, we recently developed a high-throughput experimental screen, called dsTβL, and inferred apparent insertion energies for each amino acid at dozens of positions across the bacterial plasma membrane. Here, we express these profiles as lipophilicity energy terms in Rosetta and demonstrate that the new energy function outperforms previous ones in modelling and design benchmarks. Rosetta ab initio simulations starting from an extended chain recapitulate two-thirds of the experimentally determined structures of membrane-spanning homo-oligomers with <2.5Å root-mean-square deviation within the top-predicted five models (available online: http://tmhop.weizmann.ac.il). Furthermore, in two sequence-design benchmarks, the energy function improves discrimination of stabilizing point mutations and recapitulates natural membrane-protein sequences of known structure, thereby recommending this new energy function for membrane-protein modelling and design.
Collapse
Affiliation(s)
| | - Assaf Elazar
- Department of Biomolecular Sciences, Weizmann Institute of Science, Rehovot, Israel
| | - Sarel Jacob Fleishman
- Department of Biomolecular Sciences, Weizmann Institute of Science, Rehovot, Israel
- * E-mail:
| |
Collapse
|
25
|
Lapin J, Nevzorov AA. Validation of protein backbone structures calculated from NMR angular restraints using Rosetta. JOURNAL OF BIOMOLECULAR NMR 2019; 73:229-244. [PMID: 31076969 DOI: 10.1007/s10858-019-00251-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/24/2019] [Accepted: 05/02/2019] [Indexed: 06/09/2023]
Abstract
Multidimensional solid-state NMR spectra of oriented membrane proteins can be used to infer the backbone torsion angles and hence the overall protein fold by measuring dipolar couplings and chemical shift anisotropies, which depend on the orientation of each peptide plane with respect to the external magnetic field. However, multiple peptide plane orientations can be consistent with a given set of angular restraints. This ambiguity is further exacerbated by experimental uncertainty in obtaining and interpreting such restraints. The previously developed algorithms for structure calculations using angular restraints typically involve a sequential walkthrough along the backbone to find the torsion angles between the consecutive peptide plane orientations that are consistent with the experimental data. This method is sensitive to experimental uncertainty in interpreting the peak positions of as low as ± 10 Hz, often yielding high structural RMSDs for the calculated structures. Here we present a significantly improved version of the algorithm which includes the fitting of several peptide planes at once in order to prevent propagation of error along the backbone. In addition, a protocol has been devised for filtering the structural solutions using Rosetta scoring functions in order to find the structures that both fit the spectrum and satisfy bioinformatics restraints. The robustness of the new algorithm has been tested using synthetic angular restraints generated from the known structures for two proteins: a soluble protein 2gb1 (56 residues), chosen for its diverse secondary structure elements, i.e. an alpha-helix and two beta-sheets, and a membrane protein 4a2n, from which the first two transmembrane helices (having a total of 64 residues) have been used. Extensive simulations have been performed by varying the number of fitted planes, experimental error, and the number of NMR dimensions. It has been found that simultaneously fitting two peptide planes always shifted the distribution of the calculated structures toward lower structural RMSD values as compared to fitting a single torsion-angle pair. For each protein, irrespective of the simulation parameters, Rosetta was able to distinguish the most plausible structures, often having structural RMSDs lower than 2 Å with respect to the original structure. This study establishes a framework for de-novo protein structure prediction using a combination of solid-state NMR angular restraints and bioinformatics.
Collapse
Affiliation(s)
- Joel Lapin
- Department of Chemistry, North Carolina State University, 2620 Yarbrough Drive, Raleigh, NC, 27695-8204, USA
| | - Alexander A Nevzorov
- Department of Chemistry, North Carolina State University, 2620 Yarbrough Drive, Raleigh, NC, 27695-8204, USA.
| |
Collapse
|
26
|
In silico prediction of prolactin molecules as a tool for equine genomics reproduction. Mol Divers 2019; 23:1019-1028. [PMID: 30740642 DOI: 10.1007/s11030-018-09914-3] [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: 11/06/2018] [Accepted: 12/31/2018] [Indexed: 10/27/2022]
Abstract
The prolactin hormone is involved in several biological functions, although its main role resides on reproduction. As it interferes on fertility changes, studies focused on human health have established a linkage of this hormone to fertility losses. Regarding animal research, there is still a lack of information about the structure of prolactin. In case of horse breeding, prolactin has a particular influence; once there is an individualization of these animals and equines are known for presenting several reproductive disorders. As there is no molecular structure available for the prolactin hormone and receptor, we performed several bioinformatics analyses through prediction and refinement softwares, as well as manual modifications. Aiming to elucidate the first computational structure of both molecules and analyse structural and functional aspects related to these proteins, here we provide the first known equine model for prolactin and prolactin receptor, which obtained high global quality scores in diverse software's for quality assessment. QMEAN overall score obtained for ePrl was (- 4.09) and QMEANbrane for ePrlr was (- 8.45), which proves the structures' reliability. This study will implement another tool in equine genomics in order to give light to interactions of these molecules, structural and functional alterations and therefore help diagnosing fertility problems, contributing in the selection of a high genetic herd.
Collapse
|
27
|
Abstract
The pentameric γ-aminobutyric acid type A receptors are ion channels activated by ligands, which intervene in the rapid inhibitory transmission in the mammalian CNS. Due to their rich pharmacology and therapeutic potential, it is essential to understand their structure and function thoroughly. This deep characterization was hampered by the lack of experimental structural information for many years. Thus, computational techniques have been extensively combined with experimental data, in order to undertake the study of γ-aminobutyric acid type A receptors and their interaction with drugs. Here, we review the exciting journey made to assess the structures of these receptors and outline major outcomes. Finally, we discuss the brand new structure of the α1β2γ2 subtype and the amazing advances it brings to the field.
Collapse
|
28
|
MacCarthy E, Perry D, Kc DB. Advances in Protein Super-Secondary Structure Prediction and Application to Protein Structure Prediction. Methods Mol Biol 2019; 1958:15-45. [PMID: 30945212 DOI: 10.1007/978-1-4939-9161-7_2] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Due to the advancement in various sequencing technologies, the gap between the number of protein sequences and the number of experimental protein structures is ever increasing. Community-wide initiatives like CASP have resulted in considerable efforts in the development of computational methods to accurately model protein structures from sequences. Sequence-based prediction of super-secondary structure has direct application in protein structure prediction, and there have been significant efforts in the prediction of super-secondary structure in the last decade. In this chapter, we first introduce the protein structure prediction problem and highlight some of the important progress in the field of protein structure prediction. Next, we discuss recent methods for the prediction of super-secondary structures. Finally, we discuss applications of super-secondary structure prediction in structure prediction/analysis of proteins. We also discuss prediction of protein structures that are composed of simple super-secondary structure repeats and protein structures that are composed of complex super-secondary structure repeats. Finally, we also discuss the recent trends in the field.
Collapse
Affiliation(s)
- Elijah MacCarthy
- Department of Computational Science and Engineering, North Carolina A&T State University, Greensboro, NC, USA
| | - Derrick Perry
- Department of Computational Science and Engineering, North Carolina A&T State University, Greensboro, NC, USA
| | - Dukka B Kc
- Department of Computational Science and Engineering, North Carolina A&T State University, Greensboro, NC, USA.
| |
Collapse
|
29
|
DeMarco KR, Bekker S, Vorobyov I. Challenges and advances in atomistic simulations of potassium and sodium ion channel gating and permeation. J Physiol 2018; 597:679-698. [PMID: 30471114 DOI: 10.1113/jp277088] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2018] [Accepted: 10/15/2018] [Indexed: 12/19/2022] Open
Abstract
Ion channels are implicated in many essential physiological events such as electrical signal propagation and cellular communication. The advent of K+ and Na+ ion channel structure determination has facilitated numerous investigations of molecular determinants of their behaviour. At the same time, rapid development of computer hardware and molecular simulation methodologies has made computational studies of large biological molecules in all-atom representation tractable. The concurrent evolution of experimental structural biology with biomolecular computer modelling has yielded mechanistic details of fundamental processes unavailable through experiments alone, such as ion conduction and ion channel gating. This review is a short survey of the atomistic computational investigations of K+ and Na+ ion channels, focusing on KcsA and several voltage-gated channels from the KV and NaV families, which have garnered many successes and engendered several long-standing controversies regarding the nature of their structure-function relationship. We review the latest advancements and challenges facing the field of molecular modelling and simulation regarding the structural and energetic determinants of ion channel function and their agreement with experimental observations.
Collapse
Affiliation(s)
- Kevin R DeMarco
- Department of Physiology and Membrane Biology, School of Medicine, University of California, Davis, CA, USA.,Department of Pharmacology, School of Medicine, University of California, Davis, CA, USA
| | - Slava Bekker
- Department of Physiology and Membrane Biology, School of Medicine, University of California, Davis, CA, USA.,Chemistry Department, American River College, Sacramento, CA, USA
| | - Igor Vorobyov
- Department of Physiology and Membrane Biology, School of Medicine, University of California, Davis, CA, USA.,Department of Pharmacology, School of Medicine, University of California, Davis, CA, USA
| |
Collapse
|
30
|
Wang Z, Jumper JM, Wang S, Freed KF, Sosnick TR. A Membrane Burial Potential with H-Bonds and Applications to Curved Membranes and Fast Simulations. Biophys J 2018; 115:1872-1884. [PMID: 30413241 DOI: 10.1016/j.bpj.2018.10.012] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2018] [Revised: 09/21/2018] [Accepted: 10/10/2018] [Indexed: 10/28/2022] Open
Abstract
We use the statistics of a large and curated training set of transmembrane helical proteins to develop a knowledge-based potential that accounts for the dependence on both the depth of burial of the protein in the membrane and the degree of side-chain exposure. Additionally, the statistical potential includes depth-dependent energies for unsatisfied backbone hydrogen bond donors and acceptors, which are found to be relatively small, ∼2 RT. Our potential accurately places known proteins within the bilayer. The potential is applied to the mechanosensing MscL channel in membranes of varying thickness and curvature, as well as to the prediction of protein structure. The potential is incorporated into our new Upside molecular dynamics algorithm. Notably, we account for the exchange of protein-lipid interactions for protein-protein interactions as helices contact each other, thereby avoiding overestimating the energetics of helix association within the membrane. Simulations of most multimeric complexes find that isolated monomers and the oligomers retain the same orientation in the membrane, suggesting that the assembly of prepositioned monomers presents a viable mechanism of oligomerization.
Collapse
Affiliation(s)
- Zongan Wang
- Department of Chemistry, The University of Chicago, Chicago, Illinois; James Franck Institute, The University of Chicago, Chicago, Illinois
| | - John M Jumper
- Department of Chemistry, The University of Chicago, Chicago, Illinois; James Franck Institute, The University of Chicago, Chicago, Illinois; Department of Biochemistry and Molecular Biology, The University of Chicago, Chicago, Illinois
| | - Sheng Wang
- Computational Bioscience Research Center, King Abdullah University of Science and Technology, Thuwal, Saudi Arabia; Toyota Technological Institute at Chicago, Chicago, Illinois
| | - Karl F Freed
- Department of Chemistry, The University of Chicago, Chicago, Illinois; James Franck Institute, The University of Chicago, Chicago, Illinois.
| | - Tobin R Sosnick
- Department of Biochemistry and Molecular Biology, The University of Chicago, Chicago, Illinois; Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois.
| |
Collapse
|
31
|
Nikolaev D, Shtyrov AA, Panov MS, Jamal A, Chakchir OB, Kochemirovsky VA, Olivucci M, Ryazantsev MN. A Comparative Study of Modern Homology Modeling Algorithms for Rhodopsin Structure Prediction. ACS OMEGA 2018; 3:7555-7566. [PMID: 30087916 PMCID: PMC6068592 DOI: 10.1021/acsomega.8b00721] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/13/2018] [Accepted: 06/21/2018] [Indexed: 06/08/2023]
Abstract
Rhodopsins are seven α-helical membrane proteins that are of great importance in chemistry, biology, and modern biotechnology. Any in silico study on rhodopsin properties and functioning requires a high-quality three-dimensional structure. Due to particular difficulties with obtaining membrane protein structures from the experiment, in silico prediction of the three-dimensional rhodopsin structure based only on its primary sequence is an especially important task. For the last few years, significant progress was made in the field of protein structure prediction, especially for methods based on comparative modeling. However, the majority of this progress was made for soluble proteins and further investigations are needed to achieve similar progress for membrane proteins. In this paper, we evaluate the performance of modern protein structure prediction methodologies (implemented in the Medeller, I-TASSER, and Rosetta packages) for their ability to predict rhodopsin structures. Three widely used methodologies were considered: two general methodologies that are commonly applied to soluble proteins and a methodology that uses constraints that are specific for membrane proteins. The test pool consisted of 36 target-template pairs with different sequence similarities that was constructed on the basis of 24 experimental rhodopsin structures taken from the RCSB database. As a result, we showed that all three considered methodologies allow obtaining rhodopsin structures with the quality that is close to the crystallographic one (root mean square deviation (RMSD) of the predicted structure from the corresponding X-ray structure up to 1.5 Å) if the target-template sequence identity is higher than 40%. Moreover, all considered methodologies provided structures of average quality (RMSD < 4.0 Å) if the target-template sequence identity is higher than 20%. Such structures can be subsequently used for further investigation of molecular mechanisms of protein functioning and for the development of modern protein-based biotechnologies.
Collapse
Affiliation(s)
- Dmitrii
M. Nikolaev
- Nanotechnology
Research and Education Centre RAS, Saint-Petersburg
Academic University, 8/3 Khlopina Street, St. Petersburg 194021, Russia
| | - Andrey A. Shtyrov
- Nanotechnology
Research and Education Centre RAS, Saint-Petersburg
Academic University, 8/3 Khlopina Street, St. Petersburg 194021, Russia
| | - Maxim S. Panov
- Institute
of Chemistry, Saint Petersburg State University, 7/9 Universitetskaya emb., St. Petersburg 199034, Russia
| | - Adeel Jamal
- Department
of Chemical Engineering, Massachusetts Institute
of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Oleg B. Chakchir
- Nanotechnology
Research and Education Centre RAS, Saint-Petersburg
Academic University, 8/3 Khlopina Street, St. Petersburg 194021, Russia
| | - Vladimir A. Kochemirovsky
- Institute
of Chemistry, Saint Petersburg State University, 7/9 Universitetskaya emb., St. Petersburg 199034, Russia
| | - Massimo Olivucci
- Department
of Biotechnology, Chemistry and Pharmacy, Università di Siena, via A. Moro 2, Siena I-53100, Italy
| | - Mikhail N. Ryazantsev
- Institute
of Chemistry, Saint Petersburg State University, 7/9 Universitetskaya emb., St. Petersburg 199034, Russia
- Institute
of Macromolecular Compounds of the Russian Academy of Sciences, 31 Bolshoy pr., St. Petersburg 199004, Russia
| |
Collapse
|
32
|
Dutagaci B, Heo L, Feig M. Structure refinement of membrane proteins via molecular dynamics simulations. Proteins 2018; 86:738-750. [PMID: 29675899 PMCID: PMC6013386 DOI: 10.1002/prot.25508] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2018] [Revised: 04/09/2018] [Accepted: 04/14/2018] [Indexed: 12/12/2022]
Abstract
A refinement protocol based on physics-based techniques established for water soluble proteins is tested for membrane protein structures. Initial structures were generated by homology modeling and sampled via molecular dynamics simulations in explicit lipid bilayer and aqueous solvent systems. Snapshots from the simulations were selected based on scoring with either knowledge-based or implicit membrane-based scoring functions and averaged to obtain refined models. The protocol resulted in consistent and significant refinement of the membrane protein structures similar to the performance of refinement methods for soluble proteins. Refinement success was similar between sampling in the presence of lipid bilayers and aqueous solvent but the presence of lipid bilayers may benefit the improvement of lipid-facing residues. Scoring with knowledge-based functions (DFIRE and RWplus) was found to be as good as scoring using implicit membrane-based scoring functions suggesting that differences in internal packing is more important than orientations relative to the membrane during the refinement of membrane protein homology models.
Collapse
Affiliation(s)
- Bercem Dutagaci
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI, USA
| | - Lim Heo
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI, USA
| | - Michael Feig
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI, USA
| |
Collapse
|
33
|
Corradi V, Mendez-Villuendas E, Ingólfsson HI, Gu RX, Siuda I, Melo MN, Moussatova A, DeGagné LJ, Sejdiu BI, Singh G, Wassenaar TA, Delgado Magnero K, Marrink SJ, Tieleman DP. Lipid-Protein Interactions Are Unique Fingerprints for Membrane Proteins. ACS CENTRAL SCIENCE 2018; 4:709-717. [PMID: 29974066 PMCID: PMC6028153 DOI: 10.1021/acscentsci.8b00143] [Citation(s) in RCA: 205] [Impact Index Per Article: 34.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/06/2018] [Indexed: 05/08/2023]
Abstract
Cell membranes contain hundreds of different proteins and lipids in an asymmetric arrangement. Our current understanding of the detailed organization of cell membranes remains rather elusive, because of the challenge to study fluctuating nanoscale assemblies of lipids and proteins with the required spatiotemporal resolution. Here, we use molecular dynamics simulations to characterize the lipid environment of 10 different membrane proteins. To provide a realistic lipid environment, the proteins are embedded in a model plasma membrane, where more than 60 lipid species are represented, asymmetrically distributed between the leaflets. The simulations detail how each protein modulates its local lipid environment in a unique way, through enrichment or depletion of specific lipid components, resulting in thickness and curvature gradients. Our results provide a molecular glimpse of the complexity of lipid-protein interactions, with potentially far-reaching implications for our understanding of the overall organization of real cell membranes.
Collapse
Affiliation(s)
- Valentina Corradi
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Eduardo Mendez-Villuendas
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Helgi I. Ingólfsson
- Groningen
Biomolecular Sciences and Biotechnology Institute and Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Ruo-Xu Gu
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Iwona Siuda
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Manuel N. Melo
- Groningen
Biomolecular Sciences and Biotechnology Institute and Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Anastassiia Moussatova
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Lucien J. DeGagné
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Besian I. Sejdiu
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Gurpreet Singh
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Tsjerk A. Wassenaar
- Groningen
Biomolecular Sciences and Biotechnology Institute and Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Karelia Delgado Magnero
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Siewert J. Marrink
- Groningen
Biomolecular Sciences and Biotechnology Institute and Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - D. Peter Tieleman
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
- E-mail:
| |
Collapse
|
34
|
Folding a viral peptide in different membrane environments: pathway and sampling analyses. J Biol Phys 2018; 44:195-209. [PMID: 29644513 DOI: 10.1007/s10867-018-9490-y] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2017] [Accepted: 03/16/2018] [Indexed: 10/17/2022] Open
Abstract
Flock House virus (FHV) is a well-characterized model system to study infection mechanisms in non-enveloped viruses. A key stage of the infection cycle is the disruption of the endosomal membrane by a component of the FHV capsid, the membrane active γ peptide. In this study, we perform all-atom molecular dynamics simulations of the 21 N-terminal residues of the γ peptide interacting with membranes of differing compositions. We carry out umbrella sampling calculations to study the folding of the peptide to a helical state in homogenous and heterogeneous membranes consisting of neutral and anionic lipids. From the trajectory data, we evaluate folding energetics and dissect the mechanism of folding in the different membrane environments. We conclude the study by analyzing the extent of configurational sampling by performing time-lagged independent component analysis.
Collapse
|
35
|
Xia Y, Fischer AW, Teixeira P, Weiner B, Meiler J. Integrated Structural Biology for α-Helical Membrane Protein Structure Determination. Structure 2018; 26:657-666.e2. [PMID: 29526436 PMCID: PMC5884713 DOI: 10.1016/j.str.2018.02.006] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2017] [Revised: 06/14/2017] [Accepted: 02/05/2018] [Indexed: 01/12/2023]
Abstract
While great progress has been made, only 10% of the nearly 1,000 integral, α-helical, multi-span membrane protein families are represented by at least one experimentally determined structure in the PDB. Previously, we developed the algorithm BCL::MP-Fold, which samples the large conformational space of membrane proteins de novo by assembling predicted secondary structure elements guided by knowledge-based potentials. Here, we present a case study of rhodopsin fold determination by integrating sparse and/or low-resolution restraints from multiple experimental techniques including electron microscopy, electron paramagnetic resonance spectroscopy, and nuclear magnetic resonance spectroscopy. Simultaneous incorporation of orthogonal experimental restraints not only significantly improved the sampling accuracy but also allowed identification of the correct fold, which is demonstrated by a protein size-normalized transmembrane root-mean-square deviation as low as 1.2 Å. The protocol developed in this case study can be used for the determination of unknown membrane protein folds when limited experimental restraints are available.
Collapse
Affiliation(s)
- Yan Xia
- Department of Chemistry, Vanderbilt University, Stevenson Center, Station B 351822, Room 7330, Nashville, TN 37232, USA; Center for Structural Biology, Vanderbilt University, Nashville, TN 37232, USA
| | - Axel W Fischer
- Department of Chemistry, Vanderbilt University, Stevenson Center, Station B 351822, Room 7330, Nashville, TN 37232, USA; Center for Structural Biology, Vanderbilt University, Nashville, TN 37232, USA
| | - Pedro Teixeira
- Center for Structural Biology, Vanderbilt University, Nashville, TN 37232, USA
| | - Brian Weiner
- Center for Structural Biology, Vanderbilt University, Nashville, TN 37232, USA
| | - Jens Meiler
- Department of Chemistry, Vanderbilt University, Stevenson Center, Station B 351822, Room 7330, Nashville, TN 37232, USA; Center for Structural Biology, Vanderbilt University, Nashville, TN 37232, USA.
| |
Collapse
|
36
|
Koehler Leman J, Bonneau R, Ulmschneider MB. Statistically derived asymmetric membrane potentials from α-helical and β-barrel membrane proteins. Sci Rep 2018. [PMID: 29535329 PMCID: PMC5849751 DOI: 10.1038/s41598-018-22476-6] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022] Open
Abstract
Modeling membrane protein (MP) folding, insertion, association and their interactions with other proteins, lipids, and drugs requires accurate transfer free energies (TFEs). Various TFE scales have been derived to quantify the energy required or released to insert an amino acid or protein into the membrane. Experimental measurement of TFEs is challenging, and only few scales were extended to depth-dependent energetic profiles. Statistical approaches can be used to derive such potentials; however, this requires a sufficient number of MP structures. Furthermore, MPs are tightly coupled to bilayers that are heterogeneous in terms of lipid composition, asymmetry, and protein content between organisms and organelles. Here we derived asymmetric implicit membrane potentials from β-barrel and α-helical MPs and use them to predict topology, depth and orientation of proteins in the membrane. Our data confirm the ‘charge-outside’ and ‘positive-inside’ rules for β-barrels and α-helical proteins, respectively. We find that the β-barrel profiles have greater asymmetry than the ones from α-helical proteins, as a result of the different membrane architecture of gram-negative bacterial outer membranes and the existence of lipopolysaccharide in the outer leaflet. Our data further suggest that pore-facing residues in β-barrels have a larger contribution to membrane insertion and stability than previously suggested.
Collapse
Affiliation(s)
- Julia Koehler Leman
- Center for Computational Biology, Flatiron Institute, Simons Foundation, New York, 10010 NY, USA. .,Department of Biology, Center for Genomics and Systems Biology, New York University, New York, 10003 NY, USA.
| | - Richard Bonneau
- Center for Computational Biology, Flatiron Institute, Simons Foundation, New York, 10010 NY, USA.,Department of Biology, Center for Genomics and Systems Biology, New York University, New York, 10003 NY, USA.,Department of Computer Science, New York University, New York, 10012 NY, USA
| | | |
Collapse
|
37
|
Popov P, Grudinin S. Eurecon: Equidistant uniform rigid-body ensemble constructor. J Mol Graph Model 2018; 80:313-319. [PMID: 29427936 DOI: 10.1016/j.jmgm.2018.01.015] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2017] [Revised: 01/07/2018] [Accepted: 01/23/2018] [Indexed: 12/14/2022]
Abstract
Conformational ensembles comprise one of the fundamental concepts in statistical bioinformatics and appear in a variety of applications, e.g. molecular docking, virtual screening, searching for pharmacophores, etc. High-throughput applications require billions of conformations to be considered, thus, one often uses the rigid-body representation of molecules or its fragments to cope with the computational cost. Of particular interest is generation of the near-native conformational ensembles, which consist of conformations structurally close to the biologically relevant ones. One possible way to compose such ensembles is to control the root mean square deviation (RMSD) between the original and the generated conformations. To the best of our knowledge there is no computational approach that guarantees that all the generated conformations have the desired RMSD with respect to the reference structure. In this study we presented a fast algorithm for the construction of rigid-body conformational ensembles, which possess two main properties: (i) each generated conformation has a fixed RMSD with respect to the original conformation, (ii) generated conformations are distributed uniformly over the sphere of axes corresponding to the rigid-body motions. The algorithm is very efficient, it does not require any standard RMSD computation between the conformations and has the O(N + M) complexity to generate the required rigid-body transforms, where N is the number of atoms in the system, and M is the size of the conformational ensemble. Eurecon is applicable to an arbitrary atomic system, thus, it could be used for molecular systems of various size and type. We demonstrated Eurecon application by generating near-native conformational ensembles for a ligand placed inside a binding site, a protein dimer embedded into a membrane, and a ribosomal complex. We implemented the developed algorithm in C++ and called it Eurecon, which stands for Equidistant Uniform Rigid-body Ensemble CONstructor. A user-friendly interface allows to define the desired RMSD value, the relative amplitudes for rotation and translation motions by means of the partition parameter, and the set of axes corresponding to the rigid-body motions. Eurecon is available as the SAMSON Element (https://samson-connect.net).
Collapse
Affiliation(s)
- P Popov
- Moscow Institute of Physics and Technology, Dolgoprudny, Russia.
| | | |
Collapse
|
38
|
Abstract
[Formula: see text]-Barrel membrane proteins ([Formula: see text]MPs) play important roles, but knowledge of their structures is limited. We have developed a method to predict their 3D structures. We predict strand registers and construct transmembrane (TM) domains of [Formula: see text]MPs accurately, including proteins for which no prediction has been attempted before. Our method also accurately predicts structures from protein families with a limited number of sequences and proteins with novel folds. An average main-chain rmsd of 3.48 Å is achieved between predicted and experimentally resolved structures of TM domains, which is a significant improvement ([Formula: see text]3 Å) over a recent study. For [Formula: see text]MPs with NMR structures, the deviation between predictions and experimentally solved structures is similar to the difference among the NMR structures, indicating excellent prediction accuracy. Moreover, we can now accurately model the extended [Formula: see text]-barrels and loops in non-TM domains, increasing the overall coverage of structure prediction by [Formula: see text]%. Our method is general and can be applied to genome-wide structural prediction of [Formula: see text]MPs.
Collapse
|
39
|
Koehler Leman J, Mueller BK, Gray JJ. Expanding the toolkit for membrane protein modeling in Rosetta. Bioinformatics 2017; 33:754-756. [PMID: 28011777 DOI: 10.1093/bioinformatics/btw716] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2016] [Accepted: 11/10/2016] [Indexed: 11/14/2022] Open
Abstract
Motivation A range of membrane protein modeling tools has been developed in the past 5-10 years, yet few of these tools are integrated and make use of existing functionality for soluble proteins. To extend existing methods in the Rosetta biomolecular modeling suite for membrane proteins, we recently implemented RosettaMP, a general framework for membrane protein modeling. While RosettaMP facilitates implementation of new methods, addressing real-world biological problems also requires a set of accessory tools that are used to carry out standard modeling tasks. Results Here, we present six modeling tools, including de novo prediction of single trans-membrane helices, making mutations and refining the structure with different amounts of flexibility, transforming a protein into membrane coordinates and optimizing its embedding, computing a Rosetta energy score, and visualizing the protein in the membrane bilayer. We present these methods with complete protocol captures that allow non-expert modelers to carry out the computations. Availability and Implementation The presented tools are part of the Rosetta software suite, available at www.rosettacommons.org . Contact julia.koehler.leman@gmail.com. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Julia Koehler Leman
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, MD 21218, USA.,Simons Center for Data Analysis, Simons Foundation, New York, NY 10001, USA
| | - Benjamin K Mueller
- Department of Chemistry, Vanderbilt University, Nashville, TN 37221, USA.,Center for Structural Biology, Vanderbilt University, Nashville, TN 37221, USA
| | - Jeffrey J Gray
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, MD 21218, USA
| |
Collapse
|
40
|
Hastoy B, Clark A, Rorsman P, Lang J. Fusion pore in exocytosis: More than an exit gate? A β-cell perspective. Cell Calcium 2017; 68:45-61. [PMID: 29129207 DOI: 10.1016/j.ceca.2017.10.005] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2017] [Revised: 10/17/2017] [Accepted: 10/24/2017] [Indexed: 12/14/2022]
Abstract
Secretory vesicle exocytosis is a fundamental biological event and the process by which hormones (like insulin) are released into the blood. Considerable progress has been made in understanding this precisely orchestrated sequence of events from secretory vesicle docked at the cell membrane, hemifusion, to the opening of a membrane fusion pore. The exact biophysical and physiological regulation of these events implies a close interaction between membrane proteins and lipids in a confined space and constrained geometry to ensure appropriate delivery of cargo. We consider some of the still open questions such as the nature of the initiation of the fusion pore, the structure and the role of the Soluble N-ethylmaleimide-sensitive-factor Attachment protein REceptor (SNARE) transmembrane domains and their influence on the dynamics and regulation of exocytosis. We discuss how the membrane composition and protein-lipid interactions influence the likelihood of the nascent fusion pore forming. We relate these factors to the hypothesis that fusion pore expansion could be affected in type-2 diabetes via changes in disease-related gene transcription and alterations in the circulating lipid profile. Detailed characterisation of the dynamics of the fusion pore in vitro will contribute to understanding the larger issue of insulin secretory defects in diabetes.
Collapse
Affiliation(s)
- Benoit Hastoy
- Oxford Centre for Diabetes, Endocrinology and Metabolism, University of Oxford, Churchill Hospital, Oxford OX3 7LE, UK.
| | - Anne Clark
- Oxford Centre for Diabetes, Endocrinology and Metabolism, University of Oxford, Churchill Hospital, Oxford OX3 7LE, UK
| | - Patrik Rorsman
- Oxford Centre for Diabetes, Endocrinology and Metabolism, University of Oxford, Churchill Hospital, Oxford OX3 7LE, UK; Metabolic Research, Institute of Neuroscience and Physiology, University of Goteborg, Medicinaregatan 11, S-41309 Göteborg, Sweden
| | - Jochen Lang
- Laboratoire de Chimie et Biologie des Membranes et Nano-objets (CBMN), CNRS UMR 5248, Université de Bordeaux, Allée de Geoffrey St Hilaire, 33600 Pessac, France.
| |
Collapse
|
41
|
Membrane proteins structures: A review on computational modeling tools. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2017; 1859:2021-2039. [DOI: 10.1016/j.bbamem.2017.07.008] [Citation(s) in RCA: 62] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/23/2017] [Revised: 07/04/2017] [Accepted: 07/13/2017] [Indexed: 01/02/2023]
|
42
|
Latek D. Rosetta Broker for membrane protein structure prediction: concentrative nucleoside transporter 3 and corticotropin-releasing factor receptor 1 test cases. BMC STRUCTURAL BIOLOGY 2017; 17:8. [PMID: 28774292 PMCID: PMC5543540 DOI: 10.1186/s12900-017-0078-8] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/11/2017] [Accepted: 07/26/2017] [Indexed: 02/12/2023]
Abstract
Background Membrane proteins are difficult targets for structure prediction due to the limited structural data deposited in Protein Data Bank. Most computational methods for membrane protein structure prediction are based on the comparative modeling. There are only few de novo methods targeting that distinct protein family. In this work an example of such de novo method was used to structurally and functionally characterize two representatives of distinct membrane proteins families of solute carrier transporters and G protein-coupled receptors. The well-known Rosetta program and one of its protocols named Broker was used in two test cases. The first case was de novo structure prediction of three N-terminal transmembrane helices of the human concentrative nucleoside transporter 3 (hCNT3) homotrimer belonging to the solute carrier 28 family of transporters (SLC28). The second case concerned the large scale refinement of transmembrane helices of a homology model of the corticotropin-releasing factor receptor 1 (CRFR1) belonging to the G protein-coupled receptors family. Results The inward-facing model of the hCNT3 homotrimer was used to propose the functional impact of its single nucleotide polymorphisms. Additionally, the 100 ns molecular dynamics simulation of the unliganded hCNT3 model confirmed its validity and revealed mobility of the selected binding site and homotrimer interface residues. The large scale refinement of transmembrane helices of the CRFR1 homology model resulted in the significant improvement of its accuracy with respect to the crystal structure of CRFR1, especially in the binding site area. Consequently, the antagonist CP-376395 could be docked with Autodock VINA to the CRFR1 model without any steric clashes. Conclusions The presented work demonstrated that Rosetta Broker can be a versatile tool for solving various issues referring to protein biology. Two distinct examples of de novo membrane protein structure prediction presented here provided important insights into three major areas of protein biology. Namely, the dynamics of the inward-facing hCNT3 homotrimer system, the structural changes of the CRFR1 receptor upon the antagonist binding and finally, the role of single nucleotide polymorphisms in both, hCNT3 and CRFR1 proteins, were investigated. Electronic supplementary material The online version of this article (doi:10.1186/s12900-017-0078-8) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Dorota Latek
- Faculty of Chemistry, University of Warsaw, Pasteur St. 1, 02-093, Warsaw, Poland.
| |
Collapse
|
43
|
Wang J, Li W, Wang B, Hu B, Jiang H, Lai B, Li N, Cheng M. In Silicon Approach for Discovery of Chemopreventive Agents. ACTA ACUST UNITED AC 2017. [DOI: 10.1007/s40495-017-0094-1] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
|
44
|
Hastoy B, Scotti PA, Milochau A, Fezoua-Boubegtiten Z, Rodas J, Megret R, Desbat B, Laguerre M, Castano S, Perrais D, Rorsman P, Oda R, Lang J. A Central Small Amino Acid in the VAMP2 Transmembrane Domain Regulates the Fusion Pore in Exocytosis. Sci Rep 2017; 7:2835. [PMID: 28588281 PMCID: PMC5460238 DOI: 10.1038/s41598-017-03013-3] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2016] [Accepted: 04/27/2017] [Indexed: 11/30/2022] Open
Abstract
Exocytosis depends on cytosolic domains of SNARE proteins but the function of the transmembrane domains (TMDs) in membrane fusion remains controversial. The TMD of the SNARE protein synaptobrevin2/VAMP2 contains two highly conserved small amino acids, G100 and C103, in its central portion. Substituting G100 and/or C103 with the β-branched amino acid valine impairs the structural flexibility of the TMD in terms of α-helix/β-sheet transitions in model membranes (measured by infrared reflection-absorption or evanescent wave spectroscopy) during increase in protein/lipid ratios, a parameter expected to be altered by recruitment of SNAREs at fusion sites. This structural change is accompanied by reduced membrane fluidity (measured by infrared ellipsometry). The G100V/C103V mutation nearly abolishes depolarization-evoked exocytosis (measured by membrane capacitance) and hormone secretion (measured biochemically). Single-vesicle optical (by TIRF microscopy) and biophysical measurements of ATP release indicate that G100V/C103V retards initial fusion-pore opening, hinders its expansion and leads to premature closure in most instances. We conclude that the TMD of VAMP2 plays a critical role in membrane fusion and that the structural mobility provided by the central small amino acids is crucial for exocytosis by influencing the molecular re-arrangements of the lipid membrane that are necessary for fusion pore opening and expansion.
Collapse
Affiliation(s)
- Benoît Hastoy
- Laboratory of Membrane Chemistry and Biology (CBMN), UMR CNRS 5248, Université de Bordeaux, Allée de Geoffroy St Hilaire, 33600, Pessac, France.,Université de Bordeaux, 351 Cours de la Libération, 33400, Talence, France.,Oxford Centre for Diabetes, Endocrinology and Metabolism, University of Oxford, Churchill Hospital, Oxford, OX3 7LJ, UK
| | - Pier A Scotti
- Laboratory of Membrane Chemistry and Biology (CBMN), UMR CNRS 5248, Université de Bordeaux, Allée de Geoffroy St Hilaire, 33600, Pessac, France.,Université de Bordeaux, 351 Cours de la Libération, 33400, Talence, France
| | - Alexandra Milochau
- Laboratory of Membrane Chemistry and Biology (CBMN), UMR CNRS 5248, Université de Bordeaux, Allée de Geoffroy St Hilaire, 33600, Pessac, France.,Université de Bordeaux, 351 Cours de la Libération, 33400, Talence, France
| | - Zahia Fezoua-Boubegtiten
- Laboratory of Membrane Chemistry and Biology (CBMN), UMR CNRS 5248, Université de Bordeaux, Allée de Geoffroy St Hilaire, 33600, Pessac, France.,Université de Bordeaux, 351 Cours de la Libération, 33400, Talence, France
| | - Jorge Rodas
- Université de Bordeaux, 351 Cours de la Libération, 33400, Talence, France.,Laboratoire de l'Intégration du Matériau au Système, UMR CNRS 5218, 351 Cours de la Libération, 33400 Talence, France.,Institut Polytechnique de Bordeaux, Avernue des Facultés, 33405, Talence, France
| | - Rémi Megret
- Université de Bordeaux, 351 Cours de la Libération, 33400, Talence, France.,Laboratoire de l'Intégration du Matériau au Système, UMR CNRS 5218, 351 Cours de la Libération, 33400 Talence, France.,Institut Polytechnique de Bordeaux, Avernue des Facultés, 33405, Talence, France
| | - Bernard Desbat
- Laboratory of Membrane Chemistry and Biology (CBMN), UMR CNRS 5248, Université de Bordeaux, Allée de Geoffroy St Hilaire, 33600, Pessac, France.,Université de Bordeaux, 351 Cours de la Libération, 33400, Talence, France
| | - Michel Laguerre
- Laboratory of Membrane Chemistry and Biology (CBMN), UMR CNRS 5248, Université de Bordeaux, Allée de Geoffroy St Hilaire, 33600, Pessac, France.,Université de Bordeaux, 351 Cours de la Libération, 33400, Talence, France
| | - Sabine Castano
- Laboratory of Membrane Chemistry and Biology (CBMN), UMR CNRS 5248, Université de Bordeaux, Allée de Geoffroy St Hilaire, 33600, Pessac, France.,Université de Bordeaux, 351 Cours de la Libération, 33400, Talence, France
| | - David Perrais
- Université de Bordeaux, 351 Cours de la Libération, 33400, Talence, France.,Interdisciplinary Institute for Neuroscience, UMR CNRS 5287, 146, rue Léo-Saignat, 33077, Bordeaux, France
| | - Patrik Rorsman
- Oxford Centre for Diabetes, Endocrinology and Metabolism, University of Oxford, Churchill Hospital, Oxford, OX3 7LJ, UK
| | - Reiko Oda
- Laboratory of Membrane Chemistry and Biology (CBMN), UMR CNRS 5248, Université de Bordeaux, Allée de Geoffroy St Hilaire, 33600, Pessac, France.,Université de Bordeaux, 351 Cours de la Libération, 33400, Talence, France
| | - Jochen Lang
- Laboratory of Membrane Chemistry and Biology (CBMN), UMR CNRS 5248, Université de Bordeaux, Allée de Geoffroy St Hilaire, 33600, Pessac, France. .,Université de Bordeaux, 351 Cours de la Libération, 33400, Talence, France.
| |
Collapse
|
45
|
Simkovic F, Ovchinnikov S, Baker D, Rigden DJ. Applications of contact predictions to structural biology. IUCRJ 2017; 4:291-300. [PMID: 28512576 PMCID: PMC5414403 DOI: 10.1107/s2052252517005115] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/12/2016] [Accepted: 04/03/2017] [Indexed: 06/07/2023]
Abstract
Evolutionary pressure on residue interactions, intramolecular or intermolecular, that are important for protein structure or function can lead to covariance between the two positions. Recent methodological advances allow much more accurate contact predictions to be derived from this evolutionary covariance signal. The practical application of contact predictions has largely been confined to structural bioinformatics, yet, as this work seeks to demonstrate, the data can be of enormous value to the structural biologist working in X-ray crystallo-graphy, cryo-EM or NMR. Integrative structural bioinformatics packages such as Rosetta can already exploit contact predictions in a variety of ways. The contribution of contact predictions begins at construct design, where structural domains may need to be expressed separately and contact predictions can help to predict domain limits. Structure solution by molecular replacement (MR) benefits from contact predictions in diverse ways: in difficult cases, more accurate search models can be constructed using ab initio modelling when predictions are available, while intermolecular contact predictions can allow the construction of larger, oligomeric search models. Furthermore, MR using supersecondary motifs or large-scale screens against the PDB can exploit information, such as the parallel or antiparallel nature of any β-strand pairing in the target, that can be inferred from contact predictions. Contact information will be particularly valuable in the determination of lower resolution structures by helping to assign sequence register. In large complexes, contact information may allow the identity of a protein responsible for a certain region of density to be determined and then assist in the orientation of an available model within that density. In NMR, predicted contacts can provide long-range information to extend the upper size limit of the technique in a manner analogous but complementary to experimental methods. Finally, predicted contacts can distinguish between biologically relevant interfaces and mere lattice contacts in a final crystal structure, and have potential in the identification of functionally important regions and in foreseeing the consequences of mutations.
Collapse
Affiliation(s)
- Felix Simkovic
- Institute of Integrative Biology, University of Liverpool, Liverpool L69 7ZB, England
| | - Sergey Ovchinnikov
- Department of Biochemistry, University of Washington, Seattle, WA 98195, USA
- Institute for Protein Design, University of Washington, Seattle, WA 98195, USA
- Howard Hughes Medical Institute, University of Washington, Box 357370, Seattle, WA 98195, USA
| | - David Baker
- Department of Biochemistry, University of Washington, Seattle, WA 98195, USA
- Institute for Protein Design, University of Washington, Seattle, WA 98195, USA
- Howard Hughes Medical Institute, University of Washington, Box 357370, Seattle, WA 98195, USA
| | - Daniel J. Rigden
- Institute of Integrative Biology, University of Liverpool, Liverpool L69 7ZB, England
| |
Collapse
|
46
|
Feig M. Computational protein structure refinement: Almost there, yet still so far to go. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2017; 7:e1307. [PMID: 30613211 PMCID: PMC6319934 DOI: 10.1002/wcms.1307] [Citation(s) in RCA: 46] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Protein structures are essential in modern biology yet experimental methods are far from being able to catch up with the rapid increase in available genomic data. Computational protein structure prediction methods aim to fill the gap while the role of protein structure refinement is to take approximate initial template-based models and bring them closer to the true native structure. Current methods for computational structure refinement rely on molecular dynamics simulations, related sampling methods, or iterative structure optimization protocols. The best methods are able to achieve moderate degrees of refinement but consistent refinement that can reach near-experimental accuracy remains elusive. Key issues revolve around the accuracy of the energy function, the inability to reliably rank multiple models, and the use of restraints that keep sampling close to the native state but also limit the degree of possible refinement. A different aspect is the question of what exactly the target of high-resolution refinement should be as experimental structures are affected by experimental conditions and different biological questions require varying levels of accuracy. While improvement of the global protein structure is a difficult problem, high-resolution refinement methods that improves local structural quality such as favorable stereochemistry and the avoidance of atomic clashes are much more successful.
Collapse
Affiliation(s)
- Michael Feig
- Department of Biochemistry and Molecular Biology, Michigan State University, 603 Wilson Rd., Room 218 BCH, East Lansing, MI, USA, ; 517-432-7439
| |
Collapse
|
47
|
Grassmann AA, Kremer FS, Dos Santos JC, Souza JD, Pinto LDS, McBride AJA. Discovery of Novel Leptospirosis Vaccine Candidates Using Reverse and Structural Vaccinology. Front Immunol 2017; 8:463. [PMID: 28496441 PMCID: PMC5406399 DOI: 10.3389/fimmu.2017.00463] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2016] [Accepted: 04/04/2017] [Indexed: 12/03/2022] Open
Abstract
Leptospira spp. are diderm (two membranes) bacteria that infect mammals causing leptospirosis, a public health problem with global implications. Thousands of people die every year due to leptospirosis, especially in developing countries with tropical climates. Prophylaxis is difficult due to multiple factors, including the large number of asymptomatic hosts that transmit the bacteria, poor sanitation, increasing numbers of slum dwellers, and the lack of an effective vaccine. Several leptospiral recombinant antigens were evaluated as a replacement for the inactivated (bacterin) vaccine; however, success has been limited. A prospective vaccine candidate is likely to be a surface-related protein that can stimulate the host immune response to clear leptospires from blood and organs. In this study, a comprehensive bioinformatics approach based on reverse and structural vaccinology was applied toward the discovery of novel leptospiral vaccine candidates. The Leptospira interrogans serovar Copenhageni strain L1-130 genome was mined in silico for the enhanced identification of conserved β-barrel (βb) transmembrane proteins and outer membrane (OM) lipoproteins. Orthologs of the prospective vaccine candidates were screened in the genomes of 20 additional Leptospira spp. Three-dimensional structural models, with a high degree of confidence, were created for each of the surface-exposed proteins. Major histocompatibility complex II (MHC-II) epitopes were identified, and their locations were mapped on the structural models. A total of 18 βb transmembrane proteins and 8 OM lipoproteins were identified. These proteins were conserved among the pathogenic Leptospira spp. and were predicted to have epitopes for several variants of MHC-II receptors. A structural and functional analysis of the sequence of these surface proteins demonstrated that most βb transmembrane proteins seem to be TonB-dependent receptors associated with transportation. Other proteins identified included, e.g., TolC efflux pump proteins, a BamA-like OM component of the βb transmembrane protein assembly machinery, and the LptD-like LPS assembly protein. The structural mapping of the immunodominant epitopes identified the location of conserved, surface-exposed, immunogenic regions for each vaccine candidate. The proteins identified in this study are currently being evaluated for experimental evidence for their involvement in virulence, disease pathogenesis, and physiology, in addition to vaccine development.
Collapse
Affiliation(s)
- André Alex Grassmann
- Biotechnology Unit, Technological Development Center, Federal University of Pelotas, Pelotas, Rio Grande do Sul, Brazil
| | - Frederico Schmitt Kremer
- Biotechnology Unit, Technological Development Center, Federal University of Pelotas, Pelotas, Rio Grande do Sul, Brazil
| | - Júlia Cougo Dos Santos
- Biotechnology Unit, Technological Development Center, Federal University of Pelotas, Pelotas, Rio Grande do Sul, Brazil
| | - Jéssica Dias Souza
- Biotechnology Unit, Technological Development Center, Federal University of Pelotas, Pelotas, Rio Grande do Sul, Brazil
| | - Luciano da Silva Pinto
- Biotechnology Unit, Technological Development Center, Federal University of Pelotas, Pelotas, Rio Grande do Sul, Brazil
| | - Alan John Alexander McBride
- Biotechnology Unit, Technological Development Center, Federal University of Pelotas, Pelotas, Rio Grande do Sul, Brazil.,Gonçalo Moniz Institute, Oswaldo Cruz Foundation, Ministry of Health, Salvador, Bahia, Brazil
| |
Collapse
|
48
|
Abstract
Membrane proteins play a most important part in metabolism, signaling, cell motility, transport, development, and many other biochemical and biophysical processes which constitute fundamentals of life on the molecular level. Detailed understanding of these processes is necessary for the progress of life sciences and biomedical applications. Nanodiscs provide a new and powerful tool for a broad spectrum of biochemical and biophysical studies of membrane proteins and are commonly acknowledged as an optimal membrane mimetic system that provides control over size, composition, and specific functional modifications on the nanometer scale. In this review we attempted to combine a comprehensive list of various applications of nanodisc technology with systematic analysis of the most attractive features of this system and advantages provided by nanodiscs for structural and mechanistic studies of membrane proteins.
Collapse
Affiliation(s)
- Ilia G Denisov
- Department of Biochemistry and Department of Chemistry, University of Illinois , Urbana, Illinois 61801, United States
| | - Stephen G Sligar
- Department of Biochemistry and Department of Chemistry, University of Illinois , Urbana, Illinois 61801, United States
| |
Collapse
|
49
|
Reichel K, Fisette O, Braun T, Lange OF, Hummer G, Schäfer LV. Systematic evaluation of CS-Rosetta for membrane protein structure prediction with sparse NOE restraints. Proteins 2017; 85:812-826. [PMID: 27936510 DOI: 10.1002/prot.25224] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2016] [Revised: 10/25/2016] [Accepted: 11/23/2016] [Indexed: 11/06/2022]
Abstract
We critically test and validate the CS-Rosetta methodology for de novo structure prediction of α-helical membrane proteins (MPs) from NMR data, such as chemical shifts and NOE distance restraints. By systematically reducing the number and types of NOE restraints, we focus on determining the regime in which MP structures can be reliably predicted and pinpoint the boundaries of the approach. Five MPs of known structure were used as test systems, phototaxis sensory rhodopsin II (pSRII), a subdomain of pSRII, disulfide binding protein B (DsbB), microsomal prostaglandin E2 synthase-1 (mPGES-1), and translocator protein (TSPO). For pSRII and DsbB, where NMR and X-ray structures are available, resolution-adapted structural recombination (RASREC) CS-Rosetta yields structures that are as close to the X-ray structure as the published NMR structures if all available NMR data are used to guide structure prediction. For mPGES-1 and Bacillus cereus TSPO, where only X-ray crystal structures are available, highly accurate structures are obtained using simulated NMR data. One main advantage of RASREC CS-Rosetta is its robustness with respect to even a drastic reduction of the number of NOEs. Close-to-native structures were obtained with one randomly picked long-range NOEs for every 14, 31, 38, and 8 residues for full-length pSRII, the pSRII subdomain, TSPO, and DsbB, respectively, in addition to using chemical shifts. For mPGES-1, atomically accurate structures could be predicted even from chemical shifts alone. Our results show that atomic level accuracy for helical membrane proteins is achievable with CS-Rosetta using very sparse NOE restraint sets to guide structure prediction. Proteins 2017; 85:812-826. © 2016 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Katrin Reichel
- Center for Theoretical Chemistry, Ruhr-University Bochum, Bochum, 44780, Germany.,Max Planck Institute of Biophysics, 60438, Frankfurt am Main, Germany
| | - Olivier Fisette
- Center for Theoretical Chemistry, Ruhr-University Bochum, Bochum, 44780, Germany
| | - Tatjana Braun
- ICS-6 Structural Biochemistry, Institute of Complex Systems, Forschungszentrum Jülich, Jülich, 52425, Germany
| | - Oliver F Lange
- Biomolecular NMR and Munich Center for Integrated Protein Science, Department Chemie, Technische Universität München, Garching, 85747, Germany
| | - Gerhard Hummer
- Max Planck Institute of Biophysics, 60438, Frankfurt am Main, Germany.,Institute of Biophysics, Goethe University, 60438, Frankfurt am Main, Germany
| | - Lars V Schäfer
- Center for Theoretical Chemistry, Ruhr-University Bochum, Bochum, 44780, Germany
| |
Collapse
|
50
|
Koehler Leman J, Lyskov S, Bonneau R. Computing structure-based lipid accessibility of membrane proteins with mp_lipid_acc in RosettaMP. BMC Bioinformatics 2017; 18:115. [PMID: 28219343 PMCID: PMC5319049 DOI: 10.1186/s12859-017-1541-z] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2016] [Accepted: 02/08/2017] [Indexed: 11/17/2022] Open
Abstract
BACKGROUND Membrane proteins are underrepresented in structural databases, which has led to a lack of computational tools and the corresponding inappropriate use of tools designed for soluble proteins. For membrane proteins, lipid accessibility is an essential property. Although programs are available for sequence-based prediction of lipid accessibility and structure-based identification of solvent-accessible surface area, the latter does not distinguish between water accessible and lipid accessible residues in membrane proteins. RESULTS Here we present mp_lipid_acc, the first method to identify lipid accessible residues from the protein structure, implemented in the RosettaMP framework and available as a webserver. Our method uses protein structures transformed in membrane coordinates, for instance from PDBTM or OPM databases, and a defined membrane thickness to classify lipid accessibility of residues. mp_lipid_acc is applicable to both α-helical and β-barrel membrane proteins of diverse architectures with or without water-filled pores and uses a concave hull algorithm for surface-residue classification. We further provide a manually curated benchmark dataset that can be used for further method development. CONCLUSIONS We present a novel tool to classify lipid accessibility from the protein structure, which is applicable to proteins of diverse architectures and achieves prediction accuracies of 90% on a manually curated database. mp_lipid_acc is part of the Rosetta software suite, available at www.rosettacommons.org . The webserver is available at http://rosie.graylab.jhu.edu/mp_lipid_acc/submit and the benchmark dataset is available at http://tinyurl.com/mp-lipid-acc-dataset .
Collapse
Affiliation(s)
- Julia Koehler Leman
- Center for Computational Biology, Flatiron Institute, 162 Fifth Avenue, New York, NY 10010 USA
- Departments of Biology and Computer Science, Center for Genomics and Systems Biology, New York University, New York, NY 10003 USA
| | - Sergey Lyskov
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, MD 21218 USA
| | - Richard Bonneau
- Center for Computational Biology, Flatiron Institute, 162 Fifth Avenue, New York, NY 10010 USA
- Departments of Biology and Computer Science, Center for Genomics and Systems Biology, New York University, New York, NY 10003 USA
| |
Collapse
|