201
|
Wang T, Yang Y, Zhou Y, Gong H. LRFragLib: an effective algorithm to identify fragments for de novo protein structure prediction. Bioinformatics 2017; 33:677-684. [PMID: 27797773 DOI: 10.1093/bioinformatics/btw668] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2016] [Accepted: 10/18/2016] [Indexed: 11/13/2022] Open
Abstract
Motivation The quality of fragment library determines the efficiency of fragment assembly, an approach that is widely used in most de novo protein-structure prediction algorithms. Conventional fragment libraries are constructed mainly based on the identities of amino acids, sometimes facilitated by predicted information including dihedral angles and secondary structures. However, it remains challenging to identify near-native fragment structures with low sequence homology. Results We introduce a novel fragment-library-construction algorithm, LRFragLib, to improve the detection of near-native low-homology fragments of 7-10 residues, using a multi-stage, flexible selection protocol. Based on logistic regression scoring models, LRFragLib outperforms existing techniques by achieving a significantly higher precision and a comparable coverage on recent CASP protein sets in sampling near-native structures. The method also has a comparable computational efficiency to the fastest existing techniques with substantially reduced memory usage. Availability and Implementation The source code is available for download at http://166.111.152.91/Downloads.html. Contact hgong@tsinghua.edu.cn. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Tong Wang
- MOE Key Laboratory of Bioinformatics, School of Life Sciences.,Beijing Innovation Center of Structural Biology, Tsinghua University, Beijing 100084, China
| | - Yuedong Yang
- Institute for Glycomics and School of Information and Communication Technology, Griffith University, Gold Coast, QLD 4222, Australia
| | - Yaoqi Zhou
- Institute for Glycomics and School of Information and Communication Technology, Griffith University, Gold Coast, QLD 4222, Australia
| | - Haipeng Gong
- MOE Key Laboratory of Bioinformatics, School of Life Sciences.,Beijing Innovation Center of Structural Biology, Tsinghua University, Beijing 100084, China
| |
Collapse
|
202
|
Thomas JMH, Simkovic F, Keegan R, Mayans O, Zhang C, Zhang Y, Rigden DJ. Approaches to ab initio molecular replacement of α-helical transmembrane proteins. Acta Crystallogr D Struct Biol 2017; 73:985-996. [PMID: 29199978 PMCID: PMC5713875 DOI: 10.1107/s2059798317016436] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2017] [Accepted: 11/15/2017] [Indexed: 02/06/2023] Open
Abstract
α-Helical transmembrane proteins are a ubiquitous and important class of proteins, but present difficulties for crystallographic structure solution. Here, the effectiveness of the AMPLE molecular replacement pipeline in solving α-helical transmembrane-protein structures is assessed using a small library of eight ideal helices, as well as search models derived from ab initio models generated both with and without evolutionary contact information. The ideal helices prove to be surprisingly effective at solving higher resolution structures, but ab initio-derived search models are able to solve structures that could not be solved with the ideal helices. The addition of evolutionary contact information results in a marked improvement in the modelling and makes additional solutions possible.
Collapse
Affiliation(s)
- Jens M. H. Thomas
- Institute of Integrative Biology, University of Liverpool, Liverpool L69 7ZB, England
| | - Felix Simkovic
- Institute of Integrative Biology, University of Liverpool, Liverpool L69 7ZB, England
| | - Ronan Keegan
- Research Complex at Harwell, STFC Rutherford Appleton Laboratory, Didcot OX11 0FA, England
| | - Olga Mayans
- Fachbereich Biologie, Universität Konstanz, D-78457 Konstanz, Germany
| | - Chengxin Zhang
- Department of Computational Medicine and Bioinformatics, Department of Biological Chemistry, Medical School, University of Michigan, 100 Washtenaw Avenue, Ann Arbor, MI 48109-2218, USA
| | - Yang Zhang
- Department of Computational Medicine and Bioinformatics, Department of Biological Chemistry, Medical School, University of Michigan, 100 Washtenaw Avenue, Ann Arbor, MI 48109-2218, USA
| | - Daniel J. Rigden
- Institute of Integrative Biology, University of Liverpool, Liverpool L69 7ZB, England
| |
Collapse
|
203
|
Vetrivel I, Mahajan S, Tyagi M, Hoffmann L, Sanejouand YH, Srinivasan N, de Brevern AG, Cadet F, Offmann B. Knowledge-based prediction of protein backbone conformation using a structural alphabet. PLoS One 2017; 12:e0186215. [PMID: 29161266 PMCID: PMC5697859 DOI: 10.1371/journal.pone.0186215] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2017] [Accepted: 09/27/2017] [Indexed: 01/19/2023] Open
Abstract
Libraries of structural prototypes that abstract protein local structures are known as structural alphabets and have proven to be very useful in various aspects of protein structure analyses and predictions. One such library, Protein Blocks, is composed of 16 standard 5-residues long structural prototypes. This form of analyzing proteins involves drafting its structure as a string of Protein Blocks. Predicting the local structure of a protein in terms of protein blocks is the general objective of this work. A new approach, PB-kPRED is proposed towards this aim. It involves (i) organizing the structural knowledge in the form of a database of pentapeptide fragments extracted from all protein structures in the PDB and (ii) applying a knowledge-based algorithm that does not rely on any secondary structure predictions and/or sequence alignment profiles, to scan this database and predict most probable backbone conformations for the protein local structures. Though PB-kPRED uses the structural information from homologues in preference, if available. The predictions were evaluated rigorously on 15,544 query proteins representing a non-redundant subset of the PDB filtered at 30% sequence identity cut-off. We have shown that the kPRED method was able to achieve mean accuracies ranging from 40.8% to 66.3% depending on the availability of homologues. The impact of the different strategies for scanning the database on the prediction was evaluated and is discussed. Our results highlight the usefulness of the method in the context of proteins without any known structural homologues. A scoring function that gives a good estimate of the accuracy of prediction was further developed. This score estimates very well the accuracy of the algorithm (R2 of 0.82). An online version of the tool is provided freely for non-commercial usage at http://www.bo-protscience.fr/kpred/.
Collapse
Affiliation(s)
- Iyanar Vetrivel
- Université de Nantes, Unité Fonctionnalité et Ingénierie des Protéines (UFIP), UMR 6286 CNRS, UFR Sciences et Techniques, 2, chemin de la Houssinière, France
| | - Swapnil Mahajan
- Université de Nantes, Unité Fonctionnalité et Ingénierie des Protéines (UFIP), UMR 6286 CNRS, UFR Sciences et Techniques, 2, chemin de la Houssinière, France
- DSIMB, INSERM, UMR S-1134, Laboratory of Excellence, GR-Ex, Université de La Réunion, Faculty of Sciences and Technology, Saint Denis Cedex, La Réunion, France
| | - Manoj Tyagi
- Université de La Réunion, Saint Denis Cedex, La Réunion, France
| | - Lionel Hoffmann
- Université de Nantes, Unité Fonctionnalité et Ingénierie des Protéines (UFIP), UMR 6286 CNRS, UFR Sciences et Techniques, 2, chemin de la Houssinière, France
| | - Yves-Henri Sanejouand
- Université de Nantes, Unité Fonctionnalité et Ingénierie des Protéines (UFIP), UMR 6286 CNRS, UFR Sciences et Techniques, 2, chemin de la Houssinière, France
| | | | - Alexandre G. de Brevern
- INSERM UMR_S 1134, DSIMB team, Laboratory of Excellence, GR-Ex, Univ Paris Diderot, Univ Sorbonne Paris Cité, INTS, rue Alexandre Cabanel, Paris, France
| | - Frédéric Cadet
- DSIMB, INSERM, UMR S-1134, Laboratory of Excellence, GR-Ex, Université de La Réunion, Faculty of Sciences and Technology, Saint Denis Cedex, La Réunion, France
- PEACCEL SAS, Paris, France
| | - Bernard Offmann
- Université de Nantes, Unité Fonctionnalité et Ingénierie des Protéines (UFIP), UMR 6286 CNRS, UFR Sciences et Techniques, 2, chemin de la Houssinière, France
| |
Collapse
|
204
|
Xu G, Ma T, Zang T, Wang Q, Ma J. OPUS-CSF: A C-atom-based scoring function for ranking protein structural models. Protein Sci 2017; 27:286-292. [PMID: 29047165 PMCID: PMC5734313 DOI: 10.1002/pro.3327] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2017] [Revised: 10/14/2017] [Accepted: 10/16/2017] [Indexed: 12/12/2022]
Abstract
We report a C‐atom‐based scoring function, named OPUS‐CSF, for ranking protein structural models. Rather than using traditional Boltzmann formula, we built a scoring function (CSF score) based on the native distributions (derived from the entire PDB) of coordinate components of mainchain C (carbonyl) atoms on selected residues of peptide segments of 5, 7, 9, and 11 residues in length. In testing OPUS‐CSF on decoy recognition, it maximally recognized 257 native structures out of 278 targets in 11 commonly used decoy sets, significantly outperforming other popular all‐atom empirical potentials. The average correlation coefficient with TM‐score was also comparable with those of other potentials. OPUS‐CSF is a highly coarse‐grained scoring function, which only requires input of partial mainchain information, and very fast. Thus, it is suitable for applications at early stage of structural building.
Collapse
Affiliation(s)
- Gang Xu
- School of Life Sciences, Tsinghua University, Beijing, China
| | - Tianqi Ma
- Applied Physics Program, Rice University, Houston, Texas.,Department of Bioengineering, Rice University, Houston, Texas
| | - Tianwu Zang
- Applied Physics Program, Rice University, Houston, Texas.,Department of Bioengineering, Rice University, Houston, Texas
| | - Qinghua Wang
- Verna and Marrs McLean Department of Biochemistry and Molecular Biology, Baylor College of Medicine, One Baylor Plaza, BCM-125, Houston, Texas
| | - Jianpeng Ma
- School of Life Sciences, Tsinghua University, Beijing, China.,Applied Physics Program, Rice University, Houston, Texas.,Department of Bioengineering, Rice University, Houston, Texas.,Verna and Marrs McLean Department of Biochemistry and Molecular Biology, Baylor College of Medicine, One Baylor Plaza, BCM-125, Houston, Texas
| |
Collapse
|
205
|
The soluble domains of Gpi8 and Gaa1, two subunits of glycosylphosphatidylinositol transamidase (GPI-T), assemble into a complex. Arch Biochem Biophys 2017; 633:58-67. [DOI: 10.1016/j.abb.2017.09.006] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2017] [Revised: 09/06/2017] [Accepted: 09/07/2017] [Indexed: 11/23/2022]
|
206
|
Optimal Tilt Angle and Orientation of Photovoltaic Modules Using HS Algorithm in Different Climates of China. APPLIED SCIENCES-BASEL 2017. [DOI: 10.3390/app7101028] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
207
|
Hao XH, Zhang GJ, Zhou XG. Conformational Space Sampling Method Using Multi-Subpopulation Differential Evolution for De novo Protein Structure Prediction. IEEE Trans Nanobioscience 2017; 16:618-633. [DOI: 10.1109/tnb.2017.2749243] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
208
|
Sønderby P, Rinnan Å, Madsen JJ, Harris P, Bukrinski JT, Peters GHJ. Small-Angle X-ray Scattering Data in Combination with RosettaDock Improves the Docking Energy Landscape. J Chem Inf Model 2017; 57:2463-2475. [DOI: 10.1021/acs.jcim.6b00789] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Affiliation(s)
- Pernille Sønderby
- Department
of Chemistry, Technical University of Denmark, DK-2800 Kongens
Lyngby, Denmark
| | - Åsmund Rinnan
- Department
of Food Science, Faculty of Science, University of Copenhagen, DK-1958 Frederiksberg C, Denmark
| | - Jesper J. Madsen
- Department
of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Pernille Harris
- Department
of Chemistry, Technical University of Denmark, DK-2800 Kongens
Lyngby, Denmark
| | | | - Günther H. J. Peters
- Department
of Chemistry, Technical University of Denmark, DK-2800 Kongens
Lyngby, Denmark
| |
Collapse
|
209
|
Suply T, Hannedouche S, Carte N, Li J, Grosshans B, Schaefer M, Raad L, Beck V, Vidal S, Hiou-Feige A, Beluch N, Barbieri S, Wirsching J, Lageyre N, Hillger F, Debon C, Dawson J, Smith P, Lannoy V, Detheux M, Bitsch F, Falchetto R, Bouwmeester T, Porter J, Baumgarten B, Mansfield K, Carballido JM, Seuwen K, Bassilana F. A natural ligand for the orphan receptor GPR15 modulates lymphocyte recruitment to epithelia. Sci Signal 2017; 10:10/496/eaal0180. [DOI: 10.1126/scisignal.aal0180] [Citation(s) in RCA: 46] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
|
210
|
Xu G, Ma T, Zang T, Sun W, Wang Q, Ma J. OPUS-DOSP: A Distance- and Orientation-Dependent All-Atom Potential Derived from Side-Chain Packing. J Mol Biol 2017; 429:3113-3120. [PMID: 28864201 DOI: 10.1016/j.jmb.2017.08.013] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2017] [Revised: 07/27/2017] [Accepted: 08/22/2017] [Indexed: 01/18/2023]
Abstract
We report a new distance- and orientation-dependent, all-atom statistical potential derived from side-chain packing, named OPUS-DOSP, for protein structure modeling. The framework of OPUS-DOSP is based on OPUS-PSP, previously developed by us [JMB (2008), 376, 288-301], with refinement and new features. In particular, distance or orientation contribution is considered depending on the range of contact distance. A new auxiliary function in energy function is also introduced, in addition to the traditional Boltzmann term, in order to adjust the contributions of extreme cases. OPUS-DOSP was tested on 11 decoy sets commonly used for statistical potential benchmarking. Among 278 native structures, 239 and 249 native structures were recognized by OPUS-DOSP without and with the auxiliary function, respectively. The results show that OPUS-DOSP has an increased decoy recognition capability comparing with those of other relevant potentials to date.
Collapse
Affiliation(s)
- Gang Xu
- School of Life Sciences, Tsinghua University, Beijing 100084, China
| | - Tianqi Ma
- Applied Physics Program, Rice University, Houston, TX 77005, United States; Department of Bioengineering, Rice University, Houston, TX 77005, United States
| | - Tianwu Zang
- Applied Physics Program, Rice University, Houston, TX 77005, United States; Department of Bioengineering, Rice University, Houston, TX 77005, United States
| | - Weitao Sun
- Zhou Pei-Yuan Center for Applied Mathematics, Tsinghua University, Beijing 100084, China
| | - Qinghua Wang
- Verna and Marrs Mclean Department of Biochemistry and Molecular Biology, Baylor College of Medicine, Houston, TX 77030, United States
| | - Jianpeng Ma
- School of Life Sciences, Tsinghua University, Beijing 100084, China; Applied Physics Program, Rice University, Houston, TX 77005, United States; Department of Bioengineering, Rice University, Houston, TX 77005, United States; Verna and Marrs Mclean Department of Biochemistry and Molecular Biology, Baylor College of Medicine, Houston, TX 77030, United States.
| |
Collapse
|
211
|
Alexander NS, Palczewski K. Crowd sourcing difficult problems in protein science . Protein Sci 2017; 26:2118-2125. [PMID: 28762619 DOI: 10.1002/pro.3247] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2017] [Accepted: 07/21/2017] [Indexed: 11/08/2022]
Abstract
Dedicated computing resources are expensive to develop, maintain, and administrate. Frequently, research groups require bursts of computing power, during which progress is still limited by available computing resources. One way to alleviate this bottleneck would be to use additional computing resources. Today, many computing devices remain idle most of the time. Passive volunteer computing exploits this unemployed reserve of computing power by allowing device-owners to donate computing time on their own devices. Another complementary way to alleviate bottlenecks in computing resources is to use more efficient algorithms. Engaging volunteer computing employs human intuition to help solve challenging problems for which efficient algorithms are difficult to develop or unavailable. Designing engaging volunteer computing projects is challenging but can result in high-quality solutions. Here, we highlight four examples.
Collapse
Affiliation(s)
- Nathan S Alexander
- Department of Pharmacology, School of Medicine, Case Western Reserve University, Cleveland, Ohio, 44106
| | - Krzysztof Palczewski
- Department of Pharmacology, School of Medicine, Case Western Reserve University, Cleveland, Ohio, 44106.,Cleveland Center for Membrane and Structural Biology, Case Western Reserve University, Cleveland, Ohio, 44106
| |
Collapse
|
212
|
Convergent and divergent genetic changes in the genome of Chinese and European pigs. Sci Rep 2017; 7:8662. [PMID: 28819228 PMCID: PMC5561219 DOI: 10.1038/s41598-017-09061-z] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2017] [Accepted: 07/20/2017] [Indexed: 01/17/2023] Open
Abstract
Since 10,000 BC, continuous human selection has led to intense genetic and phenotypic changes in pig (Sus scrofa) domestication. Through whole genome analysis of 257 individuals, we demonstrated artificial unidirectional and bidirectional selection as the primary force to shape the convergent and divergent changes between Chinese domestic pigs (CHD) and European domestic pigs (EUD). We identified 31 genes in unidirectional selection regions that might be related to fundamental domestication requirements in pigs. And these genes belong predominantly to categories related to the nervous system, muscle development, and especially to metabolic diseases. In addition, 35 genes, representing different breeding preference, were found under bidirectional selection for the distinct leanness and reproduction traits between CHD and EUD. The convergent genetic changes, contributing physical and morphological adaption, represent the common concerns on pig domestication. And the divergent genetic changes reflect distinct breeding goals between Chinese and European pigs. Using ITPR3, AHR and NMU as examples, we explored and validated how the genetic variations contribute to the phenotype changes.
Collapse
|
213
|
Zhu J, Zhang H, Li SC, Wang C, Kong L, Sun S, Zheng WM, Bu D. Improving protein fold recognition by extracting fold-specific features from predicted residue–residue contacts. Bioinformatics 2017; 33:3749-3757. [DOI: 10.1093/bioinformatics/btx514] [Citation(s) in RCA: 39] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2017] [Accepted: 08/09/2017] [Indexed: 01/05/2023] Open
Affiliation(s)
- Jianwei Zhu
- Key Lab of Intelligent Information Processing, Institute of Computing Technology, Chinese Academy of Sciences, Beijing, China
- University of Chinese Academy of Sciences, Beijing, China
| | - Haicang Zhang
- Key Lab of Intelligent Information Processing, Institute of Computing Technology, Chinese Academy of Sciences, Beijing, China
| | - Shuai Cheng Li
- Department of Computer Science, City University of Hong Kong, Kowloon, Hong Kong
| | - Chao Wang
- Key Lab of Intelligent Information Processing, Institute of Computing Technology, Chinese Academy of Sciences, Beijing, China
| | - Lupeng Kong
- Key Lab of Intelligent Information Processing, Institute of Computing Technology, Chinese Academy of Sciences, Beijing, China
- University of Chinese Academy of Sciences, Beijing, China
| | - Shiwei Sun
- Key Lab of Intelligent Information Processing, Institute of Computing Technology, Chinese Academy of Sciences, Beijing, China
| | - Wei-Mou Zheng
- Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing, China
| | - Dongbo Bu
- Key Lab of Intelligent Information Processing, Institute of Computing Technology, Chinese Academy of Sciences, Beijing, China
| |
Collapse
|
214
|
Farrel A, Guo JT. An efficient algorithm for improving structure-based prediction of transcription factor binding sites. BMC Bioinformatics 2017; 18:342. [PMID: 28715997 PMCID: PMC5514533 DOI: 10.1186/s12859-017-1755-0] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2017] [Accepted: 07/12/2017] [Indexed: 01/07/2023] Open
Abstract
Background Gene expression is regulated by transcription factors binding to specific target DNA sites. Understanding how and where transcription factors bind at genome scale represents an essential step toward our understanding of gene regulation networks. Previously we developed a structure-based method for prediction of transcription factor binding sites using an integrative energy function that combines a knowledge-based multibody potential and two atomic energy terms. While the method performs well, it is not computationally efficient due to the exponential increase in the number of binding sequences to be evaluated for longer binding sites. In this paper, we present an efficient pentamer algorithm by splitting DNA binding sequences into overlapping fragments along with a simplified integrative energy function for transcription factor binding site prediction. Results A DNA binding sequence is split into overlapping pentamers (5 base pairs) for calculating transcription factor-pentamer interaction energy. To combine the results from overlapping pentamer scores, we developed two methods, Kmer-Sum and PWM (Position Weight Matrix) stacking, for full-length binding motif prediction. Our results show that both Kmer-Sum and PWM stacking in the new pentamer approach along with a simplified integrative energy function improved transcription factor binding site prediction accuracy and dramatically reduced computation time, especially for longer binding sites. Conclusion Our new fragment-based pentamer algorithm and simplified energy function improve both efficiency and accuracy. To our knowledge, this is the first fragment-based method for structure-based transcription factor binding sites prediction. Electronic supplementary material The online version of this article (doi:10.1186/s12859-017-1755-0) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Alvin Farrel
- Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, 9201 University City Blvd, Charlotte, NC, 28223, USA
| | - Jun-Tao Guo
- Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, 9201 University City Blvd, Charlotte, NC, 28223, USA.
| |
Collapse
|
215
|
Dalton JAR, Pin JP, Giraldo J. Analysis of positive and negative allosteric modulation in metabotropic glutamate receptors 4 and 5 with a dual ligand. Sci Rep 2017; 7:4944. [PMID: 28694498 PMCID: PMC5504000 DOI: 10.1038/s41598-017-05095-5] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2017] [Accepted: 05/25/2017] [Indexed: 12/24/2022] Open
Abstract
As class C GPCRs and regulators of synaptic activity, human metabotropic glutamate receptors (mGluRs) 4 and 5 are prime targets for allosteric modulation, with mGlu5 inhibition or mGlu4 stimulation potentially treating conditions like chronic pain and Parkinson’s disease. As an allosteric modulator that can bind both receptors, 2-Methyl-6-(phenylethynyl)pyridine (MPEP) is able to negatively modulate mGlu5 or positively modulate mGlu4. At a structural level, how it elicits these responses and how mGluRs undergo activation is unclear. Here, we employ homology modelling and 30 µs of atomistic molecular dynamics (MD) simulations to probe allosteric conformational change in mGlu4 and mGlu5, with and without docked MPEP. Our results identify several structural differences between mGlu4 and mGlu5, as well as key differences responsible for MPEP-mediated positive and negative allosteric modulation, respectively. A novel mechanism of mGlu4 activation is revealed, which may apply to all mGluRs in general. This involves conformational changes in TM3, TM4 and TM5, separation of intracellular loop 2 (ICL2) from ICL1/ICL3, and destabilization of the ionic-lock. On the other hand, mGlu5 experiences little disturbance when MPEP binds, maintaining its inactive state with reduced conformational fluctuation. In addition, when MPEP is absent, a lipid molecule can enter the mGlu5 allosteric pocket.
Collapse
Affiliation(s)
- James A R Dalton
- Laboratory of Molecular Neuropharmacology and Bioinformatics, Institut de Neurociències and Unitat de Bioestadística, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain.,Network Biomedical Research Centre on Mental Health (CIBERSAM), Madrid, Spain
| | - Jean-Philippe Pin
- Institute of Functional Genomics, Université de Montpellier, Unité Mixte de Recherche 5302 CNRS, Montpellier, France.,Unité de recherche U1191, INSERM, Montpellier, France
| | - Jesús Giraldo
- Laboratory of Molecular Neuropharmacology and Bioinformatics, Institut de Neurociències and Unitat de Bioestadística, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain. .,Network Biomedical Research Centre on Mental Health (CIBERSAM), Madrid, Spain.
| |
Collapse
|
216
|
Li B, Mendenhall J, Nguyen ED, Weiner BE, Fischer AW, Meiler J. Improving prediction of helix-helix packing in membrane proteins using predicted contact numbers as restraints. Proteins 2017; 85:1212-1221. [PMID: 28263405 PMCID: PMC5476507 DOI: 10.1002/prot.25281] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2016] [Revised: 01/20/2017] [Accepted: 02/17/2017] [Indexed: 01/21/2023]
Abstract
One of the challenging problems in tertiary structure prediction of helical membrane proteins (HMPs) is the determination of rotation of α-helices around the helix normal. Incorrect prediction of helix rotations substantially disrupts native residue-residue contacts while inducing only a relatively small effect on the overall fold. We previously developed a method for predicting residue contact numbers (CNs), which measure the local packing density of residues within the protein tertiary structure. In this study, we tested the idea of incorporating predicted CNs as restraints to guide the sampling of helix rotation. For a benchmark set of 15 HMPs with simple to rather complicated folds, the average contact recovery (CR) of best-sampled models was improved for all targets, the likelihood of sampling models with CR greater than 20% was increased for 13 targets, and the average RMSD100 of best-sampled models was improved for 12 targets. This study demonstrated that explicit incorporation of CNs as restraints improves the prediction of helix-helix packing. Proteins 2017; 85:1212-1221. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Bian Li
- Department of Chemistry, Vanderbilt University, Nashville, TN 37232, USA
- Center for Structural Biology, Vanderbilt University, Nashville, TN 37232, USA
| | - Jeffrey Mendenhall
- Department of Chemistry, Vanderbilt University, Nashville, TN 37232, USA
- Center for Structural Biology, Vanderbilt University, Nashville, TN 37232, USA
| | | | - Brian E. Weiner
- Department of Chemistry, Vanderbilt University, Nashville, TN 37232, USA
- Center for Structural Biology, Vanderbilt University, Nashville, TN 37232, USA
| | - Axel W. Fischer
- Department of Chemistry, Vanderbilt University, Nashville, TN 37232, USA
- Center for Structural Biology, Vanderbilt University, Nashville, TN 37232, USA
| | - Jens Meiler
- Department of Chemistry, Vanderbilt University, Nashville, TN 37232, USA
- Center for Structural Biology, Vanderbilt University, Nashville, TN 37232, USA
| |
Collapse
|
217
|
Alford RF, Leaver-Fay A, Jeliazkov JR, O’Meara MJ, DiMaio FP, Park H, Shapovalov MV, Renfrew PD, Mulligan VK, Kappel K, Labonte JW, Pacella MS, Bonneau R, Bradley P, Dunbrack RL, Das R, Baker D, Kuhlman B, Kortemme T, Gray JJ. The Rosetta All-Atom Energy Function for Macromolecular Modeling and Design. J Chem Theory Comput 2017; 13:3031-3048. [PMID: 28430426 PMCID: PMC5717763 DOI: 10.1021/acs.jctc.7b00125] [Citation(s) in RCA: 795] [Impact Index Per Article: 113.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Over the past decade, the Rosetta biomolecular modeling suite has informed diverse biological questions and engineering challenges ranging from interpretation of low-resolution structural data to design of nanomaterials, protein therapeutics, and vaccines. Central to Rosetta's success is the energy function: a model parametrized from small-molecule and X-ray crystal structure data used to approximate the energy associated with each biomolecule conformation. This paper describes the mathematical models and physical concepts that underlie the latest Rosetta energy function, called the Rosetta Energy Function 2015 (REF15). Applying these concepts, we explain how to use Rosetta energies to identify and analyze the features of biomolecular models. Finally, we discuss the latest advances in the energy function that extend its capabilities from soluble proteins to also include membrane proteins, peptides containing noncanonical amino acids, small molecules, carbohydrates, nucleic acids, and other macromolecules.
Collapse
Affiliation(s)
- Rebecca F. Alford
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, 3400 North Charles Street, Baltimore, Maryland 21218, United States
| | - Andrew Leaver-Fay
- Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, 120 Mason Farm Road, Chapel Hill, North Carolina 27599, United States
| | - Jeliazko R. Jeliazkov
- Program in Molecular Biophysics, Johns Hopkins University, 3400 North Charles Street, Baltimore, Maryland 21218, United States
| | - Matthew J. O’Meara
- Department of Pharmaceutical Chemistry, University of California at San Francisco, 1700 Fourth Street, San Francisco, California 94158, United States
| | - Frank P. DiMaio
- Department of Biochemistry, University of Washington, J-Wing Health Sciences Building, Box 357350, Seattle, Washington 98195, United States
| | - Hahnbeom Park
- Department of Biochemistry, University of Washington, Molecular Engineering and Sciences, Box 357350, 4000 15 Ave NE, Seattle, Washington 98195, United States
| | - Maxim V. Shapovalov
- Institute for Cancer Research, Fox Chase Cancer Center, 333 Cottman Avenue, Philadelphia, Pennsylvania 19111, United States
| | - P. Douglas Renfrew
- Department of Biology, Center for Genomics and Systems Biology, New York University, 100 Washington Square East, New York, New York 10003
- Center for Computational Biology, Flatiron Institute, Simons Foundation, 162 5 Avenue, New York, New York 10010, United States
| | - Vikram K. Mulligan
- Department of Biochemistry, University of Washington, Molecular Engineering and Sciences, Box 357350, 4000 15 Ave NE, Seattle, Washington 98195, United States
| | - Kalli Kappel
- Biophysics Program, Stanford University, 450 Serra Mall, Stanford, California 94305, United States
| | - Jason W. Labonte
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, 3400 North Charles Street, Baltimore, Maryland 21218, United States
| | - Michael S. Pacella
- Department of Biomedical Engineering, Johns Hopkins University, 3400 North Charles Street, Baltimore, Maryland 21218, United States
| | - Richard Bonneau
- Department of Biology, Center for Genomics and Systems Biology, New York University, 100 Washington Square East, New York, New York 10003
- Center for Computational Biology, Flatiron Institute, Simons Foundation, 162 5 Avenue, New York, New York 10010, United States
| | - Philip Bradley
- Computational Biology Program, Fred Hutchinson Cancer Research Center, 1100 Fairview Avenue North, Seattle, Washington 98109, United States
| | - Roland L. Dunbrack
- Institute for Cancer Research, Fox Chase Cancer Center, 333 Cottman Avenue, Philadelphia, Pennsylvania 19111, United States
| | - Rhiju Das
- Biophysics Program, Stanford University, 450 Serra Mall, Stanford, California 94305, United States
| | - David Baker
- Department of Biochemistry, University of Washington, Molecular Engineering and Sciences, Box 357350, 4000 15 Ave NE, Seattle, Washington 98195, United States
| | - Brian Kuhlman
- Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, 120 Mason Farm Road, Chapel Hill, North Carolina 27599, United States
| | - Tanja Kortemme
- Department of Bioengineering and Therapeutic Sciences, University of California at San Francisco, San Francisco, California 94158, United States
| | - Jeffrey J. Gray
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, 3400 North Charles Street, Baltimore, Maryland 21218, United States
- Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, 120 Mason Farm Road, Chapel Hill, North Carolina 27599, United States
| |
Collapse
|
218
|
Mackenzie CO, Grigoryan G. Protein structural motifs in prediction and design. Curr Opin Struct Biol 2017; 44:161-167. [PMID: 28460216 PMCID: PMC5513761 DOI: 10.1016/j.sbi.2017.03.012] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2016] [Revised: 03/18/2017] [Accepted: 03/28/2017] [Indexed: 01/11/2023]
Abstract
The Protein Data Bank (PDB) has been an integral resource for shaping our fundamental understanding of protein structure and for the advancement of such applications as protein design and structure prediction. Over the years, information from the PDB has been used to generate models ranging from specific structural mechanisms to general statistical potentials. With accumulating structural data, it has become possible to mine for more complete and complex structural observations, deducing more accurate generalizations. Motif libraries, which capture recurring structural features along with their sequence preferences, have exposed modularity in the structural universe and found successful application in various problems of structural biology. Here we summarize recent achievements in this arena, focusing on subdomain level structural patterns and their applications to protein design and structure prediction, and suggest promising future directions as the structural database continues to grow.
Collapse
Affiliation(s)
- Craig O Mackenzie
- Institute for Quantitative Biomedical Sciences, Dartmouth College, Hanover, NH 03755, United States
| | - Gevorg Grigoryan
- Institute for Quantitative Biomedical Sciences, Dartmouth College, Hanover, NH 03755, United States; Department of Computer Science, Dartmouth College, Hanover, NH 03755, United States.
| |
Collapse
|
219
|
Sequence statistics of tertiary structural motifs reflect protein stability. PLoS One 2017; 12:e0178272. [PMID: 28552940 PMCID: PMC5446159 DOI: 10.1371/journal.pone.0178272] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2017] [Accepted: 05/10/2017] [Indexed: 11/19/2022] Open
Abstract
The Protein Data Bank (PDB) has been a key resource for learning general rules of sequence-structure relationships in proteins. Quantitative insights have been gained by defining geometric descriptors of structure (e.g., distances, dihedral angles, solvent exposure, etc.) and observing their distributions and sequence preferences. Here we argue that as the PDB continues to grow, it may become unnecessary to reduce structure into a set of elementary descriptors. Instead, it could be possible to deduce quantitative sequence-structure relationships in the context of precisely-defined complex structural motifs by mining the PDB for closely matching backbone geometries. To validate this idea, we turned to the the task of predicting changes in protein stability upon amino-acid substitution—a difficult problem of broad significance. We defined non-contiguous tertiary motifs (TERMs) around a protein site of interest and extracted sequence preferences from ensembles of closely-matching substructures in the PDB to predict mutational stability changes at the site, ΔΔGm. We demonstrate that these ensemble statistics predict ΔΔGm on par with state-of-the-art statistical and machine-learning methods on large thermodynamic datasets, and outperform these, along with a leading structure-based modeling approach, when tested in the context of unbiased diverse mutations. Further, we show that the performance of the TERM-based method is directly related to the amount of available relevant structural data, automatically improving with the growing PDB. This enables a means of estimating prediction accuracy. Our results clearly demonstrate that: 1) statistics of non-contiguous structural motifs in the PDB encode fundamental sequence-structure relationships related to protein thermodynamic stability, and 2) the PDB is now large enough that such statistics are already useful in practice, with their accuracy expected to continue increasing as the database grows. These observations suggest new ways of using structural data towards addressing problems of computational structural biology.
Collapse
|
220
|
Luiz DP, Almeida JF, Goulart LR, Nicolau-Junior N, Ueira-Vieira C. Heterologous expression of abaecin peptide from Apis mellifera in Pichia pastoris. Microb Cell Fact 2017; 16:76. [PMID: 28468662 PMCID: PMC5414229 DOI: 10.1186/s12934-017-0689-6] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2016] [Accepted: 04/21/2017] [Indexed: 02/03/2023] Open
Abstract
Background Antimicrobial peptides (AMPs) are the first line of host immune defense against pathogens. Among AMPs from the honeybee Apis mellifera, abaecin is a major broad-spectrum antibacterial proline-enriched cationic peptide. Results For heterologous expression of abaecin in Pichia pastoris, we designed an ORF with HisTag, and the codon usage was optimized. The gene was chemically synthetized and cloned in the pUC57 vector. The new ORF was sub-cloned in the pPIC9 expression vector and transformed into P. pastoris. After selection of positive clones, the expression was induced by methanol. The supernatant was analyzed at different times to determine the optimal time for the recombinant peptide expression. As a proof-of-concept, Escherichia coli was co-incubated with the recombinant peptide to verify its antimicrobial potential. Discussion Briefly, the recombinant Abaecin (rAbaecin) has efficiently decreased E. coli growth (P < 0.05) through an in vitro assay, and may be considered as a novel therapeutic agent that may complement other conventional antibiotic therapies.
Collapse
Affiliation(s)
- Denis Prudencio Luiz
- Genetics Laboratory, Institute of Genetics and Biochemistry, Federal University of Uberlândia, 1720 Pará, Uberlândia, MG, 38400902, Brazil.
| | - Juliana Franco Almeida
- Nanobiotechnology Laboratory, Institute of Genetics and Biochemistry, Federal University of Uberlândia, 1720 Pará, Uberlândia, MG, 38400902, Brazil
| | - Luiz Ricardo Goulart
- Nanobiotechnology Laboratory, Institute of Genetics and Biochemistry, Federal University of Uberlândia, 1720 Pará, Uberlândia, MG, 38400902, Brazil
| | - Nilson Nicolau-Junior
- Genetics Laboratory, Institute of Genetics and Biochemistry, Federal University of Uberlândia, 1720 Pará, Uberlândia, MG, 38400902, Brazil
| | - Carlos Ueira-Vieira
- Genetics Laboratory, Institute of Genetics and Biochemistry, Federal University of Uberlândia, 1720 Pará, Uberlândia, MG, 38400902, Brazil
| |
Collapse
|
221
|
Porto W, Pires A, Franco O. Computational tools for exploring sequence databases as a resource for antimicrobial peptides. Biotechnol Adv 2017; 35:337-349. [DOI: 10.1016/j.biotechadv.2017.02.001] [Citation(s) in RCA: 79] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2016] [Revised: 01/12/2017] [Accepted: 02/08/2017] [Indexed: 12/22/2022]
|
222
|
Computational Prediction of the Heterodimeric and Higher-Order Structure of gpE1/gpE2 Envelope Glycoproteins Encoded by Hepatitis C Virus. J Virol 2017; 91:JVI.02309-16. [PMID: 28148799 DOI: 10.1128/jvi.02309-16] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2016] [Accepted: 01/25/2017] [Indexed: 12/24/2022] Open
Abstract
Despite the recent success of newly developed direct-acting antivirals against hepatitis C, the disease continues to be a global health threat due to the lack of diagnosis of most carriers and the high cost of treatment. The heterodimer formed by glycoproteins E1 and E2 within the hepatitis C virus (HCV) lipid envelope is a potential vaccine candidate and antiviral target. While the structure of E1/E2 has not yet been resolved, partial crystal structures of the E1 and E2 ectodomains have been determined. The unresolved parts of the structure are within the realm of what can be modeled with current computational modeling tools. Furthermore, a variety of additional experimental data is available to support computational predictions of E1/E2 structure, such as data from antibody binding studies, cryo-electron microscopy (cryo-EM), mutational analyses, peptide binding analysis, linker-scanning mutagenesis, and nuclear magnetic resonance (NMR) studies. In accordance with these rich experimental data, we have built an in silico model of the full-length E1/E2 heterodimer. Our model supports that E1/E2 assembles into a trimer, which was previously suggested from a study by Falson and coworkers (P. Falson, B. Bartosch, K. Alsaleh, B. A. Tews, A. Loquet, Y. Ciczora, L. Riva, C. Montigny, C. Montpellier, G. Duverlie, E. I. Pecheur, M. le Maire, F. L. Cosset, J. Dubuisson, and F. Penin, J. Virol. 89:10333-10346, 2015, https://doi.org/10.1128/JVI.00991-15). Size exclusion chromatography and Western blotting data obtained by using purified recombinant E1/E2 support our hypothesis. Our model suggests that during virus assembly, the trimer of E1/E2 may be further assembled into a pentamer, with 12 pentamers comprising a single HCV virion. We anticipate that this new model will provide a useful framework for HCV envelope structure and the development of antiviral strategies.IMPORTANCE One hundred fifty million people have been estimated to be infected with hepatitis C virus, and many more are at risk for infection. A better understanding of the structure of the HCV envelope, which is responsible for attachment and fusion, could aid in the development of a vaccine and/or new treatments for this disease. We draw upon computational techniques to predict a full-length model of the E1/E2 heterodimer based on the partial crystal structures of the envelope glycoproteins E1 and E2. E1/E2 has been widely studied experimentally, and this provides valuable data, which has assisted us in our modeling. Our proposed structure is used to suggest the organization of the HCV envelope. We also present new experimental data from size exclusion chromatography that support our computational prediction of a trimeric oligomeric state of E1/E2.
Collapse
|
223
|
Davies HA, Rigden DJ, Phelan MM, Madine J. Probing Medin Monomer Structure and its Amyloid Nucleation Using 13C-Direct Detection NMR in Combination with Structural Bioinformatics. Sci Rep 2017; 7:45224. [PMID: 28327552 PMCID: PMC5361114 DOI: 10.1038/srep45224] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2016] [Accepted: 02/20/2017] [Indexed: 12/21/2022] Open
Abstract
Aortic medial amyloid is the most prevalent amyloid found to date, but remarkably little is known about it. It is characterised by aberrant deposition of a 5.4 kDa protein called medin within the medial layer of large arteries. Here we employ a combined approach of ab initio protein modelling and 13C-direct detection NMR to generate a model for soluble monomeric medin comprising a stable core of three β-strands and shorter more labile strands at the termini. Molecular dynamics simulations suggested that detachment of the short, C-terminal β-strand from the soluble fold exposes key amyloidogenic regions as a potential site of nucleation enabling dimerisation and subsequent fibril formation. This mechanism resembles models proposed for several other amyloidogenic proteins suggesting that despite variations in sequence and protomer structure these proteins may share a common pathway for amyloid nucleation and subsequent protofibril and fibril formation.
Collapse
Affiliation(s)
- Hannah A. Davies
- Institute of Integrative Biology, University of Liverpool, Biosciences Building, Crown Street, L69 7ZB, UK
| | - Daniel J. Rigden
- Institute of Integrative Biology, University of Liverpool, Biosciences Building, Crown Street, L69 7ZB, UK
| | - Marie M. Phelan
- Institute of Integrative Biology, University of Liverpool, Biosciences Building, Crown Street, L69 7ZB, UK
| | - Jillian Madine
- Institute of Integrative Biology, University of Liverpool, Biosciences Building, Crown Street, L69 7ZB, UK
| |
Collapse
|
224
|
Poosarla VG, Li T, Goh BC, Schulten K, Wood TK, Maranas CD. Computational de novo design of antibodies binding to a peptide with high affinity. Biotechnol Bioeng 2017; 114:1331-1342. [PMID: 28059445 DOI: 10.1002/bit.26244] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2016] [Revised: 12/23/2016] [Accepted: 01/05/2017] [Indexed: 12/24/2022]
Abstract
Antibody drugs play a critical role in infectious diseases, cancer, autoimmune diseases, and inflammation. However, experimental methods for the generation of therapeutic antibodies such as using immunized mice or directed evolution remain time consuming and cannot target a specific antigen epitope. Here, we describe the application of a computational framework called OptMAVEn combined with molecular dynamics to de novo design antibodies. Our reference system is antibody 2D10, a single-chain antibody (scFv) that recognizes the dodecapeptide DVFYPYPYASGS, a peptide mimic of mannose-containing carbohydrates. Five de novo designed scFvs sharing less than 75% sequence similarity to all existing natural antibody sequences were generated using OptMAVEn and their binding to the dodecapeptide was experimentally characterized by biolayer interferometry and isothermal titration calorimetry. Among them, three scFvs show binding affinity to the dodecapeptide at the nM level. Critically, these de novo designed scFvs exhibit considerably diverse modeled binding modes with the dodecapeptide. The results demonstrate the potential of OptMAVEn for the de novo design of thermally and conformationally stable antibodies with high binding affinity to antigens and encourage the targeting of other antigen targets in the future. Biotechnol. Bioeng. 2017;114: 1331-1342. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Venkata Giridhar Poosarla
- Department of Chemical Engineering, University Park, Pennsylvania, 16802.,Department of Biochemistry and Molecular Biology, Pennsylvania State University, University Park, Pennsylvania, 16802
| | - Tong Li
- Department of Chemical Engineering, University Park, Pennsylvania, 16802
| | - Boon Chong Goh
- Department of Physics and Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, Illinois, 61801
| | - Klaus Schulten
- Department of Physics and Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, Illinois, 61801
| | - Thomas K Wood
- Department of Chemical Engineering, University Park, Pennsylvania, 16802.,Department of Biochemistry and Molecular Biology, Pennsylvania State University, University Park, Pennsylvania, 16802
| | - Costas D Maranas
- Department of Chemical Engineering, University Park, Pennsylvania, 16802
| |
Collapse
|
225
|
Hoffmann T, Marion A, Antes I. DynaDom: structure-based prediction of T cell receptor inter-domain and T cell receptor-peptide-MHC (class I) association angles. BMC STRUCTURAL BIOLOGY 2017; 17:2. [PMID: 28148269 PMCID: PMC5289058 DOI: 10.1186/s12900-016-0071-7] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/04/2016] [Accepted: 12/29/2016] [Indexed: 11/22/2022]
Abstract
Background T cell receptor (TCR) molecules are involved in the adaptive immune response as they distinguish between self- and foreign-peptides, presented in major histocompatibility complex molecules (pMHC). Former studies showed that the association angles of the TCR variable domains (Vα/Vβ) can differ significantly and change upon binding to the pMHC complex. These changes can be described as a rotation of the domains around a general Center of Rotation, characterized by the interaction of two highly conserved glutamine residues. Methods We developed a computational method, DynaDom, for the prediction of TCR Vα/Vβ inter-domain and TCR/pMHC orientations in TCRpMHC complexes, which allows predicting the orientation of multiple protein-domains. In addition, we implemented a new approach to predict the correct orientation of the carboxamide endgroups in glutamine and asparagine residues, which can also be used as an external, independent tool. Results The approach was evaluated for the remodeling of 75 and 53 experimental structures of TCR and TCRpMHC (class I) complexes, respectively. We show that the DynaDom method predicts the correct orientation of the TCR Vα/Vβ angles in 96 and 89% of the cases, for the poses with the best RMSD and best interaction energy, respectively. For the concurrent prediction of the TCR Vα/Vβ and pMHC orientations, the respective rates reached 74 and 72%. Through an exhaustive analysis, we could show that the pMHC placement can be further improved by a straightforward, yet very time intensive extension of the current approach. Conclusions The results obtained in the present remodeling study prove the suitability of our approach for interdomain-angle optimization. In addition, the high prediction rate obtained specifically for the energetically highest ranked poses further demonstrates that our method is a powerful candidate for blind prediction. Therefore it should be well suited as part of any accurate atomistic modeling pipeline for TCRpMHC complexes and potentially other large molecular assemblies. Electronic supplementary material The online version of this article (doi:10.1186/s12900-016-0071-7) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Thomas Hoffmann
- Department of Biosciences and Center for Integrated Protein Science Munich, Technische Universität München, Emil-Erlenmeyer-Forum 8, 85354, Freising, Germany
| | - Antoine Marion
- Department of Biosciences and Center for Integrated Protein Science Munich, Technische Universität München, Emil-Erlenmeyer-Forum 8, 85354, Freising, Germany
| | - Iris Antes
- Department of Biosciences and Center for Integrated Protein Science Munich, Technische Universität München, Emil-Erlenmeyer-Forum 8, 85354, Freising, Germany.
| |
Collapse
|
226
|
Critical Features of Fragment Libraries for Protein Structure Prediction. PLoS One 2017; 12:e0170131. [PMID: 28085928 PMCID: PMC5235372 DOI: 10.1371/journal.pone.0170131] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2016] [Accepted: 12/29/2016] [Indexed: 11/19/2022] Open
Abstract
The use of fragment libraries is a popular approach among protein structure prediction methods and has proven to substantially improve the quality of predicted structures. However, some vital aspects of a fragment library that influence the accuracy of modeling a native structure remain to be determined. This study investigates some of these features. Particularly, we analyze the effect of using secondary structure prediction guiding fragments selection, different fragments sizes and the effect of structural clustering of fragments within libraries. To have a clearer view of how these factors affect protein structure prediction, we isolated the process of model building by fragment assembly from some common limitations associated with prediction methods, e.g., imprecise energy functions and optimization algorithms, by employing an exact structure-based objective function under a greedy algorithm. Our results indicate that shorter fragments reproduce the native structure more accurately than the longer. Libraries composed of multiple fragment lengths generate even better structures, where longer fragments show to be more useful at the beginning of the simulations. The use of many different fragment sizes shows little improvement when compared to predictions carried out with libraries that comprise only three different fragment sizes. Models obtained from libraries built using only sequence similarity are, on average, better than those built with a secondary structure prediction bias. However, we found that the use of secondary structure prediction allows greater reduction of the search space, which is invaluable for prediction methods. The results of this study can be critical guidelines for the use of fragment libraries in protein structure prediction.
Collapse
|
227
|
Honey Bee Deformed Wing Virus Structures Reveal that Conformational Changes Accompany Genome Release. J Virol 2017; 91:JVI.01795-16. [PMID: 27852845 DOI: 10.1128/jvi.01795-16] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2016] [Accepted: 11/02/2016] [Indexed: 11/20/2022] Open
Abstract
The picornavirus-like deformed wing virus (DWV) has been directly linked to colony collapse; however, little is known about the mechanisms of host attachment or entry for DWV or its molecular and structural details. Here we report the three-dimensional (3-D) structures of DWV capsids isolated from infected honey bees, including the immature procapsid, the genome-filled virion, the putative entry intermediate (A-particle), and the empty capsid that remains after genome release. The capsids are decorated by large spikes around the 5-fold vertices. The 5-fold spikes had an open flower-like conformation for the procapsid and genome-filled capsids, whereas the putative A-particle and empty capsids that had released the genome had a closed tube-like spike conformation. Between the two conformations, the spikes undergo a significant hinge-like movement that we predicted using a Robetta model of the structure comprising the spike. We conclude that the spike structures likely serve a function during host entry, changing conformation to release the genome, and that the genome may escape from a 5-fold vertex to initiate infection. Finally, the structures illustrate that, similarly to picornaviruses, DWV forms alternate particle conformations implicated in assembly, host attachment, and RNA release. IMPORTANCE Honey bees are critical for global agriculture, but dramatic losses of entire hives have been reported in numerous countries since 2006. Deformed wing virus (DWV) and infestation with the ectoparasitic mite Varroa destructor have been linked to colony collapse disorder. DWV was purified from infected adult worker bees to pursue biochemical and structural studies that allowed the first glimpse into the conformational changes that may be required during transmission and genome release for DWV.
Collapse
|
228
|
Annotation of Alternatively Spliced Proteins and Transcripts with Protein-Folding Algorithms and Isoform-Level Functional Networks. Methods Mol Biol 2017; 1558:415-436. [PMID: 28150250 DOI: 10.1007/978-1-4939-6783-4_20] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Tens of thousands of splice isoforms of proteins have been catalogued as predicted sequences from transcripts in humans and other species. Relatively few have been characterized biochemically or structurally. With the extensive development of protein bioinformatics, the characterization and modeling of isoform features, isoform functions, and isoform-level networks have advanced notably. Here we present applications of the I-TASSER family of algorithms for folding and functional predictions and the IsoFunc, MIsoMine, and Hisonet data resources for isoform-level analyses of network and pathway-based functional predictions and protein-protein interactions. Hopefully, predictions and insights from protein bioinformatics will stimulate many experimental validation studies.
Collapse
|
229
|
EpiSweep: Computationally Driven Reengineering of Therapeutic Proteins to Reduce Immunogenicity While Maintaining Function. Methods Mol Biol 2017; 1529:375-398. [PMID: 27914063 DOI: 10.1007/978-1-4939-6637-0_20] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Therapeutic proteins are yielding ever more advanced and efficacious new drugs, but the biological origins of these highly effective therapeutics render them subject to immune surveillance within the patient's body. When recognized by the immune system as a foreign agent, protein drugs elicit a coordinated response that can manifest a range of clinical complications including rapid drug clearance, loss of functionality and efficacy, delayed infusion-like allergic reactions, more serious anaphylactic shock, and even induced auto-immunity. It is thus often necessary to deimmunize an exogenous protein in order to enable its clinical application; critically, the deimmunization process must also maintain the desired therapeutic activity.To meet the growing need for effective, efficient, and broadly applicable protein deimmunization technologies, we have developed the EpiSweep suite of protein design algorithms. EpiSweep seamlessly integrates computational prediction of immunogenic T cell epitopes with sequence- or structure-based assessment of the impacts of mutations on protein stability and function, in order to select combinations of mutations that make Pareto optimal trade-offs between the competing goals of low immunogenicity and high-level function. The methods are applicable both to the design of individual functionally deimmunized variants as well as the design of combinatorial libraries enriched in functionally deimmunized variants. After validating EpiSweep in a series of retrospective case studies providing comparisons to conventional approaches to T cell epitope deletion, we have experimentally demonstrated it to be highly effective in prospective application to deimmunization of a number of different therapeutic candidates. We conclude that our broadly applicable computational protein design algorithms guide the engineer towards the most promising deimmunized therapeutic candidates, and thereby have the potential to accelerate development of new protein drugs by shortening time frames and improving hit rates.
Collapse
|
230
|
Priya P, Kesheri M, Sinha RP, Kanchan S. Molecular Dynamics Simulations for Biological Systems. PHARMACEUTICAL SCIENCES 2017. [DOI: 10.4018/978-1-5225-1762-7.ch040] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022] Open
Abstract
Molecular dynamics simulation is an important tool to capture the dynamicity of biological molecule and the atomistic insights. These insights are helpful to explore biological functions. Molecular dynamics simulation from femto seconds to milli seconds scale give a large ensemble of conformations that can reveal many biological mysteries. The main focus of the chapter is to throw light on theories, requirement of molecular dynamics for biological studies and application of molecular dynamics simulations. Molecular dynamics simulations are widely used to study protein-protein interaction, protein-ligand docking, effects of mutation on interactions, protein folding and flexibility of the biological molecules. This chapter also deals with various methods/algorithms of protein tertiary structure prediction, their strengths and weaknesses.
Collapse
|
231
|
Abstract
Protein-protein interactions play critical roles in essentially every cellular process. These interactions are often mediated by protein interaction domains that enable proteins to recognize their interaction partners, often by binding to short peptide motifs. For example, PDZ domains, which are among the most common protein interaction domains in the human proteome, recognize specific linear peptide sequences that are often at the C-terminus of other proteins. Determining the set of peptide sequences that a protein interaction domain binds, or it's "peptide specificity," is crucial for understanding its cellular function, and predicting how mutations impact peptide specificity is important for elucidating the mechanisms underlying human diseases. Moreover, engineering novel cellular functions for synthetic biology applications, such as the biosynthesis of biofuels or drugs, requires the design of protein interaction specificity to avoid crosstalk with native metabolic and signaling pathways. The ability to accurately predict and design protein-peptide interaction specificity is therefore critical for understanding and engineering biological function. One approach that has recently been employed toward accomplishing this goal is computational protein design. This chapter provides an overview of recent methodological advances in computational protein design and highlights examples of how these advances can enable increased accuracy in predicting and designing peptide specificity.
Collapse
Affiliation(s)
- Noah Ollikainen
- Division of Biology and Biological Engineering, California Institute of Technology, 1200 East California Blvd., Pasadena, CA, 91125, USA
| |
Collapse
|
232
|
Leelananda SP, Lindert S. Computational methods in drug discovery. Beilstein J Org Chem 2016; 12:2694-2718. [PMID: 28144341 PMCID: PMC5238551 DOI: 10.3762/bjoc.12.267] [Citation(s) in RCA: 285] [Impact Index Per Article: 35.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2016] [Accepted: 11/22/2016] [Indexed: 12/11/2022] Open
Abstract
The process for drug discovery and development is challenging, time consuming and expensive. Computer-aided drug discovery (CADD) tools can act as a virtual shortcut, assisting in the expedition of this long process and potentially reducing the cost of research and development. Today CADD has become an effective and indispensable tool in therapeutic development. The human genome project has made available a substantial amount of sequence data that can be used in various drug discovery projects. Additionally, increasing knowledge of biological structures, as well as increasing computer power have made it possible to use computational methods effectively in various phases of the drug discovery and development pipeline. The importance of in silico tools is greater than ever before and has advanced pharmaceutical research. Here we present an overview of computational methods used in different facets of drug discovery and highlight some of the recent successes. In this review, both structure-based and ligand-based drug discovery methods are discussed. Advances in virtual high-throughput screening, protein structure prediction methods, protein-ligand docking, pharmacophore modeling and QSAR techniques are reviewed.
Collapse
Affiliation(s)
- Sumudu P Leelananda
- Department of Chemistry and Biochemistry, Ohio State University, Columbus, OH 43210, USA
| | - Steffen Lindert
- Department of Chemistry and Biochemistry, Ohio State University, Columbus, OH 43210, USA
| |
Collapse
|
233
|
Structural studies of RNA-protein complexes: A hybrid approach involving hydrodynamics, scattering, and computational methods. Methods 2016; 118-119:146-162. [PMID: 27939506 DOI: 10.1016/j.ymeth.2016.12.002] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2016] [Revised: 12/01/2016] [Accepted: 12/05/2016] [Indexed: 01/01/2023] Open
Abstract
The diverse functional cellular roles played by ribonucleic acids (RNA) have emphasized the need to develop rapid and accurate methodologies to elucidate the relationship between the structure and function of RNA. Structural biology tools such as X-ray crystallography and Nuclear Magnetic Resonance are highly useful methods to obtain atomic-level resolution models of macromolecules. However, both methods have sample, time, and technical limitations that prevent their application to a number of macromolecules of interest. An emerging alternative to high-resolution structural techniques is to employ a hybrid approach that combines low-resolution shape information about macromolecules and their complexes from experimental hydrodynamic (e.g. analytical ultracentrifugation) and solution scattering measurements (e.g., solution X-ray or neutron scattering), with computational modeling to obtain atomic-level models. While promising, scattering methods rely on aggregation-free, monodispersed preparations and therefore the careful development of a quality control pipeline is fundamental to an unbiased and reliable structural determination. This review article describes hydrodynamic techniques that are highly valuable for homogeneity studies, scattering techniques useful to study the low-resolution shape, and strategies for computational modeling to obtain high-resolution 3D structural models of RNAs, proteins, and RNA-protein complexes.
Collapse
|
234
|
Cao R, Bhattacharya D, Hou J, Cheng J. DeepQA: improving the estimation of single protein model quality with deep belief networks. BMC Bioinformatics 2016; 17:495. [PMID: 27919220 PMCID: PMC5139030 DOI: 10.1186/s12859-016-1405-y] [Citation(s) in RCA: 112] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2016] [Accepted: 12/01/2016] [Indexed: 01/02/2023] Open
Abstract
BACKGROUND Protein quality assessment (QA) useful for ranking and selecting protein models has long been viewed as one of the major challenges for protein tertiary structure prediction. Especially, estimating the quality of a single protein model, which is important for selecting a few good models out of a large model pool consisting of mostly low-quality models, is still a largely unsolved problem. RESULTS We introduce a novel single-model quality assessment method DeepQA based on deep belief network that utilizes a number of selected features describing the quality of a model from different perspectives, such as energy, physio-chemical characteristics, and structural information. The deep belief network is trained on several large datasets consisting of models from the Critical Assessment of Protein Structure Prediction (CASP) experiments, several publicly available datasets, and models generated by our in-house ab initio method. Our experiments demonstrate that deep belief network has better performance compared to Support Vector Machines and Neural Networks on the protein model quality assessment problem, and our method DeepQA achieves the state-of-the-art performance on CASP11 dataset. It also outperformed two well-established methods in selecting good outlier models from a large set of models of mostly low quality generated by ab initio modeling methods. CONCLUSION DeepQA is a useful deep learning tool for protein single model quality assessment and protein structure prediction. The source code, executable, document and training/test datasets of DeepQA for Linux is freely available to non-commercial users at http://cactus.rnet.missouri.edu/DeepQA/ .
Collapse
Affiliation(s)
- Renzhi Cao
- Department of Computer Science, Pacific Lutheran University, Tacoma, WA, 98447, USA
| | - Debswapna Bhattacharya
- Department of Electrical Engineering and Computer Science, Wichita State University, Wichita, KS, 67260, USA
| | - Jie Hou
- Department of Computer Science, University of Missouri, Columbia, MO, 65211, USA
| | - Jianlin Cheng
- Department of Computer Science, University of Missouri, Columbia, MO, 65211, USA. .,Informatics Institute, University of Missouri, Columbia, MO, 65211, USA.
| |
Collapse
|
235
|
Rehman Z, Sadia H, Fahim A, Niazi UHK, Azam MZ. Mutational analysis and interactions of HBV preS1 with asialoglycoprotein receptor. Future Virol 2016. [DOI: 10.2217/fvl-2016-0083] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Background: The mutations in preS1 of a large envelop protein of HBV may have profound implications in HBV receptor binding to hepatocytes and subsequent entry of the virus into host cells. Aims: This study aimed to identify the mutations in preS1 region and the receptor binding interactions of preS1 with hepatocytes. Methods: The mutations were searched through direct sequencing of the preS1 region. Sequence analysis was done through ClustalX and Jalview. Ab initio modeling of preS1 was done through Rosetta and QUARK followed by glycosylation of best model of preS1. Finally the interactions of preS1 with ASGPR was studied using PatchDock and analysis was done using MOE and pyMol. Results: Sequence comparison revealed changes in the preS1 region. Ab initio modeling results showed that preS1 is an overall unstructured protein with the presence of three structural motifs. Docking of preS1 with asialoglycoprotein receptor showed mostly hydrophobic interactions. Conclusion: In conclusion, preS1 sequences from Pakistani isolates were found to be 90% conserved and the predicted structure of preS1 was near to native structure.
Collapse
Affiliation(s)
- Zaira Rehman
- Healthcare Biotechnology Department, Atta-ur-Rahman School of Applied Biosciences (ASAB), National University of Sciences & Technology, Islamabad, Pakistan
| | - Hajra Sadia
- Healthcare Biotechnology Department, Atta-ur-Rahman School of Applied Biosciences (ASAB), National University of Sciences & Technology, Islamabad, Pakistan
| | - Ammad Fahim
- Healthcare Biotechnology Department, Atta-ur-Rahman School of Applied Biosciences (ASAB), National University of Sciences & Technology, Islamabad, Pakistan
| | - Umer HK Niazi
- IBERS, Aberystwyth University, Edward Llwyd Building, Penglais Campus, Aberystwyth, Ceredigion, Wales SY23 3FG, UK
| | - Muhammad Z Azam
- National Institute of Liver & Gastrointestinal Diseases, Dow University of Health Sciences, Karachi, Pakistan
| |
Collapse
|
236
|
|
237
|
Abstract
Here, we systematically decompose the known protein structural universe into its basic elements, which we dub tertiary structural motifs (TERMs). A TERM is a compact backbone fragment that captures the secondary, tertiary, and quaternary environments around a given residue, comprising one or more disjoint segments (three on average). We seek the set of universal TERMs that capture all structure in the Protein Data Bank (PDB), finding remarkable degeneracy. Only ∼600 TERMs are sufficient to describe 50% of the PDB at sub-Angstrom resolution. However, more rare geometries also exist, and the overall structural coverage grows logarithmically with the number of TERMs. We go on to show that universal TERMs provide an effective mapping between sequence and structure. We demonstrate that TERM-based statistics alone are sufficient to recapitulate close-to-native sequences given either NMR or X-ray backbones. Furthermore, sequence variability predicted from TERM data agrees closely with evolutionary variation. Finally, locations of TERMs in protein chains can be predicted from sequence alone based on sequence signatures emergent from TERM instances in the PDB. For multisegment motifs, this method identifies spatially adjacent fragments that are not contiguous in sequence-a major bottleneck in structure prediction. Although all TERMs recur in diverse proteins, some appear specialized for certain functions, such as interface formation, metal coordination, or even water binding. Structural biology has benefited greatly from previously observed degeneracies in structure. The decomposition of the known structural universe into a finite set of compact TERMs offers exciting opportunities toward better understanding, design, and prediction of protein structure.
Collapse
|
238
|
Nastri F, Chino M, Maglio O, Bhagi-Damodaran A, Lu Y, Lombardi A. Design and engineering of artificial oxygen-activating metalloenzymes. Chem Soc Rev 2016; 45:5020-54. [PMID: 27341693 PMCID: PMC5021598 DOI: 10.1039/c5cs00923e] [Citation(s) in RCA: 133] [Impact Index Per Article: 16.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Many efforts are being made in the design and engineering of metalloenzymes with catalytic properties fulfilling the needs of practical applications. Progress in this field has recently been accelerated by advances in computational, molecular and structural biology. This review article focuses on the recent examples of oxygen-activating metalloenzymes, developed through the strategies of de novo design, miniaturization processes and protein redesign. Considerable progress in these diverse design approaches has produced many metal-containing biocatalysts able to adopt the functions of native enzymes or even novel functions beyond those found in Nature.
Collapse
Affiliation(s)
- Flavia Nastri
- Department of Chemical Sciences, University of Naples “Federico II”, Via Cintia, 80126 Naples, Italy
| | - Marco Chino
- Department of Chemical Sciences, University of Naples “Federico II”, Via Cintia, 80126 Naples, Italy
| | - Ornella Maglio
- Department of Chemical Sciences, University of Naples “Federico II”, Via Cintia, 80126 Naples, Italy
- IBB, CNR, Via Mezzocannone 16, 80134 Naples, Italy
| | - Ambika Bhagi-Damodaran
- Department of Chemistry, University of Illinois at Urbana-Champaign, A322 CLSL, 600 South Mathews Avenue, Urbana, IL 61801
| | - Yi Lu
- Department of Chemistry, University of Illinois at Urbana-Champaign, A322 CLSL, 600 South Mathews Avenue, Urbana, IL 61801
| | - Angela Lombardi
- Department of Chemical Sciences, University of Naples “Federico II”, Via Cintia, 80126 Naples, Italy
| |
Collapse
|
239
|
Piepenbrink KH, Lillehoj E, Harding CM, Labonte JW, Zuo X, Rapp CA, Munson RS, Goldblum SE, Feldman MF, Gray JJ, Sundberg EJ. Structural Diversity in the Type IV Pili of Multidrug-resistant Acinetobacter. J Biol Chem 2016; 291:22924-22935. [PMID: 27634041 PMCID: PMC5087714 DOI: 10.1074/jbc.m116.751099] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2016] [Indexed: 11/06/2022] Open
Abstract
Acinetobacter baumannii is a Gram-negative coccobacillus found primarily in hospital settings that has recently emerged as a source of hospital-acquired infections. A. baumannii expresses a variety of virulence factors, including type IV pili, bacterial extracellular appendages often essential for attachment to host cells. Here, we report the high resolution structures of the major pilin subunit, PilA, from three Acinetobacter strains, demonstrating that A. baumannii subsets produce morphologically distinct type IV pilin glycoproteins. We examine the consequences of this heterogeneity for protein folding and assembly as well as host-cell adhesion by Acinetobacter Comparisons of genomic and structural data with pilin proteins from other species of soil gammaproteobacteria suggest that these structural differences stem from evolutionary pressure that has resulted in three distinct classes of type IVa pilins, each found in multiple species.
Collapse
Affiliation(s)
| | | | - Christian M Harding
- Department of Molecular Microbiology, Washington University School of Medicine, St. Louis, Missouri 63110
| | - Jason W Labonte
- Department of Chemical and Biomolecular Engineering, Whiting School of Engineering, The Johns Hopkins University, Baltimore, Maryland 21218
| | - Xiaotong Zuo
- Department of Chemical and Biomolecular Engineering, Whiting School of Engineering, The Johns Hopkins University, Baltimore, Maryland 21218
| | | | - Robert S Munson
- The Center for Microbial Pathogenesis in the Research Institute at Nationwide Children's Hospital and Department of Pediatrics, The Ohio State University College of Medicine, Columbus, Ohio 43205, and
| | - Simeon E Goldblum
- Departments of Medicine.,Baltimore Veterans Affairs Medical Center, Baltimore, Maryland 21201.,Pathology, University of Maryland School of Medicine, Baltimore, Maryland 21201
| | - Mario F Feldman
- Department of Molecular Microbiology, Washington University School of Medicine, St. Louis, Missouri 63110
| | - Jeffrey J Gray
- Department of Chemical and Biomolecular Engineering, Whiting School of Engineering, The Johns Hopkins University, Baltimore, Maryland 21218
| | - Eric J Sundberg
- From the Institute of Human Virology and .,Departments of Medicine.,Microbiology and Immunology
| |
Collapse
|
240
|
The coming of age of de novo protein design. Nature 2016; 537:320-7. [DOI: 10.1038/nature19946] [Citation(s) in RCA: 803] [Impact Index Per Article: 100.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2016] [Accepted: 07/20/2016] [Indexed: 12/24/2022]
|
241
|
Kuroda D, Gray JJ. Pushing the Backbone in Protein-Protein Docking. Structure 2016; 24:1821-1829. [PMID: 27568930 DOI: 10.1016/j.str.2016.06.025] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2015] [Revised: 05/17/2016] [Accepted: 06/19/2016] [Indexed: 02/01/2023]
Abstract
Conformational changes of proteins that occur upon binding typically confound computational docking algorithms. In this study, we test computational methods to capture protein backbone conformational change related to binding. To address how well existing algorithms can sample bound-like backbones, we query seven techniques including Monte Carlo-based sampling, molecular dynamics, and normal mode analysis. All methods tested rarely sample near-bound states from the unbound conformation. Nevertheless, the direction of the predicted motions overlap with the actual conformational change. We next forced the backbone from the unbound toward the bound conformation to create a family of docking energy landscapes. Seventy percent of docking targets succeed when the unbound backbones is pushed to within 0.6 Å of the bound. Current methods can capture an average of 22% of unbound-bound transitions through conformer selection methods and another 57% through induced-fit methodologies, delineating a stubborn gap (21%) in backbone motion not covered by any current approach.
Collapse
Affiliation(s)
- Daisuke Kuroda
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, MD 21218, USA; Department of Analytical and Physical Chemistry, Showa University School of Pharmacy, Tokyo 142-8555, Japan
| | - Jeffrey J Gray
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, MD 21218, USA; Program in Molecular Biophysics, Johns Hopkins University, Baltimore, MD 21218, USA.
| |
Collapse
|
242
|
Bender BJ, Cisneros A, Duran AM, Finn JA, Fu D, Lokits AD, Mueller BK, Sangha AK, Sauer MF, Sevy AM, Sliwoski G, Sheehan JH, DiMaio F, Meiler J, Moretti R. Protocols for Molecular Modeling with Rosetta3 and RosettaScripts. Biochemistry 2016; 55:4748-63. [PMID: 27490953 PMCID: PMC5007558 DOI: 10.1021/acs.biochem.6b00444] [Citation(s) in RCA: 137] [Impact Index Per Article: 17.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
![]()
Previously, we published an article
providing an overview of the
Rosetta suite of biomacromolecular modeling software and a series
of step-by-step tutorials [Kaufmann, K. W., et al. (2010) Biochemistry 49, 2987–2998]. The overwhelming positive
response to this publication we received motivates us to here share
the next iteration of these tutorials that feature de novo folding, comparative modeling, loop construction, protein docking,
small molecule docking, and protein design. This updated and expanded
set of tutorials is needed, as since 2010 Rosetta has been fully redesigned
into an object-oriented protein modeling program Rosetta3. Notable
improvements include a substantially improved energy function, an
XML-like language termed “RosettaScripts” for flexibly
specifying modeling task, new analysis tools, the addition of the
TopologyBroker to control conformational sampling, and support for
multiple templates in comparative modeling. Rosetta’s ability
to model systems with symmetric proteins, membrane proteins, noncanonical
amino acids, and RNA has also been greatly expanded and improved.
Collapse
Affiliation(s)
- Brian J Bender
- Department of Pharmacology, Vanderbilt University , Nashville, Tennessee 37232-6600, United States.,Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States
| | - Alberto Cisneros
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Chemical and Physical Biology Program, Vanderbilt University , Nashville, Tennessee 37232-0301, United States
| | - Amanda M Duran
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Department of Chemistry, Vanderbilt University , Nashville, Tennessee 37235, United States
| | - Jessica A Finn
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Department of Pathology, Microbiology and Immunology, Vanderbilt University , Nashville, Tennessee 37232-2561, United States
| | - Darwin Fu
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Department of Chemistry, Vanderbilt University , Nashville, Tennessee 37235, United States
| | - Alyssa D Lokits
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Neuroscience Program, Vanderbilt University , Nashville, Tennessee 37235, United States
| | - Benjamin K Mueller
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Department of Chemistry, Vanderbilt University , Nashville, Tennessee 37235, United States
| | - Amandeep K Sangha
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Department of Chemistry, Vanderbilt University , Nashville, Tennessee 37235, United States
| | - Marion F Sauer
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Chemical and Physical Biology Program, Vanderbilt University , Nashville, Tennessee 37232-0301, United States
| | - Alexander M Sevy
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Chemical and Physical Biology Program, Vanderbilt University , Nashville, Tennessee 37232-0301, United States
| | - Gregory Sliwoski
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Department of Chemistry, Vanderbilt University , Nashville, Tennessee 37235, United States
| | - Jonathan H Sheehan
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States
| | - Frank DiMaio
- Department of Biochemistry, University of Washington , Seattle, Washington 98195, United States
| | - Jens Meiler
- Department of Pharmacology, Vanderbilt University , Nashville, Tennessee 37232-6600, United States.,Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Chemical and Physical Biology Program, Vanderbilt University , Nashville, Tennessee 37232-0301, United States.,Department of Chemistry, Vanderbilt University , Nashville, Tennessee 37235, United States.,Department of Pathology, Microbiology and Immunology, Vanderbilt University , Nashville, Tennessee 37232-2561, United States.,Neuroscience Program, Vanderbilt University , Nashville, Tennessee 37235, United States
| | - Rocco Moretti
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Department of Chemistry, Vanderbilt University , Nashville, Tennessee 37235, United States
| |
Collapse
|
243
|
Rakesh R, Joseph AP, Bhaskara RM, Srinivasan N. Structural and mechanistic insights into human splicing factor SF3b complex derived using an integrated approach guided by the cryo-EM density maps. RNA Biol 2016; 13:1025-1040. [PMID: 27618338 DOI: 10.1080/15476286.2016.1218590] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Abstract
Pre-mRNA splicing in eukaryotes is performed by the spliceosome, a highly complex macromolecular machine. SF3b is a multi-protein complex which recognizes the branch point adenosine of pre-mRNA as part of a larger U2 snRNP or U11/U12 di-snRNP in the dynamic spliceosome machinery. Although a cryo-EM map is available for human SF3b complex, the structure and relative spatial arrangement of all components in the complex are not yet known. We have recognized folds of domains in various proteins in the assembly and generated comparative models. Using an integrative approach involving structural and other experimental data, guided by the available cryo-EM density map, we deciphered a pseudo-atomic model of the closed form of SF3b which is found to be a "fuzzy complex" with highly flexible components and multiplicity of folds. Further, the model provides structural information for 5 proteins (SF3b10, SF3b155, SF3b145, SF3b130 and SF3b14b) and localization information for 4 proteins (SF3b10, SF3b145, SF3b130 and SF3b14b) in the assembly for the first time. Integration of this model with the available U11/U12 di-snRNP cryo-EM map enabled elucidation of an open form. This now provides new insights on the mechanistic features involved in the transition between closed and open forms pivoted by a hinge region in the SF3b155 protein that also harbors cancer causing mutations. Moreover, the open form guided model of the 5' end of U12 snRNA, which includes the branch point duplex, shows that the architecture of SF3b acts as a scaffold for U12 snRNA: pre-mRNA branch point duplex formation with potential implications for branch point adenosine recognition fidelity.
Collapse
Affiliation(s)
- Ramachandran Rakesh
- a Molecular Biophysics Unit, Indian Institute of Science , Bangalore , India
| | - Agnel Praveen Joseph
- b National Center for Biological Sciences, TIFR, GKVK Campus , Bangalore , India
| | - Ramachandra M Bhaskara
- a Molecular Biophysics Unit, Indian Institute of Science , Bangalore , India.,b National Center for Biological Sciences, TIFR, GKVK Campus , Bangalore , India
| | | |
Collapse
|
244
|
Wright KM, Du H, Dagnachew M, Massiah MA. Solution structure of the microtubule-targeting COS domain of MID1. FEBS J 2016; 283:3089-102. [PMID: 27367845 DOI: 10.1111/febs.13795] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2016] [Revised: 05/27/2016] [Accepted: 06/29/2016] [Indexed: 11/27/2022]
Abstract
UNLABELLED The human MID1 protein is required for the proper development during embryogenesis. Mutations of MID1 are associated with X-linked Opitz G syndrome, characterized by midline anomalies. MID1 associates with the microtubules and functions as an ubiquitin E3 ligase, targeting protein phosphatase 2A for ubiquitin-mediated regulation. The mechanism of microtubule association is not known. Recently, a 60-amino acid region termed the C-terminal subgroup One Signature (COS) box/domain was identified at the C-terminal end of the coiled-coil (CC) domain that facilitates microtubule localization. Insertion of the MID1 COS domain at the C-terminal end of the CC domain of a nonmicrotubule-associated TRIM protein confers microtubule localization. Here, we report the solution structure of the COS domain of MID1. The domain adopts a helix-loop-helix structure in which the N- and C-terminal ends are in close proximity. Hydrophobic residues stabilizing the interaction of the two α-helices form a central hydrophobic core. The loop separating the α-helices is structured, with two of its hydrophobic residues making contact with the central core. On the outer surface, positively charged residues form a distinct basic patch near the termini that we postulate is important for microtubule binding. A model of the structure of the preceding coiled-coil and COS domains (CC-COS) show that the COS domain forms a helical bundle at the C-terminal end of the CC domain similar to the spectrin-like fold observed with some known microtubule-binding proteins. Interestingly, the CC-COS domains bind to microtubules, demonstrating for the first time that MID1 can directly associate with the microtubules. DATABASE Structural data are available in PDB database under the accession number 5IM8.
Collapse
Affiliation(s)
- Katharine M Wright
- Department of Chemistry and Center of Biomolecular Sciences, George Washington University, DC, USA
| | - Haijuan Du
- Department of Chemistry and Center of Biomolecular Sciences, George Washington University, DC, USA
| | - Mesgana Dagnachew
- Department of Chemistry and Center of Biomolecular Sciences, George Washington University, DC, USA
| | - Michael A Massiah
- Department of Chemistry and Center of Biomolecular Sciences, George Washington University, DC, USA
| |
Collapse
|
245
|
Topham CM, Barbe S, André I. An Atomistic Statistically Effective Energy Function for Computational Protein Design. J Chem Theory Comput 2016; 12:4146-68. [PMID: 27341125 DOI: 10.1021/acs.jctc.6b00090] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Shortcomings in the definition of effective free-energy surfaces of proteins are recognized to be a major contributory factor responsible for the low success rates of existing automated methods for computational protein design (CPD). The formulation of an atomistic statistically effective energy function (SEEF) suitable for a wide range of CPD applications and its derivation from structural data extracted from protein domains and protein-ligand complexes are described here. The proposed energy function comprises nonlocal atom-based and local residue-based SEEFs, which are coupled using a novel atom connectivity number factor to scale short-range, pairwise, nonbonded atomic interaction energies and a surface-area-dependent cavity energy term. This energy function was used to derive additional SEEFs describing the unfolded-state ensemble of any given residue sequence based on computed average energies for partially or fully solvent-exposed fragments in regions of irregular structure in native proteins. Relative thermal stabilities of 97 T4 bacteriophage lysozyme mutants were predicted from calculated energy differences for folded and unfolded states with an average unsigned error (AUE) of 0.84 kcal mol(-1) when compared to experiment. To demonstrate the utility of the energy function for CPD, further validation was carried out in tests of its capacity to recover cognate protein sequences and to discriminate native and near-native protein folds, loop conformers, and small-molecule ligand binding poses from non-native benchmark decoys. Experimental ligand binding free energies for a diverse set of 80 protein complexes could be predicted with an AUE of 2.4 kcal mol(-1) using an additional energy term to account for the loss in ligand configurational entropy upon binding. The atomistic SEEF is expected to improve the accuracy of residue-based coarse-grained SEEFs currently used in CPD and to extend the range of applications of extant atom-based protein statistical potentials.
Collapse
Affiliation(s)
- Christopher M Topham
- Université de Toulouse; INSA, UPS, INP; LISBP , 135 Avenue de Rangueil, F-31077 Toulouse, France.,CNRS, UMR5504 , F-31400 Toulouse, France.,INRA, UMR792 Ingénierie des Systèmes Biologiques et des Procédés , F-31400 Toulouse, France
| | - Sophie Barbe
- Université de Toulouse; INSA, UPS, INP; LISBP , 135 Avenue de Rangueil, F-31077 Toulouse, France.,CNRS, UMR5504 , F-31400 Toulouse, France.,INRA, UMR792 Ingénierie des Systèmes Biologiques et des Procédés , F-31400 Toulouse, France
| | - Isabelle André
- Université de Toulouse; INSA, UPS, INP; LISBP , 135 Avenue de Rangueil, F-31077 Toulouse, France.,CNRS, UMR5504 , F-31400 Toulouse, France.,INRA, UMR792 Ingénierie des Systèmes Biologiques et des Procédés , F-31400 Toulouse, France
| |
Collapse
|
246
|
Multiphase Simulated Annealing Based on Boltzmann and Bose-Einstein Distribution Applied to Protein Folding Problem. Adv Bioinformatics 2016; 2016:7357123. [PMID: 27413369 PMCID: PMC4931103 DOI: 10.1155/2016/7357123] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2015] [Revised: 04/05/2016] [Accepted: 04/19/2016] [Indexed: 11/17/2022] Open
Abstract
A new hybrid Multiphase Simulated Annealing Algorithm using Boltzmann and Bose-Einstein distributions (MPSABBE) is proposed. MPSABBE was designed for solving the Protein Folding Problem (PFP) instances. This new approach has four phases: (i) Multiquenching Phase (MQP), (ii) Boltzmann Annealing Phase (BAP), (iii) Bose-Einstein Annealing Phase (BEAP), and (iv) Dynamical Equilibrium Phase (DEP). BAP and BEAP are simulated annealing searching procedures based on Boltzmann and Bose-Einstein distributions, respectively. DEP is also a simulated annealing search procedure, which is applied at the final temperature of the fourth phase, which can be seen as a second Bose-Einstein phase. MQP is a search process that ranges from extremely high to high temperatures, applying a very fast cooling process, and is not very restrictive to accept new solutions. However, BAP and BEAP range from high to low and from low to very low temperatures, respectively. They are more restrictive for accepting new solutions. DEP uses a particular heuristic to detect the stochastic equilibrium by applying a least squares method during its execution. MPSABBE parameters are tuned with an analytical method, which considers the maximal and minimal deterioration of problem instances. MPSABBE was tested with several instances of PFP, showing that the use of both distributions is better than using only the Boltzmann distribution on the classical SA.
Collapse
|
247
|
Simpkin AJ, Rigden DJ. GP0.4 from bacteriophage T7: in silico characterisation of its structure and interaction with E. coli FtsZ. BMC Res Notes 2016; 9:343. [PMID: 27411831 PMCID: PMC4944311 DOI: 10.1186/s13104-016-2149-5] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2016] [Accepted: 07/06/2016] [Indexed: 01/21/2023] Open
Abstract
Background Proteins produced by bacteriophages can have potent antimicrobial activity. The study of phage-host interactions can therefore inform small molecule drug discovery by revealing and characterising new drug targets. Here we characterise in silico the predicted interaction of gene protein 0.4 (GP0.4) from the Escherichia coli (E. coli) phage T7 with E. coli filamenting temperature-sensitive mutant Z division protein (FtsZ). FtsZ is a tubulin homolog which plays a key role in bacterial cell division and that has been proposed as a drug target. Results Using ab initio, fragment assembly structure modelling, we predicted the structure of GP0.4 with two programs. A structure similarity-based network was used to identify a U-shaped helix-turn-helix candidate fold as being favoured. ClusPro was used to dock this structure prediction to a homology model of E. coli FtsZ resulting in a favourable predicted interaction mode. Alternative docking methods supported the proposed mode which offered an immediate explanation for the anti-filamenting activity of GP0.4. Importantly, further strong support derived from a previously characterised insertion mutation, known to abolish GP0.4 activity, that is positioned in close proximity to the proposed GP0.4/FtsZ interface. Conclusions The mode of interaction predicted by bioinformatics techniques strongly suggests a mechanism through which GP0.4 inhibits FtsZ and further establishes the latter’s druggable intrafilament interface as a potential drug target. Electronic supplementary material The online version of this article (doi:10.1186/s13104-016-2149-5) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Adam J Simpkin
- Institute of Integrative Biology, University of Liverpool, Liverpool, L69 7ZB, UK
| | - Daniel J Rigden
- Institute of Integrative Biology, University of Liverpool, Liverpool, L69 7ZB, UK.
| |
Collapse
|
248
|
Bottaro S, Gil-Ley A, Bussi G. RNA folding pathways in stop motion. Nucleic Acids Res 2016; 44:5883-91. [PMID: 27091499 PMCID: PMC4937309 DOI: 10.1093/nar/gkw239] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2016] [Revised: 03/29/2016] [Accepted: 03/29/2016] [Indexed: 11/19/2022] Open
Abstract
We introduce a method for predicting RNA folding pathways, with an application to the most important RNA tetraloops. The method is based on the idea that ensembles of three-dimensional fragments extracted from high-resolution crystal structures are heterogeneous enough to describe metastable as well as intermediate states. These ensembles are first validated by performing a quantitative comparison against available solution nuclear magnetic resonance (NMR) data of a set of RNA tetranucleotides. Notably, the agreement is better with respect to the one obtained by comparing NMR with extensive all-atom molecular dynamics simulations. We then propose a procedure based on diffusion maps and Markov models that makes it possible to obtain reaction pathways and their relative probabilities from fragment ensembles. This approach is applied to study the helix-to-loop folding pathway of all the tetraloops from the GNRA and UNCG families. The results give detailed insights into the folding mechanism that are compatible with available experimental data and clarify the role of intermediate states observed in previous simulation studies. The method is computationally inexpensive and can be used to study arbitrary conformational transitions.
Collapse
Affiliation(s)
- Sandro Bottaro
- Scuola Internazionale Superiore di Studi Avanzati, International School for Advanced Studies, 265, Via Bonomea I-34136 Trieste, Italy
| | - Alejandro Gil-Ley
- Scuola Internazionale Superiore di Studi Avanzati, International School for Advanced Studies, 265, Via Bonomea I-34136 Trieste, Italy
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, International School for Advanced Studies, 265, Via Bonomea I-34136 Trieste, Italy
| |
Collapse
|
249
|
Simkovic F, Thomas JMH, Keegan RM, Winn MD, Mayans O, Rigden DJ. Residue contacts predicted by evolutionary covariance extend the application of ab initio molecular replacement to larger and more challenging protein folds. IUCRJ 2016; 3:259-70. [PMID: 27437113 PMCID: PMC4937781 DOI: 10.1107/s2052252516008113] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/10/2016] [Accepted: 05/18/2016] [Indexed: 05/05/2023]
Abstract
For many protein families, the deluge of new sequence information together with new statistical protocols now allow the accurate prediction of contacting residues from sequence information alone. This offers the possibility of more accurate ab initio (non-homology-based) structure prediction. Such models can be used in structure solution by molecular replacement (MR) where the target fold is novel or is only distantly related to known structures. Here, AMPLE, an MR pipeline that assembles search-model ensembles from ab initio structure predictions ('decoys'), is employed to assess the value of contact-assisted ab initio models to the crystallographer. It is demonstrated that evolutionary covariance-derived residue-residue contact predictions improve the quality of ab initio models and, consequently, the success rate of MR using search models derived from them. For targets containing β-structure, decoy quality and MR performance were further improved by the use of a β-strand contact-filtering protocol. Such contact-guided decoys achieved 14 structure solutions from 21 attempted protein targets, compared with nine for simple Rosetta decoys. Previously encountered limitations were superseded in two key respects. Firstly, much larger targets of up to 221 residues in length were solved, which is far larger than the previously benchmarked threshold of 120 residues. Secondly, contact-guided decoys significantly improved success with β-sheet-rich proteins. Overall, the improved performance of contact-guided decoys suggests that MR is now applicable to a significantly wider range of protein targets than were previously tractable, and points to a direct benefit to structural biology from the recent remarkable advances in sequencing.
Collapse
Affiliation(s)
- Felix Simkovic
- Institute of Integrative Biology, University of Liverpool, Liverpool L69 7ZB, England
| | - Jens M. H. Thomas
- Institute of Integrative Biology, University of Liverpool, Liverpool L69 7ZB, England
| | - Ronan M. Keegan
- Research Complex at Harwell, STFC Rutherford Appleton Laboratory, Didcot OX11 0FA, England
| | - Martyn D. Winn
- Science and Technology Facilities Council, Daresbury Laboratory, Warrington WA4 4AD, England
| | - Olga Mayans
- Institute of Integrative Biology, University of Liverpool, Liverpool L69 7ZB, England
| | - Daniel J. Rigden
- Institute of Integrative Biology, University of Liverpool, Liverpool L69 7ZB, England
| |
Collapse
|
250
|
Lindert S, McCammon JA. Improved cryoEM-Guided Iterative Molecular Dynamics--Rosetta Protein Structure Refinement Protocol for High Precision Protein Structure Prediction. J Chem Theory Comput 2016; 11:1337-46. [PMID: 25883538 PMCID: PMC4393324 DOI: 10.1021/ct500995d] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2014] [Indexed: 12/13/2022]
Abstract
![]()
Many excellent methods exist that
incorporate cryo-electron microscopy
(cryoEM) data to constrain computational protein structure prediction
and refinement. Previously, it was shown that iteration of two such
orthogonal sampling and scoring methods – Rosetta and molecular
dynamics (MD) simulations – facilitated exploration of conformational
space in principle. Here, we go beyond a proof-of-concept study and
address significant remaining limitations of the iterative MD–Rosetta
protein structure refinement protocol. Specifically, all parts of
the iterative refinement protocol are now guided by medium-resolution
cryoEM density maps, and previous knowledge about the native structure
of the protein is no longer necessary. Models are identified solely
based on score or simulation time. All four benchmark proteins showed
substantial improvement through three rounds of the iterative refinement
protocol. The best-scoring final models of two proteins had sub-Ångstrom
RMSD to the native structure over residues in secondary structure
elements. Molecular dynamics was most efficient in refining secondary
structure elements and was thus highly complementary to the Rosetta
refinement which is most powerful in refining side chains and loop
regions.
Collapse
|